CN106651628A - Regional cool and thermal power comprehensive energy optimizing configuration method and apparatus based on graph theory - Google Patents

Abstract

The invention relates to a regional cool and thermal power comprehensive energy optimizing configuration method and apparatus based on a graph theory. The technical characteristics of the invention are that the method includes the following steps: step 1, constructing an optimizing model of the lowest regional cool and thermal power comprehensive energy cost according to a target function of the lowest regional cool and thermal power comprehensive energy cost and constraint conditions set in advance; step 2, on the basis of the optimizing model of the lowest regional cool and thermal power comprehensive energy cost established in the step 1, establishing a regional cool and thermal power comprehensive energy optimizing configuration model based on the graph theory through a network topology connection of a cool and thermal power comprehensive energy system; and step 3, solving the regional cool and thermal power comprehensive energy optimizing configuration model based on the graph theory through an improved genetic algorithm to realize cool and thermal power comprehensive energy optimizing configuration. By comprehensively taking the mutual coupling characteristics of cool and thermal power hybrid energy into consideration and also taking output constraints of the cool and thermal power hybrid energy into account, the method constructs the optimizing model at the lowest regional cool and thermal power comprehensive energy cost, and provides a solution for optimizing configuration of cool and thermal power hybrid energy within regions.

Description

Region cool and thermal power comprehensive energy Optimal Configuration Method and device based on graph theory

Technical field

The invention belongs to cool and thermal power comprehensive energy optimisation technique field, especially a kind of region cool and thermal power based on graph theory is comprehensive Close energy source optimization collocation method and device.

Background technology

With the development of economic society, hot and cold, the electric equal energy source consumption amplification of user increases year by year, and total energy consumption is annual Increase, social total energy consumption has remained high；Under the call that national energy-saving reduces discharging with new forms of energy efficient utilization, user is in the energy The Optimum utilization of cool and thermal power energy mix is carried out in selection, also become user's energy-conservation, saved the important component part of expense, be also to be lifted Efficiency of energy utilization, effectively reduces a kind of effective way of social total energy consumption；But currently for the optimization of cool and thermal power energy mix Configuration lacks effective Optimal Configuration Method, lacks feasible technological means.Traditional distribution system Optimal Configuration Method, typically All be from the economic optimum of bulk power grid/power distribution network, consider one-sided electric energy etc. aspect, propose to carry out effective cold, heat and electricity three-way For the building level addressing constant volume method of unit.

In brief, traditional energy Optimal Configuration Method distributing rationally only for electric power energy, only unilaterally examines Consider the object function and constraints of electric power, do not account for the various energy resources such as hot and cold, electric and distribute rationally.

Cool and thermal power energy mix is comprehensively utilized, and is the future trend of intelligent grid/energy internet.For electric power energy Distribute rationally, some scholars have carried out research：As document (Wang Jiangjiang. building level cooling heating and power generation system optimize and many attributes it is comprehensive Close study on evaluation way [D]. North China Electric Power University, 2010, thesis for the doctorate) from the cooling and heating load of single building, by cold Thermoelectricity trilogy supply unit realizes the comprehensive utilization of the energy；Because the starting point of document research is only led from electric power energy as it Research object is wanted, when the multiple building in actual area or other energy demands are varied widely, it tends to be difficult to realize Its target.Also other documents are for the planning and designing of cold, heat and electricity triple supply unit, optimal control and wait many-side to carry out depth The research for entering, but still with certain limitation：Be first in such document only study trilogy supply unit single element or Person's system, and in actual application, the mode variation of central cooling, such as ice cold-storage, water cold storage；The side of central heating Formula also has various, for example city planting ductwork, steam power plant etc.；And these modes are interlaced with each other and merge into each other, if only studying three Alliance unit single element or system, then with certain limitation, also lack practical application meaning.

The content of the invention

It is an object of the invention to overcome the deficiencies in the prior art, there is provided one kind can consider cool and thermal power energy mix Be mutually coupled feature and take into account the units limits of cool and thermal power energy mix the region cool and thermal power comprehensive energy based on graph theory it is excellent Change collocation method and device.

The present invention solves its technical problem and takes technical scheme below to realize：

A kind of region cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory, comprises the following steps：

Step 1, the object function according to region cool and thermal power comprehensive energy operating cost minimum and the constraint bar for pre-setting Part, builds the minimum Optimized model of region cool and thermal power comprehensive energy expense；

The Optimized model of step 2, the region cool and thermal power comprehensive energy expense minimum built based on step 1, with reference to cool and thermal power The network topology connection of integrated energy system, sets up based on the region cool and thermal power comprehensive energy Optimal Allocation Model of graph theory；

Step 3, the region cool and thermal power comprehensive energy Optimal Allocation Model using improved adaptive GA-IAGA solution based on graph theory, it is real Show distributing rationally for cool and thermal power comprehensive energy.

And, the object function of the Optimized model of the region cool and thermal power comprehensive energy expense minimum of the step 1 is:

In above formula, E_iIt is i-th user power utilization expense minimum of a value of region cool and thermal power comprehensive energy expense；Q_iIt is that region is cold I-th user's heat energy expense minimum of a value of thermoelectricity comprehensive energy expense；W_iIt is i-th use of region cool and thermal power comprehensive energy expense Family cold energy expense minimum of a value；N is the total quantity of region cool and thermal power energy mix user；

Its constraints includes：

(1) region cool and thermal power comprehensive energy trend constraint

In above formula, N_sFor total node number set；G_ij、Bi_jFor the admittance coefficient between node i, j；V_i, V_jFor node i, j's Voltage magnitude；P_Gi、Q_GiActive the exerting oneself of the respectively generator of node i is exerted oneself with idle；P_Di、Q_DiRespectively node i is active And load or burden without work；θ_ijIt is the generator rotor angle between node i and node j.

