CN107947178A - A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm - Google Patents
A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm Download PDFInfo
- Publication number
- CN107947178A CN107947178A CN201711359910.9A CN201711359910A CN107947178A CN 107947178 A CN107947178 A CN 107947178A CN 201711359910 A CN201711359910 A CN 201711359910A CN 107947178 A CN107947178 A CN 107947178A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- intelligent body
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000000446 fuel Substances 0.000 claims abstract description 54
- 238000005457 optimization Methods 0.000 claims abstract description 25
- 230000005611 electricity Effects 0.000 claims abstract description 23
- 230000008901 benefit Effects 0.000 claims abstract description 22
- 238000002485 combustion reaction Methods 0.000 claims abstract description 21
- 230000007613 environmental effect Effects 0.000 claims abstract description 20
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 238000007599 discharging Methods 0.000 claims abstract description 13
- 230000009194 climbing Effects 0.000 claims abstract description 12
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 239000003795 chemical substances by application Substances 0.000 claims description 28
- 210000000746 body region Anatomy 0.000 claims description 18
- 238000012423 maintenance Methods 0.000 claims description 15
- 239000003344 environmental pollutant Substances 0.000 claims description 9
- 231100000719 pollutant Toxicity 0.000 claims description 9
- 238000006116 polymerization reaction Methods 0.000 claims description 9
- 230000006978 adaptation Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 239000002737 fuel gas Substances 0.000 claims description 3
- 230000000379 polymerizing effect Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 208000034530 PLAA-associated neurodevelopmental disease Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/02—Circuit arrangements for ac mains or ac distribution networks using a single network for simultaneous distribution of power at different frequencies; using a single network for simultaneous distribution of ac power and of dc power
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/123—DNA computing
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- H02J3/383—
-
- H02J3/386—
-
- H02J3/387—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- Y—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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- Y—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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- Y—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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Power Engineering (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Genetics & Genomics (AREA)
- Secondary Cells (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm:For including wind-power electricity generation, photovoltaic generation, storage battery, miniature combustion engine, the alternating current-direct current mixing microgrid of fuel cell and diesel-driven generator, establish multiple target, multiple constraint, nonlinear optimization operation mathematical model, object function considers power grid purchases strategies, micro- source fuel cost, environmental benefit cost, network loss and operation expense, and obey microgrid internal power balance, points of common connection transmission capacity, controllable micro- source climbing rate, unit interval accumulator cell charging and discharging bound, storage battery charge state bound, storage battery surrounding time section power-balance and the constant constraints of storage battery charge state;Using cultural gene Algorithm for Solving alternating current-direct current mixing microgrid mathematical model;Verify the correctness and validity of the alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm.The present invention can effectively solve multiple target, the nonlinear optimization objective function of multiple constraint, improve microgrid economic benefit and environmental benefit.
Description
Technical field
The present invention relates to a kind of alternating current-direct current mixing microgrid optimizing operation method.Calculated more particularly to one kind based on cultural gene
The alternating current-direct current mixing microgrid optimizing operation method of method.
Background technology
With traditional energy increasingly depleted, the new generation mode such as wind-power electricity generation, photovoltaic generation is due to its good environment
Benefit is paid close attention to be subject to domestic and foreign scholars.Microgrid is the new supply network for including the devices such as distributed energy, load, energy storage,
Since it possesses complementary the advantages that utilizing for dissolving in time, realizing various forms of micro- sources to new energy, it has also become both at home and abroad
The hot spot of concern.At present, load species is increasingly abundanter, and the power supply reliability of DC load becomes a big academic problem.
Alternating current-direct current mixing microgrid has merged AC load and DC load, it is possible to increase power quality, suitably reduces electric power
The use of electronic device is so as to reduce harmonic pollution.In order to ensure the safe and reliable operation of alternating current-direct current mixing microgrid, domestic foreign minister
Scholar is closed just to study in development and application of the Efforts To Develop for alternating current-direct current mixing microgrid operation control technology.Alternating current-direct current mixing is micro-
The operation of net economic optimization is one of research topic, and with the complicated variation of microgrid, it is better to find optimization performance
Good intelligent search algorithm is crucial.
Currently used intelligent algorithm optimizing works well, but it is slow or do not reach convergence essence to still suffer from some convergence rates
The problems such as degree requires, it is therefore desirable to attempt new intelligent algorithm.There is a kind of evolutional algorithm in intelligent algorithm, calculated similar to heredity
Method, immune algorithm etc..This kind of algorithm is based on the rule of " survival of the fittest " in evolutionism, and the excellent genes in parent are entailed down
A generation.This theory is generalized in the development of culture in relation to scholar, foreign countries are referred to as " cultural volution ".Cultural gene algorithm is just
It is to be based on this thought, algorithm iteration process is divided into local search and global search, and it is a set of individually theoretical right each to have
Individual is screened, so as to accelerate individual to approach optimal value from two aspects, improves algorithmic statement precision.
The content of the invention
The technical problem to be solved by the invention is to provide one kind can effectively solve multidimensional, multiple target, multiple constraint and non-thread
The alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm of the object function of property.
The technical solution adopted in the present invention is:A kind of alternating current-direct current mixing microgrid optimization operation based on cultural gene algorithm
Method, includes the following steps:
1) it is straight for the friendship comprising wind-power electricity generation, photovoltaic generation, storage battery, miniature combustion engine, fuel cell and diesel-driven generator
Stream mixing microgrid, establish multiple target, multiple constraint, it is nonlinear optimization operation mathematical model, object function consider power grid power purchase into
Sheet, micro- source fuel cost, environmental benefit cost, network loss and operation expense, and obey microgrid internal power and balance, is public
Tie point transmission capacity, controllable micro- source climbing rate, unit interval accumulator cell charging and discharging bound, above and below storage battery charge state
Limit, the constraints that storage battery surrounding time section power-balance and storage battery charge state are constant;
2) cultural gene Algorithm for Solving alternating current-direct current mixing microgrid mathematical model is used;
3) correctness and validity of the alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm are verified.
Object function F described in step 1) is as follows:
F=FGrid+Fec+Floss+Fom
FGrid=FACGrid+FDCGrid
Fec=FACec+FDCec
Floss=FACloss+FDCloss
Fom=FACom+FDCom
Wherein, FGridFor power grid purchases strategies;FecFor financial cost;FlossFor network loss;FomFor equipment operation maintenance cost;
FACGridFor exchanging area power grid purchases strategies;FACecFor exchanging area financial cost, by the micro- source fuel cost F in exchanging areaACfuelAnd ring
Border benefit-cost FACenAddition obtains;FAClossFor exchanging area network loss;FAComFor exchanging area equipment operation maintenance cost;FDCGridFor
Exchanging area power grid purchases strategies;FDCecFor the micro- source financial cost in exchanging area, by the micro- source fuel cost F in DC areaDCfuelAnd Environmental Effect
Beneficial cost FDCenAddition obtains;FDClossFor exchanging area network loss;FDComFor exchanging area equipment operation maintenance cost;Wherein,
(1) exchanging area power grid purchases strategies FACGridWith DC area power grid purchases strategies FDCGridRespectively:
Wherein, △ TACPeriod for from exchanging area to power grid power purchase;T is total period in one day;T is day part in one day;It is exchanging area to power grid power purchase electricity;△TDCPeriod for from DC area to power grid power purchase;It is DC area to electricity
Net purchase power consumption;For power grid sale of electricity price in the corresponding period;
(2) exchanging area financial cost FACecWith DC area financial cost FDCec:
FACec=FACfuel+FACen
FDCec=FDCfuel+FDCen
Exchanging area includes wind-power electricity generation, miniature combustion engine and diesel-driven generator, and DC area includes photovoltaic generation, fuel cell and storage
Battery, exchanging area fuel cost FACfuelWith DC area fuel cost FDCfuelIt can be calculated by the following formula:
FACfuel=FMTfuel+FDEGfuel
FDCfuel=FFCfuel
ηFC=-0.0023 × PFC+0.6735
Wherein, FMTfuelFor miniature combustion engine fuel cost;CMTFor miniature combustion engine cooler fuel price;LHV is the low heat value of fuel gas;
PMTFor miniature combustion engine output power;ηMTFor the generating efficiency of miniature combustion engine;FFCfuelFor the fuel cost of fuel cell;CFCFor fuel electricity
The cooler fuel price in pond;PFCFor the output power of fuel cell;ηFCFor the generating efficiency of fuel cell;FDEGfuelFor diesel-driven generator
Fuel cost;PDEGFor the output power of diesel-driven generator;A, b, c are diesel-driven generator power generation coefficient, by diesel engine factory
Family provides;
Exchanging area environmental benefit cost FACenWith DC area environmental benefit cost FDCenIt is calculated by the following formula:
Wherein, n1 is the micro- source number in exchanging area;M is the type of pollutant;αjFor the conversion cost of corresponding pollutant, member/
kg;EFi,jThe unit discharge of the jth kind pollutant produced for i-th of micro- source, kg/kW;PiFor the output work in i-th of micro- source
Rate;n2For the micro- source number in DC area.
(3) network loss F is exchangedAClossWith direct current network loss FDCloss:
Wherein, L1 is exchanging area branch sum;Pk、QkActive power, reactive power for branch k transmission;L2 is DC area
Branch sum;RkFor the resistance of branch k;UkFor the voltage effective value of branch k;
(4) ac operation maintenance cost FAComWith DC operation maintenance cost FDCom:
Wherein, β i are the operation expense coefficient in i-th of micro- source.
The method of weighting is taken, finally obtaining object function is:
MinF=(1- λ) × (FGrid+Fec+Fom)+λFloss
Wherein, λ is Web-based exercise coefficient.
The Web-based exercise coefficient lambda is 0.1.
Constraints described in step 1) includes:
(1) microgrid internal power equilibrium constraint
(2) points of common connection transmission capacity constraints
(3) controllable micro- source climbing rate constraints
(4) unit interval accumulator cell charging and discharging bound constraints
(5) storage battery charge state (State of Charge, SOC) bound constraints
Socmin≤Soct≤Socmax
(6) storage battery surrounding time section power-balance constraint condition
(7) storage battery charge state constraint independent of time condition
Socinitial=Socend
Wherein, N is micro- source number;Represent that t-th of period, i-th of micro- source is contributed;Represent that t-th of period is purchased from power grid
Electrical power;For t-th of period battery power variable quantity, discharge just, to be charged as bearing;T is represented respectively
A period AC and DC area payload;For the net flow power of t-th of period points of common connection;PG,maxFor points of common connection
Transmission capacity limit value;riFor the unit interval climbing rate in i-th of micro- source;For the exhausted of t-th period battery power variable quantity
To value;For the limit value of t-th of period battery power variable quantity;SoctFor the state-of-charge of storage battery in t-th of period;
Socmin、SocmaxThe respectively upper lower limit value of state-of-charge;Uu represents charge and discharge electrostrictive coefficient, is 1 during charging, and when electric discharge is -1;η is
Accumulator cell charging and discharging efficiency, takes 95% here;QESFor battery rating;Socinitial、SocendFor the initial of state-of-charge
Value and end value.
Step 2) includes:
(1) initialize
Initial individuals are produced using random generating mode, the value range of each variable is represented with varmin, varmax, its
In any one variable xmInitial value drawn by following formula:
xm=varminm+rand(0,1)*(varmaxm-varminm)
If M is all individual numbers in colony, then m=1,2 ..., M, MagentFor the number of intelligent body, MpublicTo be common
The number of individual, that is, have:
M=Magent+Mpublic
Ascending order arrangement, preceding M are carried out according to fitness value size to initial individualsagentAs intelligent body, remaining is common
Individual, wherein intelligent body are that fitness value comes preceding M in colonyagentThe individual of position, remaining individual is average individual in colony,
Colony is divided into MagentRegion, average individual belong to the region where each intelligent body, and each intelligent body initial time is possessed
The number of average individual be to be determined by each intelligent body with respect to strength, a-th intelligent body (a=1,2 ..., Magent) it is relatively real
Power is specifically calculated by following formula:
Sa=max { sb}-sa, b=1,2 ..., Magent
saFor the fitness value of a-th of intelligent body;
The strength size of a-th of intelligent body is defined as:
The average individual number that each intelligent body region is assigned to is:
M.Sa=round { Pa×Mpublic}
M.SaThe average individual number possessed by a-th of intelligent body region, so that the intelligent body after initializing
Strength is stronger, and the number for the average individual that region is assigned to is more;
(2) local search, is specifically realized in two steps:
(2.1) polymerization movement
Intelligent body of the average individual gradually to one's respective area in each region is close, and displacement distance Move is obeyed and is uniformly distributed,
Represent as follows:
Move~U (0, β × d)
Wherein, U is to be uniformly distributed symbol;β is polymerizing factor, and β takes 2;D for intelligent body in the same area and average individual it
Between distance;
(2.2) movement is changed
In order to which the speed of the average individual movement in polymerization process is carried out slows down, increase population diversity, to each common
Individual is changed at random, and the average individual number that a-th of intelligent body region is possessed is Ma.public, wherein needing to carry out
The average individual number M of changea,rpublicFor:
Ma.rpublic=round (qa×Ma.public)
Wherein, round functions is round up function, qaFor change rate, 0.3 is taken;
Inside each region after polymerization and changing movement, each individual strength value meeting great changes have taken place, it is necessary to according to
The ideal adaptation angle value rearrangement of each region, the ideal adaptation angle value to make number one is best, is new intelligent body;
(3) global search, including:
(3.1) vie each other
Each iteration is by local search afterwards, it is necessary to calculate total strength value of each intelligent body.Intelligent body in each region
Total strength value by itself strength value and together decide on the strength value of the average individual in region, be calculated by following formula:
T.SaFor total strength of a-th of intelligent body, ξ is a positive number less than 1, represents that average individual is real in the same area
Power accounts for the weight of the total strength of intelligent body, takes 0.1, wnFor the strength value of average individual;
A-th of intelligent body competition probability is expressed as:
Wherein, N.T.SaRepresent relatively total strength value of a-th of intelligent body, be defined as:
M.T.Sa=max { T.Sb}-T.Sa, b=1,2 ... Magent
Thus the competition probability of each intelligent body is calculated, if vector p is:
Introduce and vector p is with the random vector R of dimension, be expressed as:
Wherein, r~U (0,1) represents being uniformly distributed for the element obedience 0 to 1 in R;
Definition vector V is the difference of vector p and vector R:
V=p-R
The corresponding intelligent body of maximum element finally obtains the average individual competed in vectorial V;
(3.2) cooperate
Increase cooperates operation after intelligent body is vied each other, when the distance between the intelligent body in two regions is less than
During cooperation distance D, all average individuals in two intelligent bodies in the small intelligent body region of strength value return strength value big
Intelligent body region own, i.e. two intelligent bodies merge to increase strength value, so as to increase itself competitiveness;Intelligent body xc
With xdBetween cooperation distance D be defined as:
D=norm (xc-xd)×u
Wherein, c=1,2 ... Magent, d=1,2 ... Magent, for norm functions to seek Norm function, u represents cooperation coefficient,
Value is 0~1;
(4) algorithm terminates
When running to there are not having average individual in some intelligent body region, which is eliminated, this
Sample intelligent body number gradually decreases, when iterating to maximum iteration, end of run.
A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm of the present invention, can effectively solve more
Target, the nonlinear optimization objective function of multiple constraint, improve microgrid economic benefit and environmental benefit.Specifically have the following advantages that:
1st, the alternating current-direct current mixing microgrid mathematical model established of the present invention, it is contemplated that power grid purchases strategies, micro- source fuel cost,
Environmental benefit cost, network loss and operation expense, object function obey microgrid internal power balance, points of common connection transmission
Before capacity, controllable micro- source climbing rate, unit interval accumulator cell charging and discharging bound, storage battery charge state bound, storage battery
Period power-balance and the constant constraints of storage battery charge state afterwards, are applied in Practical Project easy to method;
The 2nd, object function that cultural gene algorithm is applied to microgrid cost optimization solves, and can effectively solve multidimensional, more
Target, multiple constraint and nonlinear object function, a kind of new method and new approaches are provided for the optimization operation of alternating current-direct current mixing microgrid;
3rd, using the global search strategy of competition-collaboration mode, it is contemplated that global diversity, improves algorithm low optimization accuracy.
Brief description of the drawings
Fig. 1 is a kind of flow chart of the alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm of the present invention;
Fig. 2 is the structure diagram of alternating current-direct current mixing microgrid typical case's rack in the present invention;
Fig. 3 is example load prediction curve figure of the present invention;
Fig. 4 is each micro- source output situation map in exchanging area after present invention optimization;
Fig. 5 is exchanging area operating cost result figure after present invention optimization;
Fig. 6 is each micro- source output situation map in DC area after present invention optimization;
Fig. 7 is DC area operating cost result figure after present invention optimization;
Fig. 8 is storage battery charge state curve map after present invention optimization;
Fig. 9 a are each power supply output situation maps in microgrid exchanging area one day after present invention optimization;
Fig. 9 b are each power supply output situation maps in microgrid DC area one day after present invention optimization;
Fig. 9 c are battery discharging situation maps in microgrid DC area after present invention optimization;
Fig. 9 d are microgrid DC area storage battery charge condition figures after present invention optimization;
Figure 10 is points of common connection flowing power diagram after present invention optimization.
Embodiment
It is excellent to a kind of alternating current-direct current mixing microgrid based on cultural gene algorithm of the present invention with reference to embodiment and attached drawing
Change operation method to be described in detail.
A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm of the present invention, including, including such as
Lower step:
1) it is straight for the friendship comprising wind-power electricity generation, photovoltaic generation, storage battery, miniature combustion engine, fuel cell and diesel-driven generator
Stream mixing microgrid, establish multiple target, multiple constraint, it is nonlinear optimization operation mathematical model, object function consider power grid power purchase into
Sheet, micro- source fuel cost, environmental benefit cost, network loss and operation expense, and obey microgrid internal power and balance, is public
Tie point transmission capacity, controllable micro- source climbing rate, unit interval accumulator cell charging and discharging bound, above and below storage battery charge state
Limit, the constraints that storage battery surrounding time section power-balance and storage battery charge state are constant;Wherein,
The object function F is as follows:
F=FGrid+Fec+Floss+Fom
FGrid=FACGrid+FDCGrid
Fec=FACec+FDCec
Floss=FACloss+FDCloss
Fom=FACom+FDCom
Wherein, FGridFor power grid purchases strategies;FecFor financial cost;FlossFor network loss;FomFor equipment operation maintenance cost;
FACGridFor exchanging area power grid purchases strategies;FACecFor exchanging area financial cost, by the micro- source fuel cost F in exchanging areaACfuelAnd ring
Border benefit-cost FACenAddition obtains;FAClossFor exchanging area network loss;FAComFor exchanging area equipment operation maintenance cost;FDCGridFor
Exchanging area power grid purchases strategies;FDCecFor the micro- source financial cost in exchanging area, by the micro- source fuel cost F in DC areaDCfuelAnd Environmental Effect
Beneficial cost FDCenAddition obtains;FDClossFor exchanging area network loss;FDComFor exchanging area equipment operation maintenance cost;Wherein,
(1) exchanging area power grid purchases strategies FACGridWith DC area power grid purchases strategies FDCGridRespectively:
Wherein, △ TACPeriod for from exchanging area to power grid power purchase;T is total period in one day;T is day part in one day;It is exchanging area to power grid power purchase electricity;△TDCPeriod for from DC area to power grid power purchase;It is DC area to electricity
Net purchase power consumption;For power grid sale of electricity price in the corresponding period;
(2) exchanging area financial cost FACecWith DC area financial cost FDCec:
FACec=FACfuel+FACen
FDCec=FDCfuel+FDCen
Exchanging area includes wind-power electricity generation (WT), miniature combustion engine (MT) and diesel-driven generator (DEG), and DC area includes photovoltaic generation
(PV), fuel cell (FC) and storage battery (ES), exchanging area fuel cost FACfuelWith DC area fuel cost FDCfuelCan be by
The following formula is calculated:
FACfuel=FMTfuel+FDEGfuel
FDCfuel=FFCfuel
ηFC=-0.0023 × PFC+0.6735
Wherein, FMTfuelFor miniature combustion engine fuel cost;CMTFor miniature combustion engine cooler fuel price;LHV is the low heat value of fuel gas;
PMTFor miniature combustion engine output power;ηMTFor the generating efficiency of miniature combustion engine;FFCfuelFor the fuel cost of fuel cell;CFCFor fuel electricity
The cooler fuel price in pond;PFCFor the output power of fuel cell;ηFCFor the generating efficiency of fuel cell;FDEGfuelFor diesel-driven generator
Fuel cost;PDEGFor the output power of diesel-driven generator;A, b, c are diesel-driven generator power generation coefficient, by diesel engine factory
Family provides;
Exchanging area environmental benefit cost FACenWith DC area environmental benefit cost FDCenIt is calculated by the following formula:
Wherein, n1 is the micro- source number in exchanging area;M is the type of pollutant;αjFor the conversion cost of corresponding pollutant, member/
kg;EFi,jThe unit discharge of the jth kind pollutant produced for i-th of micro- source, kg/kW;PiFor the output work in i-th of micro- source
Rate;n2For the micro- source number in DC area.
(3) network loss F is exchangedAClossWith direct current network loss FDCloss:
Wherein, L1 is exchanging area branch sum;Pk、QkActive power, reactive power for branch k transmission;L2 is DC area
Branch sum;RkFor the resistance of branch k;UkFor the voltage effective value of branch k;
(4) ac operation maintenance cost FAComWith DC operation maintenance cost FDCom:
Wherein, β i are the operation expense coefficient in i-th of micro- source.
The method of weighting is taken, finally obtaining object function is:
MinF=(1- λ) × (FGrid+Fec+Fom)+λFloss
Wherein, λ is Web-based exercise coefficient, and the Web-based exercise coefficient lambda is 0.1.
The constraints includes:
(1) microgrid internal power equilibrium constraint
(2) points of common connection transmission capacity constraints
(3) controllable micro- source climbing rate constraints
(4) unit interval accumulator cell charging and discharging bound constraints
(5) storage battery charge state (State of Charge, SOC) bound constraints
Socmin≤Soct≤Socmax
(6) storage battery surrounding time section power-balance constraint condition
(7) storage battery charge state constraint independent of time condition
Socinitial=Socend
Wherein, N is micro- source number;Represent that t-th of period, i-th of micro- source is contributed;Represent that t-th of period is purchased from power grid
Electrical power;For t-th of period battery power variable quantity, discharge just, to be charged as bearing;Represent respectively t-th
Period AC and DC area payload;For the net flow power of t-th of period points of common connection;PG,maxPassed for points of common connection
Defeated capacity limit value;riFor the unit interval climbing rate in i-th of micro- source;For the absolute of t-th period battery power variable quantity
Value;For the limit value of t-th of period battery power variable quantity;SoctFor the state-of-charge of storage battery in t-th of period;
Socmin、SocmaxThe respectively upper lower limit value of state-of-charge;Uu represents charge and discharge electrostrictive coefficient, is 1 during charging, and when electric discharge is -1;η is
Accumulator cell charging and discharging efficiency, takes 95% here;QESFor battery rating;Socinitial、SocendFor the initial of state-of-charge
Value and end value.
2) cultural gene Algorithm for Solving alternating current-direct current mixing microgrid mathematical model is used;Including:
(1) initialize
Initial individuals are produced using random generating mode, the value range of each variable is represented with varmin, varmax, its
In any one variable xmInitial value drawn by following formula:
xm=varminm+rand(0,1)*(varmaxm-varminm)
If M is all individual numbers in colony, then m=1,2 ..., M, MagentFor the number of intelligent body, MpublicTo be common
The number of individual, that is, have:
M=Magent+Mpublic
Ascending order arrangement, preceding M are carried out according to fitness value size to initial individualsagentAs intelligent body, remaining is common
Individual, wherein intelligent body are that fitness value comes preceding M in colonyagentThe individual of position, remaining individual is average individual in colony,
Colony is divided into MagentRegion, average individual belong to the region where each intelligent body, and each intelligent body initial time is possessed
The number of average individual be to be determined by each intelligent body with respect to strength, a-th intelligent body (a=1,2 ..., Magent) it is relatively real
Power is specifically calculated by following formula:
Sa=max { sb}-sa, b=1,2 ..., Magent
saFor the fitness value of a-th of intelligent body;
The strength size of a-th of intelligent body is defined as:
The average individual number that each intelligent body region is assigned to is:
M.Sa=round { Pa×Mpublic}
M.SaThe average individual number possessed by a-th of intelligent body region, so that the intelligent body after initializing
Strength is stronger, and the number for the average individual that region is assigned to is more;
(2) local search, the local searching strategy used are for optimal Selection Strategy, i.e., real in each region in each iteration
Force value is best for intelligent body.Local search is specifically realized in two steps:
(2.1) polymerization movement
In the cultural gene algorithm based on competition-collaboration mode, the average individual in each region is gradually to one's respective area
Intelligent body is close, and displacement distance Move is obeyed and is uniformly distributed, and represents as follows:
Move~U (0, β × d)
Wherein, U is to be uniformly distributed symbol;β is polymerizing factor, by many experiments when β takes 2 algorithm optimizing performance compared with
It is good;D is the distance between intelligent body and average individual in the same area;
(2.2) movement is changed
After polymerization campaign is carried out, average individual can be drawn close rapidly cultural gene algorithm to intelligent body in each region, be carried
High convergence speed of the algorithm, but be likely to result in algorithm and be absorbed in local convergence, the solution finally obtained is not optimal solution.In order to
The speed that average individual moves in polymerization process is carried out slows down, and increases population diversity, each average individual is carried out random
Change, the average individual number that a-th of intelligent body region is possessed is Ma.public, wherein needing common changed
Body number Ma,rpublicFor:
Ma.rpublic=round (qa×Ma.public)
Wherein, round functions is round up function, qaFor change rate, q is drawn after many experimentsaAlgorithm is sought when taking 0.3
Excellent better performances;
Inside each region after polymerization and changing movement, each individual strength value meeting great changes have taken place, it is necessary to according to
The ideal adaptation angle value rearrangement of each region, the ideal adaptation angle value to make number one is best, is new intelligent body;
(3) global search, including:
(3.1) vie each other
Each iteration is by local search afterwards, it is necessary to calculate total strength value of each intelligent body.Intelligent body in each region
Total strength value by itself strength value and together decide on the strength value of the average individual in region, be calculated by following formula:
T.SaFor total strength of a-th of intelligent body, ξ is a positive number less than 1, represents that average individual is real in the same area
Power accounts for the weight of the total strength of intelligent body, takes 0.1, wnFor the strength value of average individual;
Competition is obtained between being embodied in intelligent body by strength value size in cultural gene algorithm based on competition-collaboration mode
Corresponding probability is obtained to compete the average individual that strength is most weak in all individuals, so as to strengthen the strength of itself.Strength is stronger
The probability that intelligent body competes to obtain the most weak average individual of global strength is bigger, and a-th of intelligent body competition probability is expressed as:
Wherein, N.T.SaRepresent relatively total strength value of a-th of intelligent body, be defined as:
M.T.Sa=max { T.Sb}-T.Sa, b=1,2 ... Magent
Thus the competition probability of each intelligent body is calculated, if vector p is:
Introduce and vector p is with the random vector R of dimension, be expressed as:
Wherein, r~U (0,1) represents being uniformly distributed for the element obedience 0 to 1 in R;
Definition vector V is the difference of vector p and vector R:
V=p-R
The corresponding intelligent body of maximum element finally obtains the average individual competed in vectorial V;
(3.2) cooperate
Iteration of the cultural gene algorithm Jing Guo former steps, disclosure satisfy that algorithm diversity, effectively prevent algorithmic statement to office
Portion's optimal solution.But experimental data shows, in the algorithm later stage, algorithm the convergence speed is slower, it is therefore desirable to increases acceleration mechanism, carries
High algorithm the convergence speed.Increase cooperates operation after intelligent body is vied each other, when between the intelligent body in two regions
When distance is less than cooperation distance D, all average individuals in two intelligent bodies in the small intelligent body region of strength value are returned
The big intelligent body region of strength value owns, i.e., two intelligent bodies merge to increase strength value, so as to increase itself competitiveness;
Intelligent body xcWith xdBetween cooperation distance D be defined as:
D=norm (xc-xd)×u
Wherein, c=1,2 ... Magent, d=1,2 ... Magent, for norm functions to seek Norm function, u represents cooperation coefficient,
Value is 0~1, and by many experiments, algorithmic statement performance is good when u takes 0.2;
(4) algorithm terminates
When running to there are not having average individual in some intelligent body region, which is eliminated, this
Sample intelligent body number gradually decreases, when iterating to maximum iteration, end of run.
3) correctness and validity of the alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm are verified.
Example is given below:
Consider alternating current-direct current mixing microgrid as shown in Figure 2, used for the flow controller between exchanging area and DC area
Bi-directional current controller, realizes the power circulation between exchanging area and DC area, so as to reduce microgrid operating cost, improves to micro-
The utilization rate in source.The micro- source in exchanging area includes miniature combustion engine, diesel-driven generator and wind-power electricity generation, the micro- source in DC area include fuel cell,
Photovoltaic generation and the storage battery 1 that capacity is 250kWh, each micro- source parameter are as shown in table 1.Storage battery charge state changes model
Enclose for 0.3~0.9, SocinitialAnd Socend0.3 is taken, in order to give full play to the effect of storage battery peak load shifting, according to Fig. 4 institutes
The load prediction curve figure shown, determines storage battery in peak times of power consumption 11:00-13:00 and 20:00-24:00, as long as charged shape
State meets constraints, and storage battery necessarily be in discharge condition, remaining period storage battery charges, when state-of-charge reaches
When 0.9, stop charging.Table 2 converts cost and emission factor for micro- source, and table 3 is power grid purchase electricity price at different moments.
The parameter in each micro- source of table 1
Convert cost and emission factor in 2 micro- source of table
3 tou power price of table
(1) load prediction data in alternating current-direct current mixing microgrid is collected, as shown in figure 3, to alternating current-direct current mixing microgrid future one
Each micro- source, which is contributed, when 24 is small in it optimizes, to consider power grid purchases strategies, micro- source fuel cost, environmental benefit cost, net
Damage and the minimum object function of operation expense, and obey microgrid internal power balance, points of common connection transmission capacity, can
Control micro- source climbing rate, unit interval accumulator cell charging and discharging bound, storage battery charge state bound, storage battery surrounding time section
Multiple constraints such as power-balance and storage battery charge state are constant;
(2) object function is solved with cultural gene algorithm, the alternating current-direct current mixing obtained as Figure 4-Figure 7 is micro-
Each micro- source output situation curve and cost curve of exchanging area and DC area are netted, wherein financial cost is micro- source operation expense
The sum of with fuel cost.As can be seen that since the micro- source fuel cell in DC area, storage battery and photovoltaic generation environmental benefit are preferable,
So environmental protection converts cost close to zero;
(3) implementation requirements storage battery charge state of the invention meets 0.3~0.9, and whole story state is 0.3 pact
Beam.Fig. 8 is storage battery charge state change curve after optimization, meets constraints;
(4) pass through cultural gene Algorithm for Solving optimization object function, finally obtain each micro- source in alternating current-direct current mixing microgrid
Output situation and points of common connection flow power diagram as shown in Fig. 9 a- Figure 10.Wherein PLA and PLD represents exchanging area and straight respectively
The load in area is flowed, GRID represents to represent that storage battery is charged and discharged electricity respectively from power grid power purchase electricity, ESC and ESD.
The example of the present invention uses cultural gene Algorithm for Solving, show that alternating current-direct current mixing microgrid total operating cost is in example
2608.0 yuan, wherein exchanging area operating cost is 1964.1 yuan, and DC area operating cost is 643.9 yuan, and each micro- source is contributed can
Meet various constraints, storage battery reaches full hair-like state when load is larger, since the micro- source cost of electricity-generating in exchanging area is more micro- than DC area
Source higher, so the flow of power there are DC area to exchanging area.
In conclusion test result through this embodiment, illustrates a kind of friendship based on cultural gene algorithm of the present invention
Direct current mixing microgrid optimizing operation method can effectively solve multiple target, the nonlinear optimization objective function of multiple constraint, improve microgrid
Economic benefit and environmental benefit, it was confirmed that the validity and correctness of cultural gene algorithm.
Claims (5)
1. a kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm, it is characterised in that including following step
Suddenly:
1) mixed for the alternating current-direct current comprising wind-power electricity generation, photovoltaic generation, storage battery, miniature combustion engine, fuel cell and diesel-driven generator
Microgrid is closed, establishes multiple target, multiple constraint, nonlinear optimization operation mathematical model, object function considers power grid purchases strategies, micro-
Source fuel cost, environmental benefit cost, network loss and operation expense, and obey microgrid internal power and balance, is commonly connected
Point transmission capacity, controllable micro- source climbing rate, unit interval accumulator cell charging and discharging bound, storage battery charge state bound, storage
Battery surrounding time section power-balance and the constant constraints of storage battery charge state;
2) cultural gene Algorithm for Solving alternating current-direct current mixing microgrid mathematical model is used;
3) correctness and validity of the alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm are verified.
2. a kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm according to claim 1, its
It is characterized in that, the object function F described in step 1) is as follows:
F=FGrid+Fec+Floss+Fom
FGrid=FACGrid+FDCGrid
Fec=FACec+FDCec
Floss=FACloss+FDCloss
Fom=FACom+FDCom
Wherein, FGridFor power grid purchases strategies;FecFor financial cost;FlossFor network loss;FomFor equipment operation maintenance cost;
FACGridFor exchanging area power grid purchases strategies;FACecFor exchanging area financial cost, by the micro- source fuel cost F in exchanging areaACfuelAnd ring
Border benefit-cost FACenAddition obtains;FAClossFor exchanging area network loss;FAComFor exchanging area equipment operation maintenance cost;FDCGridFor
Exchanging area power grid purchases strategies;FDCecFor the micro- source financial cost in exchanging area, by the micro- source fuel cost F in DC areaDCfuelAnd Environmental Effect
Beneficial cost FDCenAddition obtains;FDClossFor exchanging area network loss;FDComFor exchanging area equipment operation maintenance cost;Wherein,
(1) exchanging area power grid purchases strategies FACGridWith DC area power grid purchases strategies FDCGridRespectively:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>G</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>G</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>c</mi>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>G</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>G</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>c</mi>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
Wherein, △ TACPeriod for from exchanging area to power grid power purchase;T is total period in one day;T is day part in one day;
It is exchanging area to power grid power purchase electricity;△TDCPeriod for from DC area to power grid power purchase;It is DC area to power grid power purchase
Electricity;For power grid sale of electricity price in the corresponding period;
(2) exchanging area financial cost FACecWith DC area financial cost FDCec:
FACec=FACfuel+FACen
FDCec=FDCfuel+FDCen
Exchanging area includes wind-power electricity generation, miniature combustion engine and diesel-driven generator, and DC area includes photovoltaic generation, fuel cell and electric power storage
Pond, exchanging area fuel cost FACfuelWith DC area fuel cost FDCfuelIt can be calculated by the following formula:
FACfuel=FMTfuel+FDEGfuel
FDCfuel=FFCfuel
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mrow>
<mi>L</mi>
<mi>H</mi>
<mi>V</mi>
</mrow>
</mfrac>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>G</mi>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>+</mo>
<msub>
<mi>bP</mi>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>G</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>cP</mi>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>F</mi>
<mi>C</mi>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>F</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mi>L</mi>
<mi>H</mi>
<mi>V</mi>
</mrow>
</mfrac>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>F</mi>
<mi>C</mi>
</mrow>
</msub>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>F</mi>
<mi>C</mi>
</mrow>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.0753</mn>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mn>65</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>-</mo>
<mn>0.3095</mn>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mn>65</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>0.4174</mn>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mn>65</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>0.1068</mn>
</mrow>
ηFC=-0.0023 × PFC+0.6735
Wherein, FMTfuelFor miniature combustion engine fuel cost;CMTFor miniature combustion engine cooler fuel price;LHV is the low heat value of fuel gas;PMTFor
Miniature combustion engine output power;ηMTFor the generating efficiency of miniature combustion engine;FFCfuelFor the fuel cost of fuel cell;CFCFor fuel cell
Cooler fuel price;PFCFor the output power of fuel cell;ηFCFor the generating efficiency of fuel cell;FDEGfuelFor the combustion of diesel-driven generator
Expect cost;PDEGFor the output power of diesel-driven generator;A, b, c are diesel-driven generator power generation coefficient, and by diesel engine, manufacturer gives
Go out;
Exchanging area environmental benefit cost FACenWith DC area environmental benefit cost FDCenIt is calculated by the following formula:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>EF</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>EF</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
Wherein, n1 is the micro- source number in exchanging area;M is the type of pollutant;αjFor the conversion cost of corresponding pollutant, member/kg;
EFi,jThe unit discharge of the jth kind pollutant produced for i-th of micro- source, kg/kW;PiFor the output power in i-th of micro- source;n2
For the micro- source number in DC area.
(3) network loss F is exchangedAClossWith direct current network loss FDCloss:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mn>1</mn>
</mrow>
</munderover>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mrow>
<msubsup>
<mi>U</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mn>2</mn>
</mrow>
</munderover>
<mfrac>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>U</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
</mrow>
Wherein, L1 is exchanging area branch sum;Pk、QkActive power, reactive power for branch k transmission;L2 is DC area branch
Sum;RkFor the resistance of branch k;UkFor the voltage effective value of branch k;
(4) ac operation maintenance cost FAComWith DC operation maintenance cost FDCom:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>o</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>o</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</munderover>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
Wherein, β i are the operation expense coefficient in i-th of micro- source.
The method of weighting is taken, finally obtaining object function is:
MinF=(1- λ) × (FGrid+Fec+Fom)+λFloss
Wherein, λ is Web-based exercise coefficient.
3. a kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm according to claim 2, its
It is characterized in that, the Web-based exercise coefficient lambda is 0.1.
4. a kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm according to claim 1, its
It is characterized in that, the constraints described in step 1) includes:
(1) microgrid internal power equilibrium constraint
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mi>G</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mi>L</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mi>L</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
(2) points of common connection transmission capacity constraints
<mrow>
<mo>|</mo>
<msubsup>
<mi>P</mi>
<mi>G</mi>
<mi>t</mi>
</msubsup>
<mo>|</mo>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
(3) controllable micro- source climbing rate constraints
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
(4) unit interval accumulator cell charging and discharging bound constraints
<mrow>
<mo>|</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>|</mo>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
(5) storage battery charge state (State of Charge, SOC) bound constraints
Socmin≤Soct≤Socmax
(6) storage battery surrounding time section power-balance constraint condition
<mrow>
<msup>
<mi>Soc</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>Soc</mi>
<mi>t</mi>
</msup>
<mo>+</mo>
<mi>u</mi>
<mi>u</mi>
<mo>&times;</mo>
<mi>&eta;</mi>
<mo>&times;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>/</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
(7) storage battery charge state constraint independent of time condition
Socinitial=Socend
Wherein, N is micro- source number;Represent that t-th of period, i-th of micro- source is contributed;Represent t-th of period from power grid power purchase work(
Rate;For t-th of period battery power variable quantity, discharge just, to be charged as bearing;When representing respectively t-th
Section AC and DC area payload;For the net flow power of t-th of period points of common connection;PG,maxTransmitted for points of common connection
Capacity limit value;riFor the unit interval climbing rate in i-th of micro- source;For the absolute of t-th period battery power variable quantity
Value;For the limit value of t-th of period battery power variable quantity;SoctFor the state-of-charge of storage battery in t-th of period;
Socmin、SocmaxThe respectively upper lower limit value of state-of-charge;Uu represents charge and discharge electrostrictive coefficient, is 1 during charging, and when electric discharge is -1;η is
Accumulator cell charging and discharging efficiency, takes 95% here;QESFor battery rating;Socinitial、SocendFor the initial of state-of-charge
Value and end value.
5. a kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm according to claim 1, its
It is characterized in that, step 2) includes:
(1) initialize
Initial individuals are produced using random generating mode, the value range of each variable are represented with varmin, varmax, wherein appointing
Anticipate a variable xmInitial value drawn by following formula:
xm=varminm+rand(0,1)*(varmaxm-varminm)
If M is all individual numbers in colony, then m=1,2 ..., M, MagentFor the number of intelligent body, MpublicFor average individual
Number, that is, have:
M=Magent+Mpublic
Ascending order arrangement, preceding M are carried out according to fitness value size to initial individualsagentAs intelligent body, remaining is average individual,
Wherein intelligent body is that fitness value comes preceding M in colonyagentThe individual of position, remaining individual is average individual in colony, by group
Body is divided into MagentRegion, average individual belong to the region where each intelligent body, and each intelligent body initial time is possessed general
The number for leading to individual is to be determined by each intelligent body with respect to strength, a-th intelligent body (a=1,2 ..., Magent) opposite strength tool
Body is calculated by following formula:
Sa=max { sb}-sa, b=1,2 ..., Magent
saFor the fitness value of a-th of intelligent body;
The strength size of a-th of intelligent body is defined as:
<mrow>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
</msubsup>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
</mrow>
The average individual number that each intelligent body region is assigned to is:
M.Sa=round { Pa×Mpublic}
M.SaThe average individual number possessed by a-th of intelligent body region, so that the intelligent body strength after initializing is got over
By force, the number for the average individual that region is assigned to is more;
(2) local search, is specifically realized in two steps:
(2.1) polymerization movement
Intelligent body of the average individual gradually to one's respective area in each region is close, and displacement distance Move is obeyed and is uniformly distributed, and represents
It is as follows:
Move~U (0, β × d)
Wherein, U is to be uniformly distributed symbol;β is polymerizing factor, and β takes 2;D is in the same area between intelligent body and average individual
Distance;
(2.2) movement is changed
In order to which the speed of the average individual movement in polymerization process is carried out slows down, increase population diversity, to each average individual
Changed at random, the average individual number that a-th of intelligent body region is possessed is Ma.public, wherein needing to be changed
Average individual number Ma,rpublicFor:
Ma.rpublic=round (qa×Ma.public)
Wherein, round functions is round up function, qaFor change rate, 0.3 is taken;
Inside each region after polymerizeing and changing movement, great changes have taken place, it is necessary to according to each for each individual strength value meeting
Region ideal adaptation angle value rearrangement, the ideal adaptation angle value to make number one is best, is new intelligent body;
(3) global search, including:
(3.1) vie each other
Each iteration is by local search afterwards, it is necessary to calculate total strength value of each intelligent body.Intelligent body is always real in each region
Force value by itself strength value and together decide on the strength value of the average individual in region, be calculated by following formula:
<mrow>
<mi>T</mi>
<mo>.</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msub>
<mi>s</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mi>&xi;</mi>
<mo>&times;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>.</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
</mrow>
</msubsup>
<msub>
<mi>w</mi>
<mi>b</mi>
</msub>
</mrow>
<mrow>
<mi>M</mi>
<mo>.</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
</mrow>
T.SaFor total strength of a-th of intelligent body, ξ is a positive number less than 1, represents that average individual strength accounts in the same area
The weight of the total strength of intelligent body, takes 0.1, wnFor the strength value of average individual;
A-th of intelligent body competition probability is expressed as:
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>P</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>|</mo>
<mfrac>
<mrow>
<mi>M</mi>
<mo>.</mo>
<mi>T</mi>
<mo>.</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
</msubsup>
<mi>M</mi>
<mo>.</mo>
<mi>T</mi>
<mo>.</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
Wherein, N.T.SaRepresent relatively total strength value of a-th of intelligent body, be defined as:
M.T.Sa=max { T.Sb}-T.Sa, b=1,2 ... Magent
Thus the competition probability of each intelligent body is calculated, if vector p is:
<mrow>
<mi>p</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>p</mi>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</msub>
</mtd>
<mtd>
<msub>
<mi>p</mi>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>p</mi>
<msub>
<mi>P</mi>
<msub>
<mi>M</mi>
<mrow>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
</msub>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
</mrow>
Introduce and vector p is with the random vector R of dimension, be expressed as:
<mrow>
<mi>R</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<msub>
<mi>M</mi>
<mrow>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>r</mi>
<mo>~</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, r~U (0,1) represents being uniformly distributed for the element obedience 0 to 1 in R;
Definition vector V is the difference of vector p and vector R:
V=p-R
The corresponding intelligent body of maximum element finally obtains the average individual competed in vectorial V;
(3.2) cooperate
Increase cooperates operation after intelligent body is vied each other, when the distance between the intelligent body in two regions is less than cooperation
During distance D, all average individuals in two intelligent bodies in the small intelligent body region of strength value return the big intelligence of strength value
Energy body region owns, i.e., two intelligent bodies merge to increase strength value, so as to increase itself competitiveness;Intelligent body xcWith xd
Between cooperation distance D be defined as:
D=norm (xc-xd)×u
Wherein, c=1,2 ... Magent, d=1,2 ... Magent, for norm functions to seek Norm function, u represents cooperation coefficient, value
For 0~1;
(4) algorithm terminates
When running to there are not having average individual in some intelligent body region, which is eliminated, such intelligence
Energy body number gradually decreases, when iterating to maximum iteration, end of run.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711359910.9A CN107947178B (en) | 2017-12-15 | 2017-12-15 | A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711359910.9A CN107947178B (en) | 2017-12-15 | 2017-12-15 | A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107947178A true CN107947178A (en) | 2018-04-20 |
CN107947178B CN107947178B (en) | 2019-03-05 |
Family
ID=61944545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711359910.9A Active CN107947178B (en) | 2017-12-15 | 2017-12-15 | A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107947178B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109768566A (en) * | 2018-11-30 | 2019-05-17 | 天合光能股份有限公司 | A kind of novel friendship directly-heated is electrically coupled mixing microgrid main circuit |
CN111293719A (en) * | 2020-02-29 | 2020-06-16 | 华北电力大学(保定) | Alternating current-direct current hybrid micro-grid optimized operation method based on multi-factor evolution algorithm |
CN111697635A (en) * | 2020-06-03 | 2020-09-22 | 华北电力大学(保定) | Alternating current-direct current hybrid micro-grid optimized operation method considering random fuzzy double uncertainty |
US11742667B2 (en) | 2018-05-03 | 2023-08-29 | Vestas Wind Systems A/S | Integrated hybrid power plants for off-grid systems |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166877A (en) * | 2014-05-31 | 2014-11-26 | 徐多 | Microgrid optimization operation method based on improved binary system particle swarm optimization algorithm |
CN105977991A (en) * | 2016-05-10 | 2016-09-28 | 浙江工业大学 | Independent micro grid optimization configuration method considering price-type demand response |
CN106100002A (en) * | 2016-07-28 | 2016-11-09 | 华北电力大学(保定) | A kind of optimizing operation method of alternating current-direct current mixing microgrid |
US9564757B2 (en) * | 2013-07-08 | 2017-02-07 | Eaton Corporation | Method and apparatus for optimizing a hybrid power system with respect to long-term characteristics by online optimization, and real-time forecasts, prediction or processing |
CN106611379A (en) * | 2016-02-22 | 2017-05-03 | 四川用联信息技术有限公司 | Improved culture gene algorithm for solving multi-objective flexible job shop scheduling problem |
-
2017
- 2017-12-15 CN CN201711359910.9A patent/CN107947178B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9564757B2 (en) * | 2013-07-08 | 2017-02-07 | Eaton Corporation | Method and apparatus for optimizing a hybrid power system with respect to long-term characteristics by online optimization, and real-time forecasts, prediction or processing |
CN104166877A (en) * | 2014-05-31 | 2014-11-26 | 徐多 | Microgrid optimization operation method based on improved binary system particle swarm optimization algorithm |
CN106611379A (en) * | 2016-02-22 | 2017-05-03 | 四川用联信息技术有限公司 | Improved culture gene algorithm for solving multi-objective flexible job shop scheduling problem |
CN105977991A (en) * | 2016-05-10 | 2016-09-28 | 浙江工业大学 | Independent micro grid optimization configuration method considering price-type demand response |
CN106100002A (en) * | 2016-07-28 | 2016-11-09 | 华北电力大学(保定) | A kind of optimizing operation method of alternating current-direct current mixing microgrid |
Non-Patent Citations (1)
Title |
---|
刘漫丹: "文化基因算法研究进展", 《自动化技术与应用》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11742667B2 (en) | 2018-05-03 | 2023-08-29 | Vestas Wind Systems A/S | Integrated hybrid power plants for off-grid systems |
CN109768566A (en) * | 2018-11-30 | 2019-05-17 | 天合光能股份有限公司 | A kind of novel friendship directly-heated is electrically coupled mixing microgrid main circuit |
CN111293719A (en) * | 2020-02-29 | 2020-06-16 | 华北电力大学(保定) | Alternating current-direct current hybrid micro-grid optimized operation method based on multi-factor evolution algorithm |
CN111293719B (en) * | 2020-02-29 | 2023-06-27 | 华北电力大学(保定) | AC/DC hybrid micro-grid optimized operation method based on multi-factor evolution algorithm |
CN111697635A (en) * | 2020-06-03 | 2020-09-22 | 华北电力大学(保定) | Alternating current-direct current hybrid micro-grid optimized operation method considering random fuzzy double uncertainty |
CN111697635B (en) * | 2020-06-03 | 2023-07-25 | 华北电力大学(保定) | AC/DC hybrid micro-grid optimization operation method considering random fuzzy double uncertainty |
Also Published As
Publication number | Publication date |
---|---|
CN107947178B (en) | 2019-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104701871B (en) | One kind is containing the honourable complementary microgrid hybrid energy-storing capacity optimum proportioning method of water multi-source | |
Mesbahi et al. | Combined optimal sizing and control of Li-ion battery/supercapacitor embedded power supply using hybrid particle Swarm–Nelder–Mead algorithm | |
CN107947178B (en) | A kind of alternating current-direct current mixing microgrid optimizing operation method based on cultural gene algorithm | |
Chen et al. | Optimal sizing and sitting of EVCS in the distribution system using metaheuristics: A case study | |
CN104362677B (en) | A kind of active distribution network distributes structure and its collocation method rationally | |
CN110070292B (en) | Micro-grid economic dispatching method based on cross variation whale optimization algorithm | |
CN106026152A (en) | Charging and discharging scheduling method for electric vehicles connected to micro-grid | |
CN108090277A (en) | A kind of electric vehicle microgrid dual-layer optimization dispatching method for considering satisfaction and dispatching | |
CN106549419B (en) | Independent microgrid system design method based on universal gravitation algorithm | |
CN106058855A (en) | Active power distribution network multi-target optimization scheduling method of coordinating stored energy and flexible load | |
CN108964098A (en) | A kind of independent micro-grid system capacity configuration optimizing method | |
CN105226688A (en) | Based on the polymorphic type energy storage system capacity configuration optimizing method of Chance-constrained Model | |
CN106712114A (en) | Layered energy optimization scheduling method of active power distribution network considering environmental benefits | |
CN110084443A (en) | A kind of electrical changing station optimal operation model analysis method based on QPSO optimization algorithm | |
CN110363363A (en) | A kind of dual-stage scheduling model optimization method for integrated energy system | |
CN106026169A (en) | Decomposition-coordination optimization method based on multi-microgrid merging into power distribution network | |
CN107546773A (en) | A kind of more micro-capacitance sensor dynamic networking methods in region based on graph theory | |
CN110071496A (en) | A kind of configuration of direct-current grid power optimization and operation method based on wave-activated power generation | |
CN102856899A (en) | Method of reducing network loss of micro power grid | |
CN107482660A (en) | A kind of active distribution network energy storage configuration method based on double-energy storage system | |
CN107273968A (en) | A kind of Multiobjective Scheduling method and device based on dynamic fuzzy Chaos-Particle Swarm Optimization | |
CN107732945A (en) | A kind of energy-storage units optimization method based on simulated annealing particle cluster algorithm | |
CN108039741A (en) | The alternating current-direct current mixing microgrid optimizing operation method of electricity online more than meter and micro- source | |
CN108512259A (en) | A kind of alternating current-direct current mixing microgrid dual blank-holder based on Demand Side Response | |
CN109494813A (en) | A kind of power dispatching method, electronic equipment and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |