Nothing Special   »   [go: up one dir, main page]

CN113469430B - Multi-energy complementary capacity configuration method for comprehensive energy park - Google Patents

Multi-energy complementary capacity configuration method for comprehensive energy park Download PDF

Info

Publication number
CN113469430B
CN113469430B CN202110725103.4A CN202110725103A CN113469430B CN 113469430 B CN113469430 B CN 113469430B CN 202110725103 A CN202110725103 A CN 202110725103A CN 113469430 B CN113469430 B CN 113469430B
Authority
CN
China
Prior art keywords
power
electric
cold
natural gas
constraint
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.)
Active
Application number
CN202110725103.4A
Other languages
Chinese (zh)
Other versions
CN113469430A (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

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Photovoltaic Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

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

Description

Multi-energy complementary capacity configuration method for comprehensive energy park
Technical Field
The invention relates to the technical field of power economy, in particular to a method for configuring a multifunctional complementary capacity of an integrated energy park.
Background
With the continuous increase of the development proportion of new energy, the basic characteristics of the traditional energy system are gradually changed, the limit of the original different energy systems is broken through by the comprehensive energy utilization park, and the integration of multiple types of energy is realized. The characteristics of randomness and indirection of resources such as wind, light, water, fire and the like, time-space complementation and the like bring challenges to the safe and stable operation of the power system, and the reasonable configuration of the capacity of each energy source is the key for ensuring the feasibility of the comprehensive energy utilization park, thereby being beneficial to improving the power supply reliability of the comprehensive energy source system. The reasonable distribution of the capacities of all power supplies in the multi-energy complementary power generation system is fully utilized, the advantages of all power supplies are brought into play, 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 allocation of the domestic comprehensive energy utilization park is only reserved in the multi-energy complementary capacity allocation of primary energy sources such as wind, light, water, fire and the like, the multi-energy complementary capacity allocation of secondary energy sources such as cold, heat, electricity and the like is not considered, and the allocation of electric loads, heat loads, cold loads and the like at the user side in different seasons is not related.
Therefore, it is necessary to study a method for configuring the multi-energy complementary capacity of the comprehensive energy park, which can optimize the multi-energy complementary capacity configuration of the comprehensive energy park, integrate the multi-energy sources, and improve the economical efficiency and reliability of the system of the comprehensive energy park.
Disclosure of Invention
The invention aims to provide a multi-energy complementary capacity configuration method for a comprehensive energy utilization park, which is used for modeling the distributed wind power, distributed photovoltaic, cold-heat-electricity triple supply system and the output of electric heating and electric refrigerating equipment of the comprehensive energy system and analyzing the typical energy utilization scene and load of the comprehensive energy utilization park, so that the multi-energy complementary capacity configuration optimization is carried out on the comprehensive energy utilization park, the integration of multiple energy sources is realized, and the economical efficiency and the reliability of the comprehensive energy utilization park system are improved.
The invention provides a method for configuring the multi-energy complementary capacity of a comprehensive energy utilization park, which comprises the following steps:
s1, building an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple supply system and electric heating and refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy consumption scenes and loads of the comprehensive energy utilization park, wherein the typical energy consumption scenes and loads comprise electric loads, heat loads and cold loads in a heating season, a cooling season and a transition season;
s3, taking the lowest total running cost of the system as an optimization target, taking electric power balance constraint of a power supply system, natural gas cold-hot-motor set constraint, electricity purchasing power constraint, wind power photovoltaic output constraint, heating system output constraint, start-stop time constraint, conveying channel capacity constraint and standby capacity constraint as constraint conditions, establishing a comprehensive energy 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:
(1) the fan output is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
wherein,maximum output electric power of the wind turbine generator set in period t, v t Is the wind speed of period t, v in Is cut-in wind speed v rated Is rated wind speed v out Is the cut-out wind speed, P S WT The rated output electric power of the wind turbine generator is used;
(2) the principle of photovoltaic power generation is to convert solar energy into electric energy by using a photovoltaic module according to the photovoltaic effect. The output force of the photovoltaic unit is related to the intensity and the temperature of solar radiation, and the output force model is as follows:
wherein,maximum output electric power of the photovoltaic unit in time period t, < >>Is the output electric power of the photovoltaic unit under the standard condition, I stc Is standard condition solar radiation intensity, generally 1000W/m 2 ,T stc Is at standard temperature, typically 25 ℃, I t Is the solar radiation intensity of period T, T t Is the temperature of period t, alpha is the power temperature coefficient of the photovoltaic cell, and is typically 0.0039 DEG C -1
(3) The natural gas cooling-heating-electricity triple supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, and can simultaneously supply three energy sources of cooling, heating and electricity, wherein the output model is as follows:
wherein,is the residual heat power of the natural gas cold-heat-electricity triple supply system in the period t, +.>Is the output electric power of the natural gas cold-heat-electricity triple supply system in the period t, eta MT,e Is the power generation efficiency eta HL Is the heat dissipation loss coefficient, generally 0.03,is the output heat power, eta of the natural gas cold-heat-electricity triple supply system in the period t WHR Is waste heat recovery efficiency, O H Heating coefficient is 1.20 # -, generally>Is the output cold power of the natural gas cold-heat-electricity triple supply system in the period t, O C Is a refrigeration coefficient, generally 0.95, V NG Is the consumption of natural gas, q NG Is the low heating value of natural gas;
(4) the electric energy belongs to high-quality energy, the heat energy can be produced by electric heating equipment such as a heat pump, and the cold energy can be produced by electric refrigerating equipment, and the output model is as follows:
wherein,is the output thermal power of the time period t electric heating equipment, eta EH Is the electric heat conversion efficiency->Is the input power rate of the time period t electric heating equipment; />Is the output cold power of the time period t electric refrigerating equipment, eta EC Is the conversion efficiency of electric cooling,is the input power rate of the time period t-electric refrigeration equipment.
Further, the step S2 includes:
(1) the electric load, the heat load and the cold load in the heating season are calculated according to the following calculation formula:
wherein,is the electrical load for heating Ji Shiduan t, +.>Is the heat load of the heat supply Ji Shiduan t, +.>Is the cold load of the heat supply Ji Shiduan t;
(2) the electrical, thermal and cold loads for the cold season were calculated as follows:
wherein,is the electrical load for the cold season period t, < ->Is the heat load for the cold season period t, < ->Is the cold load for the cold supply season period t;
(3) the electrical, thermal and cold loads in transition season were calculated as follows:
wherein,is the electrical load of the transition season period t, +.>Is the thermal load of the transition season t, +.>Is the cold load of the transition season period t.
Further, the comprehensive energy park multi-energy complementary capacity configuration optimization model in the step S3 is as follows:
taking the minimum total running cost of the system as an optimization target, and establishing an objective function:
wherein F is the total cost of system operation, including the power generation cost of the natural gas cold-hot-motor group and the deviation planned output penalty cost, the electricity purchase cost and the deviation planned electricity purchase penalty cost and the waste wind and waste light cost, NrMT is the number of r-zone natural gas cold-hot-motor units,is the electric force of the i-th natural gas cold-hot-motor group in the r region t moment under the scene s,/->For the planned power output of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is the heat output of the ith natural gas cold-heat-motor unit in the r region t moment under the scene s,/->Is the electricity purchasing power of the ith natural gas cold-hot-motor unit in the r region t moment under the scene s, < +.>For the planned power supply of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is a cost function of a natural gas cold-hot-motor unit, alpha 0,r,i Initial investment cost for natural gas cooling-heating-motor set, alpha 1,r,i 、α 3,r,i The primary term coefficient and the secondary term coefficient of the electric output in the function of the natural gas cold-hot-motor composition are respectively alpha 2,r,i 、α 4,r,i The primary term coefficient and the secondary term coefficient of the thermal output of the natural gas cold-hot-motor composition function are respectively alpha 5,r,i For the cross term coefficient, alpha is calculated from historical data and experience, and +>Is the wind-abandoning and light-abandoning power, < >>Is the power of wind power, which is->Is the power of the photovoltaic, in +.>Is the Internet power of wind power, < >>Is the internet power of the photovoltaic, +.>Is the heat supply power of wind power, < >>Is the heating power of the photovoltaic, +.>Is the purchase price of electricity, ρ 1 Is punishment price of the natural gas cold-hot-motor set deviating from planned output, ρ 2 Is punishment price of electricity purchasing power deviating from planned electricity purchasing power, ρ 3 Is the punishment price of wind and light abandoning.
Further, the constraint conditions in the step S3 are:
power supply system electric power balance constraint:
wherein L is r Is the set of links of the r region with other regions,is the transmission power of the link l, +.>Indicating at time t that tie l inputs power into the area,/>Indicating that at time t the tie line l is delivering power outside the area,/>Is the electrical load power, ">Is the power of the electric refrigerator;
thermoelectric ratio constraint, upper and lower output limit constraint and climbing constraint of a natural gas cold-hot-motor unit:
wherein K is hp,r,t Is the thermoelectric ratio of the ith natural gas cold-hot-motor unit in the r region,is the lower limit of the electric output of the ith natural gas cold-hot-motor unit in the r area,/->Is the electricity of the ith natural gas cold-hot-motor group in the r areaThe upper limit of the output force is set,the power is maximally adjusted upwards in unit time of the ith natural gas cooling-heating-motor unit in the r area;
and (3) the power purchasing constraint, wherein each functional area only allows power purchasing and gas purchasing, and does not allow power selling and gas selling:
wherein,is the capacity of a grid-connected point transformer substation, which is%>Is the outsourcing gas power, < >>Is the upper limit of the conveying capacity of the natural gas pipeline;
wind-powered photovoltaic force constraint:
and the constraint of a heating system is that an electric boiler in an r area only consumes wind power and photovoltaic power in the r area, so that the constraint of an electric-heating conversion relation of the electric boiler and the constraint of the output of the electric boiler are met:
wherein,is the thermal output, eta of the electric boiler at the moment t in the r region under the scene s EB,r Is the electrical conversion efficiency of the r-zone electric boiler, < >>Is the upper limit of the output of the electric boiler in the r area;
equipment start-up time constraint:
wherein,is the duration of the natural gas cold-hot-electric machine group,/->Is the continuous down time of the natural gas cold-hot-electric motor group, < >>Is the duration of the electric boiler plant, +.>Is the continuous down time of the electric boiler plant,
is the minimum start-up time of the natural gas cold-hot-electricity, +.>Is the minimum down time of natural gas cold-hot-electricity, +.>Respectively the minimum start-up time of the electric boiler plant,/->Minimum downtime of electric boiler plants, u MT,r,t-1 Is a natural gas cold-hot-motor group at t-1The operating state of the scale is 0-1 variable, u EB,r,t-1 The operation state of the electric boiler equipment at the time t-1 is a variable of 0-1;
the delivery channel capacity satisfies the constraint:
wherein,is the line capacity between nodes i and j at the initial moment, is +.>Is the newly added line capacity between nodes i and j at time t,/>Is the maximum capacity of the line allowed to be erected by the transmission channel between the nodes i and j, M is the penalty coefficient, so thatWhen (I)>0, when->When large enough to have enough room to select a new line, +>Is a boolean variable, a continuously constrained "switch," producing greater than or equal to the maximum capacity allowed by the line upgrade.
Spare capacity constraint:
wherein r is 1 Is the upper standby coefficient of the load, r 2 Is the upper standby coefficient of wind power, r 3 Is the upper standby coefficient of the photovoltaic, r 4 Is the lower standby coefficient of wind power, r 5 Is the lower standby coefficient of 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 load, heat load, cold load and the like at the user side in different seasons, fully utilizes the reasonable distribution of the capacities of all power sources in the multi-energy complementary power generation system, exerts the advantages of all power sources, realizes the integration of multiple energy sources, can increase the consumption of clean energy sources, reduces the total investment of the system, and improves the economy and the reliability of the system in the comprehensive energy park.
Drawings
FIG. 1 is a flow chart of a method for configuring multiple complementary capacities of an integrated energy park provided in an embodiment of the present invention;
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention become more apparent, 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 some, but not all, embodiments of the invention. The embodiments described below by referring to the drawings are illustrative and intended to explain the present invention and should not be construed as limiting the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the 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 park according to the embodiment includes the following steps:
s1, building an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple supply system and electric heating and refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy consumption scenes and loads of the comprehensive energy utilization park, wherein the typical energy consumption scenes and loads comprise electric loads, heat loads and cold loads in a heating season, a cooling season and a transition season;
s3, taking the lowest total running cost of the system as an optimization target, taking electric power balance constraint of a power supply system, natural gas cold-hot-motor set constraint, electricity purchasing power constraint, wind power photovoltaic output constraint, heating system output constraint, start-stop time constraint, conveying channel capacity constraint and standby capacity constraint as constraint conditions, establishing a comprehensive energy 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:
(1) the fan output is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
wherein,maximum output electric power of the wind turbine generator set in period t, v t Is the wind speed of period t, v in Is cut-in wind speed v rated Is rated wind speed v out Is the cut-out wind speed, P S WT The rated output electric power of the wind turbine generator is used;
(2) the principle of photovoltaic power generation is to convert solar energy into electric energy by using a photovoltaic module according to the photovoltaic effect. The output force of the photovoltaic unit is related to the intensity and the temperature of solar radiation, and the output force model is as follows:
wherein,maximum output electric power of the photovoltaic unit in time period t, < >>Is the output electric power of the photovoltaic unit under the standard condition, I stc Is standard condition solar radiation intensity, generally 1000W/m 2 ,T stc Is at standard temperature, typically 25 ℃, I t Is the solar radiation intensity of period T, T t Is the temperature of period t, alpha is the power temperature coefficient of the photovoltaic cell, and is typically 0.0039 DEG C -1
(3) The natural gas cooling-heating-electricity triple supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, and can simultaneously supply three energy sources of cooling, heating and electricity, wherein the output model is as follows:
wherein,is the residual heat power of the natural gas cold-heat-electricity triple supply system in the period t, +.>Is the output electric power of the natural gas cold-heat-electricity triple supply system in the period t, eta MT,e Is the power generation efficiency eta HL Is the heat dissipation loss coefficient, generally 0.03,is the output heat power, eta of the natural gas cold-heat-electricity triple supply system in the period t WHR Is waste heat recovery efficiency, O H Heating coefficient is 1.20 # -, generally>Is the output cold power of the natural gas cold-heat-electricity triple supply system in the period t, O C Is a refrigeration coefficient, generally 0.95, V NG Is the consumption of natural gas, q NG Is the low heating value of natural gas;
(4) the electric energy belongs to high-quality energy, the heat energy can be produced by electric heating equipment such as a heat pump, and the cold energy can be produced by electric refrigerating equipment, and the output model is as follows:
wherein,is the output thermal power of the time period t electric heating equipment, eta EH Is the electric heat conversion efficiency->Is the input power rate of the time period t electric heating equipment; />Is the output cold power of the time period t electric refrigerating equipment, eta EC Is the conversion efficiency of electric cooling,is the input power rate of the time period t-electric refrigeration equipment.
Further, the step S2 includes:
(1) the electric load, the heat load and the cold load in the heating season are calculated according to the following calculation formula:
wherein,is the electrical load for heating Ji Shiduan t, +.>Is the heat load of the heat supply Ji Shiduan t, +.>Is the cold load of the heat supply Ji Shiduan t;
(2) the electrical, thermal and cold loads for the cold season were calculated as follows:
wherein,is the electrical load for the cold season period t, < ->Is cold Ji Shiduant heat load, < >>Is the cold load for the cold supply season period t;
(3) the electrical, thermal and cold loads in transition season were calculated as follows:
wherein,is the electrical load of the transition season period t, +.>Is the thermal load of the transition season t, +.>Is the cold load of the transition season period t.
Further, the comprehensive energy park multi-energy complementary capacity configuration optimization model in the step S3 is as follows:
taking the minimum total running cost of the system as an optimization target, and establishing an objective function:
wherein F is the total cost of system operation, including the natural gas cold-hot-electric machine group power generation cost and the deviation plan output punishment cost, the electricity purchase cost and the deviation plan electricity purchase punishmentThe cost and the wind and light discarding cost are calculated, NrMT is the number of r-zone natural gas cold-hot-motor units,is the electric force of the i-th natural gas cold-hot-motor group in the r region t moment under the scene s,/->For the planned power output of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is the heat output of the ith natural gas cold-heat-motor unit in the r region t moment under the scene s,/->Is the electricity purchasing power of the ith natural gas cold-hot-motor unit in the r region t moment under the scene s, < +.>For the planned power supply of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is a cost function of a natural gas cold-hot-motor unit, alpha 0,r,i Initial investment cost for natural gas cooling-heating-motor set, alpha 1,r,i 、α 3,r,i The primary term coefficient and the secondary term coefficient of the electric output in the function of the natural gas cold-hot-motor composition are respectively alpha 2,r,i 、α 4,r,i The primary term coefficient and the secondary term coefficient of the thermal output of the natural gas cold-hot-motor composition function are respectively alpha 5,r,i For the cross term coefficient, alpha is calculated from historical data and experience, and +>Is the wind-abandoning and light-abandoning power, < >>Is the power of wind power, which is->Is the power of the photovoltaic, in +.>Is the Internet power of wind power, < >>Is the internet power of the photovoltaic, +.>Is the heat supply power of wind power, < >>Is the heating power of the photovoltaic, +.>Is the purchase price of electricity, ρ 1 Is punishment price of the natural gas cold-hot-motor set deviating from planned output, ρ 2 Is punishment price of electricity purchasing power deviating from planned electricity purchasing power, ρ 3 Is the punishment price of wind and light abandoning.
Further, the constraint conditions in the step S3 are:
power supply system electric power balance constraint:
wherein L is r Is the set of links of the r region with other regions,is the transmission power of the link l, +.>Indicating at time t that tie l inputs power into the area,/>Indicating that at time t the tie line l is delivering power outside the area,/>Is the electrical load power, ">Is the power of the electric refrigerator;
thermoelectric ratio constraint, upper and lower output limit constraint and climbing constraint of a natural gas cold-hot-motor unit:
wherein K is hp,r,t Is the thermoelectric ratio of the ith natural gas cold-hot-motor unit in the r region,is the lower limit of the electric output of the ith natural gas cold-hot-motor unit in the r area,/->Is the upper limit of the electric output of the ith natural gas cold-hot-motor group in the r region,the power is maximally adjusted upwards in unit time of the ith natural gas cooling-heating-motor unit in the r area;
and (3) the power purchasing constraint, wherein each functional area only allows power purchasing and gas purchasing, and does not allow power selling and gas selling:
wherein,is the capacity of a grid-connected point transformer substation, which is%>Is the outsourcing gas power, < >>Is the upper limit of the conveying capacity of the natural gas pipeline; wind-powered photovoltaic force constraint:
and the constraint of a heating system is that an electric boiler in an r area only consumes wind power and photovoltaic power in the r area, so that the constraint of an electric-heating conversion relation of the electric boiler and the constraint of the output of the electric boiler are met:
wherein,is the thermal output, eta of the electric boiler at the moment t in the r region under the scene s EB,r Is the electrical conversion efficiency of the r-zone electric boiler, < >>Is the upper limit of the output of the electric boiler in the r area;
equipment start-up time constraint:
wherein,is the duration of the natural gas cold-hot-electric machine group,/->Is the continuous down time of the natural gas cold-hot-electric motor group, < >>Is the duration of the electric boiler plant, +.>Is the continuous down time of the electric boiler plant,
is the minimum start-up time of the natural gas cold-hot-electricity, +.>Is the minimum down time of natural gas cold-hot-electricity, +.>Respectively the minimum start-up time of the electric boiler plant,/->Minimum downtime of electric boiler plants, u MT,r,t-1 Is the operation state of the natural gas cold-hot-motor group at the time t-1, is 0-1 variable, u EB,r,t-1 The operation state of the electric boiler equipment at the time t-1 is a variable of 0-1;
the delivery channel capacity satisfies the constraint:
wherein,is the line capacity between nodes i and j at the initial moment, is +.>Is the newly added line capacity between nodes i and j at time t,/>Is the maximum capacity of the line allowed to be erected by the transmission channel between the nodes i and j, M is the penalty coefficient, so thatWhen (I)>0, when->When large enough to have enough room to select a new line, +>Is a boolean variable, a continuously constrained "switch," producing greater than or equal to the maximum capacity allowed by the line upgrade.
Spare capacity constraint:
wherein r is 1 Is the upper standby coefficient of the load, r 2 Is the upper standby coefficient of wind power, r 3 Is the upper standby coefficient of the photovoltaic, r 4 Is the lower standby coefficient of wind power, r 5 Is the lower standby coefficient of the photovoltaic.
The invention is further described in connection with specific embodiments below:
an industrial park power supply system in a certain area is taken as a simulation example. The park comprises photovoltaic, wind power and natural gas, and the annual average radiation intensity of the park is 1200W/m 2 Average wind speed 4.7m/s, park area about 64000m 2 Annual electricity consumption is approximately 2120MWh.
The price of electricity purchased from the park to the power grid adopts time-sharing electricity price, the price of natural gas purchased from the park to the natural gas grid also adopts time-sharing price, and the heat value of the natural gas is 9.7kWh/m 3 And performing power conversion. The upper limit of the outsourcing electricity is 50MW and the upper limit of the outsourcing gas is 80m 3 . The specific prices are shown in Table 1.
TABLE 1 energy trading prices for parks
The solution algorithm proposed in this embodiment is adopted to obtain the optimal configuration result shown in table 2:
table 2 park device configuration results
Apparatus and method for controlling the operation of a device Configuration capacity/kW Apparatus and method for controlling the operation of a device Configuration capacity/kW
Photovoltaic deviceGenerating 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 electrical load on typical days in the cold season is 2095kW, the cold load is 4413kW, the heat load is 716kW, and the system efficiency and economic index under the typical days in the cold season are solved and analyzed:
TABLE 3 energy interaction costs for the campus
Cost of energy interaction energy/kWh Cost/element
Purchase electricity 14017.6 8273.1
Purchasing gas 71175.9 22013.2
System operation / 30302.3
As can be seen from table 3, the main cost of the system is on purchasing natural gas with the system fully absorbing renewable energy. Table 4 shows the energy, efficiency and required operating costs of each plant on a typical day.
TABLE 4 Garden plant operation costs
Apparatus and method for controlling the operation of a device energy/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 shows the ratio of the energy of 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 lower, and the efficiency of the gas turbine and the refrigeration plant is higher in the cooling season, corresponding to the large cooling load demand in the cooling season. Table 5 is a park energy cost savings.
Table 5 park energy cost savings
Pre-optimization/meta-optimization Optimized post/meta Cost savings/metaments
123191.4 116598.5 6592.9(5.4%)
Therefore, the energy consumption cost can be reduced by the multi-energy complementary capacity optimization configuration, additional economic benefits are obtained, and the energy consumption cost is saved by 5.4% after the optimization configuration.
Finally, it should be pointed out that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting. Although the 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 scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (4)

1. The method for configuring the multi-energy complementary capacity of the comprehensive energy park is characterized by comprising the following steps of:
s1, building an output model of a distributed wind power, distributed photovoltaic, cold-heat-electricity triple supply system and electric heating and refrigerating equipment in comprehensive energy supply;
s2, analyzing typical energy consumption scenes and loads of the comprehensive energy utilization park, wherein the typical energy consumption scenes and loads comprise electric loads, heat loads and cold loads in a heating season, a cooling season and a transition season;
s3, taking the lowest total running cost of the system as an optimization target, taking electric power balance constraint of a power supply system, natural gas cold-hot-motor set constraint, electricity purchasing power constraint, wind power photovoltaic output constraint, heating system output constraint, start-stop time constraint, conveying channel capacity constraint and standby capacity constraint as constraint conditions, establishing a comprehensive energy park multi-energy complementary capacity configuration optimization model, and solving the model by adopting a genetic algorithm;
the comprehensive energy park multi-energy complementary capacity configuration optimization model in the step S3 is as follows:
taking the minimum total running cost of the system as an optimization target, and establishing an objective function:
wherein F is the total cost of system operation, including the power generation cost of the natural gas cold-hot-motor group and the deviation planned output penalty cost, the electricity purchase cost and the deviation planned electricity purchase penalty cost and the waste wind and waste light cost, NrMT is the number of r-zone natural gas cold-hot-motor units,is the ith natural station in r region t moment under scene sElectric force of air-cooled-heated-motor group, < >>For the planned power output of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is the heat output of the ith natural gas cold-heat-motor unit in the r region t moment under the scene s,/->Is the electricity purchasing power of the ith natural gas cold-hot-motor unit in the r region t moment under the scene s,for the planned power supply of the ith natural gas cooling-heating-motor unit in the r region t time,/>Is a cost function of a natural gas cold-hot-motor unit, alpha 0,r,i Initial investment cost for natural gas cooling-heating-motor set, alpha 1,r,i 、α 3,r,i The primary term coefficient and the secondary term coefficient of the electric output in the function of the natural gas cold-hot-motor composition are respectively alpha 2,r,i 、α 4,r,i The primary term coefficient and the secondary term coefficient of the thermal output of the natural gas cold-hot-motor composition function are respectively alpha 5,r,i For cross term coefficients, ++>Is the wind-abandoning and light-abandoning power, < >>Is the power of wind power, which is->Is the power of the photovoltaic, in +.>Is the Internet power of wind power, < >>Is the internet power of the photovoltaic, +.>Is the heat supply power of wind power, < >>Is the heating power of the photovoltaic, +.>Is the purchase price of electricity, ρ 1 Is punishment price of the natural gas cold-hot-motor set deviating from planned output, ρ 2 Is punishment price of electricity purchasing power deviating from planned electricity purchasing power, ρ 3 Is the punishment price of wind and light abandoning.
2. The method for configuring a multi-energy complementary capacity of a comprehensive energy park according to claim 1, wherein the step S1 includes the steps of:
(1) the fan output is related to wind speed, cut-in wind speed, rated wind speed and cut-out wind speed, and the output model is as follows:
wherein,maximum output electric power of the wind turbine generator set in period t, v t Is the wind speed of period t, v in Is cut-in wind speed v rated Is rated wind speed v out Is the cut-out wind speed, P S WT Rated output power of wind turbine generator systemA power;
(2) 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 force of a photovoltaic unit is related to the intensity and the temperature of solar radiation, and the output force model is as follows:
wherein,maximum output electric power of the photovoltaic unit in time period t, < >>Is the output electric power of the photovoltaic unit under the standard condition, I stc Is the standard condition of solar radiation intensity, T stc Is the standard condition temperature, I t Is the solar radiation intensity of period T, T t Is the temperature of period t, α is the power temperature coefficient of the photovoltaic cell;
(3) the natural gas cooling-heating-electricity triple supply system generally comprises a gas turbine, an absorption refrigerator, a waste heat boiler and other equipment, and can simultaneously supply three energy sources of cooling, heating and electricity, wherein the output model is as follows:
wherein,is the residual heat power of the natural gas cold-heat-electricity triple supply system in the period t, +.>Is the output electric power of the natural gas cold-heat-electricity triple supply system in the period t, eta MT,e Is the power generation efficiency eta HL Is the heat dissipation loss coefficient->Is the output heat power, eta of the natural gas cold-heat-electricity triple supply system in the period t WHR Is waste heat recovery efficiency, O H Is the heating coefficient>Is the output cold power of the natural gas cold-heat-electricity triple supply system in the period t, O C Is the refrigeration coefficient, V NG Is the consumption of natural gas, q NG Is the low heating value of natural gas;
(4) the electric energy belongs to high-quality energy, the heat energy can be produced by electric heating equipment such as a heat pump, and the cold energy can be produced by electric refrigerating equipment, and the output model is as follows:
wherein,is time period t electric heatingOutput thermal power, eta of the device EH Is the electric heat conversion efficiency->Is the input power rate of the time period t electric heating equipment; />Is the output cold power of the time period t electric refrigerating equipment, eta EC Is the conversion efficiency of electric cooling,is the input power rate of the time period t-electric refrigeration equipment.
3. The method for configuring a multi-energy complementary capacity of a comprehensive energy park according to claim 1, wherein the step S2 includes:
(1) the electric load, the heat load and the cold load in the heating season are calculated according to the following calculation formula:
wherein,is the electrical load for heating Ji Shiduan t, +.>Is the heat load of the heat supply Ji Shiduan t, +.>Is the cold load of the heat supply Ji Shiduan t;
(2) the electrical, thermal and cold loads for the cold season were calculated as follows:
wherein,is the electrical load for the cold season period t, < ->Is the heat load for the cold season period t, < ->Is the cold load for the cold supply season period t;
(3) the electrical, thermal and cold loads in transition season were calculated as follows:
wherein,is the electrical load of the transition season period t, +.>Is the thermal load of the transition season t, +.>Is the cold load of the transition season period t.
4. The method for configuring a multi-energy complementary capacity of a comprehensive energy park according to claim 1, wherein the constraint condition in the step S3 is:
power supply system electric power balance constraint:
where Lr is the set of tie lines of the r region with other regions,is the transmission power of the link l, +.>Indicating at time t that tie l inputs power into the area,/>Indicating that at time t the tie line l is delivering power outside the area,/>Is the electrical load power, ">Is the power of the electric refrigerator;
thermoelectric ratio constraint, upper and lower output limit constraint and climbing constraint of a natural gas cold-hot-motor unit:
wherein K is hp,r,t Is the thermoelectric ratio of the ith natural gas cold-hot-motor unit in the r region,is the lower limit of the electric output of the ith natural gas cold-hot-motor unit in the r area,/->Is the upper limit of the electric output of the ith natural gas cold-hot-motor group in the r area,/->The power is maximally adjusted upwards in unit time of the ith natural gas cooling-heating-motor unit in the r area;
and (3) the power purchasing constraint, wherein each functional area only allows power purchasing and gas purchasing, and does not allow power selling and gas selling:
wherein,is the capacity of a grid-connected point transformer substation, which is%>Is the outsourcing gas power, < >>Is the upper limit of the conveying capacity of the natural gas pipeline;
wind-powered photovoltaic force constraint:
and the constraint of a heating system is that an electric boiler in an r area only consumes wind power and photovoltaic power in the r area, so that the constraint of an electric-heating conversion relation of the electric boiler and the constraint of the output of the electric boiler are met:
wherein,is the thermal output, eta of the electric boiler at the moment t in the r region under the scene s EB,r Is the electrical conversion efficiency of the r-zone electric boiler,is the upper limit of the output of the electric boiler in the r area;
equipment start-up time constraint:
wherein,is the duration of the natural gas cold-hot-electric machine group,/->Is the continuous down time of the natural gas cold-hot-electric motor group, < >>Is the duration of the electric boiler plant, +.>Is the continuous down time of the electric boiler plant,is natural gas cold-hot-minimum start-up time of electricity, +.>Is the minimum down time of natural gas cold-hot-electricity, +.>Respectively the minimum start-up time of the electric boiler plant,/->Minimum downtime of electric boiler plants, u MT,r,t-1 Is the operation state of the natural gas cold-hot-motor group at the time t-1, u EB,r,t-1 The operation state of the electric boiler equipment at the time t-1;
the delivery channel capacity satisfies the constraint:
wherein,is the line capacity between nodes i and j at the initial moment, is +.>Is the newly added line capacity between nodes i and j at time t,/>Is the maximum capacity of the line allowed to be erected by the transmission channel between the nodes i and j, M is the penalty coefficient, so thatWhen (I)>0, when->When large enough to have enough room to select a new line, +>Is a boolean variable, continuously constrained "switch" producing greater than or equal to the maximum capacity allowed by the line upgrade;
spare capacity constraint:
wherein r is 1 Is the upper standby coefficient of the load, r 2 Is the upper standby coefficient of wind power, r 3 Is the upper standby coefficient of the photovoltaic, r 4 Is the lower standby coefficient of wind power, r 5 Is the lower standby coefficient of 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 CN113469430A (en) 2021-10-01
CN113469430B true 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)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114186875A (en) * 2021-12-14 2022-03-15 广东电网有限责任公司 Park multi-energy optimization scheduling control method with high new energy ratio
CN114819508B (en) * 2022-03-28 2024-03-29 上海交通大学 Comprehensive energy system distributed photovoltaic maximum admittance capacity calculation method and system
CN115619145B (en) * 2022-10-14 2024-03-19 国网江苏省电力有限公司电力科学研究院 Cooperative control method and device for comprehensive energy system and computer equipment
CN116611199B (en) * 2023-02-24 2024-06-25 中国市政工程华北设计研究总院有限公司 Genetic algorithm-based capacity optimization configuration method for multi-energy complementary heating system

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015021603A1 (en) * 2013-08-13 2015-02-19 Accenture Global Services Limited System, method and apparatus for integrated multi-energy scheduling in a micro-grid and a tangible computer readable medium

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
"Research on Multi-Time Scale Optimization Strategy of Cold-Thermal-Electric Integrated Energy System Considering Feasible Interval of System Load Rate";Bin Ouyang 等;《Energies》;第12卷(第17期);1-27 *
基于冷热电联供的多园区博弈优化策略;吴福保 等;《电力系统自动化》(第13期);74-81 *
考虑多功能区差异性以及冷热电负荷平衡的综合能源系统规划;李柳松 等;《能源工程》(第03期);5-11 *
考虑电热多种负荷综合需求响应的园区微网综合能源系统优化运行;方绍凤 等;《电力系统及其自动化学报》;第32卷(第01期);50-57 *