(2) region cool and thermal power comprehensive energy node power output constraint

P_imin≤P_i≤P_imax

In above formula, P_iminAnd P_imaxIt is respectively cool and thermal power comprehensive energy node power output lower limit and the upper limit；

(3) cool and thermal power comprehensive energy in region intercouples constraint

P_inmax≤P_c≤P_outmax

In above formula, P_inmaxAnd P_outmaxIt is respectively the lower limit and the upper limit to bulk power grid sale of electricity power；

And, the concrete grammar of the step 2 is：Consider the network topology connection of actual cool and thermal power comprehensive energy, will Region cool and thermal power comprehensive energy is distributed rationally and is equivalent to a undirected authorized graph G=(V, E, W), finds undirected authorized graph G most Little weights W, so as to distribute cool and thermal power comprehensive energy rationally；

Wherein, the element in set V is the fixed point or node or point of undirected authorized graph G, represents reality in regional extent One cool and thermal power power supply, the finite non-NULL node set that it is made up of region cool and thermal power power supply, according to the number of cool and thermal power power supply Amount order starts number consecutively from 1, until all of cool and thermal power power supply is numbered and finished；

Wherein, the element of set E is side or the line of undirected authorized graph G, represent cool and thermal power power supply between interconnection switch or Person's PCC switch sets, can use e_ijRepresent, e_ijValue be 1 and 0；Wherein, 1 represent cool and thermal power power supply i and cool and thermal power power supply j it Between exist contact, 0 represent without contact；E represents side collection in V；

Wherein, the element of set W is the active power cross-over value between any two node, the power of referred to as undirected authorized graph G Value, W_ijNode i is represented to the active power value exchanged between node j, value is positive number when active power is flowed into, otherwise when stream Value is negative when going out active power.

And, the concrete steps of the step 3 include：

(1) according to the region cool and thermal power energy mix Optimal Allocation Model based on graph theory constructed by step 2, engineering is taken into account Real data, the force data that goes out of the cool and thermal power network parameter in region and each cool and thermal power power supply is initialized.

(2) binary coding scheme is adopted, the fortune of grid-connected PCC switches, single cool and thermal power power supply to each cool and thermal power power supply Three coding decision variables of row expense and capacity are encoded；

(3) whether to reach the iterations of maximum as the foundation for terminating, judge whether to meet end condition；If reaching It is i.e. out of service, obtain final result；Otherwise enter (4th) step；

(4) Population Size is set, genetic algorithm inherent parameters is determined by the operation method for selecting, intersect and make a variation；

Wherein, adaptive crossover operator function is as follows：

In above formula, f_avgIt is the average fitness of per generation colony；f_maxFor fitness maximum in the individuality to be intersected；F is Larger fitness in two individualities to be intersected；P_c1Value 0.9, P_c2Value 0.6；

Wherein, TSP question rate function is as follows：

P_m=P_m1-P_m1×i/N

In above formula, P_m1For the initial value of aberration rate, value 0.08；I is current iteration number of times；N is iteration total degree；

(5) the cool and thermal power configuration of respective nodes is changed, and it is minimum according to cool and thermal power comprehensive energy expense in region described in step 1 Bound for objective function, the configuration to cool and thermal power power supply judges one by one according to constraints, if being satisfied by condition, Into (6th) step；Otherwise, iterations N adds 1, and returns (4th) step；

(6) optimization aim and its constraint bar of the comprehensive energy expense minimum of cool and thermal power energy mix in region are considered Part, builds fitness function as follows；

In above formula, C_maxIt is a given constant, f (x) is the object function after being normalized；

(7) optimal value that fitness function is produced is replaced, and iterations N is added 1, return (4th) step, carry out maximum Iterations terminates judging.

And, the concrete grammar of (2nd) step of the step 3 is：First, the chromosome that coding total length is 13 is chosen, Described three coding decision variables are each constituted into chromosome word string, then above-mentioned chromosome word string is linked to be into a complete dye Colour solid；Front two represents that the grid-connected PCC of cool and thermal power power supply is switched, middle five operating costs for representing single cool and thermal power power supply, Last six capacity for representing single cool and thermal power power supply.

Device is distributed rationally based on the region cool and thermal power comprehensive energy of graph theory, including：

Optimized model builds module, for according to the minimum object function of region cool and thermal power comprehensive energy operating cost and in advance The constraints for first arranging, builds the minimum Optimized model of region cool and thermal power comprehensive energy expense；

Optimal Allocation Model builds module, for building the region cool and thermal power synthesis of module construction based on the optimization module The minimum Optimized model of energy expenditure, connects with reference to the network topology of cool and thermal power integrated energy system, sets up based on the area of graph theory Domain cool and thermal power comprehensive energy Optimal Allocation Model；

Module is distributed rationally, for solving based on the region cool and thermal power comprehensive energy optimization of graph theory using improved adaptive GA-IAGA Allocation models, realizes distributing rationally for cool and thermal power comprehensive energy.

Advantages of the present invention and good effect are：

1st, the present invention proposes a kind of minimum Optimized model of cool and thermal power comprehensive energy expense, and introduces a kind of think of of graph theory Want to process cool and thermal power energy mix Optimal Allocation Model, and a kind of cool and thermal power of graph theory is solved using improved adaptive GA-IAGA Energy mix Optimal Allocation Model.Compared with existing one-sided electric power energy is distributed rationally, institute's extracting method of the present invention considers cold The units limits for being mutually coupled feature, taking into account cool and thermal power energy mix of thermoelectricity energy mix, construct comprehensive with region cool and thermal power The minimum Optimized model of energy expenditure is closed, effectively cool and thermal power energy mix in regional extent is distributed rationally there is provided solution party Case.

