CN113469430A - Multi-energy complementary capacity configuration method of integrated energy park - Google Patents

Multi-energy complementary capacity configuration method of integrated energy park Download PDF

Info

Publication number
CN113469430A
CN113469430A CN202110725103.4A CN202110725103A CN113469430A CN 113469430 A CN113469430 A CN 113469430A CN 202110725103 A CN202110725103 A CN 202110725103A CN 113469430 A CN113469430 A CN 113469430A
Authority
CN
China
Prior art keywords
power
electric
natural gas
output
heat
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
Application number
CN202110725103.4A
Other languages
Chinese (zh)
Other versions
CN113469430B (en
Inventor
周保中
范炜
刘敦楠
张继广
李忆
毕圣
吴思翰
周畅游
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Huadian Electric Power Research Institute Co Ltd
Original Assignee
North China Electric Power University
Huadian Electric Power Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University, Huadian Electric Power Research Institute Co Ltd filed Critical North China Electric Power University
Priority to CN202110725103.4A priority Critical patent/CN113469430B/en
Publication of CN113469430A publication Critical patent/CN113469430A/en
Application granted granted Critical
Publication of CN113469430B publication Critical patent/CN113469430B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/067—Enterprise or organisation modelling
    • 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

Abstract

The invention provides a multi-energy complementary capacity configuration method of a comprehensive energy park, which comprises the following steps: establishing output models of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple power supply system and electric heating and refrigerating equipment in comprehensive energy supply; analyzing typical energy utilization scenes and loads of the comprehensive energy utilization park, wherein the typical energy utilization scenes and loads comprise electric loads, heat loads and cold loads in hot seasons, cold seasons and transition seasons; the method comprises the steps of establishing a comprehensive energy park multi-energy complementary capacity configuration optimization model by taking the lowest total operation cost of a system as an optimization target and taking power supply system electric power balance constraint, natural gas cooling-heat-electric generator set constraint, electricity purchasing power constraint, wind power photovoltaic output constraint, heat supply system output constraint, start-stop time constraint, conveying channel capacity constraint and reserve capacity constraint as constraint conditions, and solving the model by adopting a genetic algorithm. The invention can optimize the multi-energy complementary capacity configuration of the comprehensive energy utilization park, realize the integration of multiple energy sources and improve the economy and reliability of the park system.

Description

Multi-energy complementary capacity configuration method of integrated energy park
Technical Field
The invention relates to the technical field of power economy, in particular to a multi-energy complementary capacity configuration method of an integrated energy park.
Background
With the increasing of the development proportion of new energy, the basic characteristics of the traditional energy system are gradually changed, the comprehensive energy utilization park breaks the boundary of the original different energy systems, and the integration of various types of energy is realized. The randomness, the indirection, the space-time complementarity and other characteristics of resources such as wind, light, water, fire and the like bring challenges to the safe and stable operation of a power system, and the reasonable configuration of the capacity of each energy source is the key for ensuring the feasibility of the comprehensive energy park, thereby being beneficial to improving the power supply reliability of the comprehensive energy system. The reasonable distribution of the capacity of each power supply in the multi-energy complementary power generation system is fully utilized, the advantages of each power supply are exerted, the consumption of clean energy can be increased, the total investment of the system is reduced, and the economic benefit is improved. At present, the capacity configuration of the comprehensive energy utilization park in China only stays in the multi-energy complementary capacity configuration of primary energy sources such as wind, light, water, fire and the like, the multi-energy complementary capacity configuration of secondary energy sources such as cold, heat, electricity and the like is not considered, and the configuration of electric loads, heat loads, cold loads and the like on the user side in different seasons is not involved.
Therefore, there is a need to develop a method for configuring and optimizing the multi-energy complementary capacity of the integrated energy park, so as to integrate the multi-energy sources and improve the economic efficiency and reliability of the system of the integrated energy park.
Disclosure of Invention
The invention aims to provide a multi-energy complementary capacity configuration method for an integrated energy utilization park, which is used for modeling the output of distributed wind power, distributed photovoltaic, cold-heat-electricity triple supply systems and electric heating and electric refrigerating equipment of an integrated energy system and analyzing typical energy utilization scenes and loads of the integrated energy utilization park, so that multi-energy complementary capacity configuration optimization is performed on the integrated energy utilization park, integration of multiple energy sources is realized, and the economy and reliability of the integrated energy utilization park system are improved.
The invention provides a multi-energy complementary capacity configuration method of a comprehensive energy utilization park, which comprises the following steps:
s1, establishing an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple power supply system and electric heating and electric refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy utilization scenes and loads of the integrated energy utilization park, including electric loads, heat loads and cold loads in hot seasons, cold seasons and transition seasons;
s3, with the lowest total operation cost of the system as an optimization target, and with the constraint of electric power balance of a power supply system, the constraint of a natural gas cooling-heat-electric generator set, the constraint of electricity purchasing power, the constraint of wind power photovoltaic output, the constraint of heat supply system output, the constraint of start-stop time, the constraint of capacity of a conveying channel and the constraint of reserve capacity as constraint conditions, establishing a comprehensive energy use park multi-energy complementary capacity configuration optimization model, and solving the model by adopting a genetic algorithm.
Further, the step S1 includes the following steps:
the output of the fan is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
Figure BDA0003138290860000021
wherein the content of the first and second substances,
Figure BDA0003138290860000022
is the maximum output electric power v of the wind turbine generator at time ttIs the wind speed, v, over a period of time tinIs the cut-in wind speed, vratedIs rated wind speed, voutIs the cut-out wind speed, PS WTOutputting rated output electric power for the wind turbine generator;
secondly, the principle of photovoltaic power generation is that according to the photovoltaic effect, solar energy is converted into electric energy by utilizing a photovoltaic module. The output of the photovoltaic unit is related to the solar radiation intensity and the temperature, and the output model is as follows:
Figure BDA0003138290860000023
wherein the content of the first and second substances,
Figure BDA0003138290860000024
is the maximum output electric power of the photovoltaic generator set in the time period t,
Figure BDA0003138290860000025
is the output electric power of the photovoltaic unit under standard conditions, IstcIs the standard condition solar radiation intensity, which is generally 1000W/m2，TstcIs a standard conditioned temperature, typically 25 ℃, ItIs the intensity of solar radiation, T, of the time period TtIs the temperature of the time period t, and alpha is the power temperature coefficient of the photovoltaic cell, generally 0.0039 DEG C-1；
The natural gas cold-heat-electricity triple power supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, can supply cold, heat and electricity simultaneously, and has the following output model:
Figure BDA0003138290860000026
Figure BDA0003138290860000027
Figure BDA0003138290860000028
Figure BDA0003138290860000029
wherein the content of the first and second substances,
Figure BDA00031382908600000210
is the residual heat power of the natural gas cold-heat-electricity triple supply system in the time period t,
Figure BDA00031382908600000211
is the output electric power, eta, of the natural gas cold-heat-electricity combined supply system in the time period tMT,eIs the efficiency of electric power generation, ηHLIs a heat dissipation loss coefficient, generally 0.03,
Figure BDA00031382908600000212
is the output thermal power, eta, of the natural gas cooling-heating-power combined supply system in the time period tWHRIs the waste heat recovery efficiency, OHIs a heat generation coefficient, generally 1.20,
Figure BDA0003138290860000031
is the output cold power O of the natural gas cold-heat-electricity triple supply system in the time period tCIs a refrigeration coefficient, generally 0.95, VNGIs the consumption of natural gas, qNGIs the low heating value of natural gas;
electric energy belongs to high-quality energy, heat energy can be produced by electric heating equipment such as a heat pump and the like, cold energy can also be produced by the electric refrigerating equipment, and the output model is as follows:
Figure BDA0003138290860000032
Figure BDA0003138290860000033
wherein the content of the first and second substances,
Figure BDA0003138290860000034
is the output thermal power, η, of the electric heating apparatus over a period of time tEHIt is the efficiency of the electric-to-heat conversion,
Figure BDA0003138290860000035
is the input power rate of the electric heating device for time period t;
Figure BDA0003138290860000036
is the output cold power, eta, of the electric refrigerating apparatus in the time interval tECIt is the electric-to-cold conversion efficiency,
Figure BDA0003138290860000037
is the input power rate of the electric refrigeration appliance for time period t.
Further, the step S2 includes:
calculating the electric load, the heat load and the cold load in the hot season, wherein the calculation formula is as follows:
Figure BDA0003138290860000038
wherein the content of the first and second substances,
Figure BDA0003138290860000039
is the electrical load for the heating season t,
Figure BDA00031382908600000310
is the heat load for the heating season t,
Figure BDA00031382908600000311
is the cooling load of the heating season t;
secondly, calculating the electric load, the heat load and the cold load in the cooling season, wherein the calculation formula is as follows:
Figure BDA00031382908600000312
wherein the content of the first and second substances,
Figure BDA00031382908600000313
is the electrical load for the cooling season period t,
Figure BDA00031382908600000314
is the thermal load for the cooling season period t,
Figure BDA00031382908600000315
is the cooling load for the cooling season period t;
thirdly, calculating the electric load, the heat load and the cold load in the transition season, wherein the calculation formula is as follows:
Figure BDA00031382908600000316
wherein the content of the first and second substances,
Figure BDA00031382908600000317
is the electrical load for the transitional season period t,
Figure BDA00031382908600000318
is the thermal load for the transitional season period t,
Figure BDA00031382908600000319
is the cooling load for the transitional season period t.
Further, the comprehensive energy utilization park multipotency complementary capacity allocation optimization model in the step S3 is as follows:
and (3) establishing an objective function by taking the lowest total operation cost of the system as an optimization objective:
Figure BDA00031382908600000320
Figure BDA0003138290860000041
Figure BDA0003138290860000042
wherein F is the total system operation cost which comprises the natural gas cooling-heating-electric machine set power generation cost, the penalty cost of deviating from the planned output, the power purchasing cost, the penalty cost of deviating from the planned power purchasing and the wind and light abandoning cost,NrMTis the number of natural gas cooling-heating-electric units in the r region,
Figure BDA0003138290860000043
is the electric output of the ith natural gas cooling-heating-electric machine set at the moment t of an r area under a scene s,
Figure BDA0003138290860000044
the planned electric output of the ith natural gas cooling-heating-electric machine set at the moment t of the r area,
Figure BDA0003138290860000045
is the thermal output of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure BDA0003138290860000046
is the electricity purchasing power of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure BDA0003138290860000047
the planned purchase power of the ith natural gas cooling-heating-electric machine set at the moment t in the r area,
Figure BDA0003138290860000048
is a cost function of the natural gas cooling-heating-electric machine set, alpha0,r,iFor the initial investment cost of the natural gas cooling-heating-electric machine, alpha1,r,i、α3,r,iRespectively a first term coefficient and a second term coefficient of the electric output in the cost function of the natural gas cold-heat-electric machine set, alpha2,r,i、α4,r,iPrimary and secondary coefficient of thermal output in the cost function of natural gas cold-heat-electric machine, alpha5,r,iIn order to be a cross-term coefficient,alpha is measured by historical data and experience,
Figure BDA0003138290860000049
the method is characterized in that the wind and the light power are abandoned,
Figure BDA00031382908600000410
is the power that can be generated by the wind power,
Figure BDA00031382908600000411
is the amount of power that can be generated by the photovoltaic,
Figure BDA00031382908600000412
is the power of the wind power on the internet,
Figure BDA00031382908600000413
is the power of the photovoltaic grid-connected system,
Figure BDA00031382908600000414
is the heat supply power of the wind power,
Figure BDA00031382908600000415
is the heating power of the photovoltaic system,
Figure BDA00031382908600000416
is the price of electricity purchase, ρ1Is the punishment price of the deviation of the natural gas cooling-heating-electric machine group from the planned output, rho2Is the penalty price, rho, of the deviation of the power purchase from the planned power purchase3Is the punishment price of wind and light abandonment.
Further, the constraint conditions in step S3 are:
electric power balance constraint of the power supply system:
Figure BDA00031382908600000417
wherein L isrIs a set of links between the r region and other regions,
Figure BDA00031382908600000418
is the transmission power of the tie-line l,
Figure BDA00031382908600000419
indicating that the tie i is inputting power into the area at time t,
Figure BDA00031382908600000420
indicating that at time t the tie i delivers power outside this area,
Figure BDA00031382908600000421
is the power of the electrical load and,
Figure BDA00031382908600000422
is the electric refrigerator power;
thermoelectric ratio constraint, upper and lower output limits constraint and climbing constraint of the natural gas cooling-heating-electric machine set:
Figure BDA0003138290860000051
Figure BDA0003138290860000052
Figure BDA0003138290860000053
wherein, Khp,r,tIs the thermoelectric ratio of the natural gas cooling-heating-electric machine set of the ith station in the r area,
Figure BDA0003138290860000054
is the lower limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure BDA0003138290860000055
is the upper limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure BDA0003138290860000056
the maximum up-regulated power in unit time of the ith natural gas cooling-heating-electric machine set in the r area;
the power of purchasing electricity is restricted, each functional area only allows purchasing electricity and gas, but not selling electricity and gas:
Figure BDA0003138290860000057
Figure BDA0003138290860000058
wherein the content of the first and second substances,
Figure BDA0003138290860000059
is the capacity of the substation at the point of connection,
Figure BDA00031382908600000510
is the power of the external gas purchase,
Figure BDA00031382908600000511
is the upper limit of the natural gas pipeline delivery capacity;
wind power photovoltaic output restraint:
Figure BDA00031382908600000512
Figure BDA00031382908600000513
Figure BDA00031382908600000514
Figure BDA00031382908600000515
Figure BDA00031382908600000516
Figure BDA00031382908600000517
the heat supply system is restricted, the electric boiler in the r area only consumes the wind power and the photovoltaic power in the r area, and the electric heat conversion relation restriction of the electric boiler and the output restriction of the electric boiler are met:
Figure BDA00031382908600000518
Figure BDA00031382908600000519
Figure BDA00031382908600000520
wherein the content of the first and second substances,
Figure BDA00031382908600000521
is the thermal output, eta, of the electric boiler at the time of r region t under the scene sEB,rIs the electrical conversion efficiency of the r-zone electric boiler,
Figure BDA00031382908600000522
is the upper limit of the output of the r area electric boiler;
and (3) constraint of start-stop time of the equipment:
Figure BDA0003138290860000061
Figure BDA0003138290860000062
Figure BDA0003138290860000063
Figure BDA0003138290860000064
wherein the content of the first and second substances,
Figure BDA0003138290860000065
is the continuous operation time of the natural gas cooling-heating-electric generating set,
Figure BDA0003138290860000066
is the continuous downtime of the natural gas cooling-heating-electric machine set,
Figure BDA0003138290860000067
is the duration of the operation of the electric boiler plant,
Figure BDA0003138290860000068
is the continuous down time of the electric boiler plant,
Figure BDA0003138290860000069
is the minimum start-up time of natural gas cooling-heating-power,
Figure BDA00031382908600000610
is the natural gas cold-heat-electricity minimum down time,
Figure BDA00031382908600000611
respectively the minimum start-up time of the electric boiler plant,
Figure BDA00031382908600000612
minimum down time of the electric boiler plant, uMT,r,t-1Is the operation state of the natural gas cooling-heating-electric machine set at the time t-1, and is a variable of 0-1 uEB,r,t-1The operation state of the electric boiler equipment at the time t-1 is a variable of 0-1;
the transfer channel capacity satisfies the constraint:
Figure BDA00031382908600000613
Figure BDA00031382908600000614
wherein the content of the first and second substances,
Figure BDA00031382908600000615
is the line capacity between nodes i and j at the initial time,
Figure BDA00031382908600000616
is the newly added line capacity between nodes i and j at time t,
Figure BDA00031382908600000617
is the maximum capacity of the line allowed to be erected in the power transmission channel between the nodes i and j, and M is a penalty coefficient, so that
Figure BDA00031382908600000618
When the temperature of the water is higher than the set temperature,
Figure BDA00031382908600000619
is 0, when
Figure BDA00031382908600000620
Large enough to have enough room to select a new line,
Figure BDA00031382908600000621
is a boolean variable, continuously constrained "switch", that produces a maximum capacity greater than or equal to that allowed by the line upgrade.
And (4) constraint of spare capacity:
Figure BDA00031382908600000622
Figure BDA00031382908600000623
wherein r is1Is the upper spare coefficient of the load, r2Is the upper spare coefficient of wind power, r3Is the upper standby coefficient of the photovoltaic, r4Is the lower standby coefficient of wind power, r5Is the lower standby factor for the photovoltaic.
The invention has the advantages and positive effects that:
compared with the prior art, the invention fully considers the multi-energy complementary capacity configuration of secondary energy sources such as cold, heat, electricity and the like, also relates to the configuration of electric loads, heat loads, cold loads and the like at the user side in different seasons, fully utilizes the reasonable distribution of the capacity of each power supply in the multi-energy complementary power generation system, exerts the advantages of each power supply, realizes the integration of multi-energy sources, can increase the consumption of clean energy sources, reduces the total investment of the system and improves the economy and reliability of the comprehensive energy park system.
Drawings
Fig. 1 is a flowchart of a method for configuring a multi-energy complementary capacity of an integrated energy campus, according to an embodiment of the present invention;
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in more detail below with reference to the accompanying drawings in the embodiments of the present invention. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are only some, but not all embodiments of the invention. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
As shown in fig. 1, the method for configuring the multi-energy complementary capacity of the integrated energy consumption park provided in this embodiment includes the following steps:
s1, establishing an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple power supply system and electric heating and electric refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy utilization scenes and loads of the integrated energy utilization park, including electric loads, heat loads and cold loads in hot seasons, cold seasons and transition seasons;
s3, with the lowest total operation cost of the system as an optimization target, and with the constraint of electric power balance of a power supply system, the constraint of a natural gas cooling-heat-electric generator set, the constraint of electricity purchasing power, the constraint of wind power photovoltaic output, the constraint of heat supply system output, the constraint of start-stop time, the constraint of capacity of a conveying channel and the constraint of reserve capacity as constraint conditions, establishing a comprehensive energy use park multi-energy complementary capacity configuration optimization model, and solving the model by adopting a genetic algorithm.
Further, the step S1 includes the following steps:
the output of the fan is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
Figure BDA0003138290860000071
wherein the content of the first and second substances,
Figure BDA0003138290860000081
is the maximum output electric power v of the wind turbine generator at time ttIs the wind speed, v, over a period of time tinIs the cut-in wind speed, vratedIs rated wind speed, voutIs the cut-out wind speed, PS WTOutputting rated output electric power for the wind turbine generator;
secondly, the principle of photovoltaic power generation is that according to the photovoltaic effect, solar energy is converted into electric energy by utilizing a photovoltaic module. The output of the photovoltaic unit is related to the solar radiation intensity and the temperature, and the output model is as follows:
Figure BDA0003138290860000082
wherein the content of the first and second substances,
Figure BDA0003138290860000083
is the maximum output electric power of the photovoltaic generator set in the time period t,
Figure BDA0003138290860000084
is the output electric power of the photovoltaic unit under standard conditions, IstcIs the standard condition solar radiation intensity, which is generally 1000W/m2，TstcIs a standard conditioned temperature, typically 25 ℃, ItIs the intensity of solar radiation, T, of the time period TtIs the temperature of the time period t, and alpha is the power temperature coefficient of the photovoltaic cell, generally 0.0039 DEG C-1；
The natural gas cold-heat-electricity triple power supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, can supply cold, heat and electricity simultaneously, and has the following output model:
Figure BDA0003138290860000085
Figure BDA0003138290860000086
Figure BDA0003138290860000087
Figure BDA0003138290860000088
wherein the content of the first and second substances,
Figure BDA0003138290860000089
is time period t natural gas cold-heat-electricityThe waste heat power of the combined supply system,
Figure BDA00031382908600000810
is the output electric power, eta, of the natural gas cold-heat-electricity combined supply system in the time period tMT,eIs the efficiency of electric power generation, ηHLIs a heat dissipation loss coefficient, generally 0.03,
Figure BDA00031382908600000811
is the output thermal power, eta, of the natural gas cooling-heating-power combined supply system in the time period tWHRIs the waste heat recovery efficiency, OHIs a heat generation coefficient, generally 1.20,
Figure BDA00031382908600000812
is the output cold power O of the natural gas cold-heat-electricity triple supply system in the time period tCIs a refrigeration coefficient, generally 0.95, VNGIs the consumption of natural gas, qNGIs the low heating value of natural gas;
electric energy belongs to high-quality energy, heat energy can be produced by electric heating equipment such as a heat pump and the like, cold energy can also be produced by the electric refrigerating equipment, and the output model is as follows:
Figure BDA00031382908600000813
Figure BDA00031382908600000814
wherein the content of the first and second substances,
Figure BDA0003138290860000091
is the output thermal power, η, of the electric heating apparatus over a period of time tEHIt is the efficiency of the electric-to-heat conversion,
Figure BDA0003138290860000092
is the input power rate of the electric heating device for time period t;
Figure BDA0003138290860000093
is the output cold power, eta, of the electric refrigerating apparatus in the time interval tECIt is the electric-to-cold conversion efficiency,
Figure BDA0003138290860000094
is the input power rate of the electric refrigeration appliance for time period t.
Further, the step S2 includes:
calculating the electric load, the heat load and the cold load in the hot season, wherein the calculation formula is as follows:
Figure BDA0003138290860000095
wherein the content of the first and second substances,
Figure BDA0003138290860000096
is the electrical load for the heating season t,
Figure BDA0003138290860000097
is the heat load for the heating season t,
Figure BDA0003138290860000098
is the cooling load of the heating season t;
secondly, calculating the electric load, the heat load and the cold load in the cooling season, wherein the calculation formula is as follows:
Figure BDA0003138290860000099
wherein the content of the first and second substances,
Figure BDA00031382908600000910
is the electrical load for the cooling season period t,
Figure BDA00031382908600000911
is the thermal load for the cooling season period t,
Figure BDA00031382908600000912
is the cooling load for the cooling season period t;
thirdly, calculating the electric load, the heat load and the cold load in the transition season, wherein the calculation formula is as follows:
Figure BDA00031382908600000913
wherein the content of the first and second substances,
Figure BDA00031382908600000914
is the electrical load for the transitional season period t,
Figure BDA00031382908600000915
is the thermal load for the transitional season period t,
Figure BDA00031382908600000916
is the cooling load for the transitional season period t.
Further, the comprehensive energy utilization park multipotency complementary capacity allocation optimization model in the step S3 is as follows:
and (3) establishing an objective function by taking the lowest total operation cost of the system as an optimization objective:
Figure BDA00031382908600000917
Figure BDA00031382908600000918
Figure BDA00031382908600000919
wherein F is the total system operation cost which comprises the natural gas cooling-heating-electric machine set power generation cost, the penalty cost of deviating from the planned output, the power purchasing cost, the penalty cost of deviating from the planned power purchasing and the wind and light abandoning cost,NrMTis the number of natural gas cooling-heating-electric units in the r region,
Figure BDA00031382908600000920
is the electric output of the ith natural gas cooling-heating-electric machine set at the moment t of an r area under a scene s,
Figure BDA00031382908600000921
the planned electric output of the ith natural gas cooling-heating-electric machine set at the moment t of the r area,
Figure BDA0003138290860000101
is the thermal output of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure BDA0003138290860000102
is the electricity purchasing power of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure BDA0003138290860000103
the planned purchase power of the ith natural gas cooling-heating-electric machine set at the moment t in the r area,
Figure BDA0003138290860000104
is a cost function of the natural gas cooling-heating-electric machine set, alpha0,r,iFor the initial investment cost of the natural gas cooling-heating-electric machine, alpha1,r,i、α3,r,iRespectively a first term coefficient and a second term coefficient of the electric output in the cost function of the natural gas cold-heat-electric machine set, alpha2,r,i、α4,r,iPrimary and secondary coefficient of thermal output in the cost function of natural gas cold-heat-electric machine, alpha5,r,iAlpha is a cross term coefficient and is obtained by historical data and empirical calculation,
Figure BDA0003138290860000105
the method is characterized in that the wind and the light power are abandoned,
Figure BDA0003138290860000106
is the power that can be generated by the wind power,
Figure BDA0003138290860000107
is the amount of power that can be generated by the photovoltaic,
Figure BDA0003138290860000108
is the power of the wind power on the internet,
Figure BDA0003138290860000109
is the power of the photovoltaic grid-connected system,
Figure BDA00031382908600001010
is the heat supply power of the wind power,
Figure BDA00031382908600001011
is the heating power of the photovoltaic system,
Figure BDA00031382908600001012
is the price of electricity purchase, ρ1Is the punishment price of the deviation of the natural gas cooling-heating-electric machine group from the planned output, rho2Is the penalty price, rho, of the deviation of the power purchase from the planned power purchase3Is the punishment price of wind and light abandonment.
Further, the constraint conditions in step S3 are:
electric power balance constraint of the power supply system:
Figure BDA00031382908600001013
wherein L isrIs a set of links between the r region and other regions,
Figure BDA00031382908600001014
is the transmission power of the tie-line l,
Figure BDA00031382908600001015
indicating that the tie i is inputting power into the area at time t,
Figure BDA00031382908600001016
indicating that at time t the tie i delivers power outside this area,
Figure BDA00031382908600001017
is the power of the electrical load and,
Figure BDA00031382908600001018
is the electric refrigerator power;
thermoelectric ratio constraint, upper and lower output limits constraint and climbing constraint of the natural gas cooling-heating-electric machine set:
Figure BDA00031382908600001019
Figure BDA00031382908600001020
Figure BDA00031382908600001021
wherein, Khp,r,tIs the thermoelectric ratio of the natural gas cooling-heating-electric machine set of the ith station in the r area,
Figure BDA00031382908600001022
is the lower limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure BDA00031382908600001023
is the upper limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure BDA00031382908600001024
the maximum up-regulated power in unit time of the ith natural gas cooling-heating-electric machine set in the r area;
the power of purchasing electricity is restricted, each functional area only allows purchasing electricity and gas, but not selling electricity and gas:
Figure BDA0003138290860000111
Figure BDA0003138290860000112
wherein the content of the first and second substances,
Figure BDA0003138290860000113
is the capacity of the substation at the point of connection,
Figure BDA0003138290860000114
is the power of the external gas purchase,
Figure BDA0003138290860000115
is the upper limit of the natural gas pipeline delivery capacity; wind power photovoltaic output restraint:
Figure BDA0003138290860000116
Figure BDA0003138290860000117
Figure BDA0003138290860000118
Figure BDA0003138290860000119
Figure BDA00031382908600001110
Figure BDA00031382908600001111
the heat supply system is restricted, the electric boiler in the r area only consumes the wind power and the photovoltaic power in the r area, and the electric heat conversion relation restriction of the electric boiler and the output restriction of the electric boiler are met:
Figure BDA00031382908600001112
Figure BDA00031382908600001113
Figure BDA00031382908600001114
wherein the content of the first and second substances,
Figure BDA00031382908600001115
is the thermal output, eta, of the electric boiler at the time of r region t under the scene sEB,rIs the electrical conversion efficiency of the r-zone electric boiler,
Figure BDA00031382908600001116
is the upper limit of the output of the r area electric boiler;
and (3) constraint of start-stop time of the equipment:
Figure BDA00031382908600001117
Figure BDA00031382908600001118
Figure BDA00031382908600001119
Figure BDA00031382908600001120
wherein the content of the first and second substances,
Figure BDA00031382908600001121
is the continuous operation time of the natural gas cooling-heating-electric generating set,
Figure BDA00031382908600001122
is the continuous downtime of the natural gas cooling-heating-electric machine set,
Figure BDA00031382908600001123
is the duration of the operation of the electric boiler plant,
Figure BDA00031382908600001124
is the continuous down time of the electric boiler plant,
Figure BDA0003138290860000121
is the minimum start-up time of natural gas cooling-heating-power,
Figure BDA0003138290860000122
is the natural gas cold-heat-electricity minimum down time,
Figure BDA0003138290860000123
respectively the minimum start-up time of the electric boiler plant,
Figure BDA0003138290860000124
minimum down time of the electric boiler plant, uMT,r,t-1Is the operation state of the natural gas cooling-heating-electric machine set at the time t-1, and is a variable of 0-1 uEB,r,t-1The operation state of the electric boiler equipment at the time t-1 is a variable of 0-1;
the transfer channel capacity satisfies the constraint:
Figure BDA0003138290860000125
Figure BDA0003138290860000126
wherein the content of the first and second substances,
Figure BDA0003138290860000127
is the line capacity between nodes i and j at the initial time,
Figure BDA0003138290860000128
is the newly added line capacity between nodes i and j at time t,
Figure BDA0003138290860000129
is the maximum capacity of the line allowed to be erected in the power transmission channel between the nodes i and j, and M is a penalty coefficient, so that
Figure BDA00031382908600001210
When the temperature of the water is higher than the set temperature,
Figure BDA00031382908600001211
is 0, when
Figure BDA00031382908600001212
Large enough to have enough room to select a new line,
Figure BDA00031382908600001213
is a boolean variable, continuously constrained "switch", that produces a maximum capacity greater than or equal to that allowed by the line upgrade.
And (4) constraint of spare capacity:
Figure BDA00031382908600001214
Figure BDA00031382908600001215
wherein r is1Is the upper spare coefficient of the load, r2Is the upper spare coefficient of wind power, r3Is the upper standby coefficient of the photovoltaic, r4Is the lower standby coefficient of wind power, r5Is the lower standby factor for the photovoltaic.
The invention is further described below with reference to specific embodiments:
will be in a certain areaThe industrial park power supply system serves as a simulation example. The garden comprises photovoltaic power, wind power and natural gas, and the annual average radiation intensity of the region is 1200W/m2Average wind speed 4.7m/s and area of 64000m2Annual electricity consumption is about 2120 MWh.
The electricity purchasing price from the park to the power grid adopts a time-of-use electricity price, the price for purchasing natural gas from the park to the natural gas grid also adopts a time-of-use price, and the heat value of the natural gas is 9.7kWh/m3And performing power conversion. The upper limit of the external gas is 50MW and the upper limit of the external gas is 80m3. Specific prices are shown in table 1.
TABLE 1 park energy trade price
Figure BDA00031382908600001216
By adopting the solving algorithm proposed in this embodiment, the optimal configuration result shown in table 2 is obtained:
TABLE 2 configuration results of park devices
Device Configured capacity/kW Device Configured capacity/kW
Photovoltaic generator set 2160 Waste heat boiler 1566
Wind generating set 1830 Gas boiler 1409
Gas turbine 13062 Electric refrigerator 701
Absorption refrigerator 2000 Electric heating machine 179
The maximum electric load of a typical day of a cooling season is 2095kW, the cooling load is 4413kW, the heat load is 716kW, and the system efficiency and economic indexes under the condition of the typical day of the cooling season are solved and analyzed:
TABLE 3 park energy interaction cost
Cost of energy interaction energy/kWh Cost/dollar
Electricity purchasing device 14017.6 8273.1
Gas purchase 71175.9 22013.2
Operation of the system / 30302.3
As can be seen from table 3, the main cost of the system is on the purchase of natural gas in the case where the system is consuming renewable energy in full. Table 4 shows the energy supply, efficiency and operating costs required for each plant on a typical day.
TABLE 4 park Equipment operating costs
Device Energy supply/kWh Efficiency of Running cost
Photovoltaic generator set 3610.6 7.5％ 38.96
Wind generating set 6215.32 14.6％ 385.52
Gas turbine 213520.8 68.0％ 1259.81
Absorption refrigerator 31902.1 66.5％ 41.5
Waste heat boiler 9360.23 55.9％ 280.1
Gas boiler 2125.5 12.1％ 91.4
Electric refrigerator 5105.3 33.3％ 23.7
Electric heating machine 345.1 8.0％ 12.8
The efficiency in table 4 represents the ratio of the energy supplied by the plant to the energy of the plant operating at maximum power throughout the day, and it can be seen that the operating efficiency of the new energy unit is low, and in the cooling season, the efficiency of the gas turbine and the refrigeration plant is high, corresponding to the large cooling load demand in the cooling season. Table 5 energy use cost savings for the campus.
TABLE 5 energy usage cost savings in parks
Before/after optimization After/after optimization Cost saving/element
123191.4 116598.5 6592.9(5.4％)
Therefore, the energy consumption cost can be reduced by the optimal configuration of the multi-energy complementary capacity, additional economic benefits are obtained, and 5.4% of energy consumption cost is saved after the optimal configuration.
Finally, it should be pointed out that: the above examples are only for illustrating the technical solutions of the present invention, and are not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (5)

1. The method for configuring the multi-energy complementary capacity of the integrated energy park is characterized by comprising the following steps of:
s1, establishing an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple power supply system and electric heating and electric refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy utilization scenes and loads of the integrated energy utilization park, including electric loads, heat loads and cold loads in hot seasons, cold seasons and transition seasons;
s3, with the lowest total operation cost of the system as an optimization target, and with the constraint of electric power balance of a power supply system, the constraint of a natural gas cooling-heat-electric generator set, the constraint of electricity purchasing power, the constraint of wind power photovoltaic output, the constraint of heat supply system output, the constraint of start-stop time, the constraint of capacity of a conveying channel and the constraint of reserve capacity as constraint conditions, establishing a comprehensive energy use park multi-energy complementary capacity configuration optimization model, and solving the model by adopting a genetic algorithm.
2. The method of configuring multipotent complementary capacity of integrated energy use park according to claim 1, wherein said step S1 comprises the following steps:
the output of the fan is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
Figure FDA0003138290850000011
wherein the content of the first and second substances,
Figure FDA0003138290850000012
is the maximum output electric power v of the wind turbine generator at time ttIs the wind speed, v, over a period of time tinIs the cut-in wind speed, vratedIs rated wind speed, voutIs the cut-out wind speed, PS WTOutputting rated output electric power for the wind turbine generator;
secondly, the principle of photovoltaic power generation is that according to the photovoltaic effect, solar energy is converted into electric energy by utilizing a photovoltaic module. The output of the photovoltaic unit is related to the solar radiation intensity and the temperature, and the output model is as follows:
Figure FDA0003138290850000013
wherein the content of the first and second substances,
Figure FDA0003138290850000014
is the maximum output electric power of the photovoltaic generator set in the time period t,
Figure FDA0003138290850000015
is the output electric power of the photovoltaic unit under standard conditions, IstcIs the standard condition solar radiation intensity, TstcIs the standard condition temperature, ItIs the intensity of solar radiation, T, of the time period TtIs the temperature of time period t, α is the power temperature coefficient of the photovoltaic cell;
the natural gas cold-heat-electricity triple power supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, can supply cold, heat and electricity simultaneously, and has the following output model:
Figure FDA0003138290850000016
Figure FDA0003138290850000021
Figure FDA0003138290850000022
Figure FDA0003138290850000023
wherein the content of the first and second substances,
Figure FDA0003138290850000024
is the residual heat power of the natural gas cold-heat-electricity triple supply system in the time period t,
Figure FDA0003138290850000025
is the output electric power, eta, of the natural gas cold-heat-electricity combined supply system in the time period tMT,eIs the efficiency of electric power generation, ηHLIs the coefficient of heat dissipation loss, and,
Figure FDA0003138290850000026
is the output thermal power, eta, of the natural gas cooling-heating-power combined supply system in the time period tWHRIs the waste heat recovery efficiency, OHIs the coefficient of heat generation, and,
Figure FDA0003138290850000027
is the output cold power O of the natural gas cold-heat-electricity triple supply system in the time period tCIs the coefficient of refrigeration, VNGIs the consumption of natural gas, qNGIs the low heating value of natural gas;
electric energy belongs to high-quality energy, heat energy can be produced by electric heating equipment such as a heat pump and the like, cold energy can also be produced by the electric refrigerating equipment, and the output model is as follows:
Figure FDA0003138290850000028
Figure FDA0003138290850000029
wherein the content of the first and second substances,
Figure FDA00031382908500000210
is the output thermal power, η, of the electric heating apparatus over a period of time tEHIt is the efficiency of the electric-to-heat conversion,
Figure FDA00031382908500000211
is the input power rate of the electric heating device for time period t;
Figure FDA00031382908500000212
is the output cold power, eta, of the electric refrigerating apparatus in the time interval tECIt is the electric-to-cold conversion efficiency,
Figure FDA00031382908500000213
is the input power rate of the electric refrigeration appliance for time period t.
3. The method of configuring multipotent complementary capacity of integrated energy use park according to claim 1, wherein said step S2 comprises:
calculating the electric load, the heat load and the cold load in the hot season, wherein the calculation formula is as follows:
Figure FDA00031382908500000214
wherein the content of the first and second substances,
Figure FDA00031382908500000215
is the electrical load for the heating season t,
Figure FDA00031382908500000216
is the heat load for the heating season t,
Figure FDA00031382908500000217
is the cooling load of the heating season t;
secondly, calculating the electric load, the heat load and the cold load in the cooling season, wherein the calculation formula is as follows:
Figure FDA00031382908500000218
wherein the content of the first and second substances,
Figure FDA00031382908500000219
is the electrical load for the cooling season period t,
Figure FDA00031382908500000220
is the thermal load for the cooling season period t,
Figure FDA00031382908500000221
is the cooling load for the cooling season period t;
thirdly, calculating the electric load, the heat load and the cold load in the transition season, wherein the calculation formula is as follows:
Figure FDA0003138290850000031
wherein the content of the first and second substances,
Figure FDA0003138290850000032
is the electrical load for the transitional season period t,
Figure FDA0003138290850000033
is the thermal load for the transitional season period t,
Figure FDA0003138290850000034
is the cooling load for the transitional season period t.
4. The method of claim 1, wherein the optimization model of the multipotent complementary capacity configuration of the integrated energy park in step S3 is:
and (3) establishing an objective function by taking the lowest total operation cost of the system as an optimization objective:
Figure FDA0003138290850000035
Figure FDA0003138290850000036
Figure FDA0003138290850000037
wherein F is the total system operation cost which comprises the power generation cost of the natural gas cold-heat-electric generator set, the penalty cost of deviating from the planned output, the power purchase cost, the penalty cost of deviating from the planned power purchase and the cost of abandoned wind and abandoned light，NrMTIs the number of natural gas cooling-heating-electric units in the r region,
Figure FDA0003138290850000038
is the electric output of the ith natural gas cooling-heating-electric machine set at the moment t of an r area under a scene s,
Figure FDA0003138290850000039
the planned electric output of the ith natural gas cooling-heating-electric machine set at the moment t of the r area,
Figure FDA00031382908500000310
is the thermal output of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure FDA00031382908500000311
is the electricity purchasing power of the ith natural gas cooling-heating-electric machine set at the time t of an r area under a scene s,
Figure FDA00031382908500000312
the planned purchase power of the ith natural gas cooling-heating-electric machine set at the moment t in the r area,
Figure FDA00031382908500000313
is a cost function of the natural gas cooling-heating-electric machine set, alpha0,r,iFor the initial investment cost of the natural gas cooling-heating-electric machine, alpha1,r,i、α3,r,iRespectively a first term coefficient and a second term coefficient of the electric output in the cost function of the natural gas cold-heat-electric machine set, alpha2,r,i、α4,r,iPrimary and secondary coefficient of thermal output in the cost function of natural gas cold-heat-electric machine, alpha5,r,iIn order to be a cross-term coefficient,
Figure FDA00031382908500000314
the method is characterized in that the wind and the light power are abandoned,
Figure FDA00031382908500000315
is the power that can be generated by the wind power,
Figure FDA00031382908500000316
is the amount of power that can be generated by the photovoltaic,
Figure FDA00031382908500000317
is the power of the wind power on the internet,
Figure FDA00031382908500000318
is the power of the photovoltaic grid-connected system,
Figure FDA00031382908500000319
is the heat supply power of the wind power,
Figure FDA00031382908500000320
is the heating power of the photovoltaic system,
Figure FDA00031382908500000321
is the price of electricity purchase, ρ1Is the punishment price of the deviation of the natural gas cooling-heating-electric machine group from the planned output, rho2Is the penalty price, rho, of the deviation of the power purchase from the planned power purchase3Is the punishment price of wind and light abandonment.
5. The method of configuring multipotent complementary capacity of integrated energy usage park according to claim 4, wherein the constraint conditions in the step S3 are:
electric power balance constraint of the power supply system:
Figure FDA0003138290850000041
wherein L isrIs a set of links between the r region and other regions,
Figure FDA0003138290850000042
is the transmission power of the tie-line l,
Figure FDA0003138290850000043
indicating that the tie i is inputting power into the area at time t,
Figure FDA0003138290850000044
indicating that at time t the tie i delivers power outside this area,
Figure FDA0003138290850000045
is the power of the electrical load and,
Figure FDA0003138290850000046
is the electric refrigerator power;
thermoelectric ratio constraint, upper and lower output limits constraint and climbing constraint of the natural gas cooling-heating-electric machine set:
Figure FDA0003138290850000047
Figure FDA0003138290850000048
Figure FDA0003138290850000049
wherein, Khp,r,tIs the thermoelectric ratio of the natural gas cooling-heating-electric machine set of the ith station in the r area,
Figure FDA00031382908500000410
is the lower limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure FDA00031382908500000411
is the upper limit of the electric output of the ith natural gas cooling-heating-electric machine set in the r area,
Figure FDA00031382908500000412
the maximum up-regulated power in unit time of the ith natural gas cooling-heating-electric machine set in the r area;
the power of purchasing electricity is restricted, each functional area only allows purchasing electricity and gas, but not selling electricity and gas:
Figure FDA00031382908500000413
Figure FDA00031382908500000414
wherein the content of the first and second substances,
Figure FDA00031382908500000415
is the capacity of the substation at the point of connection,
Figure FDA00031382908500000416
is the power of the external gas purchase,
Figure FDA00031382908500000417
is the upper limit of the natural gas pipeline delivery capacity;
wind power photovoltaic output restraint:
Figure FDA00031382908500000418
Figure FDA00031382908500000419
Figure FDA00031382908500000420
Figure FDA0003138290850000051
Figure FDA0003138290850000052
Figure FDA0003138290850000053
the heat supply system is restricted, the electric boiler in the r area only consumes the wind power and the photovoltaic power in the r area, and the electric heat conversion relation restriction of the electric boiler and the output restriction of the electric boiler are met:
Figure FDA0003138290850000054
Figure FDA0003138290850000055
Figure FDA0003138290850000056
wherein the content of the first and second substances,
Figure FDA0003138290850000057
is the thermal output, eta, of the electric boiler at the time of r region t under the scene sEB,rIs the electrical conversion efficiency of the r-zone electric boiler,
Figure FDA0003138290850000058
is the upper limit of the output of the r area electric boiler;
and (3) constraint of start-stop time of the equipment:
Figure FDA0003138290850000059
Figure FDA00031382908500000510
Figure FDA00031382908500000511
Figure FDA00031382908500000512
wherein the content of the first and second substances,
Figure FDA00031382908500000513
is the continuous operation time of the natural gas cooling-heating-electric generating set,
Figure FDA00031382908500000514
is the continuous downtime of the natural gas cooling-heating-electric machine set,
Figure FDA00031382908500000515
is the duration of the operation of the electric boiler plant,
Figure FDA00031382908500000516
is the continuous down time of the electric boiler plant,
Figure FDA00031382908500000517
is the minimum start-up time of natural gas cooling-heating-power,
Figure FDA00031382908500000518
is the natural gas cold-heat-electricity minimum down time,
Figure FDA00031382908500000519
are respectively provided withFor the minimum start-up time of the electric boiler plant,
Figure FDA00031382908500000520
minimum down time of the electric boiler plant, uMT,r,t-1Is the operation state of the natural gas cooling-heating-electric machine set at the moment t-1, uEB,r,t-1Is the operation state of the electric boiler equipment at the time t-1;
the transfer channel capacity satisfies the constraint:
Figure FDA00031382908500000521
Figure FDA00031382908500000522
wherein the content of the first and second substances,
Figure FDA0003138290850000061
is the line capacity between nodes i and j at the initial time,
Figure FDA0003138290850000062
is the newly added line capacity between nodes i and j at time t,
Figure FDA0003138290850000063
is the maximum capacity of the line allowed to be erected in the power transmission channel between the nodes i and j, and M is a penalty coefficient, so that
Figure FDA0003138290850000064
When the temperature of the water is higher than the set temperature,
Figure FDA0003138290850000065
is 0, when
Figure FDA0003138290850000066
Large enough to have enough room to select a new line,
Figure FDA0003138290850000067
is a boolean variable, continuously constrained "switch" that produces a maximum capacity greater than or equal to that allowed for line upgrades;
and (4) constraint of spare capacity:
Figure FDA0003138290850000068
Figure FDA0003138290850000069
wherein r is1Is the upper spare coefficient of the load, r2Is the upper spare coefficient of wind power, r3Is the upper standby coefficient of the photovoltaic, r4Is the lower standby coefficient of wind power, r5Is the lower standby factor for the photovoltaic.
CN202110725103.4A 2021-06-29 2021-06-29 Multi-energy complementary capacity configuration method for comprehensive energy park Active CN113469430B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110725103.4A CN113469430B (en) 2021-06-29 2021-06-29 Multi-energy complementary capacity configuration method for comprehensive energy park

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110725103.4A CN113469430B (en) 2021-06-29 2021-06-29 Multi-energy complementary capacity configuration method for comprehensive energy park

Publications (2)

Publication Number Publication Date
CN113469430A true CN113469430A (en) 2021-10-01
CN113469430B CN113469430B (en) 2024-02-09

Family

ID=77873921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110725103.4A Active CN113469430B (en) 2021-06-29 2021-06-29 Multi-energy complementary capacity configuration method for comprehensive energy park

Country Status (1)

Country Link
CN (1) CN113469430B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114819508A (en) * 2022-03-28 2022-07-29 上海交通大学 Method and system for calculating distributed photovoltaic maximum access capacity of comprehensive energy system
CN115619145A (en) * 2022-10-14 2023-01-17 国网江苏省电力有限公司电力科学研究院 Cooperative control method and device for comprehensive energy system and computer equipment

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160172854A1 (en) * 2013-08-13 2016-06-16 Accenture Global Services Limited System for integrated multi-energy scheduling control in a micro-grid
CN108183500A (en) * 2017-11-24 2018-06-19 国网甘肃省电力公司电力科学研究院 A kind of rural area provided multiple forms of energy to complement each other is micro- can net capacity configuration optimizing method and device
CN108537409A (en) * 2018-03-15 2018-09-14 广东电网有限责任公司电网规划研究中心 A kind of industrial park power distribution network collaborative planning method considering multiple-energy-source coupled characteristic
CN108665188A (en) * 2018-05-23 2018-10-16 国网天津市电力公司电力科学研究院 A kind of garden multiple-energy-source main body synthesis matching method based on Optimized model
CN109784564A (en) * 2019-01-22 2019-05-21 天津大学 Consider the garden integrated energy system energy source station planing method of renewable energy access
WO2019196375A1 (en) * 2018-04-13 2019-10-17 华南理工大学 Demand side response-based microgrid optimal unit and time-of-use electricity price optimization method
CN110661254A (en) * 2019-09-26 2020-01-07 广东电网有限责任公司 Method for quantifying cooling, heating and power complementary benefits of regional comprehensive energy system
CN111445090A (en) * 2020-04-21 2020-07-24 清华大学 Double-layer planning method for off-grid type comprehensive energy system
CN112803499A (en) * 2021-03-17 2021-05-14 河海大学 Wind, light and water multi-energy complementary capacity optimal configuration method with power/electric quantity compensation cooperation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160172854A1 (en) * 2013-08-13 2016-06-16 Accenture Global Services Limited System for integrated multi-energy scheduling control in a micro-grid
CN108183500A (en) * 2017-11-24 2018-06-19 国网甘肃省电力公司电力科学研究院 A kind of rural area provided multiple forms of energy to complement each other is micro- can net capacity configuration optimizing method and device
CN108537409A (en) * 2018-03-15 2018-09-14 广东电网有限责任公司电网规划研究中心 A kind of industrial park power distribution network collaborative planning method considering multiple-energy-source coupled characteristic
WO2019196375A1 (en) * 2018-04-13 2019-10-17 华南理工大学 Demand side response-based microgrid optimal unit and time-of-use electricity price optimization method
CN108665188A (en) * 2018-05-23 2018-10-16 国网天津市电力公司电力科学研究院 A kind of garden multiple-energy-source main body synthesis matching method based on Optimized model
CN109784564A (en) * 2019-01-22 2019-05-21 天津大学 Consider the garden integrated energy system energy source station planing method of renewable energy access
CN110661254A (en) * 2019-09-26 2020-01-07 广东电网有限责任公司 Method for quantifying cooling, heating and power complementary benefits of regional comprehensive energy system
CN111445090A (en) * 2020-04-21 2020-07-24 清华大学 Double-layer planning method for off-grid type comprehensive energy system
CN112803499A (en) * 2021-03-17 2021-05-14 河海大学 Wind, light and water multi-energy complementary capacity optimal configuration method with power/electric quantity compensation cooperation

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BIN OUYANG 等: ""Research on Multi-Time Scale Optimization Strategy of Cold-Thermal-Electric Integrated Energy System Considering Feasible Interval of System Load Rate"", 《ENERGIES》, vol. 12, no. 17, pages 1 - 27 *
吴福保 等: "基于冷热电联供的多园区博弈优化策略", 《电力系统自动化》, no. 13, pages 74 - 81 *
方绍凤 等: "考虑电热多种负荷综合需求响应的园区微网综合能源系统优化运行", 《电力系统及其自动化学报》, vol. 32, no. 01, pages 50 - 57 *
李柳松 等: "考虑多功能区差异性以及冷热电负荷平衡的综合能源系统规划", 《能源工程》, no. 03, pages 5 - 11 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114819508A (en) * 2022-03-28 2022-07-29 上海交通大学 Method and system for calculating distributed photovoltaic maximum access capacity of comprehensive energy system
CN114819508B (en) * 2022-03-28 2024-03-29 上海交通大学 Comprehensive energy system distributed photovoltaic maximum admittance capacity calculation method and system
CN115619145A (en) * 2022-10-14 2023-01-17 国网江苏省电力有限公司电力科学研究院 Cooperative control method and device for comprehensive energy system and computer equipment
CN115619145B (en) * 2022-10-14 2024-03-19 国网江苏省电力有限公司电力科学研究院 Cooperative control method and device for comprehensive energy system and computer equipment

Also Published As

Publication number Publication date
CN113469430B (en) 2024-02-09

Similar Documents

Publication Publication Date Title
Li et al. Two-stage optimal operation of integrated energy system considering multiple uncertainties and integrated demand response
Zhang et al. A two-stage operation optimization method of integrated energy systems with demand response and energy storage
CN109919478B (en) Comprehensive energy microgrid planning method considering comprehensive energy supply reliability
Cao et al. An efficient and economical storage and energy sharing model for multiple multi-energy microgrids
CN111950807B (en) Comprehensive energy system optimization operation method considering uncertainty and demand response
CN112529244B (en) Comprehensive energy system collaborative optimization operation method considering electric load demand response
CN107358345B (en) Distributed combined cooling heating and power system optimization operation method considering demand side management
CN113256045B (en) Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
CN113469430B (en) Multi-energy complementary capacity configuration method for comprehensive energy park
CN109543889A (en) A kind of regional complex energy resource system cooperates with optimizing operation method a few days ago
CN112836882B (en) Regional comprehensive energy system operation optimization method considering equipment load rate change
CN112446546B (en) Comprehensive energy system two-stage optimal configuration method considering energy reliability
CN114707289A (en) Opportunity constraint-based multi-objective optimization method for electrothermal coupling comprehensive energy system
CN110620403A (en) Day-ahead scheduling method and system for collaborative operation of energy system considering renewable energy
CN114844124B (en) Operation control method of comprehensive energy system based on target optimization
CN109615125B (en) Multi-region random production simulation method considering extra-high voltage peak regulation and application
CN113806952A (en) Source-load-storage-considered cooling, heating and power comprehensive energy system and optimized operation method thereof
Cheng et al. A stochastic optimal model of micro energy internet contains rooftop PV and CCHP system
CN109193666B (en) General energy bus-based comprehensive energy system time sequence energy flow calculation method
CN111724026A (en) Optimization method for coupling operation of multi-energy network and water distribution network
CN107528352B (en) power distribution network active power optimization method based on renewable energy high permeability
CN115659585A (en) Micro-energy network low-carbon cooperative scheduling method and device considering demand response, memory and equipment
CN111262240B (en) Optimized operation method and system for comprehensive energy system
Meng et al. Economic optimization operation approach of integrated energy system considering wind power consumption and flexible load regulation
CN112436514A (en) Multi-microgrid interconnection optimization method considering photovoltaic uncertainty

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