Also Published As

Publication number Publication date
CN113469430A (en) 2021-10-01

Similar Documents

Publication Publication Date Title
CN113469430B (en) Multi-energy complementary capacity configuration method for comprehensive energy park
Zhang et al. A two-stage operation optimization method of integrated energy systems with demand response and energy storage
CN111950807B (en) Comprehensive energy system optimization operation method considering uncertainty and demand response
CN106447529A (en) Distributed energy system modeling and running optimization method considering hot water pipe network
CN113256045B (en) Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
CN107358345B (en) Distributed combined cooling heating and power system optimization operation method considering demand side management
CN114707289B (en) Multi-objective optimization method of electrothermal coupling comprehensive energy system based on opportunity constraint
CN112418479B (en) Optimal configuration method for park comprehensive energy system
CN109543889A (en) A kind of regional complex energy resource system cooperates with optimizing operation method a few days ago
CN102509175A (en) Reliability optimization method of distributed power supply system
CN113592200B (en) Low-carbon optimized operation method for regional comprehensive energy system of water-containing source heat pump
CN110807588A (en) Optimized scheduling method of multi-energy coupling comprehensive energy system
CN112434915B (en) Regional comprehensive energy system flexibility optimal configuration method for abandoned wind digestion
CN111244939A (en) Two-stage optimization design method for multi-energy complementary system considering demand side response
CN112446546B (en) Comprehensive energy system two-stage optimal configuration method considering energy reliability
CN113850409A (en) Comprehensive energy system optimized operation method considering renewable energy access
Meng et al. Economic optimization operation approach of integrated energy system considering wind power consumption and flexible load regulation
CN114844124B (en) Operation control method of comprehensive energy system based on target optimization
CN105447586A (en) Micro-grid energy optimization method based on versatile flow model
CN114362152B (en) Multi-time scale scheduling method for comprehensive energy system
CN113806952B (en) Cold-hot electricity comprehensive energy system considering source-charge-storage and optimal operation method thereof
CN114757388A (en) Regional comprehensive energy system equipment capacity optimization method based on improved NSGA-III
CN117035202B (en) Double-layer collaborative expansion planning method for electric heating comprehensive energy system considering demand response
CN110472364B (en) Optimization method of off-grid type combined heat and power generation system considering renewable energy sources
Cheng et al. A stochastic optimal model of micro energy internet contains rooftop PV and CCHP system

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