2nd, the present invention from a kind of Optimized model for proposing that cool and thermal power comprehensive energy expense is minimum on the whole of regional extent from And the starting point of part research is made up only from electric power energy as its main study subject, when the multiple building in actual area Or other energy demands are when varying widely, it tends to be difficult to realize the defect of its target.

Description of the drawings

Fig. 1 is new Tianjin ecological city in the region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory of the invention The schematic diagram of real system；

Fig. 2 is the improved adaptive GA-IAGA stream of the region cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory of the present invention Cheng Tu.

Specific embodiment

The embodiment of the present invention is described in further detail below in conjunction with accompanying drawing：

What the present invention considered cool and thermal power energy mix is mutually coupled feature, takes into account exerting oneself about for cool and thermal power energy mix Beam, constructs with the minimum Optimized model of region cool and thermal power comprehensive energy expense, it is proposed that a kind of cool and thermal power based on graph theory is mixed Energy source optimization collocation method is closed, effectively cool and thermal power energy mix in regional extent is distributed rationally there is provided solution.

A kind of region cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory, as shown in Fig. 2 comprising the following steps：

Step 1, the object function according to region cool and thermal power comprehensive energy operating cost minimum and the constraint bar for pre-setting Part, builds the minimum Optimized model of region cool and thermal power comprehensive energy expense；

The object function of the minimum Optimized model of the region cool and thermal power comprehensive energy expense is:

In above formula, E_iIt is i-th user power utilization expense minimum of a value of region cool and thermal power comprehensive energy expense；Q_iIt is that region is cold I-th user's heat energy expense minimum of a value of thermoelectricity comprehensive energy expense；W_iIt is i-th use of region cool and thermal power comprehensive energy expense Family cold energy expense minimum of a value；N is the total quantity of region cool and thermal power energy mix user.

Its constraints includes：

(1) region cool and thermal power comprehensive energy trend constraint

Cool and thermal power trend constraint is the active power of each node and reactive power equilibrium constraint, i.e. the trend constraint side of system Cheng Wei：

In above formula, N_sFor total node number set；G_ij、B_ijFor the admittance coefficient between node i, j；V_i, V_jFor node i, j's Voltage magnitude；P_Gi、Q_GiActive the exerting oneself of the respectively generator of node i is exerted oneself with idle；P_Di、Q_DiRespectively node i is active And load or burden without work；θ_ijIt is the generator rotor angle between node i and node j.

(2) region cool and thermal power comprehensive energy node power output constraint

P_imin≤P_i≤P_imax

In above formula, P_iminAnd P_imaxIt is respectively cool and thermal power comprehensive energy node power output lower limit and the upper limit；

Its power output should be between its peak power output and minimum output power.

(3) cool and thermal power comprehensive energy in region intercouples constraint

The constraint for intercoupling refers to that the maximum capacity interacted at cool and thermal power points of common connection must is fulfilled for the thing of connecting line The supply and demand agreement that transmission capacity is limited or they are reached is managed, its capacity-constrained is：

P_inmax≤P_c≤P_outmax

In above formula, P_inmaxAnd P_outmaxIt is respectively the lower limit and the upper limit to bulk power grid sale of electricity power；

The Optimized model of step 2, the region cool and thermal power comprehensive energy expense minimum built based on step 1, with reference to cool and thermal power The network topology connection of integrated energy system, sets up based on the region cool and thermal power comprehensive energy Optimal Allocation Model of graph theory；

The network topology connection of actual cool and thermal power comprehensive energy is considered, by any one cold and hot power supply according to active power A node is equivalent to, according to practical topology connection, will be connected and had between the node that the energy is exchanged with direct physical topology Carry out corresponding line, it is as shown in Figure 1 in new Tianjin ecological city cool and thermal power comprehensive energy topological diagram, by region cool and thermal power synthesis Energy source optimization configuration is equivalent to a undirected authorized graph G=(V, E, W), and the region cool and thermal power comprehensive energy set up based on graph theory is excellent Change allocation models.

The concrete grammar of the step 2 is：Region cool and thermal power energy mix is distributed rationally and is equivalent to one as shown in Figure 1 Undirected authorized graph G=(V, E, W), according to more than thinking, distributed rationally based on the region cool and thermal power comprehensive energy of graph theory Process, is exactly the process of the W weights minimum for finding figure G.

Therefore, the minimum weights W of undirected authorized graph G is found, so as to distribute cool and thermal power comprehensive energy rationally.

Wherein, the element in set V is the fixed point or node or point of undirected authorized graph G, represents reality in regional extent One cool and thermal power power supply, the finite non-NULL node set that it is made up of region cool and thermal power power supply, according to the number of cool and thermal power power supply Amount order starts number consecutively from 1, until all of cool and thermal power power supply is numbered and finished；

Wherein, the element of set E is side or the line of undirected authorized graph G, represent cool and thermal power power supply between interconnection switch or Person's PCC switch sets, can use e_ijRepresent, e_ijValue be 1 and 0；Wherein, 1 represent cool and thermal power power supply i and cool and thermal power power supply j it Between exist contact, 0 represent without contact；E represents side collection in V；

Wherein, the element of set W is the active power cross-over value between any two node, referred to as schemes the power of G, W_ijRepresent section The active power value exchanged between point i to node j, when active power is flowed into, value is positive number, otherwise when outflow active power When value be negative.

In the range of the cold and hot electric energy energy supply in region, cool and thermal power power supply on the one hand itself can independent islet operation, itself The optimum management of energy is carried out, on the other hand due to factor restrictions such as technology, cost, engineering reality, institute in region is extremely difficult to Having has actual physical connection switch PCC between cool and thermal power power supply, and after PCC switch connections, if can be formed active The interaction of power, i.e. W_ijIt is not equal to zero, various conditions such as energy storage, load for also being connected by cool and thermal power power supply are limited.

