CN114529153A - Multi-region comprehensive energy system scheduling method considering comprehensive demand response - Google Patents

Multi-region comprehensive energy system scheduling method considering comprehensive demand response Download PDF

Info

Publication number
CN114529153A
CN114529153A CN202210039372.XA CN202210039372A CN114529153A CN 114529153 A CN114529153 A CN 114529153A CN 202210039372 A CN202210039372 A CN 202210039372A CN 114529153 A CN114529153 A CN 114529153A
Authority
CN
China
Prior art keywords
energy
load
comprehensive
power
energy system
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.)
Pending
Application number
CN202210039372.XA
Other languages
Chinese (zh)
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.)
Jiangmen Minghao Electric Power Engineering Supervision Co ltd
South China University of Technology SCUT
Original Assignee
Jiangmen Minghao Electric Power Engineering Supervision Co ltd
South China University of Technology SCUT
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 Jiangmen Minghao Electric Power Engineering Supervision Co ltd, South China University of Technology SCUT filed Critical Jiangmen Minghao Electric Power Engineering Supervision Co ltd
Priority to CN202210039372.XA priority Critical patent/CN114529153A/en
Publication of CN114529153A publication Critical patent/CN114529153A/en
Pending legal-status Critical Current

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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313—Resource planning in a project environment
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

The invention discloses a multi-region comprehensive energy system scheduling method considering comprehensive demand response. The method comprises the following steps: dividing a region to be researched into a plurality of regions according to the geographic position, and acquiring wind-solar output data and cold-hot electrical load data of each region; establishing a comprehensive demand response load model and a cost and income model; constructing an energy flow topology and an energy coupling matrix of a multi-region comprehensive energy system; and establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering the inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target. The method provided by the invention enhances the multi-energy coupling conversion capability of the multi-region comprehensive energy system, and further can effectively improve the overall economy, flexibility and reliability of the system.

Description

