Summary of the invention
It is an object of the invention to overcome the problems of the above-mentioned prior art and provide it is a kind of quickly, global optimum, can
The integrated energy system energy transition equipment configuration method and configuration system of scalability.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of integrated energy system energy transition equipment configuration method, this method are based on total system transfer efficiency matrix and set
Matrix building integrated energy system multipotency stream joint equality constraint is changed, objective function is minimised as with year investment and operating cost,
Realize the configuration optimization to energy transition equipment.
Further, described based on total system transfer efficiency matrix and permutation matrix building integrated energy system multipotency stream connection
Close equality constraint specifically:
Each network of integrated energy system can be flowed and energy transition equipment carries out joint modeling, network coupling unit passes through more
Vector transfer efficiency matrix and permutation matrix are characterized.
Further, the transformational relation of the transfer efficiency matrix description different type energy, permutation matrix describe the energy
Node locating for conversion equipment is in whole energy network mapping relations corresponding with each separate network.
Further, rate of load condensate is introduced in the objective function, specifically:
In formula, NconFor the quantity of conversion equipment,For the initial outlay expense of conversion equipment i,WithThe electric power and heating power power that respectively equipment i is provided in load peak, ∨ indicate operator "or", ηloadfactor
For rate of load condensate, CO&M is operation and maintenance cost, CeFor power purchase price,For the superior power grid power purchase in load peak
Power,For the gas consumption rate total in load peak, ξeFor electric power carbon intensity, ξgFor natural gas carbon row
Put intensity, CcarbonFor carbon price.
Further, constraint condition when carrying out the configuration optimization further includes inequality constraints, the inequality constraints
Energy transition equipment capacity bound constraint, power network variable bound constraint, heat supply network variable bound constraint and
The bound of fuel gas network variable constrains.
Further, the configuration optimization to energy transition equipment includes the type for configuring energy transition equipment and each energy
The capacity of source conversion equipment.
Further, the energy transition equipment includes CHP, heat pump, gas fired-boiler, gas turbine and boilers heated electrically.
The present invention also provides a kind of integrated energy system energy transition equipments to configure system, including memory and processor,
The memory is stored with computer program, and the processor calls the computer program to execute the configuration method.
Compared with prior art, the invention has the following beneficial effects:
1) transfer efficiency matrix of the invention and permutation matrix method carry out total system Holistic modeling, are retouched with unified mode
The energy conversion and flowing of the different type energy are stated, energy transition equipment, which is systematically associated with respective electricity/heat/gas net, to flow
Equation forms whole energy resource system physical equation model, and scalability is high, is efficiently integrated and coordinates all kinds of energy transition equipments.
2) present invention introduces rate of load condensate in objective function and parallel computation has greatly accelerated planning problem globally optimal solution
Convergence rate.
3) according to conjunctive model multipotency stream calculation, the unified collocation of integrated energy system conversion equipment and energy storage is realized,
Be conducive to the reciprocation of state variable and each via net loss between network modelling and meter and heterogeneous networks.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with technical solution of the present invention
Premised on implemented, the detailed implementation method and specific operation process are given, but protection scope of the present invention is not limited to
Following embodiments.
1, integrated energy system joint modeling
Integrated energy system is that electricity, heat, gas energy mix system are established natural by the analogy with Power System Analysis
The separate networks such as gas net, heat supply network can flow analysis model.Electric power system tide is calculated for analyzing each node of power circuit
Voltage, power, to be powered ability check, line loss analyzing etc..Hot-fluid calculate for analyze heat-transfer medium temperature, flow,
Pressure, state change situation, to analyze heat-energy losses, flowing pressure loss.Airometer point counting analyses combustion gas in transmission & distribution feed channel
Interior flow, pressure, state change situation, to analyze gas flow pressure loss.
Each network of integrated energy system can be flowed and energy transition equipment carries out joint modeling analysis, consider CHP (thermoelectricity connection
Produce), heat pump, the connection such as gas fired-boiler power network, heat supply network and fuel gas network coupling element, wherein network coupling unit passes through more
Vector transfer efficiency matrix and permutation matrix are characterized.
The transfer efficiency matrix of all energy transition equipments is constructed according to input data, describes the conversion of the different type energy
Relationship, as shown in table 1.
The sample of 1 transfer efficiency matrix of table
With transfer efficiency matrixSystematization Holistic modeling is carried out, whole energy resource system physical equation equality constraint is formed:
In formula, i is the number of nodes of conversion equipment, and j is column index;Matrix EoutExpression conversion equipment use demand (for example,
Electrical power and thermal power), it is made of user input data,For matrix EoutElement;Matrix EinIndicate that conversion equipment is defeated
Enter electrical power and gas discharge,For matrix EinElement;For matrixElement;MatrixEinWith EoutHave
Identical dimension, such as Ncon×4;ElementEqual to elementDivided by the corresponding element of transfer efficiency matrix
In example as shown in Figure 1, forColumn include ge, gh, eh, eg, are that gas is arrived to electricity, gas to heat, electricity respectively
Heat, electricity arrive gas;Matrix EinColumn include gas, gas, electricity, electricity;Matrix EoutColumn include electricity, heat, heat, gas.
Permutation matrix is constructed, node locating for equipment is described and is closed in the mapping corresponding with each separate network of whole energy network
System, as shown in Figure 1.Assuming that the conversion equipment that global number is 1,2,3 and 4 be respectively gas turbine, cogeneration of heat and power CHP, heat pump,
Gas fired-boiler.The conversion equipment that one global number is 2 corresponds to the node 5, the node 3 in heating power network and combustion of electric power networks
Node 4 in gas network.
The vector of device node number and the mapping process that each network is locally numbered are as shown in Figure 2.For example, global node is compiled
Numbers 2 conversion equipment variate-valueIt is respectively mapped to the electrical power of grid nodes 5Heat supply network node 3
Thermal powerThe gas discharge of gentle net node 4
The mapping relations of energy transition equipment and power network corresponding node form electric power permutation matrix Me, wherein (Isosorbide-5-Nitrae),
The element of (2,5), (3,3) position is set to 1, and other all elements are set to 0.Therefore, couple energy transition equipment
With the permutation matrix M of power grideIt indicates are as follows:
Corresponding with Fig. 2, the energy transition equipment with coupled relation is indicated in the vector power of power network are as follows:
Given that it is known that the output (vector) of each energy device obtains the input of each energy device according to transfer efficiency matrix
(vector).By by the input vector and output vector of device numbering sequence, it is separated into electricity/independent vector of 3 groups of heat/gas.Pass through
Electricity/heat/gas vector of energy device is mapped to the vector power of all nodes of each network by permutation matrix by node.Each node
Electrical power vector, for calculating electric power system tide, each node thermal power and natural gas flow, for carrying out the heating power stream of heat supply network
It calculates and is calculated with the natural gas flow of fuel gas network.Conversion equipment model is mapped to the more of electricity/heat/gas network equation (such as formula (4))
The flow chart that conversion process can be flowed is as shown in Figure 3.
2, the energy transition equipment configuration of integrated energy system
The present invention comprehensively considers cold heat/electricity/gas energy distribution situation and load prediction situation in selected areas, with rule
The minimum target of year gross investment operating cost in the phase of drawing obtains the allocation optimum of the energy transition equipment of multipotency collaboration, including
Configure the type of energy transition equipment and the capacity of each energy transition equipment.Cost of investment passes through total equipment investment depreciation to its fortune
It is often determined over one year in the departure date.In the solution procedure of allocation optimum, by the way that the analysis of multipotency flow point is configured to the flow-optimized mould of multipotency
The equality constraint of type calls the interior point method solver of optimization algorithm to solve the problems, such as that multipotency is flow-optimized.
The present invention also provides a kind of integrated energy system energy transition equipments to configure system, including memory and processor,
The memory is stored with computer program, and the processor calls the computer program to execute configuration method.The calculating
Machine program includes input module, configuration module and output module.As shown in figure 4, input module acquisition planning region in it is cold/
Heat/electricity/gas net status parameter, Future New Energy Source and load power output, fuel price and carbon emission settlement price, network equipment parameter,
Economic parameters etc.;Configuration module carries out energy transition equipment configuration optimization according to the input data, and the constraint of configuration optimization includes
Equality constraint and inequality constraints;Output module obtain capacity, investment and operating cost etc. of energy transition equipment access as a result,
Output module is also exportable to evaluation informations such as the operational efficiency, reliability, economy of configuration result.
2.1, objective function
The method of the present invention obtains the integrated energy system conversion equipment of multipotency collaboration by minimizing year Integrated Cost
Allocation optimum realizes that entirety energy resource system investment cost CAPEX and operating cost OPEX is minimum, optimizes each network energy conversion
The production and conversion of electric power of equipment and heating power.The objective function of this method is that yearization investment is minimum with operating cost, and expression formula is such as
Under:
Min Ctotal=LFn·CCapex+COpex (4)
In formula, Ctotal- annual total cost;COPEX- annual equipment operating cost, includes superior power grid power purchase expense and combustion
Expect cost, carbon emission expense, operation and maintenance expense;CCapexThe initial total investment expenses of-all devices;LFn- recovery of the capital system
Number,D is discount rate, it is assumed that d=6%, n are the service life time limits of equipment.
The capacity of conversion equipment is equal to the electric power or heating power power that equipment i is provided in load peak, and cost of investment is (i.e.
The initial total investment expenses of all devices) as follows:
In formula, NconThe quantity of-conversion equipment;The initial outlay expense of-conversion equipment i;WithThe electric power and heating power power that-equipment i is provided in load peak, ∨ indicate operator "or".
Operating cost includes fuel cost, operation and maintenance (O&M) expense and carbon emission expense:
In formula, CO&M- operation and maintenance cost does not include fuel cost, Ce- power purchase price (£/MWh), Pimport- to
Higher level's power grid power purchase power (MW), Cg- gas price (£/MWh), vgtotal- total gas consumption rate (m3/ h), ξe- electricity
Power carbon intensity (gCO2/ kWh), ξg- natural gas carbon intensity (gCO2/ kWh), Ccarbon- carbon price (£/ton), T-
1 year hourage (T=8760).
It based on annual each hour load optimal operating cost, then adds up and obtains annual operating cost, such method calculates
Time is too slow, therefore using year electric power and heating power peak load and introduces rate of load condensate.In this case, operating cost is calculated
Using only peak load, rather than 8760 hours loads.The method significantly reduces the complex nature of the problem and calculates the time,
If operation points are 8760 × 1 hours, runing time shortens 8760 times.Therefore, equation is simplified as:
In formula,Systematic running cost when-peak load is used, ηloadfactor- rate of load condensate.
In conclusion complete objective function expression formula are as follows:
2.2, constraint condition
2.2.1 equality constraint
Equality constraint includes the Nonlinear Equality Constrained of integrated energy system.
The coupled relation of comprehensive energy network passes through decision variableThe connection of connection, these variables is logical
Cross transfer efficiency matrixIt is described, as shown in the equality constraint (9) of energy transition equipment.
Integrated energy system electric power, heating power, combustion gas and the multipotency stream power equation constraint of balancing the load are comprehensive energy
System combined equation.Pass through transfer efficiency matrixAnd permutation matrix, the conversion being connected in multipotency stream joint equation are set
Standby model (P, Φ, vq), realize the coupling of electricity/heat/gas network.
In formula: P is electric power active power vector, and Q is electric power reactive power vector, and Φ is consumed or mentioned in each heating power node
The thermal power vector of confession, vqIt is node air pressure;K is the resistance coefficient of each pipeline, and A is network associate matrix, and B is loop association square
Battle array, subscript h indicate that heat supply network, subscript g indicate fuel gas network;C is coefficient matrix and b is the solution formed by heating network thermodynamic model
Column vector;Subscript s indicates to supply water, and subscript r indicates return water.
State variable and control variable are expressed as:
P in formulaconWithIt is the electric power active power and heating power power (MW) of conversion equipment.
Shown in transfer efficiency matrix and permutation matrix equation such as formula (1), (2).
2.2.2 inequality constraints
1) conversion equipment
The bound constraint representation of each energy transition equipment capacity are as follows:
In formula, subscript lb and ub indicate the bound of place capacity, and subscript con indicates conversion equipment.
2) power network
The upper voltage limit and lower limit of each node
Control the upper and lower bound of variable: the active and reactive power output of generator
The upper and lower bound of each branch electrical power
In formula, NeIt is interstitial content, NgeIt is generator number, PlIt is the effective power flow of each route, NleIt is number of lines.
3) heat supply network
The bound constraint representation of each pipeline flow are as follows:
The bound constraint representation of each node significance for water supply are as follows:
Ts_min≤Ts≤Ts_max (13)
Tr_min≤Tr≤Tr_max
4) fuel gas network
The gaseous-pressure bound constraint representation of each node are as follows:
pgmin≤pg≤pgmax (14)
In formula, pgIt (bar) is gaseous-pressure.
The gas flow bound constraint representation of each pipeline are as follows:
vg_min≤vg≤vg_max (15)
In formula, vg(m3/ h) be pipeline gas flow rate.