Step 3, the region cool and thermal power comprehensive energy Optimal Allocation Model using improved adaptive GA-IAGA solution based on graph theory, it is real Show distributing rationally for cool and thermal power comprehensive energy.

Because the region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory is a nonlinear multi-constrained problem Solve problems, adoptable heuristic improved adaptive GA-IAGA is solved, so as to obtain optimal solution.

The concrete steps of the step 3 are as shown in Fig. 2 including step (1) as described below to step (7).

(1) according to the region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory constructed by step 2, engineering is taken into account Real data, the force data that goes out of the cool and thermal power network parameter in region and each cool and thermal power power supply is initialized.

(2) binary coding scheme is adopted, the fortune of grid-connected PCC switches, single cool and thermal power power supply to each cool and thermal power power supply Three coding decision variables of row expense and capacity are encoded；

Coding is a most important link in genetic manipulation, and the quality of coding method directly affects algorithm and calculates speed Degree, result of calculation correctness etc..In traditional binary coding scheme, each node with a binary representation, but this Kind of encoding scheme can increase that code length is long with node, so as to the speed for affecting to calculate.

The concrete grammar of (2nd) step of the step 3 is：First, the chromosome that coding total length is 13 is chosen, will be described Three coding decision variables each constitute chromosome word string, then above-mentioned chromosome word string are linked to be into a complete chromosome； Front two represents that the grid-connected PCC of cool and thermal power power supply is switched, middle five operating costs for representing single cool and thermal power power supply, last six Position represents the capacity of single cool and thermal power power supply.

In the present embodiment, using binary coding scheme, the decision variable of coding is opened for the grid-connected PCC of cool and thermal power power supply Pass, the operating cost of single cool and thermal power power supply and capacity.Due to optimization problem in contain three decision variable (cool and thermal power power supplys Grid-connected PCC switches, the operating cost of single cool and thermal power power supply and capacity) participate in optimization, belong to multi-parameter coding.Using think of Think it is each decision variable first respective composition chromosome word string, then again these substrings are linked to be into a complete chromosome.Choosing dye Colour solid coding total length is 13.Front two represents the grid-connected PCC switches of cool and thermal power power supply, and middle five represent single cool and thermal power electricity The operating cost in source, last six capacity for representing single cool and thermal power power supply.This kind of coded system code length is short, and encodes Length will not become long with increasing for node, and calculating speed is fast, and convergence is good.

(3) whether to reach the iterations of maximum as the foundation for terminating, judge whether to meet end condition；If reaching It is i.e. out of service, obtain final result；Otherwise enter (4th) step；

(4) Population Size n is set, genetic algorithm inherent parameters is determined by the operation method for selecting, intersect and make a variation；

The selection of Population Size is critically important in genetic algorithm calculating.If it is little that Population Size n is selected, hereditary calculation can be improved The arithmetic speed of method, but but reduce the diversity of population, it is possible to cause the precocious phenomenon of genetic algorithm；And n it is too big when The efficiency of algorithm can be caused again to be reduced.Generally, n generally takes 20 to 100.

Selecting operation：Selection operation is set up on the basis evaluated individual fitness, and fitness is higher The probability that individuality is genetic to population of future generation is larger, otherwise then less.

Crossing-over rate computing：Crossover operation is used to produce new individuality, and it plays the role of a nucleus in genetic manipulation.It is typically based on The individual different fitness to be evolved select different crossover probabilities.Adaptive crossover operator function is as follows：

In above formula, f_avgIt is the average fitness of per generation colony；f_maxFor fitness maximum in the individuality to be intersected；F is Larger fitness in two individualities to be intersected；P_c1Value 0.9, P_c2Value 0.6；

Aberration rate computing：Variation belongs to complementary operation herein using TSP question rate in genetic manipulation, is evolving Starting stage, adopt higher aberration rate, with algebraically increase, aberration rate is steadily decreasing, to ensure the diversity of population And Optimality.TSP question rate function is as follows：

P_m=P_m1-P_m1×i/N

In above formula, P_m1For the initial value of aberration rate, 0.08 can be set to；I is current iteration number of times；N is iteration total degree.

(5) the cool and thermal power configuration of respective nodes is changed, and it is minimum according to cool and thermal power comprehensive energy expense in region described in step 1 Bound for objective function, the configuration to cool and thermal power power supply judges one by one according to constraints, if being satisfied by condition, Into (6th) step；Otherwise, iterations N adds 1, and returns (4th) step；

(6) optimization aim and its constraint bar of the comprehensive energy expense minimum of cool and thermal power energy mix in region are considered Part, builds fitness function as follows；

In above formula, C_maxIt is a given constant, f (x) is the object function after being normalized.

(7) optimal value that fitness function is produced is replaced, and iterations N is added 1, return (4th) step, carry out maximum Iterations terminates judging.

A kind of region cool and thermal power comprehensive energy based on graph theory distributes device rationally, including：Optimized model builds module, excellent Change allocation models to build module and distribute module rationally.

Optimized model builds module, for according to the minimum object function of region cool and thermal power comprehensive energy operating cost and in advance The constraints for first arranging, builds the minimum Optimized model of region cool and thermal power comprehensive energy expense；

The object function of the minimum Optimized model of the region cool and thermal power comprehensive energy expense is:

In above formula, E_iIt is i-th user power utilization expense minimum of a value of region cool and thermal power comprehensive energy expense；Q_iIt is that region is cold I-th user's heat energy expense minimum of a value of thermoelectricity comprehensive energy expense；W_iIt is i-th use of region cool and thermal power comprehensive energy expense Family cold energy expense minimum of a value；N is the total quantity of region cool and thermal power energy mix user；