Multi-region comprehensive energy system scheduling method considering comprehensive demand response
Technical Field
The invention relates to the field of comprehensive energy system research, in particular to a multi-region comprehensive energy system scheduling method considering comprehensive demand response.
Background
With the increasing exhaustion of traditional fossil energy and the rapid development of energy internet related technologies, the coupling complementation and integration optimization technology of various energy sources has become one of the research hot problems in the field of comprehensive energy systems at home and abroad. The comprehensive energy system particularly refers to an energy production, supply and marketing integrated system formed by organically coordinating, scientifically scheduling and cascade utilizing links such as energy generation, transmission and distribution, conversion, storage, consumption and the like by exerting the multi-energy complementary characteristics of various types of energy. When a plurality of regional comprehensive energy systems exist in the same power distribution region, a multi-region comprehensive energy system or a multi-region electric cold and heat combined supply system is formed.
In order to reduce energy waste and realize optimal configuration of energy, the traditional energy system is upgraded and transformed in recent years, and the call for forming a regional comprehensive energy system cluster which can meet the complementary operation of various energy sources is more and more loud. The invention patent with the publication number of CN105183991A provides a regional comprehensive energy system planning and designing method, which establishes a time sequence model of electric load, heat load and cold load of a certain region, but the time sequence model does not consider a flexible and changeable electric-heat load adjustable model and is difficult to adapt to the actual situation; the invention patent with the publication number of CN110266004A provides a standardized construction method of an energy hub model of an integrated energy system, which considers the condition of energy circulation by adding a limited equation in the energy hub model of the integrated energy system, and avoids the condition that the coupling matrixes of the energy hub models of different simple integrated energy systems are not matched by adding a virtual energy conversion device; the invention patent with publication number CN111313429B provides a reliability evaluation method and system for an integrated energy system, which considers the differences of network topology, transmission delay, terminal thermal inertia and user reliability requirements, performs load flow calculation on the integrated energy system, determines the reliability index of the integrated energy system by using the load flow calculation result and combining a sequential monte carlo method, and improves the accuracy of the reliability evaluation of the integrated energy system.
In the prior art, most researches on a multi-region comprehensive energy system are in a starting stage, the load types considered by demand response are single, the coverage area is small, the complementary characteristics of electricity, gas, cold and heat on time and space are difficult to embody, and the multi-region comprehensive energy system is limited in time scale and capacity scale. Under the large background of rapid development of the energy interconnection technology, a new technical means is urgently needed to realize that a plurality of regional comprehensive energy systems cooperatively provide comprehensive demand response and perform multi-energy coupling optimization, so that the coupling characteristic of the interconnection system is enhanced, the overall response capability of the interconnection system is improved, and the economy, flexibility and reliability of the interconnection system are improved. Accordingly, the prior art is yet to be improved and developed.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a multi-region comprehensive energy system scheduling method considering comprehensive demand response. The method not only reduces the daily operation cost of the system, reduces the external energy purchase and the network side electricity purchase, but also obtains the peak clipping response compensation income by participating in the demand response market transaction, obviously reduces the total scheduling cost of the multi-region comprehensive energy system, and improves the overall economy of the system; the peak clipping and valley filling are realized by the participation of multi-energy adjustable load in the comprehensive demand response, and the overall energy utilization pressure is relieved; the load complementary characteristics among different regions and in each region are fully utilized, the electric energy mutual-aid capability and the multi-energy coupling conversion capability of the multi-region comprehensive energy system are improved, the supply and demand of each region are balanced, the reasonable allocation of resources is realized, the overall flexibility and reliability of the system are improved, and the system has good application value.
The purpose of the invention is realized by at least one of the following technical solutions.
A method for multi-zone integrated energy system dispatch taking into account integrated demand response, the method comprising the steps of:
s1: dividing a region to be researched into a plurality of regions according to the geographic position, and acquiring wind-solar output data and cold-hot electrical load data of each region;
s2: establishing a comprehensive demand response load model and a cost and income model;
s3: constructing an energy flow topology and an energy coupling matrix of a multi-region comprehensive energy system;
s4: and establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering the inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target.
Further, in step S1, dividing the entire area to be studied into a plurality of areas according to an administrative region location division method or a natural geographical location division method, where each area has a corresponding area-level integrated energy system;
and respectively acquiring typical solar wind luminous output active power statistical data and cold and hot electrical load data of each region arranged according to a time sequence, wherein the minimum adjacent data time interval is 1 hour.
Further, in step S2, a comprehensive demand response load model and a cost and revenue model considering various adjustable loads are established;
in order to adapt to the overall operation requirements of the system, the comprehensive demand response load models are distinguished according to different user types, wherein for pure electric users, the comprehensive demand response load models comprise transferable load models and reducible load models, and for comprehensive energy users, the comprehensive demand response load models comprise transferable load models, reducible load models and flexible load models;
the transferable load model and the reducible load model are specifically as follows:
Figure BDA0003469572300000021
wherein the content of the first and second substances,
Figure BDA0003469572300000022
representing the ith area-level comprehensive energy system, and in a time period t, adjusting the power of a load y of an energy x, wherein x belongs to { e g h c }, e, g, h and c are electric energy, gas energy, heat energy and cold energy respectively, y belongs to { mov cut }, and mov and cut are transferable loads and reducible loads respectively;
Figure BDA0003469572300000031
the maximum power adjustment proportionality coefficient of the load y of the energy x is expressed by the ith area level comprehensive energy system in the period t;
Figure BDA0003469572300000032
representing the load demand of the ith regional comprehensive energy system on the energy x in the period t;
the flexible load model is specifically as follows:
Figure BDA0003469572300000033
wherein the content of the first and second substances,
Figure BDA0003469572300000034
and
Figure BDA0003469572300000035
respectively representing the flexible power adjustment amount and the flexible power adjustment upper and lower limits of the ith area level comprehensive energy system to the energy x in a time period t; a. thex,flexThe fixed correlation quantity of the energy source x is a fixed value and is determined by the characteristics of different energy sources;
Figure BDA0003469572300000036
and
Figure BDA0003469572300000037
respectively, the i-th regional comprehensive energy system is characterized in that the variation related quantity of the energy x in the t period is an indeterminate value and is determined by the characteristics of different energy sources.
Further, the cost model is specifically as follows:
Figure BDA0003469572300000038
wherein, I is the total number of the regional comprehensive energy systems, and delta t is the duration of an optimized dispatching cycle; fCOMCompensating costs for the user side of the integrated demand response; ce,cut、Ce,mov、Ch,flexAnd Cc,flexThe unit power compensation cost coefficients of the electric load, the transferable electric load, the flexible heat load and the flexible cold load can be reduced respectively; the cost model indicates that when the multi-region comprehensive energy system implements comprehensive demand response, users participating in the comprehensive demand response make contribution to system operation, and operators should give certain economic compensation to the users participating in the comprehensive demand response according to the adjustment amount of the load;
the revenue model is specifically as follows:
Figure BDA0003469572300000039
wherein, FIDRCompensating for revenue for peak clipping response of the integrated demand response;
Figure BDA00034695723000000310
a response state quantity of a t period, wherein the value of the response state quantity is 1, the period is a response period, and the value of the response state quantity is 0, the period is a non-response period; deltae,IDR,cpA compensation price coefficient for the power peak clipping demand response; revenue model representation as a multi-regional integrated energy sourceWhen the system implements comprehensive demand response, the operator acts scattered power users in each region to participate in power peak clipping demand response, so that certain response compensation is obtained.
Further, in step S3, constructing an energy flow topology and an energy coupling matrix of the multi-region integrated energy system having an energy flow direction of 'source side output-input end energy node-energy hub-output end energy node-terminal load';
the multi-region comprehensive energy system energy flow topology comprises an inter-region comprehensive energy system connection mode topology and a region-level comprehensive energy system connection mode topology;
in the inter-regional comprehensive energy system connection mode topology, every two regional comprehensive energy systems are connected through an electric power connecting line, and each regional comprehensive energy system has electric energy interaction with an external power grid; fan photovoltaic installations with certain capacity are arranged in each regional comprehensive energy system to meet the basic electric load demand, and the surplus wind and light are sold to an external power grid according to the local wind and light internet price; meanwhile, each regional comprehensive energy system can purchase energy from an external natural gas source and a heat source to meet the multi-energy load requirement, and the energy purchase prices are all unified pricing;
in the connection mode topology of the regional-level comprehensive energy system, an energy flow topology of the regional-level comprehensive energy system is constructed according to the energy flow direction of 'source side output-input end energy node-energy concentrator-output end energy node-terminal load'; electric energy, gas energy and heat energy obtained by each regional comprehensive energy system from the source side are firstly concentrated and summarized at an input end energy node, then are subjected to multi-energy coupling conversion through an energy concentrator, and are conveyed to a terminal load after being concentrated and summarized at an output end energy node, so that the energy demand diversity of the load side is met.
The energy coupling matrix CmnThe general form of (a) is specifically as follows:
Figure BDA0003469572300000041
energy coupling matrix CmnThe transformation relation among four energy forms of electricity, gas, cold and heat in the energy concentrator of the comprehensive energy system of each regional level is generally described; wherein the content of the first and second substances,
Figure BDA0003469572300000042
and
Figure BDA0003469572300000043
supplying electric energy power, cold energy power, heat energy power and natural gas power of terminal loads of the ith regional comprehensive energy system in a time period t;
Figure BDA0003469572300000044
and
Figure BDA0003469572300000045
electric energy power, natural gas power and heat energy power obtained by the ith regional comprehensive energy system from the source side of the ith regional comprehensive energy system in a time period t; cmnFor coupling matrix C to energy sourcemnThe coupling factor in (1) describes the conversion relation of converting m energy into n energy in the ith area-level comprehensive energy system, and m and n are E { e g h c };
wherein the terminal load matrix
Figure BDA0003469572300000051
Source side output matrix
Figure BDA0003469572300000052
Pout＝CmnPin。
Further, in step S4, establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering inter-region electric energy interaction with the minimum total system scheduling cost as an optimization objective;
the multi-region integrated energy system represents a combination of a plurality of regional integrated energy systems; the optimization target of the multi-region comprehensive energy system collaborative optimization scheduling model is that the total scheduling cost of the multi-region comprehensive energy system is the minimum, and the constraint conditions of the multi-region comprehensive energy system collaborative optimization scheduling model comprise power balance constraint, unit operation constraint, tie line interaction power constraint and source side capacity constraint;
the objective function of the multi-region comprehensive energy system collaborative optimization scheduling model is specifically as follows:
Figure BDA0003469572300000059
wherein, FMRIESTotal dispatch cost for multi-zone integrated energy system, FOPEThe daily operating cost of the multi-region comprehensive energy system; fPROCost for external energy purchase; fPURThe cost of purchasing electricity for the network side; fPUNPenalizing costs for load shedding.
Further, the power balance constraint comprises a source side power balance constraint, a load side power balance constraint and an energy coupling link power balance constraint;
the source side power balance constraint considers a source side output end, including wind-solar power generation, network side electricity purchasing and selling, other area electricity supply and utilization, heat source heat production and gas source gas production, and specifically comprises the following steps:
Figure BDA0003469572300000053
wherein the content of the first and second substances,
Figure BDA0003469572300000054
and
Figure BDA0003469572300000055
respectively representing wind power and photoelectricity consumed by the ith regional comprehensive energy system in the time period t;
Figure BDA0003469572300000056
net purchased electric power at the grid side for the ith zone during the time period t,
Figure BDA0003469572300000057
meaning the purchase of electricity from the electricity grid,
Figure BDA0003469572300000058
indicating selling electricity to the grid;
because the user participates in the IDR, the terminal load not only considers the electric cold and hot basic load, but also considers the adjustable load including reducible and transferable electric load and flexible cold and hot load; thereby, the terminal load matrix P in the formula (5) is formedoutExpressed as the load-side power balance constraint, the following is specific:
Figure BDA0003469572300000061
wherein the content of the first and second substances,
Figure BDA0003469572300000062
and
Figure BDA0003469572300000063
respectively providing a basic cold load, a basic heat load and a basic gas load of the ith area-level comprehensive energy system in a t period;
Figure BDA0003469572300000064
and
Figure BDA0003469572300000065
flexible cold load and flexible heat load of the ith area level comprehensive energy system in a time period t are respectively;
the energy coupling link power balance constraint comprises an energy coupling link electric energy balance constraint, an energy coupling link heat energy balance constraint, an energy coupling link cold energy balance constraint and an energy coupling ring solar energy balance constraint;
the electric energy balance constraint of the energy coupling link is as follows:
Figure BDA0003469572300000066
the energy coupling link has the following heat energy balance constraints:
Figure BDA0003469572300000067
the cold energy balance constraint of the energy coupling link is as follows:
Figure BDA0003469572300000068
the energy coupling ring solar energy balance constraint is as follows:
Figure BDA0003469572300000069
wherein LHV is the low heating value of natural gas, also called net heating value, which represents the heat released by complete combustion of unit volume of natural gas without calculating the heat of vaporization, and has a value of 9.7 x 10-3MWh/Nm 3.
Further, the unit operation constraints include gas turbine operation constraints and gas boiler operation constraints;
the gas turbine operating constraints are specifically as follows:
Figure BDA00034695723000000610
wherein the content of the first and second substances,
Figure BDA00034695723000000611
and
Figure BDA00034695723000000612
respectively representing the upper limit and the lower limit of the power generation power of the gas turbine of the ith regional comprehensive energy system;
Figure BDA00034695723000000613
the starting and stopping variable of the gas turbine of the ith area-level comprehensive energy system in the t period is represented, the starting state is represented when the value of the starting and stopping variable is 1, and the stopping state is represented when the value of the starting and stopping variable is 0; r isMT,uAnd rMT,dAs a gasUpward and downward climbing coefficients of the turbine;
Figure BDA00034695723000000614
rated power of the gas turbine of the ith regional level integrated energy system;
the operation constraints of the gas boiler are as follows:
Figure BDA0003469572300000071
wherein the content of the first and second substances,
Figure BDA0003469572300000072
respectively representing the upper and lower limits of the heat production power of the gas-fired boiler of the ith regional comprehensive energy system;
Figure BDA0003469572300000073
the starting and stopping variable of the gas turbine of the ith area-level comprehensive energy system in the t period is represented, the starting state is represented when the value of the starting and stopping variable is 1, and the stopping state is represented when the value of the starting and stopping variable is 0; r isGB,uAnd rGB,dUpward and downward climbing coefficients for the gas turbine;
Figure BDA0003469572300000074
the rated power of the gas turbine of the ith regional level integrated energy system.
Further, the tie line interaction power constraint is specifically as follows:
Figure BDA0003469572300000075
wherein the content of the first and second substances,
Figure BDA0003469572300000076
and
Figure BDA0003469572300000077
respectively representing the upper and lower limits of the interactive power of the i-th regional comprehensive energy system and a tie line of a large power grid in a t period;
Figure BDA0003469572300000078
and
Figure BDA0003469572300000079
and respectively representing the upper and lower limits of the junctor interaction power of the ith area-level comprehensive energy system with other areas in the t period.
Further, the source side capacity constraint is specifically as follows:
Figure BDA00034695723000000710
wherein the content of the first and second substances,
Figure BDA00034695723000000711
and
Figure BDA00034695723000000712
respectively representing the upper limit and the lower limit of heat source heat production power of the ith regional comprehensive energy system in a time period t;
Figure BDA00034695723000000713
and
Figure BDA00034695723000000714
respectively representing the upper and lower limits of the gas source gas production rate of the ith area in the t period.
Compared with the prior art, the invention has the advantages that:
the invention considers various load types of demand response, has a large coverage area, can effectively embody the complementary characteristics of electricity, gas, cold and heat in space and time, breaks the limitation of the traditional method in both time scale and capacity scale, enhances the coupling characteristic of the interconnected system, improves the integral response capability of the interconnected system, and improves the economy, flexibility and reliability of the interconnected system.
Drawings
FIG. 1 is a schematic flow chart illustrating a method for scheduling a multi-zone integrated energy system with consideration of integrated demand response according to an embodiment of the present invention;
FIG. 2 is a topological diagram of connection modes of inter-regional integrated energy systems according to an embodiment of the present invention;
fig. 3 is a topological diagram of a connection mode of the regional-level integrated energy system according to an embodiment of the present invention.
Detailed Description
The technical scheme of the application is further described in detail through specific embodiments by combining the drawings in the specification.
Example 1:
a method for scheduling a multi-region integrated energy system considering integrated demand response, as shown in fig. 1, includes the following steps:
s1: dividing a region to be researched into a plurality of regions according to the geographic position, and acquiring wind-solar output data and cold-hot electrical load data of each region;
dividing the whole area to be researched into a plurality of areas according to an administrative region position division method or a natural geographical position division method, wherein each area is provided with a corresponding area-level comprehensive energy system;
and respectively acquiring typical solar wind luminous output active power statistical data and cold and hot electrical load data of each region arranged according to a time sequence, wherein the minimum adjacent data time interval is 1 hour.
S2: establishing a comprehensive demand response load model and a cost and income model;
establishing a comprehensive demand response load model and a cost and income model considering various adjustable loads;
in order to adapt to the overall operation requirements of the system, the comprehensive demand response load models are distinguished according to different user types, wherein for pure electric users, the comprehensive demand response load models comprise transferable load models and reducible load models, and for comprehensive energy users, the comprehensive demand response load models comprise transferable load models, reducible load models and flexible load models;
the transferable load model and the reducible load model are specifically as follows:
Figure BDA0003469572300000081
wherein the content of the first and second substances,
Figure BDA0003469572300000082
representing the ith area-level comprehensive energy system, and in a time period t, adjusting the power of a load y of an energy x, wherein x belongs to { e g h c }, e, g, h and c are electric energy, gas energy, heat energy and cold energy respectively, y belongs to { mov cut }, and mov and cut are transferable loads and reducible loads respectively;
Figure BDA0003469572300000083
the maximum power adjustment proportionality coefficient of the load y of the energy x is expressed in the ith region level comprehensive energy system in the t period;
Figure BDA0003469572300000084
representing the load demand of the ith regional comprehensive energy system on the energy x in the period t;
the flexible load model is specifically as follows:
Figure BDA0003469572300000085
wherein the content of the first and second substances,
Figure BDA0003469572300000086
and
Figure BDA0003469572300000087
respectively representing the flexible power adjustment amount and the flexible power adjustment upper and lower limits of the ith area level comprehensive energy system to the energy x in a time period t; a. thex,flexThe fixed related quantity of the energy x is a fixed value and is determined by the characteristics of different energy sources, for example, the fixed related quantity of the hot water load is the specific heat capacity and density of water, and the fixed related quantity of the air conditioning load is the thermal resistance of a house;
Figure BDA0003469572300000088
and
Figure BDA0003469572300000091
and respectively representing the ith regional comprehensive energy system, wherein the variation related quantity of the energy x is an uncertain value in a period t and is determined by the characteristics of different energy sources, for example, the variation related quantity of the hot water load is the volume of hot water required to be heated in different periods of different regions, and the variation related quantity of the air conditioning load is the number of cooling rooms and the outdoor temperature in different periods of different regions.
The cost model is specifically as follows:
Figure BDA0003469572300000092
wherein, I is the total number of the regional comprehensive energy systems, and delta t is the duration of an optimized dispatching cycle; fCOMCompensating costs for the user side of the integrated demand response; ce,cut、Ce,mov、Ch,flexAnd Cc,flexThe unit power compensation cost coefficients of the electric load, the transferable electric load, the flexible heat load and the flexible cold load can be reduced respectively; the cost model indicates that when the multi-region comprehensive energy system implements comprehensive demand response, users participating in the comprehensive demand response make contribution to system operation, and operators should give certain economic compensation to the users participating in the comprehensive demand response according to the adjustment amount of the load;
the revenue model is specifically as follows:
Figure BDA0003469572300000093
wherein, FIDRCompensating for revenue for peak clipping response of the integrated demand response;
Figure BDA0003469572300000094
a response state quantity of a t period, wherein the value of the response state quantity is 1, the period is a response period, and the value of the response state quantity is 0, the period is a non-response period; deltae,IDR,cpCompensating prices for power peak clipping demand responseA coefficient; the income model shows that when the multi-region comprehensive energy system implements comprehensive demand response, the operator acts scattered power users in each region to participate in power peak clipping demand response, so that certain response compensation is obtained.
S3: constructing an energy flow topology and an energy coupling matrix of a multi-region comprehensive energy system;
constructing a multi-region comprehensive energy system energy flow topology and an energy coupling matrix thereof, wherein the energy flow direction is 'source side output-input end energy node-energy concentrator-output end energy node-terminal load';
the multi-region comprehensive energy system energy flow topology comprises an inter-region comprehensive energy system connection mode topology and a region-level comprehensive energy system connection mode topology;
as shown in fig. 2, in the inter-regional integrated energy system connection topology, the regional integrated energy systems are connected in pairs by power connecting lines, and each regional integrated energy system has electric energy interaction with an external power grid; fan photovoltaic installations with certain capacity are arranged in each regional comprehensive energy system to meet the basic electric load demand, and the surplus wind and light are sold to an external power grid according to the local wind and light internet price; meanwhile, each regional comprehensive energy system can purchase energy from an external natural gas source and a heat source to meet the multi-energy load requirement, and the energy purchase prices are all unified pricing;
as shown in fig. 3, in the topology of the connection mode of the regional-level integrated energy system, an energy flow topology of the regional-level integrated energy system is constructed according to an energy flow direction of 'source side output-input end energy node-energy concentrator-output end energy node-terminal load'; electric energy, gas energy and heat energy obtained by each regional comprehensive energy system from the source side are firstly concentrated and summarized at an input end energy node, then are subjected to multi-energy coupling conversion through an energy concentrator, and are conveyed to a terminal load after being concentrated and summarized at an output end energy node, so that the energy demand diversity of the load side is met; in the context of figure 3, it is shown,
Figure BDA0003469572300000101
and
Figure BDA0003469572300000102
directly supplying electric energy and gas energy of a terminal load to the ith regional comprehensive energy system without an energy coupling link in a period t;
Figure BDA0003469572300000103
Figure BDA0003469572300000104
and
Figure BDA0003469572300000105
the heating power of an electric heating mechanism, the refrigerating power of an electric refrigerator, the gas making power of an electric gas conversion device, the gas making power of a gas turbine, the heating power of a gas boiler and the refrigerating power of an endothermic refrigerator of the ith area-level comprehensive energy system in the period t are respectively.
The energy coupling matrix CmnThe general form of (a) is specifically as follows:
Figure BDA0003469572300000106
energy coupling matrix CmnThe transformation relation among four energy forms of electricity, gas, cold and heat in the energy concentrator of the comprehensive energy system of each regional level is generally described; wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003469572300000107
and
Figure BDA0003469572300000108
supplying electric energy power, cold energy power, heat energy power and natural gas power of terminal loads of the ith regional comprehensive energy system in a t period;
Figure BDA0003469572300000109
and
Figure BDA00034695723000001010
are respectively asThe electric energy power, the natural gas power and the heat energy power are obtained by the ith area level comprehensive energy system from the source side of the ith area level comprehensive energy system in a t period; cmnFor coupling matrix C to energy sourcemnThe coupling factor in (1) describes the conversion relation of m energy in the ith area level comprehensive energy system to n energy, and m and n are belonged to { e g h c };
wherein the terminal load matrix
Figure BDA00034695723000001011
Source-side output matrix
Figure BDA00034695723000001012
Pout＝CmnPin。
In this embodiment, the specific form of the energy coupling matrix is as follows:
Figure BDA0003469572300000111
energy coupling matrix CmnThe specific form of (A) is established based on the general form thereof in cooperation with the connection mode topology of the regional-level comprehensive energy system shown in FIG. 3; aiming at the problem of conversion among different energy forms in the energy concentrator of the comprehensive energy system of each regional level, 8 energy distribution coefficients which change along with time are introduced: the natural gas heat recovery system comprises a first electric energy distribution coefficient a1, a second electric energy distribution coefficient a2, a third electric energy distribution coefficient a3 and a fourth electric energy distribution coefficient a4, a first natural gas distribution coefficient b1, a second natural gas distribution coefficient b2, a first heat energy distribution coefficient c1 and a second heat energy distribution coefficient c2, wherein the sum of the energy distribution coefficients is 1; the electric heating machine, the electric refrigerating machine and the electric gas conversion equipment can respectively convert partial electric energy of an electric energy node at an input end into heat energy, cold energy and gas energy; the micro gas turbine and the gas boiler can respectively convert part of gas energy of the heat energy node at the input end into electric energy and heat energy; the heat absorption refrigerating machine can convert the surplus heat energy gathered by the waste heat boiler into cold energy; wherein μ in the formula (6)EH、μEC、μP2G、μMT、μGBAnd muARRespectively for electric heatingMachine heating efficiency, electric refrigerator refrigeration efficiency, electric gas conversion equipment gas production efficiency, gas turbine power generation efficiency, gas boiler heating efficiency and heat absorption refrigerator refrigeration efficiency.
S4: establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target;
establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target;
the multi-region integrated energy system represents a combination of a plurality of regional integrated energy systems; the optimization target of the multi-region comprehensive energy system collaborative optimization scheduling model is that the total scheduling cost of the multi-region comprehensive energy system is the minimum, and the constraint conditions of the multi-region comprehensive energy system collaborative optimization scheduling model comprise power balance constraint, unit operation constraint, tie line interaction power constraint and source side capacity constraint;
the objective function of the multi-region comprehensive energy system collaborative optimization scheduling model is as follows:
Figure BDA0003469572300000112
wherein, FMRIESTotal dispatch cost for multi-zone integrated energy system, FOPEThe daily operating cost of the multi-region comprehensive energy system; fPROCost for external energy purchase; fPURThe cost of purchasing electricity for the network side; fPUNPenalizing costs for load shedding.
Further, the power balance constraint comprises a source side power balance constraint, a load side power balance constraint and an energy coupling link power balance constraint;
the source side power balance constraint considers a source side output end, including wind-solar power generation, network side electricity purchasing and selling, other area electricity supply and utilization, heat source heat production and gas source gas production, and specifically comprises the following steps:
Figure BDA0003469572300000121
wherein the content of the first and second substances,
Figure BDA0003469572300000122
and
Figure BDA0003469572300000123
respectively representing wind power and photoelectricity consumed by the ith regional comprehensive energy system in the time period t;
Figure BDA0003469572300000124
net purchased electric power at the grid side for the ith zone during the time period t,
Figure BDA0003469572300000125
meaning the purchase of electricity from the electricity grid,
Figure BDA0003469572300000126
indicating selling electricity to the grid;
because the user participates in the comprehensive demand side response IDR, the terminal load in FIG. 2 not only considers the electric cold and hot basic load, but also considers the adjustable load including reducible and transferable electric load and flexible cold and hot load; thereby, the terminal load matrix P in the formula (5) is formedoutExpressed as the load-side power balance constraint, the following is specific:
Figure BDA0003469572300000127
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003469572300000128
and
Figure BDA0003469572300000129
respectively providing a basic cold load, a basic heat load and a basic gas load of the ith area-level comprehensive energy system in a t period;
Figure BDA00034695723000001210
and
Figure BDA00034695723000001211
flexible cold load and flexible heat load of the ith area level comprehensive energy system in a time period t are respectively;
the energy coupling link power balance constraint comprises an energy coupling link electric energy balance constraint, an energy coupling link heat energy balance constraint, an energy coupling link cold energy balance constraint and an energy coupling ring solar energy balance constraint;
the electric energy balance constraint of the energy coupling link is as follows:
Figure BDA00034695723000001212
the energy coupling link has the following heat energy balance constraints:
Figure BDA00034695723000001213
the cold energy balance constraint of the energy coupling link is as follows:
Figure BDA00034695723000001214
the energy coupling ring solar energy balance constraint is as follows:
Figure BDA00034695723000001215
wherein LHV is the low heating value of natural gas, also called net heating value, which represents the heat released by complete combustion of unit volume of natural gas without calculating the heat of vaporization, and has a value of 9.7 x 10-3MWh/Nm 3.
The unit operation constraints comprise gas turbine operation constraints and gas boiler operation constraints;
the gas turbine operating constraints are specifically as follows:
Figure BDA0003469572300000131
wherein the content of the first and second substances,
Figure BDA0003469572300000132
and
Figure BDA0003469572300000133
respectively representing the upper limit and the lower limit of the power generation power of the gas turbine of the ith regional comprehensive energy system;
Figure BDA0003469572300000134
the starting and stopping variable of the gas turbine of the ith area-level comprehensive energy system in the t period is represented, the starting state is represented when the value of the starting and stopping variable is 1, and the stopping state is represented when the value of the starting and stopping variable is 0; r isMT,uAnd rMT,dUpward and downward climbing coefficients for the gas turbine;
Figure BDA0003469572300000135
rated power of the gas turbine of the ith regional level integrated energy system;
the operation constraints of the gas boiler are as follows:
Figure BDA0003469572300000136
wherein the content of the first and second substances,
Figure BDA0003469572300000137
respectively representing the upper and lower limits of the heat production power of the gas-fired boiler of the ith regional comprehensive energy system;
Figure BDA0003469572300000138
the starting and stopping variables of the gas turbine of the ith regional comprehensive energy system in the t period are represented, the starting state is represented when the value of the starting and stopping variables is 1, and the stopping state is represented when the value of the starting and stopping variables is 0; r isGB,uAnd rGB,dUpward and downward climbing coefficients for the gas turbine;
Figure BDA0003469572300000139
the rated power of the gas turbine of the ith regional level integrated energy system.
The tie line interaction power constraint is specifically as follows:
Figure BDA00034695723000001310
wherein the content of the first and second substances,
Figure BDA00034695723000001311
and
Figure BDA00034695723000001312
respectively representing the upper and lower limits of the interactive power of the i-th regional comprehensive energy system and a tie line of a large power grid in the time period t;
Figure BDA00034695723000001313
and
Figure BDA00034695723000001314
respectively representing the upper and lower limits of the connecting line interactive power of the ith area-level comprehensive energy system with other areas in the t period;
the source side capacity constraint is as follows:
Figure BDA00034695723000001315
wherein the content of the first and second substances,
Figure BDA00034695723000001316
and
Figure BDA00034695723000001317
respectively representing the upper and lower limits of heat source heat production power of the i-th regional comprehensive energy system in a t period;
Figure BDA00034695723000001318
and
Figure BDA00034695723000001319
respectively representing the upper and lower limits of the gas source gas production rate of the ith zone in the t period.
In this embodiment, to verify the above described collaborative optimization scheduling model of multi-region integrated energy system, 3 region-level integrated energy systems belonging to the same power distribution region in a certain area in shanxi are connected by a power line to form a multi-region integrated energy system as a research object.
Three areas are arranged, the first area 1, the second area 2 and the third area 3 are mutually connected in pairs, and the cooperation of the electric energy between the areas can be realized through an electric power connecting line. Energy hubs are arranged in each region, and when supply and demand on two sides of source loads in each region are not matched, multi-energy coupling cooperative optimization can be performed through the energy hubs. In order to make the calculation result universal, aiming at the whole-year cold and warm supply arrangement of Shanxi province (the summer cooling period is from 3 months and 15 days to 11 months and 15 days in the next year, and the winter heating period is from 11 months and 15 days to 3 months and 15 days in the next year), the article plans and selects 2 types of load typical days which are respectively the summer cooling period and the winter heating period. And selecting an optimized scheduling cycle T as 1 day, wherein the operation time interval delta T of one optimized scheduling cycle is 1 h. And calling a Cplex solver by combining matlabR2016b software and a Yalmip plug-in to solve.
In this embodiment, in order to respectively measure the influence of each regional user participating in the comprehensive demand response, the multi-energy coupling conversion in each region, and the coordination and mutual coordination of the electric energy between the regions on the total scheduling cost, embodiment 1 is set as shown in table 1.
TABLE 1
Figure BDA0003469572300000141
By using the example 1 in table 1, the optimal scheduling is performed on the multi-region integrated energy system in typical days of winter and summer, and the results are shown in table 2.
Table 2 optimization results of example 1
Figure BDA0003469572300000142
Note: all data in table 2 are the sum of the optimization results for 3 regions.
Embodiment 1 reconfigures resources by encouraging users in each area to participate in comprehensive demand response, multi-energy coupling conversion in each area, and coordination of electric energy among areas, so that supply and demand on two sides of electric, gas, cold and heat source loads of each area are matched. From the optimization result, the embodiment 1 generally reduces the daily operation cost and the load abandonment punishment cost of the system, reduces the external energy purchase and the network side electricity purchase, and obtains a certain peak clipping response compensation income by participating in the demand response market, thereby obviously reducing the total scheduling cost of the multi-region comprehensive energy system and effectively improving the overall economy of the system. In addition, peak clipping and valley filling are realized by the participation of multi-energy adjustable load in comprehensive demand response, and the overall energy utilization pressure is relieved; the load complementary characteristics among different regions and in each region are fully utilized, the electric energy mutual-aid capability and the multi-energy coupling conversion capability of the multi-region comprehensive energy system are improved, the supply and demand of each region are balanced, the reasonable allocation of resources is realized, and the overall flexibility and reliability of the system are improved.
Example 2:
the present embodiment is different from embodiment 1 in that each regional user does not participate in the integrated demand response, as shown in table 3.
TABLE 3
Figure BDA0003469572300000151
By using example 2 of table 3, optimal scheduling is performed on the multi-region integrated energy system in typical days of winter and summer, and the results are shown in table 4.
Table 4 optimization results of example 2
Figure BDA0003469572300000152
Note: all data in table 2 are the sum of the optimization results for 3 regions.
Comparing example 1 and example 2 in table 2 and table 4 at typical days of winter and summer, it can be seen whether each regional user participates in the cost impact of the overall demand response. Compared with the embodiment 2, the operator agent can transfer the electric load and can reduce the electric load to participate in external power peak clipping demand response by encouraging the user to participate in IDR, and guide the flexible cold and hot loads to be redistributed in the comfortable adjustment range of the user, thereby not only reducing the external electricity purchasing quantity, but also obviously reducing the electricity purchasing cost on the network side (136.1 ten thousand yuan is reduced in summer and 137.3 ten thousand yuan is reduced in winter); also, a large amount of response compensation revenue is obtained by participating in the demand response market; meanwhile, the 0 load abandoning rate is realized, and the load abandoning penalty cost is effectively reduced (68.3 ten thousand yuan is reduced in summer and 20.4 ten thousand yuan is reduced in winter); although additional user-side compensation costs need to be paid, embodiment 1 is significantly superior to embodiment 2 in terms of total scheduling costs, which are a drop of 1434.8 ten thousand yuan in summer and a drop of 1236.4 ten thousand yuan in winter.
Example 3:
the difference between this embodiment and embodiment 1 is that the multi-regional integrated energy system does not achieve coordination and coordination of electric energy between the regions, and the power links between the regions are substantially disconnected, as shown in table 5.
TABLE 5
Figure BDA0003469572300000161
By using example 3 of table 5, the result of performing optimal scheduling on the multi-region integrated energy system in typical days of winter and summer is shown in table 6.
Table 6 optimization results of example 3
Figure BDA0003469572300000162
Note: all data in table 2 are the sum of the optimization results for 3 regions.
Comparing example 1 and example 3 in typical days of winter and summer in table 2 and table 6, the influence on cost of whether or not electric energy coordination and complementation can be performed between the regions can be obtained. Compared with embodiment 3, in embodiment 1, when the supply and demand of electric power at two sides of the source load are not matched in a certain period of time in a certain region, surplus electric energy is preferentially called from other regions for support, so that the supply and demand sides are balanced, the abandoned load amount is reduced, and the abandoned load punishment cost is remarkably reduced (321.6 ten thousand yuan is reduced in summer and 120.1 ten thousand yuan is reduced in winter); the adjustment amount of the adjustable load is reduced, and the compensation cost of the user side is reduced (6.4 ten thousand yuan is reduced in summer and 54.6 ten thousand yuan is reduced in winter); although the cost of purchasing and selling electricity on the network side is increased (70.3 ten thousand yuan is increased in summer and 9.6 ten thousand yuan is increased in winter), 287.7 ten thousand yuan is reduced in summer and 68.1 ten thousand yuan is reduced in winter from the perspective of the total scheduling cost, and the embodiment 1 is obviously superior to the embodiment 3.

