CN115099007A - Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve - Google Patents
Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve Download PDFInfo
- Publication number
- CN115099007A CN115099007A CN202210608295.5A CN202210608295A CN115099007A CN 115099007 A CN115099007 A CN 115099007A CN 202210608295 A CN202210608295 A CN 202210608295A CN 115099007 A CN115099007 A CN 115099007A
- Authority
- CN
- China
- Prior art keywords
- energy
- comprehensive
- cost
- energy consumption
- gas
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- 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
Abstract
An integrated energy system optimization operation method based on an integrated cost-energy consumption curve comprises the following steps: establishing a park comprehensive energy system integrated model considering variable working condition characteristics of energy conversion equipment; constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park; calculating an optimized operation scheme of the comprehensive energy system based on the comprehensive cost-energy consumption curve; the invention establishes the integrated model of the comprehensive energy system considering the variable working condition running characteristics of the multi-type energy conversion equipment, and effectively reduces the energy deviation between the optimized scheduling scheme and the actual load requirement; a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park is constructed, and guidance can be provided for scheduling personnel in actual engineering application through analysis and calculation of a curve image guidance optimization operation strategy; the method can effectively give consideration to model precision and solving efficiency, greatly improves the accuracy of the optimized operation scheme, and ensures shorter optimized calculation time.
Description
Technical Field
The invention relates to an optimized operation method of a comprehensive energy system. In particular to an optimized operation method of a comprehensive energy system based on a comprehensive cost-energy consumption curve.
Background
The comprehensive energy system can fully exert the comprehensive advantages of complementation and mutual assistance of various heterogeneous energy sources, and the development of the comprehensive energy system has great significance for building a clean, low-carbon, safe and efficient modern energy system. The park comprehensive energy system serves as a terminal for energy interconnection, and can provide reliable and economic energy supply for users in the park by coordinating various energy devices.
Because the collaborative optimization operation model of the park comprehensive energy system is extremely complex, the current research is usually simplified to a certain extent, and most typically, the efficiency of the energy conversion equipment is set to be a constant. In practice, the efficiency of the energy conversion device may vary with load factor, temperature, air pressure, etc. However, the overall consideration of the variable working condition operation characteristics of various types of energy conversion equipment by domestic and foreign researches is not enough, and the solution of the optimized operation strategy faces the contradiction between the requirement of solution precision and the complexity of the model.
The cooperative optimization operation problem of the park comprehensive energy system considering the variable working condition characteristics of the equipment is a typical non-convex optimization problem. The common solution thought of the non-convex optimization problem in engineering can be divided into piecewise linear function approximation and heuristic algorithm. However, error accumulation is easy to form by utilizing a piecewise linear function, and the solving precision is limited by the number of segments; the solution is unstable by using a heuristic algorithm, and the approximation degree of the obtained solution and the optimal solution is difficult to explain. In view of this, it is necessary to design an optimal operation scheme calculation strategy that takes model precision and solution efficiency into consideration.
Disclosure of Invention
The invention aims to solve the technical problem of providing an optimized operation method of a comprehensive energy system based on a comprehensive cost-energy consumption curve, which can take model precision and solving efficiency into consideration, in order to overcome the defects of the prior art.
The technical scheme adopted by the invention is as follows: an integrated energy system optimization operation method based on an integrated cost-energy consumption curve comprises the following steps:
1) establishing a park comprehensive energy system integrated model considering variable working condition characteristics of energy conversion equipment; comprises that
(1.1) establishing a variable working condition operation model of energy conversion equipment of a combined heat and power generation unit, a gas boiler, an absorption refrigerator and an electric refrigerator, wherein the variable working condition operation model of the energy conversion equipment is as follows:
wherein, P out For outputting energy, P, to energy-converting apparatus in Energy is input to the energy conversion device, eta is the energy conversion efficiency of the energy conversion device, eta R For the energy conversion efficiency of the energy conversion device under rated operating conditions,
the efficiency correction coefficient is the efficiency correction coefficient of a cogeneration unit or a gas boiler or an absorption refrigerator or an electric refrigerator;
(1.2) establishing an integrated model of the comprehensive energy system;
2) constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park; the method comprises the following steps:
(2.1) integrating an objective function of the optimized operation of the comprehensive energy system;
(2.2) constructing a net cost-output curve of the energy conversion equipment;
(2.3) constructing a comprehensive cost-energy consumption curve facing the multi-element energy consumption requirement of the park;
and (2.4) carrying out regional division on the comprehensive cost-energy consumption curve of the diversified energy consumption requirements of the facing park.
3) Calculating an optimized operation scheme of the comprehensive energy system based on the comprehensive cost-energy consumption curve; comprises that
(3.1) calculating an optimized operation scheme of the comprehensive energy system without the energy storage equipment;
and (3.2) calculating an optimized operation scheme of the integrated energy system comprising the energy storage equipment.
The efficiency correction coefficients of the cogeneration unit, the gas boiler, the absorption chiller and the electric chiller in the step 1) and the step 1.1 are respectively as follows:
efficiency correction coefficient of cogeneration unit:
wherein the content of the first and second substances,
the efficiency correction coefficient of the cogeneration unit; k is a radical of CHP,n An n-order fitting coefficient of the electric efficiency of the cogeneration unit; n is a radical of CHP The electric load rate of the cogeneration unit;
the correction coefficient is the heat-electricity ratio of the cogeneration unit; k is a radical of α,n Fitting coefficients of an nth order of a thermoelectric ratio of the cogeneration unit;
efficiency correction factor of gas boiler:
wherein, the first and the second end of the pipe are connected with each other,
correcting the coefficient for the efficiency of the gas boiler; k is a radical of GB,n An n-order fitting coefficient of the gas boiler efficiency; n is a radical of GB Is the load factor of the gas boiler;
efficiency correction coefficient of absorption chiller:
wherein the content of the first and second substances,
correcting the coefficient for the efficiency of the absorption chiller; k is a radical of AC,n Is an n-order fitting coefficient of the efficiency of the absorption refrigerator; n is a radical of AC Is the load factor of the absorption chiller;
efficiency correction coefficient of electric refrigerator:
wherein the content of the first and second substances,
the efficiency correction factor of the electric refrigerator is obtained; k is a radical of EC,n Fitting coefficients of the n orders of the efficiency of the electric refrigerator; n is a radical of EC Is the load factor of the electric refrigerator.
Step 1) the integrated model of the comprehensive energy system in the step (1.2) is as follows:
wherein f is 1 ,f 2 ,…,f n An objective function for optimizing operation of the integrated energy system; g () is a supply and demand balance expression of the comprehensive energy system; g () is a device constraint expression of the comprehensive energy system; i represents the ith energy in the comprehensive energy system, C is cold, H is hot, E isElectricity, G is gas; x is the number of i A parameter representing the ith class of energy; t is a time scale; xi is an uncertain variable; l is the load demand; p is an output energy value of the energy supply equipment or the energy conversion equipment or the distribution network, namely an output value; h () is a coupling transformation matrix; η is the energy conversion efficiency of the energy conversion device.
The objective function of the optimized operation of the integrated comprehensive energy system in the step 2) and the step 2.1 is as follows:
f is the comprehensive cost of the optimized operation of the comprehensive energy system after multi-objective integration by using a membership function; f. of 1 ,f 2 ,…,f n Optimizing an objective function for operation of the integrated energy system; u is the energy conversion device type; f u Represents the net cost of the energy conversion device; s i Representing the energy cost of the ith type of energy; l is i Load demand for class i energy; f i The comprehensive cost of the i-th energy is obtained; delta is a repeated calculation part caused by the energy conversion equipment simultaneously generating multiple energies; omega is an energy conversion equipment set; i represents the ith type of energy in the comprehensive energy system of cold, heat, electricity and gas, C is cold, H is heat, E is electricity and G is gas.
Step 2) the net cost-output curve of the energy conversion equipment in the step (2.2) is a curve drawn according to the net cost-output function of the energy conversion equipment of the cogeneration unit, the gas boiler, the absorption refrigerator and the electric refrigerator; the net cost-out function of the energy conversion device is:
wherein, F u Represents the net cost of the energy conversion device;
for the output value of the energy conversion device u, the superscript ij represents the energyThe quantity conversion equipment u converts the ith type energy into the jth type energy; s i Representing the energy cost of the ith type of energy; s j Representing the energy cost of the j-th energy;
the energy conversion efficiency of the energy conversion device u; d u Loss of the energy conversion device u; i and j respectively represent the i and j energy in the cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is hot, E is electricity and G is gas.
Step 2) constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park in the step (2.3), wherein the curve is drawn according to a comprehensive cost-energy consumption function; the integrated cost-energy consumption function is:
wherein, F i The comprehensive cost of the i-th energy is obtained; f u Net cost for energy conversion equipment;
is the output value of the energy conversion device u; s i Energy cost for class i energy; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
for the actual energy demand of the ith type of energy in the park, the numerical value satisfies the following constraints:
wherein L is i Load demand for class i energySolving;
the energy conversion efficiency of the energy conversion device u;
is the output value of the energy supply device v; v is a set of energy supply devices;
outputting an energy value upper limit for the distribution network with the ith type of energy;
is the upper limit of the output value of the energy conversion device u;
is the upper limit of the output value of the energy supply device v.
Step 2), performing area division on the comprehensive cost-energy consumption curve of the diversified energy consumption demand of the facing park in the step (2.4), wherein the area division includes the division into an optimizable area and a conservation and supply area, and the optimizable area comprises: the energy supply system comprises a self-energy supply area, an energy purchasing area and an energy conversion area I, wherein a supply area is an energy conversion area II; wherein:
(2.4.1) self-energizing zone: the energy consumption of the park is generated by wind and light renewable energy sources in the park, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.2) energy purchase zone: energy consumption of the park is met by purchasing energy through an external network;
(2.4.3) energy conversion zone I: the energy consumption of the garden is met by energy conversion equipment, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.4) energy conversion zone II: the energy consumption of the park is met by energy conversion equipment, and the unit energy consumption cost is higher than the energy purchasing cost.
Step 3) the calculation of the optimized operation scheme of the comprehensive energy system without the energy storage device in the step (3.1) comprises the following steps:
(3.1.1) respectively constructing 24 groups of corresponding comprehensive cost-energy consumption curves 24 hours a day according to the comprehensive cost-energy consumption function in the step 2 and the step 2.3 in consideration of time-of-use electricity price, wherein each group of comprehensive cost-energy consumption curves comprises 4 curves of a cold energy comprehensive cost-energy consumption curve, a heat energy comprehensive cost-energy consumption curve, an electric energy comprehensive cost-energy consumption curve and a gas energy comprehensive cost-energy consumption curve;
(3.1.2) establishing 24 groups of energy distribution data tables containing output and purchased energy of park energy equipment according to 24 groups of comprehensive cost-energy consumption curves, wherein each group of energy distribution data tables comprises 4 tables of cold energy distribution, heat energy distribution, electric energy distribution and gas energy distribution;
(3.1.3) according to the comprehensive cost-energy consumption curve of the cold energy and by combining a cold energy distribution data table, obtaining the output conditions of the electric refrigerator and the absorption refrigerator;
(3.1.4) superposing the electric energy demand generated by the electric refrigerator to the total electric energy consumption amount, and obtaining the electricity purchasing quantity, the wind power and photovoltaic power generation quantity and the cogeneration output condition of the park according to the electric energy comprehensive cost-energy consumption curve and in combination with an electric energy distribution data table;
(3.1.5) superposing the heat energy demand generated by the absorption refrigerator to the total heat energy consumption, and combining a heat energy distribution data table according to a heat energy comprehensive cost-energy consumption curve to obtain the heat purchasing quantity of the park, the cogeneration output condition and the gas boiler output condition;
(3.1.6) judging whether the cogeneration output in the step (3.1.4) and the step (3.1.5) is the same, if not, taking a value with large cogeneration output according to a net cost-output curve of cogeneration if the cogeneration output is in a profit region; if the boiler is in a loss area and the gas boiler is not fully loaded, taking a value with small cogeneration output under the premise of meeting the equipment constraint of the comprehensive energy system; then, setting the cogeneration output as a fixed value, obtaining the park electricity purchasing quantity, the wind power generation quantity and the photovoltaic power generation quantity again according to the electric energy comprehensive cost-energy consumption curve and the electric energy distribution data table, and obtaining the park heat purchasing quantity and the gas boiler output condition according to the heat energy comprehensive cost-energy consumption curve and the heat energy distribution data table;
(3.1.7) superposing the gas energy requirements generated by cogeneration and a gas boiler to the total gas energy consumption, and combining a gas energy distribution data table according to a gas energy comprehensive cost-energy consumption curve to obtain the gas purchasing quantity of the park;
and (3.1.8) calculating the comprehensive cost of the optimized operation of the comprehensive energy system by combining 24 groups of comprehensive cost-energy consumption curves to obtain the optimized operation scheme of the comprehensive energy system.
Step 3) in the step (3.2), the energy storage using modes are firstly divided into an energy storage using mode 1 and an energy storage using mode 2, and the meanings are as follows:
energy storage usage mode 1: storing the energy which can meet the energy consumption requirement in the period a and can be supplied by the optimized area to replace the period b;
energy storage use mode 2: and storing the low-cost external network energy purchasing energy after the c-period optimizable area meets the energy consumption requirement, and replacing the high-cost external network energy purchasing energy needed by the d-period optimizable area.
Step 3) the calculation of the optimal operation scheme of the integrated energy system including the energy storage device in the step (3.2) comprises the following steps:
(3.2.1) calculating the residual energy value of each energy optimizing area and the using energy value of the reserve supply area in 24 periods of cold, heat, electricity and gas according to the optimized operation scheme of the comprehensive energy system obtained in the step (3.1), wherein the calculation formula is as follows:
wherein the content of the first and second substances,
optimizing the zone residual energy value for the ith type energy t period;
the energy using value of the supply area is guaranteed for the ith type energy t period;
outputting an energy value upper limit for the distribution network with the ith type of energy;
is the upper limit of the output value of the energy conversion device w; w is an energy conversion equipment set of the energy conversion I area;
is the upper limit of the output value of the energy supply device v; v is a set of energy supply devices;
the actual energy demand of the ith type of energy in the park; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
(3.2.2) calculating an energy storage use scheme under the energy storage use mode 1, wherein the calculation formula of the energy storage use scheme is as follows:
wherein, y i For the optimization target of the ith type of energy in the energy storage use mode 1, i represents the energy generated by various types of energy sources such as cold, heat, electricity and gas of the comprehensive energy system; t is a scheduling period of 24 h;
the output energy value of the energy storage device t time period of the ith type of energy;
the upper limit of the energy storage/release power of the energy storage device for the ith type of energy in unit time;
judging the energy storage working mode for the ith type energy t time period, if so
Then
Otherwise, the reverse is carried out
M is a set constant;
storing an energy value for the energy storage device for the ith type energy t time period;
storing an energy value for the energy storage device in the ith type energy t-1 time period;
the self-discharge rate of the energy storage device for the ith type of energy;
and
the upper limit and the lower limit of the stored energy value of the energy storage device of the ith type of energy respectively;
(3.2.3) superposing the energy storage use scheme obtained in the step (3.2.2) to the total electric energy consumption, the total heat energy consumption, the total gas energy consumption and the total cold energy consumption in the step (3.1) in the step (3), and recalculating the optimized operation scheme of the comprehensive energy system;
(3.2.4) calculating the secondary residual energy value of each energy optimizing area and the secondary using energy value of the supply area in 24 time intervals, wherein the calculation formula is as follows:
wherein the content of the first and second substances,
the secondary residual energy value of the optimization area is the ith type energy t time period;
ensuring the secondary use energy value of the supply area for the ith type energy t period;
(3.2.5) respectively judging whether the sum of the secondary use energy values of 24 time-keeping supply areas of cold, heat, electricity and gas is 0, if so, performing the step (3.2.6); if not, performing the step (3.2.7);
(3.2.6) calculating the energy storage usage scheme in the energy storage usage mode 2 according to the following formula:
wherein, Δ y i An optimization target of the ith type of energy in the energy storage use mode 2 is obtained;
the correction quantity of the energy value is output for the ith type energy storage device in the t period;
outputting an energy value for the distribution network in the ith type energy t period;
the energy consumption unit cost of the energy storage device for the ith type of energy;
energy cost of ith type of energy in t period;
(3.2.7) the energy storage use scheme in the energy storage use mode 1, the energy storage use scheme in the energy storage use mode 2 and the recalculated comprehensive energy system optimized operation scheme jointly form an optimized operation scheme of the comprehensive energy system comprising the energy storage equipment; at the moment, the comprehensive energy cost calculation formula of the comprehensive energy system is as follows:
wherein F is the comprehensive cost of the optimized operation of the comprehensive energy system; f i The comprehensive cost of the i-th energy is obtained; delta is the iterative part of the calculation resulting from the simultaneous generation of multiple energies by the energy conversion device.
The comprehensive energy system optimization operation method based on the comprehensive cost-energy consumption curve has the following advantages:
1. the invention establishes the integrated model of the comprehensive energy system considering the variable working condition operation characteristics of the multi-type energy conversion equipment, and effectively reduces the energy deviation between the optimized scheduling scheme and the actual load requirement.
2. The comprehensive cost-energy consumption curve oriented to the multi-energy consumption requirement of the park is constructed, and guidance can be provided for scheduling personnel in actual engineering application through analysis and calculation of the curve image guidance optimization operation strategy.
3. The method can effectively give consideration to model precision and solving efficiency, greatly improves the accuracy of the optimization operation scheme, and ensures shorter optimization calculation time.
Drawings
FIG. 1 is a graphical representation of the overall cost-energy curve of the present invention;
FIG. 2 is a schematic diagram of an example of the integrated energy system of the intelligent town park of Xiandan according to the present invention.
Detailed Description
The method for optimizing the operation of the integrated energy system based on the integrated cost-energy consumption curve according to the present invention is described in detail below with reference to the following embodiments and the accompanying drawings.
The invention relates to an optimized operation method of a comprehensive energy system based on a comprehensive cost-energy consumption curve, which comprises the following steps:
1) establishing a park comprehensive energy system integrated model considering variable working condition characteristics of energy conversion equipment; comprises that
(1.1) establishing a variable working condition operation model of energy conversion equipment of a cogeneration unit, a gas boiler, an absorption refrigerator and an electric refrigerator, wherein the variable working condition operation model of the energy conversion equipment comprises the following steps:
wherein, P out For outputting energy, P, to energy-converting apparatus in Energy is input to the energy conversion device, eta is the energy conversion efficiency of the energy conversion device, eta R For the energy conversion efficiency of the energy conversion device under rated operating conditions,
the efficiency correction coefficient is the efficiency correction coefficient of a cogeneration unit or a gas boiler or an absorption refrigerator or an electric refrigerator; wherein, the efficiency correction coefficients of the cogeneration unit, the gas boiler, the absorption refrigerator and the electric refrigerator are respectively as follows:
efficiency correction coefficient of cogeneration unit:
wherein the content of the first and second substances,
the efficiency correction coefficient of the cogeneration unit; k is a radical of CHP,n The fitting coefficient is the n-order fitting coefficient of the electric efficiency of the cogeneration unit; n is a radical of CHP The electric load rate of the cogeneration unit;
the correction coefficient is the heat-electricity ratio of the cogeneration unit; k is a radical of formula α,n Fitting coefficients of an nth order of a thermoelectric ratio of the cogeneration unit;
efficiency correction factor of gas boiler:
wherein, the first and the second end of the pipe are connected with each other,
the efficiency correction coefficient of the gas boiler; k is a radical of GB,n An n-order fitting coefficient for the gas boiler efficiency; n is a radical of GB Is the load factor of the gas boiler;
efficiency correction coefficient of absorption chiller:
wherein, the first and the second end of the pipe are connected with each other,
the efficiency correction factor of the absorption refrigerator; k is a radical of AC,n Fitting coefficients of order n for the efficiency of the absorption chiller; n is a radical of AC Is the load factor of the absorption chiller;
efficiency correction coefficient of electric refrigerator:
wherein the content of the first and second substances,
the efficiency correction factor of the electric refrigerator is obtained; k is a radical of formula EC,n Fitting coefficients of the n orders of the efficiency of the electric refrigerator; n is a radical of EC Is the load factor of the electric refrigerator.
(1.2) establishing an integrated model of the comprehensive energy system as follows:
wherein f is 1 ,f 2 ,…,f n Optimizing an objective function for operation of the integrated energy system; g () is a supply and demand balance expression of the comprehensive energy system; g () is a device constraint expression of the comprehensive energy system; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas; x is the number of i A parameter representing the energy of class i; t is a time scale; xi is an uncertain variable; l is the load demand; p is an output energy value of the energy supply equipment or the energy conversion equipment or the distribution network, namely an output value; h () is a coupling conversion matrix; η is the energy conversion efficiency of the energy conversion device.
2) Constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park; the method comprises the following steps:
(2.1) integrating an objective function of the optimized operation of the comprehensive energy system;
the objective function of the optimization operation of the integrated comprehensive energy system is as follows:
f is the comprehensive cost of the optimized operation of the comprehensive energy system after multi-objective integration by using a membership function; f. of 1 ,f 2 ,…,f n Optimizing an objective function for operation of the integrated energy system; u is the energy conversion device type; f u Represents the net cost of the energy conversion device; s i Representing the energy cost of the ith type of energy; l is a radical of an alcohol i Load demand for class i energy; f i The comprehensive cost of the ith type of energy; delta is a repeated calculation part caused by the energy conversion equipment simultaneously generating multiple energies; omega is an energy conversion equipment set; i represents the ith type of energy in the cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas.
(2.2) constructing a net cost-output curve of the energy conversion equipment;
the net cost-output curve of the energy conversion equipment is a curve drawn according to net cost-output functions of the energy conversion equipment of a cogeneration unit, a gas boiler, an absorption refrigerator and an electric refrigerator; the net cost-out function of the energy conversion device is:
wherein, F u Represents the net cost of the energy conversion device;
for the output value of the energy conversion device u, the superscript ij represents that the energy conversion device u converts the ith type energy into the jth type energy; s i Representing the energy cost of the ith type of energy; s. the j Representing the energy cost of the j-th energy;
the energy conversion efficiency of the energy conversion device u; d u Loss of the energy conversion device u; i and j respectively represent the i and j energy in the cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is hot, E is electricity and G is gas.
(2.3) constructing a comprehensive cost-energy consumption curve facing the multi-element energy consumption requirement of the park;
the comprehensive cost-energy consumption curve for building the diversified energy consumption demand of the park is a curve drawn according to a comprehensive cost-energy consumption function; the overall cost-energy consumption function is:
wherein, F i The comprehensive cost of the i-th energy is obtained; f u Net cost for energy conversion equipment;
is the output value of the energy conversion device u; s i Energy cost for class i energy; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
for the actual energy demand of class i energy in the campus, the numerical values satisfy the following constraints:
wherein L is i Load demand for class i energy;
the energy conversion efficiency of the energy conversion device u;
is the output value of the energy supply device v; v is a set of energy supply devices;
outputting an energy value upper limit for the distribution network with the ith type of energy;
the upper limit of the output value of the energy conversion device u;
is the upper limit of the output value of the energy supply device v.
And (2.4) carrying out regional division on a comprehensive cost-energy consumption curve of the diversified energy consumption requirements of the facing park. The method comprises the following steps of dividing an optimized area and a guarantee supply area, wherein the optimized area comprises: a self-energy supply area, an energy purchasing area and an energy conversion area I, wherein the energy supply area is an energy conversion area II, as shown in figure 1; wherein:
(2.4.1) self-energizing zone: the energy consumption of the park is generated by wind and light renewable energy sources in the park, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.2) energy purchase zone: the energy consumption of the park is met by purchasing energy through an external network;
(2.4.3) energy conversion region I: the energy consumption of the garden is met by energy conversion equipment, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.4) energy conversion zone II: the energy consumption of the park is met by energy conversion equipment, and the unit energy consumption cost is higher than the energy purchasing cost.
3) Calculating an optimized operation scheme of the comprehensive energy system based on the comprehensive cost-energy consumption curve; comprises that
(3.1) calculating an optimized operation scheme of the comprehensive energy system without the energy storage equipment; the method comprises the following steps:
(3.1.1) considering the time-of-use electricity price, respectively constructing 24 groups of comprehensive cost-energy consumption curves corresponding to 24 hours a day according to the comprehensive cost-energy consumption function in the step (2.3) in the step 2), wherein each group of comprehensive cost-energy consumption curves comprises 4 curves of a cold energy comprehensive cost-energy consumption curve, a heat energy comprehensive cost-energy consumption curve, an electric energy comprehensive cost-energy consumption curve and a gas energy comprehensive cost-energy consumption curve;
(3.1.2) according to the 24 groups of comprehensive cost-energy consumption curves, establishing 24 groups of energy distribution data tables containing output and purchased energy of park energy equipment, wherein each group of energy distribution data tables comprises 4 tables of cold energy distribution, heat energy distribution, electric energy distribution and gas energy distribution;
(3.1.3) according to the comprehensive cost-energy consumption curve of the cold energy and by combining a cold energy distribution data table, obtaining the output conditions of the electric refrigerator and the absorption refrigerator;
(3.1.4) superposing the electric energy demand generated by the electric refrigerator to the total electric energy consumption amount, and obtaining the electricity purchasing quantity, the wind power and photovoltaic power generation quantity and the cogeneration output condition of the park according to the electric energy comprehensive cost-energy consumption curve and in combination with an electric energy distribution data table;
(3.1.5) superposing the heat energy demand generated by the absorption refrigerator to the total heat energy consumption, and combining a heat energy distribution data table according to a heat energy comprehensive cost-energy consumption curve to obtain the heat purchasing quantity of the park, the cogeneration output condition and the gas boiler output condition;
(3.1.6) judging whether the cogeneration output in the step (3.1.4) and the step (3.1.5) is the same, if not, taking a value with large cogeneration output according to a net cost-output curve of cogeneration if the cogeneration output is in a profit region; if the gas boiler is in the loss area and the gas boiler is not fully loaded, taking a value with small cogeneration output under the premise of meeting the equipment constraint of the comprehensive energy system; then, setting the cogeneration output as a fixed value, obtaining the park electricity purchasing quantity, the wind power generation quantity and the photovoltaic power generation quantity again according to the electric energy comprehensive cost-energy consumption curve and the electric energy distribution data table, and obtaining the park heat purchasing quantity and the gas boiler output condition according to the heat energy comprehensive cost-energy consumption curve and the heat energy distribution data table;
(3.1.7) superposing the gas energy requirements generated by cogeneration and a gas boiler to the total gas energy consumption, and combining a gas energy distribution data table according to a gas energy comprehensive cost-energy consumption curve to obtain the gas purchasing quantity of the park;
and (3.1.8) calculating the comprehensive cost of the optimized operation of the comprehensive energy system by combining 24 groups of comprehensive cost-energy consumption curves to obtain the optimized operation scheme of the comprehensive energy system.
(3.2) calculating an optimized operation scheme of the comprehensive energy system comprising the energy storage equipment;
firstly, the energy storage using modes are divided into an energy storage using mode 1 and an energy storage using mode 2, and the meanings are as follows:
energy storage usage mode 1: storing the energy which can meet the energy consumption requirement in the period a and can be supplied by the optimized area to replace the period b;
energy storage use mode 2: and storing the low-cost external network energy purchasing energy after the c-period optimizable area meets the energy consumption requirement, and replacing the high-cost external network energy purchasing energy needed by the d-period optimizable area.
The calculation of the optimized operation scheme of the comprehensive energy system comprising the energy storage device specifically comprises the following steps:
(3.2.1) calculating the residual energy value of each energy optimizing area and the using energy value of the reserve supply area in 24 periods of cold, heat, electricity and gas according to the optimized operation scheme of the comprehensive energy system obtained in the step (3.1), wherein the calculation formula is as follows:
wherein the content of the first and second substances,
optimizing the zone residual energy value for the ith type energy t period;
the energy using value of the supply area is guaranteed for the ith type energy t period;
outputting an energy value upper limit for the distribution network with the ith type of energy;
is the upper limit of the output value of the energy conversion device w; w is an energy conversion equipment set of the energy conversion I area;
an upper limit of the output value for the energy supply device v; v is a set of energy supply devices;
the actual energy demand of the ith type of energy in the park; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
(3.2.2) calculating an energy storage use scheme under the energy storage use mode 1, wherein the calculation formula of the energy storage use scheme is as follows:
wherein, y i Aiming at the optimization target of the ith type energy in the energy storage use mode 1 and replacing the energy provided by the reserve supply area as much as possibleWhile using as little energy storage capacity as possible; i represents the ith energy of cold, heat, electricity and gas of the comprehensive energy system; t is a scheduling period of 24 h;
the output energy value of the energy storage device t period of the ith type energy (if the absorbed energy is negative) is obtained;
an upper limit of energy storage/release power of the energy storage device for the ith type of energy in unit time;
judging the energy storage working mode for the ith type energy t time period, if so
Then
Otherwise, the reverse is carried out
M is a set constant;
storing an energy value for the energy storage device for the ith type energy t time period;
storing an energy value for the energy storage device in the ith type energy t-1 time period;
the self-discharging rate of the energy storage device for the ith type of energy;
and with
Stored energy value of energy storage device respectively being i-th energyUpper/lower limits;
(3.2.3) superposing the energy storage use scheme obtained in the step (3.2.2) to the total electric energy consumption, the total heat energy consumption, the total gas energy consumption and the total cold energy consumption in the step (3.1) in the step (3), and recalculating the optimized operation scheme of the comprehensive energy system;
(3.2.4) calculating the secondary residual energy value of each energy optimizing area and the secondary using energy value of the supply area in 24 periods of cold, heat, electricity and gas, wherein the calculation formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,
the secondary residual energy value of the optimization area is the ith type energy t period;
ensuring the secondary use energy value of the supply area for the ith type energy t period;
(3.2.5) respectively judging whether the sum of the secondary use energy values of 24 time-keeping supply areas of cold, heat, electricity and gas is 0, if so, performing the step (3.2.6); if not, performing the step (3.2.7);
(3.2.6) calculating the energy storage usage scheme in the energy storage usage mode 2 according to the following formula:
wherein, Δ y i The optimization target of the ith type energy in the energy storage use mode 2 is to make use of time-of-use energy price to make profit under the condition of ensuring the energy supply of the guaranteed supply area;
outputting the correction quantity of the energy value for the ith type energy storage device in the t period (increasing the released energy to be positive, and increasing the absorbed energy to be negative);
outputting an energy value for the distribution network in the ith type energy t period;
the energy consumption unit cost of the energy storage device for the ith type of energy;
energy cost of ith type of energy in t period;
(3.2.7) the energy storage use scheme in the energy storage use mode 1, the energy storage use scheme in the energy storage use mode 2 and the recalculated comprehensive energy system optimized operation scheme jointly form a comprehensive energy system optimized operation scheme comprising energy storage equipment; at the moment, the comprehensive energy cost calculation formula of the comprehensive energy system is as follows:
wherein F is the comprehensive cost of the optimized operation of the comprehensive energy system; f i The comprehensive cost of the i-th energy is obtained; Δ is the iterative part of the calculation resulting from the simultaneous generation of multiple energies by the energy conversion device.
Examples are given below:
for example, a structure diagram of a small town with a peace of mind is shown in FIG. 2.4 scenes are designed for simulation and comparative analysis. The scene details are as follows:
scene 1: the variable working condition operation characteristics of the multi-type energy equipment are not considered, and a solver is used for solving a CIES optimization operation scheme;
scene 2: the variable working condition operation characteristics of the multi-type energy equipment are considered, the multi-type energy equipment is subjected to piecewise linearization (divided into 3 sections), and a solver is further utilized to solve the collaborative optimization operation scheme of the park comprehensive energy system;
scene 3: the variable working condition operation characteristics of the multi-type energy equipment are considered, the multi-type energy equipment is subjected to piecewise linearization (divided into 15 sections), and a solver is further utilized to solve the collaborative optimization operation scheme of the park comprehensive energy system;
scene 4: and (4) considering the variable working condition operation characteristics of the multi-type energy equipment, and solving the collaborative optimization operation scheme of the park comprehensive energy system by using the method provided by the invention.
Table 1 shows the comparison of simulation results in different scenarios of the campus integrated energy system without energy storage devices; table 2 shows the comparison of simulation results in different scenarios for the campus integrated energy system including the energy storage device. The total energy deviation in the table refers to the sum of the all-day energy differences between the optimized operating schedule and the actual load demand.
TABLE 1
TABLE 2
According to the data in table 1 and table 2, the effect of the proposed method can be clearly shown:
comparing scenes 1-3, it can be seen that the total amount of energy deviation can be remarkably reduced by considering the variable working condition operation characteristics of the multi-type energy conversion equipment, and the accuracy of the collaborative optimization operation scheme of the park comprehensive energy system is improved. The more accurate the variable operating characteristic characterization is, the higher the accuracy of the optimized operating scheme is, but at the same time, the optimization calculation time is increased. The deviation sum of various energies is different, the energy deviation of electric energy and heat energy is caused by the variable working condition operation of refrigeration equipment (an electric refrigerator and an absorption refrigerator), and the gas energy deviation is caused by the variable working condition operation of power generation and heating equipment (CHP and a gas boiler). The daily operation cost of the park comprehensive energy system is reduced after the variable working condition operation characteristics of the equipment are considered, mainly because the energy conversion efficiency of part of energy conversion equipment is higher under the non-rated working condition, such as the refrigeration efficiency of an electric refrigerator.
Comparing table 1 with table 2, it can be seen that the energy storage device can effectively reduce the operating cost of the system. But the cost reduction is not significant due to the smaller energy storage capacity in this example. After the energy storage equipment is considered, the optimization calculation time is obviously increased mainly because the energy storage equipment cannot be optimized through a single time section during operation, and all the time sections need to be planned overall, so that the number of the optimized variables is increased by 1 order of magnitude, and the more the optimized variables are, the slower the solving speed is. The energy deviation sum of the gas energy in the table 2 is smaller than that in the table 1, because the electricity energy storage and heat energy storage devices are used, the electricity generation and the heat generation of the CHP and the gas boiler are reduced, and further, the influence of the variable working condition operation characteristics of the two devices on the total energy consumption is reduced.
(3) Comparing the scene 4 with other scenes, it can be seen that the method provided by the invention can effectively overcome the contradiction between the calculation accuracy and the calculation efficiency of the collaborative optimization operation scheme of the park comprehensive energy system, and ensure shorter optimization calculation time while remarkably reducing the total amount of energy deviation. The method provided by the invention mainly considers the variable working condition operating characteristics of the equipment, and the calculation process only adopts logic judgment and data table query (small-scale linear programming is added when the energy storage equipment is considered). The method of the invention results in slightly higher cost (the daily operation cost is about 0.006 percent), because the heat energy required by refrigeration in actual conditions can be provided by a gas boiler on the basis of the energy purchase price of electric energy and heat energy during cold energy distribution, so that the optimal operation scheme is not globally optimal. In general occasions, the daily operation cost of the scheme is closer to the actual optimal scheme, so that the scheme can be directly used; on occasions with special requirements, the operation scheme can provide an optimized operation reference value with high reliability for other optimized operation strategies, so that the optimization accuracy of the other optimized operation strategies is improved, and the optimization time is shortened.
Claims (10)
1. An integrated energy system optimization operation method based on an integrated cost-energy consumption curve is characterized by comprising the following steps:
1) establishing a park comprehensive energy system integrated model considering variable working condition characteristics of energy conversion equipment; comprises that
(1.1) establishing a variable working condition operation model of energy conversion equipment of a combined heat and power generation unit, a gas boiler, an absorption refrigerator and an electric refrigerator, wherein the variable working condition operation model of the energy conversion equipment is as follows:
wherein, P out For outputting energy, P, to energy-converting apparatus in Energy is input to the energy conversion device, eta is the energy conversion efficiency of the energy conversion device, eta R For the energy conversion efficiency of the energy conversion device under rated operating conditions,
the efficiency correction coefficient is the efficiency correction coefficient of a cogeneration unit or a gas boiler or an absorption refrigerator or an electric refrigerator;
(1.2) establishing an integrated model of the comprehensive energy system;
2) constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park; the method comprises the following steps:
(2.1) integrating an objective function of the optimized operation of the comprehensive energy system;
(2.2) constructing a net cost-output curve of the energy conversion equipment;
(2.3) constructing a comprehensive cost-energy consumption curve facing the multi-energy consumption requirement of the park;
(2.4) carrying out regional division on a comprehensive cost-energy consumption curve of the multi-energy consumption requirement of the facing park;
3) calculating an optimized operation scheme of the comprehensive energy system based on the comprehensive cost-energy consumption curve; comprises that
(3.1) calculating an optimized operation scheme of the comprehensive energy system without the energy storage equipment;
and (3.2) calculating an optimized operation scheme of the integrated energy system comprising the energy storage equipment.
2. The comprehensive energy system optimal operation method based on the comprehensive cost-energy consumption curve according to claim 1, wherein the efficiency correction coefficients of the cogeneration unit, the gas boiler, the absorption chiller and the electric chiller in the step 1) and the step (1.1) are respectively as follows:
efficiency correction coefficient of cogeneration unit:
wherein the content of the first and second substances,
the efficiency correction coefficient of the cogeneration unit; k is a radical of CHP,n An n-order fitting coefficient of the electric efficiency of the cogeneration unit; n is a radical of CHP The electric load rate of the cogeneration unit;
the correction coefficient is the heat-electricity ratio of the cogeneration unit; k is a radical of α,n Fitting coefficients of an nth order of a thermoelectric ratio of the cogeneration unit;
efficiency correction factor of gas boiler:
wherein the content of the first and second substances,
the efficiency correction coefficient of the gas boiler; k is a radical of GB,n An n-order fitting coefficient of the gas boiler efficiency; n is a radical of GB Is the load factor of the gas boiler;
efficiency correction coefficient of absorption chiller:
wherein the content of the first and second substances,
correcting the coefficient for the efficiency of the absorption chiller; k is a radical of AC,n Is an n-order fitting coefficient of the efficiency of the absorption refrigerator; n is a radical of AC Is the load factor of the absorption chiller;
efficiency correction coefficient of electric refrigerator:
wherein the content of the first and second substances,
the efficiency correction factor of the electric refrigerator is obtained; k is a radical of EC,n Fitting coefficients of the electric refrigerator efficiency in the order of n; n is a radical of EC Is the load factor of the electric refrigerator.
3. The method for optimizing the operation of the integrated energy system based on the integrated cost-energy consumption curve according to claim 1, wherein the integrated model of the integrated energy system in the step 1) and the step (1.2) is as follows:
wherein f is 1 ,f 2 ,…,f n An objective function for optimizing operation of the integrated energy system; g () is a supply and demand balance expression of the comprehensive energy system; g () is a device constraint expression of the comprehensive energy system; i represents the ith energy in the cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is hot, E is electricity and G is gas; x is the number of i A parameter representing the energy of class i; t is a time scale; xi is uncertainDetermining variables; l is the load demand; p is an output energy value of the energy supply equipment or the energy conversion equipment or the distribution network, namely an output value; h () is a coupling conversion matrix; η is the energy conversion efficiency of the energy conversion device.
4. The method for optimizing the operation of the integrated energy system based on the integrated cost-energy consumption curve according to claim 1, wherein the objective function of the optimizing operation of the integrated energy system in the step 2) at the step (2.1) is:
f is the comprehensive cost of the optimized operation of the comprehensive energy system after multi-objective integration by using a membership function; f. of 1 ,f 2 ,…,f n An objective function for optimizing operation of the integrated energy system; u is the energy conversion device type; f u Represents the net cost of the energy conversion device; s i Representing the energy cost of the ith type of energy; l is i Load demand for class i energy; f i The comprehensive cost of the i-th energy is obtained; delta is a repeated calculation part caused by the energy conversion equipment simultaneously generating multiple energies; omega is an energy conversion equipment set; i represents the ith type of energy in the comprehensive energy system of cold, heat, electricity and gas, C is cold, H is heat, E is electricity and G is gas.
5. The method according to claim 1, wherein the net cost-output curve of the energy conversion device in step 2) (2.2) is a curve plotted according to a net cost-output function of the energy conversion device for the cogeneration unit, the gas boiler, the absorption chiller, and the electric chiller; the net cost-out function of the energy conversion device is:
wherein, F u Represents the net cost of the energy conversion device;
for the output value of the energy conversion device u, the superscript ij represents that the energy conversion device u converts the ith type energy into the jth type energy; s i Representing the energy cost of the ith type of energy; s j Representing the energy cost of the j-th energy;
the energy conversion efficiency of the energy conversion device u; d u Loss of the energy conversion device u; i and j respectively represent the i and j energy in the cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is hot, E is electricity and G is gas.
6. The integrated energy system optimized operation method based on integrated cost-energy consumption curve according to claim 1, characterized in that the integrated cost-energy consumption curve for building the diversified energy consumption requirements of the park in the step 2) (2.3) is a curve drawn according to the integrated cost-energy consumption function; the overall cost-energy consumption function is:
wherein, F i The comprehensive cost of the i-th energy is obtained; f u Net cost for energy conversion equipment;
is the output value of the energy conversion device u; s i Energy cost for class i energy; i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
for class i energy in the parkThe practical energy demand and the numerical value satisfy the following constraints:
wherein L is i Load demand for class i energy;
the energy conversion efficiency of the energy conversion device u;
is the output value of the energy supply device v; v is a set of energy supply devices;
outputting an energy value upper limit for the distribution network with the ith type of energy;
is the upper limit of the output value of the energy conversion device u;
is the upper limit of the output value of the energy supply device v.
7. The integrated energy system optimized operation method based on integrated cost-energy consumption curve according to claim 1, characterized in that the integrated cost-energy consumption curve of the multi-energy consumption demand of the facing park in the step 2) in the step (2.4) is divided into regions, including an optimizable region and a reserve region, wherein the optimizable region comprises: the energy supply system comprises a self-energy supply area, an energy purchasing area and an energy conversion area I, wherein a supply area is an energy conversion area II; wherein:
(2.4.1) self-energizing zone: the energy consumption of the park is generated by wind and light renewable energy sources in the park, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.2) energy purchase zone: energy consumption of the park is met by purchasing energy through an external network;
(2.4.3) energy conversion region I: the energy consumption of the garden is met by energy conversion equipment, and the unit energy consumption cost is lower than that of energy purchasing;
(2.4.4) energy conversion zone II: the energy consumption of the park is met by energy conversion equipment, and the unit energy consumption cost is higher than the energy purchasing cost.
8. The method for optimized operation of an integrated energy system based on integrated cost-energy consumption curve of claim 1, wherein the step 3) of calculating the optimized operation scheme of the integrated energy system without energy storage device in step (3.1) comprises:
(3.1.1) considering the time-of-use electricity price, respectively constructing 24 groups of comprehensive cost-energy consumption curves corresponding to 24 hours a day according to the comprehensive cost-energy consumption function in the step (2.3) in the step 2), wherein each group of comprehensive cost-energy consumption curves comprises 4 curves of a cold energy comprehensive cost-energy consumption curve, a heat energy comprehensive cost-energy consumption curve, an electric energy comprehensive cost-energy consumption curve and a gas energy comprehensive cost-energy consumption curve;
(3.1.2) according to the 24 groups of comprehensive cost-energy consumption curves, establishing 24 groups of energy distribution data tables containing output and purchased energy of park energy equipment, wherein each group of energy distribution data tables comprises 4 tables of cold energy distribution, heat energy distribution, electric energy distribution and gas energy distribution;
(3.1.3) obtaining the output conditions of the electric refrigerator and the absorption refrigerator according to the comprehensive cold energy cost-energy consumption curve and in combination with a cold energy distribution data table;
(3.1.4) superposing the electric energy demand generated by the electric refrigerator to the total electric energy consumption amount, and obtaining the electricity purchasing quantity, the wind power and photovoltaic power generation quantity and the cogeneration output condition of the park according to the electric energy comprehensive cost-energy consumption curve and in combination with an electric energy distribution data table;
(3.1.5) superposing the heat energy demand generated by the absorption refrigerator to the total heat energy consumption, and obtaining the district heat purchase amount, the cogeneration output condition and the gas boiler output condition according to the heat energy comprehensive cost-energy consumption curve and by combining a heat energy distribution data table;
(3.1.6) judging whether the cogeneration output in the step (3.1.4) and the step (3.1.5) is the same, if not, taking a large value of the cogeneration output according to a net cost-output curve of the cogeneration if the cogeneration output is in a profit region; if the gas boiler is in the loss area and the gas boiler is not fully loaded, taking a value with small cogeneration output under the premise of meeting the equipment constraint of the comprehensive energy system; then, setting the cogeneration output as a fixed value, obtaining the park electricity purchasing quantity, the wind power generation quantity and the photovoltaic power generation quantity according to the electric energy comprehensive cost-energy consumption curve and the electric energy distribution data table again, and obtaining the park heat purchasing quantity and the gas boiler output condition according to the heat energy comprehensive cost-energy consumption curve and the heat energy distribution data table;
(3.1.7) superposing the gas energy requirements generated by cogeneration and a gas boiler to the total gas energy consumption, and combining a gas energy distribution data table according to a gas energy comprehensive cost-energy consumption curve to obtain the gas purchasing quantity of the park;
and (3.1.8) calculating the comprehensive cost of the optimized operation of the comprehensive energy system by combining 24 groups of comprehensive cost-energy consumption curves to obtain the optimized operation scheme of the comprehensive energy system.
9. The comprehensive energy system optimized operation method based on the comprehensive cost-energy consumption curve as claimed in claim 1, wherein the step 3) and the step (3.2) firstly divide the energy storage usage modes into two types, namely an energy storage usage mode 1 and an energy storage usage mode 2, and the meanings are as follows:
energy storage usage pattern 1: storing the energy which can meet the energy consumption requirement in the period a and can be supplied by the optimized area to replace the period b;
energy storage use mode 2: and after the c-period optimized area is stored, the low-cost external network energy purchasing energy required by the energy consumption requirement is met, and the high-cost external network energy purchasing energy required by the d-period optimized area is replaced.
10. The integrated energy system optimal operation method based on the integrated cost-energy consumption curve according to claim 9, wherein the step 3) of calculating the integrated energy system optimal operation scheme comprising the energy storage device in the step (3.2) comprises:
(3.2.1) calculating the residual energy value of each energy optimizing area and the using energy value of the reserve supply area in 24 periods of cold, heat, electricity and gas according to the optimized operation scheme of the comprehensive energy system obtained in the step (3.1), wherein the calculation formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,
optimizing the zone residual energy value for the ith type energy t period;
the energy using value of the supply area is guaranteed for the ith type energy t period;
outputting an energy value upper limit for the distribution network with the ith type of energy;
is the upper limit of the output value of the energy conversion device w; w is an energy conversion equipment set of an energy conversion I area;
an upper limit of the output value for the energy supply device v; v is a set of energy supply devices;
for the actual energy demand of class i energy in the park(ii) a i represents the ith type of energy in cold, heat, electricity and gas of the comprehensive energy system, C is cold, H is heat, E is electricity and G is gas;
(3.2.2) calculating an energy storage use scheme under the energy storage use mode 1, wherein the calculation formula of the energy storage use scheme is as follows:
wherein, y i For the optimization target of the ith type of energy in the energy storage use mode 1, i represents the energy generated by various types of energy such as cold, heat, electricity and gas of the comprehensive energy system; t is a scheduling period of 24 h;
the output energy value of the energy storage device t time period of the ith type of energy;
the upper limit of the energy storage/release power of the energy storage device for the ith type of energy in unit time;
judging the energy storage working mode for the ith type energy t time period, if so
Then
Otherwise, the reverse is carried out
M is a set constant;
storing an energy value for the energy storage device for the ith type energy t time period;
storing an energy value for the energy storage device in the ith type energy t-1 time period;
the self-discharging rate of the energy storage device for the ith type of energy;
and
the upper limit and the lower limit of the stored energy value of the energy storage device for the ith type of energy are respectively set;
(3.2.3) superposing the energy storage use scheme obtained in the step (3.2.2) to the total electric energy consumption amount, the total heat energy consumption amount, the total gas energy consumption amount and the total cold energy consumption amount in the step (3.1) in the step (3), and recalculating the optimized operation scheme of the comprehensive energy system;
(3.2.4) calculating the secondary residual energy value of each energy optimizing area and the secondary using energy value of the supply area in 24 periods of cold, heat, electricity and gas, wherein the calculation formula is as follows:
wherein the content of the first and second substances,
the secondary residual energy value of the optimization area is the ith type energy t period;
ensuring the secondary use energy value of the supply area for the ith type energy t period;
(3.2.5) respectively judging whether the sum of the secondary use energy values of 24 time-keeping supply areas of cold, heat, electricity and gas is 0, if so, performing the step (3.2.6); if not, performing the step (3.2.7);
(3.2.6) calculating the energy storage usage scheme in the energy storage usage mode 2 according to the following formula:
wherein, Δ y i An optimization target of the ith type of energy in the energy storage use mode 2 is obtained;
the correction quantity of the energy value is output for the ith type energy storage device in the t period;
outputting an energy value for the distribution network in the ith type energy t period;
the energy consumption unit cost of the energy storage device for the ith type of energy;
energy cost of ith type of energy in t period;
(3.2.7) the energy storage use scheme in the energy storage use mode 1, the energy storage use scheme in the energy storage use mode 2 and the recalculated comprehensive energy system optimized operation scheme jointly form an optimized operation scheme of the comprehensive energy system comprising the energy storage equipment; at the moment, the comprehensive energy cost calculation formula of the comprehensive energy system is as follows:
wherein F is the comprehensive cost of the optimized operation of the comprehensive energy system; f i The comprehensive cost of the i-th energy is obtained; delta is energy conversion equipmentA plurality of energy-induced iterative computation portions are generated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210608295.5A CN115099007B (en) | 2022-05-31 | 2022-05-31 | Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210608295.5A CN115099007B (en) | 2022-05-31 | 2022-05-31 | Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115099007A true CN115099007A (en) | 2022-09-23 |
| CN115099007B CN115099007B (en) | 2023-04-18 |
Family
ID=83288673
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210608295.5A Active CN115099007B (en) | 2022-05-31 | 2022-05-31 | Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115099007B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115599963A (en) * | 2022-11-29 | 2023-01-13 | 深圳大学（Cn） | Shallow geothermal efficient utilization and storage system and method based on deep learning optimization |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020237700A1 (en) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | Operation scheduling method for multiple energy systems |
| CN113112087A (en) * | 2021-04-23 | 2021-07-13 | 国网宁夏电力有限公司经济技术研究院 | Comprehensive energy system operation cost optimization method considering electric heating load demand response |
| CN113487188A (en) * | 2021-07-08 | 2021-10-08 | 重庆理工大学 | Comprehensive energy system optimal scheduling method considering electric and gas joint price guide mechanism |
-
2022
- 2022-05-31 CN CN202210608295.5A patent/CN115099007B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020237700A1 (en) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | Operation scheduling method for multiple energy systems |
| CN113112087A (en) * | 2021-04-23 | 2021-07-13 | 国网宁夏电力有限公司经济技术研究院 | Comprehensive energy system operation cost optimization method considering electric heating load demand response |
| CN113487188A (en) * | 2021-07-08 | 2021-10-08 | 重庆理工大学 | Comprehensive energy system optimal scheduling method considering electric and gas joint price guide mechanism |
Non-Patent Citations (4)
| Title |
|---|
| 李建林等: "考虑变工况特性的微能源系统优化规划 (一)基本模型和分析" * |
| 王成山等: "园区型综合能源系统多时间尺度模型预测优化调度" * |
| 田立亭等: "考虑变工况特性的微能源系统优化规划 (二)基本模型和分析" * |
| 马国真等: "电力系统保护与控制" * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115599963A (en) * | 2022-11-29 | 2023-01-13 | 深圳大学（Cn） | Shallow geothermal efficient utilization and storage system and method based on deep learning optimization |
| CN115599963B (en) * | 2022-11-29 | 2023-04-07 | 深圳大学 | Shallow geothermal efficient utilization and storage system and method based on deep learning optimization |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115099007B (en) | 2023-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108229025B (en) | Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system | |
| CN109919478B (en) | Comprehensive energy microgrid planning method considering comprehensive energy supply reliability | |
| CN113344736B (en) | Park-level comprehensive energy system and control method thereof | |
| CN111950807B (en) | Comprehensive energy system optimization operation method considering uncertainty and demand response | |
| CN113255992A (en) | Optimal scheduling method and device for park-level multi-energy complementary system | |
| CN110661254A (en) | Method for quantifying cooling, heating and power complementary benefits of regional comprehensive energy system | |
| CN111737884B (en) | Multi-target random planning method for micro-energy network containing multiple clean energy sources | |
| CN113779783B (en) | Multi-uncertainty-considered planning and operation joint optimization method for regional comprehensive energy system | |
| CN113987734A (en) | Multi-objective optimization scheduling method for park comprehensive energy system under opportunity constraint condition | |
| CN116542370A (en) | Park low-carbon economic operation method considering carbon capture and carbon transaction | |
| CN111724045B (en) | Comprehensive energy system energy efficiency evaluation and improvement method based on data driving | |
| CN111639819B (en) | Multi-stage optimization control method for comprehensive energy park | |
| Fan et al. | Two-layer collaborative optimization for a renewable energy system combining electricity storage, hydrogen storage, and heat storage | |
| CN110994606A (en) | Multi-energy power supply capacity configuration method based on complex adaptive system theory | |
| CN112836882A (en) | Regional comprehensive energy system operation optimization method considering equipment load rate change | |
| Song et al. | A fuzzy‐based multi‐objective robust optimization model for a regional hybrid energy system considering uncertainty | |
| CN115170343A (en) | Distributed resource and energy storage collaborative planning method for regional comprehensive energy system | |
| CN115099007B (en) | Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve | |
| CN113158547B (en) | Regional comprehensive energy system optimal configuration method considering economy and reliability | |
| CN111126675A (en) | Multi-energy complementary microgrid system optimization method | |
| CN112883630B (en) | Multi-microgrid system day-ahead optimization economic dispatching method for wind power consumption | |
| CN117081143A (en) | Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion | |
| CN114139837B (en) | Regional multi-system double-layer distributed optimization scheduling method considering double-layer carbon emission optimization distribution model | |
| CN112950410B (en) | Energy hub system planning method considering wind-solar correlation and preventive maintenance | |
| CN213783243U (en) | Comprehensive energy system operation optimizing device for industrial park |
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 |