Its constraints includes：

(1) region cool and thermal power comprehensive energy trend constraint

In above formula, N_sFor total node number set；G_ij、B_ijFor the admittance coefficient between node i, j；V_i, V_jFor node i, j's Voltage magnitude；P_Gi、Q_GiActive the exerting oneself of the respectively generator of node i is exerted oneself with idle；P_Di、Q_DiRespectively node i is active And load or burden without work；θ_ijIt is the generator rotor angle between node i and node j；

(2) region cool and thermal power comprehensive energy node power output constraint

P_imin≤P_i≤P_imax

In above formula, P_iminAnd P_imaxIt is respectively cool and thermal power comprehensive energy node power output lower limit and the upper limit；

(3) cool and thermal power comprehensive energy in region intercouples constraint

P_inmax≤P_c≤P_outmax

In above formula, P_inmaxAnd P_outmaxIt is respectively the lower limit and the upper limit to bulk power grid sale of electricity power.

Optimal Allocation Model builds module, for building the region cool and thermal power synthesis of module construction based on the optimization module The minimum Optimized model of energy expenditure, connects with reference to the network topology of cool and thermal power integrated energy system, sets up based on the area of graph theory Domain cool and thermal power comprehensive energy Optimal Allocation Model；

The Optimal Allocation Model build module specifically for：The network topology for considering actual cool and thermal power comprehensive energy connects Connect, region cool and thermal power comprehensive energy is distributed rationally and is equivalent to a undirected authorized graph G=(V, E, W), find undirected authorized graph G Minimum weights W, so as to distribute cool and thermal power comprehensive energy rationally；

Wherein, the element in set V is the fixed point or node or point of undirected authorized graph G, represents reality in regional extent One cool and thermal power power supply, the finite non-NULL node set that it is made up of region cool and thermal power power supply, according to the number of cool and thermal power power supply Amount order starts number consecutively from 1, until all of cool and thermal power power supply is numbered and finished；

Wherein, the element of set E is side or the line of undirected authorized graph G, represent cool and thermal power power supply between interconnection switch or Person's PCC switch sets, can use e_ijRepresent, e_ijValue be 1 and 0；Wherein, 1 represent cool and thermal power power supply i and cool and thermal power power supply j it Between exist contact, 0 represent without contact；E represents side collection in V；

Wherein, the element of set W is the active power cross-over value between any two node, the power of referred to as undirected authorized graph G Value, W_ijNode i is represented to the active power value exchanged between node j, value is positive number when active power is flowed into, otherwise when stream Value is negative when going out active power.

Module is distributed rationally, for solving based on the region cool and thermal power comprehensive energy optimization of graph theory using improved adaptive GA-IAGA Allocation models, realizes distributing rationally for cool and thermal power comprehensive energy.

The module of distributing rationally is specifically for performing：

(1), according to the constructed region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory, the reality of engineering is taken into account Border data, initialize to the force data that goes out of the cool and thermal power network parameter in region and each cool and thermal power power supply；

(2) binary coding scheme is adopted, the fortune of grid-connected PCC switches, single cool and thermal power power supply to each cool and thermal power power supply Three coding decision variables of row expense and capacity are encoded；

First, the chromosome that coding total length is 13 is chosen, described three coding decision variables is each constituted into chromosome Word string, is then linked to be a complete chromosome by above-mentioned chromosome word string；Front two represents that the grid-connected PCC of cool and thermal power power supply is opened Close, middle five operating costs for representing single cool and thermal power power supply, last six capacity for representing single cool and thermal power power supply；

(3) whether to reach the iterations of maximum as the foundation for terminating, judge whether to meet end condition；If reaching It is i.e. out of service, obtain final result；Otherwise trigger parameter determination sub-module；

(4) Population Size is set, genetic algorithm inherent parameters is determined by the operation method for selecting, intersect and make a variation；

Wherein, adaptive crossover operator function is as follows：

In above formula, f_avgIt is the average fitness of per generation colony；f_maxFor fitness maximum in the individuality to be intersected；F is Larger fitness in two individualities to be intersected；P_c1Value 0.9, P_c2Value 0.6；

Wherein, TSP question rate function is as follows：

P_m=P_m1-P_m1×i/N

In above formula, P_m1For the initial value of aberration rate, value 0.08；I is current iteration number of times；N is iteration total degree；

(5) the cool and thermal power configuration of respective nodes is changed, and it is minimum according to cool and thermal power comprehensive energy expense in region described in step 1 Bound for objective function, the configuration to cool and thermal power power supply judges one by one according to constraints, if being satisfied by condition, Into (6th) step；Otherwise, iterations N adds 1, and returns (4th) step；

(6) optimization aim and its constraint bar of the comprehensive energy expense minimum of cool and thermal power energy mix in region are considered Part, builds fitness function as follows；

In above formula, C_maxIt is a given constant, f (x) is the object function after being normalized；

(7) optimal value that fitness function is produced is replaced, and iterations N is added 1, return (4th) step, carry out maximum Iterations terminates judging.

In the present embodiment, the real data of new Tianjin ecological city Itellectualized uptown in being based on, is based on to one kind of the present invention The region cool and thermal power comprehensive energy Optimal Configuration Method of graph theory is practiced, to verify the feasibility of the inventive method and have Beneficial effect.

The actual cool and thermal power load data of one Itellectualized uptown of ecological city, electric power installed capacity and city planting ductwork in Tianjin Heat supply situation, and the economic data of correlation is as shown in table 1 below.

(table 1)：The overview of Itellectualized uptown