Claims (10)

1. A multi-region comprehensive energy system scheduling method considering comprehensive demand response is characterized by comprising the following steps:
s1: dividing a region to be researched into a plurality of regions according to the geographic position, and acquiring wind-solar output data and cold-hot electrical load data of each region;
s2: establishing a comprehensive demand response load model and a cost and income model;
s3: constructing an energy flow topology and an energy coupling matrix of a multi-region comprehensive energy system;
s4: and establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering the inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target.
2. The method of claim 1, wherein in step S1, the whole area to be studied is divided into a plurality of areas according to an administrative division or a natural geographic division, and each area has a corresponding area-level integrated energy system;
and respectively acquiring typical solar wind luminous output active power statistical data and cold and hot electrical load data of each region arranged according to a time sequence, wherein the minimum adjacent data time interval is 1 hour.
3. The method for dispatching a multi-regional integrated energy system considering integrated demand response according to claim 1, wherein in step S2, an integrated demand response load model considering a plurality of adjustable loads and a cost and income model are established;
in order to adapt to the overall operation requirements of the system, the comprehensive demand response load models are distinguished according to different user types, wherein for pure electric users, the comprehensive demand response load models comprise transferable load models and reducible load models, and for comprehensive energy users, the comprehensive demand response load models comprise transferable load models, reducible load models and flexible load models;
the transferable load model and the reducible load model are specifically as follows:
Figure FDA0003469572290000011
wherein the content of the first and second substances,
Figure FDA0003469572290000012
representing the ith area-level comprehensive energy system, and in a time period t, adjusting the power of a load y of an energy x, wherein x belongs to { e g h c }, e, g, h and c are electric energy, gas energy, heat energy and cold energy respectively, y belongs to { mov cut }, and mov and cut are transferable loads and reducible loads respectively;
Figure FDA0003469572290000013
the maximum power adjustment proportionality coefficient of the load y of the energy x is expressed by the ith area level comprehensive energy system in the period t;
Figure FDA0003469572290000014
representing the load demand of the ith regional comprehensive energy system on the energy x in the period t;
the flexible load model is specifically as follows:
Figure FDA0003469572290000015
wherein the content of the first and second substances,
Figure FDA0003469572290000021
and
Figure FDA0003469572290000022
respectively representing the flexible power adjustment amount and the flexible power adjustment upper and lower limits of the ith area level comprehensive energy system to the energy x in a time period t; a. thex,flexThe fixed correlation quantity of the energy source x is a fixed value and is determined by the characteristics of different energy sources;
Figure FDA0003469572290000023
and
Figure FDA0003469572290000024
respectively, the i-th regional comprehensive energy system is characterized in that the variation related quantity of the energy x in the t period is an indeterminate value and is determined by the characteristics of different energy sources.
4. The method according to claim 3, wherein the cost model is as follows:
Figure FDA0003469572290000025
wherein, I is the total number of the regional comprehensive energy systems, and delta t is the duration of an optimized dispatching cycle; fCOMCompensating costs for the user side of the integrated demand response; ce,cut、Ce,mov、Ch,flexAnd Cc,flexCompensation for unit power of reducible electric load, transferable electric load, flexible heat load and flexible cold loadThe coefficient;
the revenue model is specifically as follows:
Figure FDA0003469572290000026
wherein, FIDRCompensating for revenue for peak clipping response of the integrated demand response;
Figure FDA0003469572290000027
a response state quantity of a t period, wherein the value of the response state quantity is 1, the period is a response period, and the value of the response state quantity is 0, the period is a non-response period; deltae,IDR,cpAnd compensating the price coefficient for the power peak clipping demand response.
5. The method for dispatching multi-region integrated energy system considering integrated demand response of claim 1, wherein in step S3, energy flow topology and energy coupling matrix of multi-region integrated energy system with energy flow direction being 'source side output-input energy node-energy hub-output energy node-terminal load' are constructed;
the multi-region comprehensive energy system energy flow topology comprises an inter-region comprehensive energy system connection mode topology and a region-level comprehensive energy system connection mode topology;
in the inter-regional comprehensive energy system connection mode topology, every two regional comprehensive energy systems are connected through an electric power connecting line, and each regional comprehensive energy system has electric energy interaction with an external power grid; fan photovoltaic installations are arranged inside each regional comprehensive energy system to meet the basic electric load requirements, and surplus wind and light are sold to an external power grid according to local wind and light internet prices; meanwhile, each regional comprehensive energy system can purchase energy from an external natural gas source and a heat source to meet the multi-energy load requirement, and the energy purchase prices are all unified pricing;
in the connection mode topology of the regional-level comprehensive energy system, an energy flow topology of the regional-level comprehensive energy system is constructed according to the energy flow direction of 'source side output-input end energy node-energy concentrator-output end energy node-terminal load'; electric energy, gas energy and heat energy obtained by each regional comprehensive energy system from the source side are firstly concentrated and summarized at an input end energy node, then are subjected to multi-energy coupling conversion through an energy concentrator, and are conveyed to a terminal load after being concentrated and summarized at an output end energy node, so that the energy demand diversity of the load side is met;
the energy coupling matrix CmnThe general form of (a) is specifically as follows:
Figure FDA0003469572290000031
energy coupling matrix CmnThe transformation relation among four energy forms of electricity, gas, cold and heat in the energy concentrator of the comprehensive energy system of each regional level is generally described; wherein the content of the first and second substances,
Figure FDA0003469572290000032
and
Figure FDA0003469572290000033
supplying electric energy power, cold energy power, heat energy power and natural gas power of terminal loads of the ith regional comprehensive energy system in a t period;
Figure FDA0003469572290000034
and
Figure FDA0003469572290000035
electric energy power, natural gas power and heat energy power obtained by the ith regional comprehensive energy system from the source side of the ith regional comprehensive energy system in a time period t; cmnFor coupling matrix C to energy sourcemnThe coupling factor in (1) describes the conversion relation of m energy in the ith area level comprehensive energy system to n energy, and m and n are belonged to { e g h c };
wherein the terminal load matrix
Figure FDA0003469572290000036
Source side output matrix
Figure FDA0003469572290000037
Pout＝CmnPin。
6. The multi-region integrated energy system dispatching method considering integrated demand response according to any one of claims 1 to 5, wherein: in the step S4, establishing a multi-region comprehensive energy system collaborative optimization scheduling model considering inter-region electric energy interaction by taking the minimum total scheduling cost of the system as an optimization target;
the multi-region integrated energy system represents a combination of a plurality of region-level integrated energy systems; the optimization target of the multi-region comprehensive energy system collaborative optimization scheduling model is that the total scheduling cost of the multi-region comprehensive energy system is the minimum, and the constraint conditions of the multi-region comprehensive energy system collaborative optimization scheduling model comprise power balance constraint, unit operation constraint, tie line interaction power constraint and source side capacity constraint;
the objective function of the multi-region comprehensive energy system collaborative optimization scheduling model is specifically as follows:
Figure FDA0003469572290000041
wherein, FMRIESTotal dispatch cost for multi-zone integrated energy system, FOPEThe daily operating cost of the multi-region comprehensive energy system; fPROCost for external energy purchase; fPURThe cost of purchasing electricity for the network side; fPUNPenalizing costs for load shedding.
7. The method of claim 6, wherein the method comprises: the power balance constraint comprises a source side power balance constraint, a load side power balance constraint and an energy coupling link power balance constraint;
the source side power balance constraint considers a source side output end, including wind-solar power generation, network side electricity purchasing and selling, other area electricity supply and utilization, heat source heat production and gas source gas production, and specifically comprises the following steps:
Figure FDA0003469572290000042
wherein the content of the first and second substances,
Figure FDA0003469572290000043
and
Figure FDA0003469572290000044
respectively representing wind power and photoelectricity consumed by the ith regional comprehensive energy system in the time period t;
Figure FDA0003469572290000045
net purchased electric power at the grid side for the ith zone during the time period t,
Figure FDA0003469572290000046
indicating the purchase of electricity from the power grid,
Figure FDA0003469572290000047
indicating selling electricity to the grid;
because the user participates in the IDR, the terminal load not only considers the electric cold and hot basic load, but also considers the adjustable load including reducible and transferable electric load and flexible cold and hot load; thereby, the terminal load matrix P in the formula (5) is formedoutExpressed as the load-side power balance constraint, the following is specific:
Figure FDA0003469572290000048
wherein the content of the first and second substances,
Figure FDA0003469572290000049
and
Figure FDA00034695722900000410
respectively providing a basic cold load, a basic heat load and a basic gas load of the ith area-level comprehensive energy system in a t period;
Figure FDA00034695722900000411
and
Figure FDA00034695722900000412
flexible cold load and flexible heat load of the ith area level comprehensive energy system in a time period t are respectively;
the energy coupling link power balance constraint comprises an energy coupling link electric energy balance constraint, an energy coupling link heat energy balance constraint, an energy coupling link cold energy balance constraint and an energy coupling ring solar energy balance constraint;
the electric energy balance constraint of the energy coupling link is as follows:
Figure FDA0003469572290000051
the energy coupling link has the following heat energy balance constraints:
Figure FDA0003469572290000052
the cold energy balance constraint of the energy coupling link is as follows:
Figure FDA0003469572290000053
the energy coupling ring solar energy balance constraint is as follows:
Figure FDA0003469572290000054
wherein LHV is the low heating value of natural gas, also called net heating value, which represents the heat released by complete combustion of unit volume of natural gas without calculating the heat of vaporization, and has a value of 9.7 x 10-3MWh/Nm 3.
8. The method of claim 6, wherein the method comprises: the unit operation constraint comprises gas turbine operation constraint and gas boiler operation constraint;
the gas turbine operating constraints are specifically as follows:
Figure FDA0003469572290000055
wherein the content of the first and second substances,
Figure FDA0003469572290000056
and
Figure FDA0003469572290000057
respectively representing the upper limit and the lower limit of the power generation power of the gas turbine of the ith regional comprehensive energy system;
Figure FDA0003469572290000058
the starting and stopping variable of the gas turbine of the ith area-level comprehensive energy system in the t period is represented, the starting state is represented when the value of the starting and stopping variable is 1, and the stopping state is represented when the value of the starting and stopping variable is 0; r isMT,uAnd rMT,dUpward and downward climbing coefficients for the gas turbine;
Figure FDA0003469572290000059
rated power of the gas turbine of the ith regional level integrated energy system;
the operation constraints of the gas boiler are as follows:
Figure FDA00034695722900000510
wherein the content of the first and second substances,
Figure FDA00034695722900000511
respectively representing the upper and lower limits of the heat production power of the gas-fired boiler of the ith regional comprehensive energy system;
Figure FDA00034695722900000512
the starting and stopping variable of the gas turbine of the ith area-level comprehensive energy system in the t period is represented, the starting state is represented when the value of the starting and stopping variable is 1, and the stopping state is represented when the value of the starting and stopping variable is 0; r is a radical of hydrogenGB,uAnd rGB,dUpward and downward climbing coefficients for the gas turbine;
Figure FDA00034695722900000513
the rated power of the gas turbine of the ith regional level integrated energy system.
9. The method of claim 6, wherein the method comprises: the tie line interaction power constraint is specifically as follows:
Figure FDA0003469572290000061
wherein the content of the first and second substances,
Figure FDA0003469572290000062
and
Figure FDA0003469572290000063
respectively representing the upper and lower limits of the interactive power of the i-th regional comprehensive energy system and a tie line of a large power grid in a t period;
Figure FDA0003469572290000064
and
Figure FDA0003469572290000065
and respectively representing the upper and lower limits of the junctor interaction power of the ith area-level comprehensive energy system with other areas in the t period.
10. The method of claim 6, wherein the method comprises: the source side capacity constraint is as follows:
Figure FDA0003469572290000066
wherein the content of the first and second substances,
Figure FDA0003469572290000067
and
Figure FDA0003469572290000068
respectively representing the upper and lower limits of heat source heat production power of the i-th regional comprehensive energy system in a t period;
Figure FDA0003469572290000069
and
Figure FDA00034695722900000610
respectively representing the upper and lower limits of the gas source gas production rate of the ith area in the t period.
CN202210039372.XA 2022-01-13 2022-01-13 Multi-region comprehensive energy system scheduling method considering comprehensive demand response Pending CN114529153A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210039372.XA CN114529153A (en) 2022-01-13 2022-01-13 Multi-region comprehensive energy system scheduling method considering comprehensive demand response

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210039372.XA CN114529153A (en) 2022-01-13 2022-01-13 Multi-region comprehensive energy system scheduling method considering comprehensive demand response