Project	Capacity	Node location	Running status	Economic data
					Thermic load	100M	1,3,7	Annual April November to next year	2.3 yuan/Kwh
Refrigeration duty	120M	1,3,8,9	Annual May is to September	1.8 yuan/Kwh
					Electric load	300M	1～9	Annual January is to December	0.6 yuan/Kwh
Bulk power grid is powered	300M	1～9	Annual January is to December	0.6 yuan/Kwh
					Municipal heat supply	100M	1,3,7	Annual April November to next year	2.3 yuan/Kwh

With the actual motion energy expenditure data of 2015 through adjusting, year, cold and hot electric energy total cost reached 31,470,000 Unit, monthly cool and thermal power energy expenditure reaches 262.25 ten thousand, wherein electricity expense with 28,210,000, municipal heat cost 3,260,000.With the country Same type is compared, and the unit-economy data of mainly cooling and heating load are higher, and real network topology connects with reference to shown in the Fig. 1 for being formed Connect, using the region cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory of the present invention, first build cool and thermal power comprehensive energy The minimum Optimized model of energy operation total cost, the secondly Optimized model of construction based on graph theory, finally using improved genetic algorithms Method is solved, and is optimized calculating.

It is optimized with postponing, in node 3 trilogy supply unit is newly increased, distributed power source is newly increased in node 8, in node 9 increase cold-storage unit on the rocks newly, and its situation is as shown in table 2：

(table 2)：Itellectualized uptown optimization planning is with postponing situation

From the correction data of table 2 and table 1 it is seen that：After distributing rationally, the economic number of the cold and hot electric load of Itellectualized uptown According to there is reduction.Thermic load unit-economy data are reduced to 1.9 yuan/Kwh from 2.3 yuan/Kwh, refrigeration duty unit-economy data from 1.8 yuan/Kwh is reduced to 1.5 yuan/Kwh；Electric load unit-economy number is reduced to 0.5 yuan/Kwh from 0.6 yuan/Kwh.

It is three from the position for distributing node rationally it is found that trilogy supply crew qiting has the node of common factor in cool and thermal power The optimum point of alliance unit.

It is emphasized that embodiment of the present invention is illustrative, rather than it is determinate, therefore present invention bag The embodiment for being not limited to described in specific embodiment is included, it is every by those skilled in the art's technology according to the present invention scheme The other embodiment for drawing, also belongs to the scope of protection of the invention.

Claims (10)

1. the region cool and thermal power comprehensive energy Optimal Configuration Method of graph theory is based on, it is characterised in that：Comprise the following steps：

Step 1, the object function according to region cool and thermal power comprehensive energy operating cost minimum and the constraints for pre-setting, structure Build the minimum Optimized model of region cool and thermal power comprehensive energy expense；

The Optimized model of step 2, the region cool and thermal power comprehensive energy expense minimum built based on step 1, with reference to cool and thermal power synthesis The network topology connection of energy resource system, sets up based on the region cool and thermal power comprehensive energy Optimal Allocation Model of graph theory；

Step 3, the region cool and thermal power comprehensive energy Optimal Allocation Model using improved adaptive GA-IAGA solution based on graph theory, realize cold Thermoelectricity comprehensive energy is distributed rationally.

2. the cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory according to claim 1, it is characterised in that：It is described The object function of the minimum Optimized model of the region cool and thermal power comprehensive energy expense of step 1 is:

$\begin{array}{cc}min& (\underset{i=1}{\overset{n}{\Σ}}{E}_i+Q_i+W_i)\end{array}$

In above formula, E_iIt is i-th user power utilization expense minimum of a value of region cool and thermal power comprehensive energy expense；Q_iIt is region cool and thermal power I-th user's heat energy expense minimum of a value of comprehensive energy expense；W_iBe region cool and thermal power comprehensive energy expense i-th user it is cold Can expense minimum of a value；N is the total quantity of region cool and thermal power energy mix user；

Its constraints includes：

(1) region cool and thermal power comprehensive energy trend constraint

$g\left(x\right)=\left\{\begin{array}{cc}{P}_{Gi}-P_{Di}-V_i\underset{j=1}{\overset{N_s}{\Σ}}{V}_j(G_{ij}{\mathrm{cos\θ}}_{ij}+B_{ij}{\mathrm{sin\θ}}_{ij})=0& i\∈N_s\\ Q_{Gi}-Q_{Di}-V_i\underset{j=1}{\overset{N_s}{\Σ}}{V}_j(G_{ij}{\mathrm{sin\θ}}_{ij}-B_{ij}{\mathrm{cos\θ}}_{ij})=0& i\∈N_s\end{array}\right.$

In above formula, N_sFor total node number set；G_ij、B_ijFor the admittance coefficient between node i, j；V_i, V_jFor node i, the voltage of j Amplitude；P_Gi、Q_GiActive the exerting oneself of the respectively generator of node i is exerted oneself with idle；P_Di、Q_DiThe respectively active and nothing of node i Workload；θ_ijIt is the generator rotor angle between node i and node j；

(2) region cool and thermal power comprehensive energy node power output constraint

P_imin≤P_i≤P_imax

In above formula, P_iminAnd P_imaxIt is respectively cool and thermal power comprehensive energy node power output lower limit and the upper limit；

(3) cool and thermal power comprehensive energy in region intercouples constraint

P_inmax≤P_c≤P_outmax

In above formula, P_inmaxAnd P_outmaxIt is respectively the lower limit and the upper limit to bulk power grid sale of electricity power.

3. the cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory according to claim 1 and 2, it is characterised in that： The concrete grammar of the step 2 is：The network topology connection of actual cool and thermal power comprehensive energy is considered, by region cool and thermal power synthesis Energy source optimization configuration is equivalent to a undirected authorized graph G=(V, E, W), the minimum weights W of undirected authorized graph G is found, so as to excellent Change configuration cool and thermal power comprehensive energy；

Wherein, the element in set V is the fixed point or node or point of undirected authorized graph G, represents one actual in regional extent Cool and thermal power power supply, the finite non-NULL node set that it is made up of region cool and thermal power power supply, the quantity according to cool and thermal power power supply is suitable Sequence starts number consecutively from 1, until all of cool and thermal power power supply is numbered and finished；

Wherein, the element of set E is side or the line of undirected authorized graph G, represents interconnection switch or PCC between cool and thermal power power supply Switch set, can use e_ijRepresent, e_ijValue be 1 and 0；Wherein, 1 represent and deposited between cool and thermal power power supply i and cool and thermal power power supply j Represent without contact in contact, 0；E represents side collection in V；

Wherein, the element of set W is the active power cross-over value between any two node, the weights of referred to as undirected authorized graph G, W_ij Node i is represented to the active power value exchanged between node j, value is positive number when active power is flowed into, otherwise when outflow has Value is negative during work(power.

4. the region cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory according to claim 1 and 2, its feature exists In：The concrete steps of the step 3 include：

(1) according to the region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory constructed by step 2, the reality of engineering is taken into account Border data, initialize to the force data that goes out of the cool and thermal power network parameter in region and each cool and thermal power power supply；

(2) binary coding scheme is adopted, the running cost of grid-connected PCC switches, single cool and thermal power power supply to each cool and thermal power power supply Encoded with three coding decision variables with capacity；

(3) whether to reach the iterations of maximum as the foundation for terminating, judge whether to meet end condition；If reach moving back Go out operation, obtain final result；Otherwise enter (4th) step；

(4) Population Size is set, genetic algorithm inherent parameters is determined by the operation method for selecting, intersect and make a variation；

Wherein, adaptive crossover operator function is as follows：

$P_c=\left\{\begin{array}{cc}{P}_{c1}-\frac{(P_{c1}-P_{c2})\&times;(f_{max}-f)}{f_{\max }-f_{avg}},& f\&GreaterEqual;f_{avg}\\ P_{c1},& f<f_{avg}\end{array}\right.$

In above formula, f_avgIt is the average fitness of per generation colony；f_maxFor fitness maximum in the individuality to be intersected；F is to hand over Larger fitness in two individualities of fork；P_c1Value 0.9, P_c2Value 0.6；

Wherein, TSP question rate function is as follows：

P_m=P_m1-P_m1×i/N

In above formula, P_m1For the initial value of aberration rate, value 0.08；I is current iteration number of times；N is iteration total degree；

(5) the cool and thermal power configuration of respective nodes is changed, and according to the minimum mesh of cool and thermal power comprehensive energy expense in region described in step 1 The constraints of scalar functions, the configuration to cool and thermal power power supply judges one by one according to constraints, if being satisfied by condition, enters (6th) step；Otherwise, iterations N adds 1, and returns (4th) step；

(6) optimization aim and its constraints of the comprehensive energy expense minimum of cool and thermal power energy mix in region are considered, Build fitness function as follows；

In above formula, C_maxIt is a given constant, f (x) is the object function after being normalized；

(7) optimal value that fitness function is produced is replaced, and iterations N is added 1, return (4th) step, carry out the iteration of maximum Number of times terminates judging.

5. the cool and thermal power comprehensive energy Optimal Configuration Method based on graph theory according to claim 4, it is characterised in that：It is described The concrete grammar of (2nd) step of step 3 is：First, the chromosome that coding total length is 13 is chosen, by described three decision-making is encoded Variable each constitutes chromosome word string, then above-mentioned chromosome word string is linked to be into a complete chromosome；Front two represents cold The grid-connected PCC switches of thermoelectric power source, middle five operating costs for representing single cool and thermal power power supply, last six expressions are single cold The capacity of thermoelectric power source.

6. device is distributed rationally based on the region cool and thermal power comprehensive energy of graph theory, it is characterised in that：Including：

Optimized model builds module, for setting according to the minimum object function of region cool and thermal power comprehensive energy operating cost and in advance The constraints put, builds the minimum Optimized model of region cool and thermal power comprehensive energy expense；

Optimal Allocation Model builds module, for building the region cool and thermal power comprehensive energy of module construction based on the optimization module The minimum Optimized model of expense, connects with reference to the network topology of cool and thermal power integrated energy system, and the region set up based on graph theory is cold Thermoelectricity comprehensive energy Optimal Allocation Model；

Distribute module rationally, distributed rationally based on the region cool and thermal power comprehensive energy of graph theory for being solved using improved adaptive GA-IAGA Model, realizes distributing rationally for cool and thermal power comprehensive energy.

7. the cool and thermal power comprehensive energy based on graph theory according to claim 6 distributes device rationally, it is characterised in that：It is described The object function of the minimum Optimized model of region cool and thermal power comprehensive energy expense is:

$\begin{array}{cc}min& (\underset{i=1}{\overset{n}{\&Sigma;}}{E}_i+Q_i+W_i)\end{array}$

In above formula, E_iIt is i-th user power utilization expense minimum of a value of region cool and thermal power comprehensive energy expense；Q_iIt is region cool and thermal power I-th user's heat energy expense minimum of a value of comprehensive energy expense；W_iBe region cool and thermal power comprehensive energy expense i-th user it is cold Can expense minimum of a value；N is the total quantity of region cool and thermal power energy mix user；

Its constraints includes：

(1) region cool and thermal power comprehensive energy trend constraint

$g\left(x\right)=\left\{\begin{array}{cc}{P}_{Gi}-P_{Di}-V_i\underset{j=1}{\overset{N_s}{\&Sigma;}}{V}_j(G_{ij}{\mathrm{cos\&theta;}}_{ij}+B_{ij}{\mathrm{sin\&theta;}}_{ij})=0& i\&Element;N_s\\ Q_{Gi}-Q_{Di}-V_i\underset{j=1}{\overset{N_s}{\&Sigma;}}{V}_j(G_{ij}{\mathrm{sin\&theta;}}_{ij}-B_{ij}{\mathrm{cos\&theta;}}_{ij})=0& i\&Element;N_s\end{array}\right.$

In above formula, N_sFor total node number set；G_ij、B_ijFor the admittance coefficient between node i, j；V_i, V_jFor node i, the voltage of j Amplitude；P_Gi、Q_GiActive the exerting oneself of the respectively generator of node i is exerted oneself with idle；P_Di、Q_DiThe respectively active and nothing of node i Workload；θ_ijIt is the generator rotor angle between node i and node j；

(2) region cool and thermal power comprehensive energy node power output constraint

P_imin≤P_i≤P_imax

In above formula, P_iminAnd P_imaxIt is respectively cool and thermal power comprehensive energy node power output lower limit and the upper limit；

(3) cool and thermal power comprehensive energy in region intercouples constraint

P_inmax≤P_c≤P_outmax

In above formula, P_inmaxAnd P_outmaxIt is respectively the lower limit and the upper limit to bulk power grid sale of electricity power.

8. device is distributed rationally based on the cool and thermal power comprehensive energy of graph theory according to claim 6 or 7, it is characterised in that： The Optimal Allocation Model build module specifically for：The network topology connection of actual cool and thermal power comprehensive energy is considered, by area Domain cool and thermal power comprehensive energy is distributed rationally and is equivalent to a undirected authorized graph G=(V, E, W), finds the minimum of undirected authorized graph G Weights W, so as to distribute cool and thermal power comprehensive energy rationally；

Wherein, the element in set V is the fixed point or node or point of undirected authorized graph G, represents one actual in regional extent Cool and thermal power power supply, the finite non-NULL node set that it is made up of region cool and thermal power power supply, the quantity according to cool and thermal power power supply is suitable Sequence starts number consecutively from 1, until all of cool and thermal power power supply is numbered and finished；

Wherein, the element of set E is side or the line of undirected authorized graph G, represents interconnection switch or PCC between cool and thermal power power supply Switch set, can use e_ijRepresent, e_ijValue be 1 and 0；Wherein, 1 represent and deposited between cool and thermal power power supply i and cool and thermal power power supply j Represent without contact in contact, 0；E represents side collection in V；

Wherein, the element of set W is the active power cross-over value between any two node, the weights of referred to as undirected authorized graph G, W_ij Node i is represented to the active power value exchanged between node j, value is positive number when active power is flowed into, otherwise when outflow has Value is negative during work(power.

9. device is distributed rationally based on the region cool and thermal power comprehensive energy of graph theory according to claim 6 or 7, its feature exists In：Module is distributed rationally specifically for performing：

(1), according to the constructed region cool and thermal power comprehensive energy Optimal Allocation Model based on graph theory, the actual number of engineering is taken into account According to initializing to the force data that goes out of the cool and thermal power network parameter in region and each cool and thermal power power supply；

(2) binary coding scheme is adopted, the running cost of grid-connected PCC switches, single cool and thermal power power supply to each cool and thermal power power supply Encoded with three coding decision variables with capacity；

(3) whether to reach the iterations of maximum as the foundation for terminating, judge whether to meet end condition；If reach moving back Go out operation, obtain final result；Otherwise trigger parameter determination sub-module；

(4) Population Size is set, genetic algorithm inherent parameters is determined by the operation method for selecting, intersect and make a variation；

Wherein, adaptive crossover operator function is as follows：

$P_c=\left\{\begin{array}{cc}{P}_{c1}-\frac{(P_{c1}-P_{c2})\&times;(f_{max}-f)}{f_{\max }-f_{avg}},& f\&GreaterEqual;f_{avg}\\ P_{c1},& f<f_{avg}\end{array}\right.$

In above formula, f_avgIt is the average fitness of per generation colony；f_maxFor fitness maximum in the individuality to be intersected；F is to hand over Larger fitness in two individualities of fork；P_c1Value 0.9, P_c2Value 0.6；

Wherein, TSP question rate function is as follows：

P_m=P_m1-P_m1×i/N

In above formula, P_m1For the initial value of aberration rate, value 0.08；I is current iteration number of times；N is iteration total degree；

(5) the cool and thermal power configuration of respective nodes is changed, and according to the minimum mesh of cool and thermal power comprehensive energy expense in region described in step 1 The constraints of scalar functions, the configuration to cool and thermal power power supply judges one by one according to constraints, if being satisfied by condition, enters (6th) step；Otherwise, iterations N adds 1, and returns (4th) step；

(6) optimization aim and its constraints of the comprehensive energy expense minimum of cool and thermal power energy mix in region are considered, Build fitness function as follows；

In above formula, C_maxIt is a given constant, f (x) is the object function after being normalized；

(7) optimal value that fitness function is produced is replaced, and iterations N is added 1, return (4th) step, carry out the iteration of maximum Number of times terminates judging.

10. the cool and thermal power comprehensive energy based on graph theory according to claim 9 distributes device rationally, it is characterised in that：It is excellent Change configuration module specifically for：First, the chromosome that coding total length is 13 is chosen, by described three coding decision variables each Composition chromosome word string, is then linked to be a complete chromosome by above-mentioned chromosome word string；Front two represents cool and thermal power power supply Grid-connected PCC switches, middle five operating costs for representing single cool and thermal power power supply, last six represent single cool and thermal power power supply Capacity.