Publications (1)

Publication Number Publication Date
CN114529153A true CN114529153A (en) 2022-05-24

Family

ID=81619952

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210039372.XA Pending CN114529153A (en) 2022-01-13 2022-01-13 Multi-region comprehensive energy system scheduling method considering comprehensive demand response

Country Status (1)

Country Link
CN (1) CN114529153A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115423252A (en) * 2022-07-29 2022-12-02 交通运输部规划研究院 Wind and light resource supply potential evaluation method and device and computer equipment

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115423252A (en) * 2022-07-29 2022-12-02 交通运输部规划研究院 Wind and light resource supply potential evaluation method and device and computer equipment

Similar Documents

Publication Publication Date Title
CN111738502B (en) Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption
CN108229025B (en) Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system
Gu et al. Residential CCHP microgrid with load aggregator: Operation mode, pricing strategy, and optimal dispatch
Anvari-Moghaddam et al. Cost-effective and comfort-aware residential energy management under different pricing schemes and weather conditions
CN110163411A (en) A kind of regional complex energy system operation optimization method
CN111860965B (en) User comprehensive energy system optimization scheduling method considering energy storage multi-type service
CN110889600A (en) Regional comprehensive energy system optimization scheduling method considering flexible thermal load
CN110378729B (en) Comprehensive demand response method based on dynamic energy price strategy
CN113361875B (en) Optimization scheduling method for multi-microgrid comprehensive energy system considering demand side response and shared energy storage
CN107330549A (en) Multi-energy complementary micro-grid energy management method containing multiagent
CN113705906A (en) Energy coordination optimization operation method and system for comprehensive energy park
CN109389248B (en) Thermoelectric power coordinated scheduling method of comprehensive energy system based on regional heat supply network
CN110535128A (en) Based on the multizone integrated energy system coordinated dispatching method with energy comfort level
CN111969657B (en) Multi-energy station system economic optimization operation method based on negotiation game
CN109740824A (en) A kind of multi-energy system Method for optimized planning considering heating network and thermic load
CN111985105B (en) Multi-micro-energy-source network system reliability assessment method considering thermal dynamic characteristics
CN110649594A (en) Industrial park comprehensive demand response scheduling method based on multi-energy cooperation
CN113988435B (en) Comprehensive energy system source-load collaborative optimization method based on service provider guidance
Zhang et al. Bi-layer economic scheduling for integrated energy system based on source-load coordinated carbon reduction
CN113537618B (en) Comprehensive energy system optimization scheduling method considering resident user demand response
CN114529153A (en) Multi-region comprehensive energy system scheduling method considering comprehensive demand response
CN117193188A (en) Generalized energy storage optimization scheduling method for multi-park comprehensive energy system
CN111724026A (en) Optimization method for coupling operation of multi-energy network and water distribution network
CN111523697A (en) Comprehensive energy service cost allocation and pricing calculation method
CN115906401A (en) Comprehensive energy system operation optimization method considering auxiliary service carbon transaction

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