Load side demand response method and system of low-carbon urban comprehensive energy system

Abstract

The invention discloses a load side demand response method and system of a low-carbon city integrated energy system. The method comprises the steps of establishing a load side demand response model of the low-carbon city integrated energy system; air conditioner function permeability division is conducted on the load side of the low-carbon city comprehensive energy system, and a function permeability division result is determined; establishing a low-carbon city integrated energy system air conditioner participation peak regulation model; and demand response is carried out based on the load side demand response model, the function permeability division result and the air conditioner participation peak regulation model, and the load side demand response total load variation is determined. According to the method, load demand response is divided into price-sensitive demand response and energy conversion demand response according to load-side energy consumption characteristic division of the low-carbon urban comprehensive energy system; a low-carbon city comprehensive energy system load side air conditioner participates in power grid peak regulation, an air conditioner function permeability quantitative calculation method is provided, and reference can be provided for a low-carbon operation optimization strategy of the low-carbon city comprehensive energy system.

Description

Technical Field

[0001] The present invention relates to the technical field of low-carbon dispatching and operation of power systems, and more specifically, to a load-side demand response method and system for a low-carbon urban integrated energy system. Background Art

[0002] As climate issues become increasingly prominent, 178 parties around the world have signed the Paris Agreement, committed to making unified arrangements for global climate change actions. Due to the diversified development of urban energy, the electricity demand and electricity consumption forms of system users have become more diverse, resulting in more sources of carbon emissions. Therefore, carbon emission management, environmental requirements, energy supply pressure, etc. will pose huge challenges to the rational operation of urban integrated energy systems.

[0003] The urban integrated energy system is developing rapidly, and the energy supply mode has undergone tremendous changes. Traditional urban energy supply is mainly based on coal, natural gas and power grid power supply with synchronous generators as the main body. However, with the low-carbon development of the energy structure, clean energy represented by wind power and photovoltaics will participate in energy supply and gradually participate in energy supply in a distributed grid-connected manner. Thermal storage equipment, cogeneration equipment, energy storage systems, etc. will also participate in the charging and discharging process of system energy on a large scale.

[0004] Therefore, a load-side demand response method for a low-carbon urban integrated energy system that takes into account the participation of electric energy flow in grid peak regulation is needed. Summary of the invention

[0005] The present invention proposes a load-side demand response method and system for a low-carbon city integrated energy system to solve the problem of lack of low-carbon operation scheduling in modern low-carbon city integrated energy systems.

[0006] In order to solve the above problems, according to one aspect of the present invention, a load-side demand response method for a low-carbon urban integrated energy system is provided, the method comprising:

[0007] Establish a load-side demand response model for low-carbon urban integrated energy systems;

[0008] Carry out air conditioning function penetration rate division on the load side of the low-carbon city comprehensive energy system and determine the function penetration rate division results;

[0009] Establish a peak load regulation model for air conditioning in low-carbon urban integrated energy systems;

[0010] Based on the load-side demand response model, the functional penetration rate division result and the air-conditioning participation peak load regulation model, demand response is performed to determine the total load change of the load-side demand response.

[0011] Preferably, the load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model;

[0012] The load reduction model includes:

[0013]

[0014] The adjustable load model comprises:

[0015]

[0016] The energy conversion demand response sub-model includes:

[0017]

[0018] in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL (t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

[0019] Preferably, the step of dividing the air conditioning function penetration rate on the load side of the low-carbon city comprehensive energy system and determining the function penetration rate division result comprises:

[0020] Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system

[0021] Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ;

[0022] use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon citiesAC,PS :

[0023] use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

[0024] Preferably, the establishment of a low-carbon city comprehensive energy system air conditioning participation peak load regulation model includes:

[0025] Establish the air conditioning thermodynamic model, including:

[0026]

[0027] Build an indoor air temperature model, including:

[0028]

[0029] Among them, the polynomial f 1 、f 2 、f 3 and f 4 satisfy:

[0030]

[0031] The power model of air conditioner participating in power grid peak load regulation at time t is established as follows:

[0032]

[0033] Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W·V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C·W is the water vapor temperature coefficient, is the operating status of the RAC, T 1 and T 0 The operating status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, Stop and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C w are the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

[0034] Preferably, the step of performing demand response based on the load-side demand response model, the functional penetration rate division result and the air-conditioning peak load regulation model to determine the total load change of the load-side demand response includes:

[0035]

[0036] in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; BCL , B AL , BR 1 , BR 2 , B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

[0037] According to another aspect of the present invention, a load-side demand response system for a low-carbon urban integrated energy system is provided, the system comprising:

[0038] Demand response model building unit, used to build a load-side demand response model for a low-carbon city integrated energy system;

[0039] The functional penetration rate division unit is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city comprehensive energy system and determine the functional penetration rate division result;

[0040] The peak load regulation model establishment unit is used to establish a peak load regulation model for air conditioners in a low-carbon city comprehensive energy system;

[0041] The total load change determination unit is used to perform demand response based on the load side demand response model, the functional penetration rate division result and the air conditioning participation peak load regulation model, and determine the total load change of the load side demand response.

[0042] Preferably, the load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model;

[0043] The load reduction model includes:

[0044]

[0045] The adjustable load model comprises:

[0046]

[0047] The energy conversion demand response sub-model includes:

[0048]

[0049] in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ kis the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL (t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

[0050] Preferably, the functional penetration rate division unit is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city comprehensive energy system, and determine the functional penetration rate division result, including:

[0051] Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system

[0052] Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ;

[0053] use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS :

[0054] use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

[0055] Preferably, the peak load regulation model establishing unit establishes a peak load regulation model for air conditioners in a low-carbon city comprehensive energy system, including:

[0056] Establish the air conditioning thermodynamic model, including:

[0057]

[0058] Build an indoor air temperature model, including:

[0059]

[0060] Among them, the polynomial f 1 、f 2 、f 3 and f 4 satisfy:

[0061]

[0062] The power model of air conditioner participating in power grid peak load regulation at time t is established as follows:

[0063]

[0064] Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W·V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C·W is the water vapor temperature coefficient, is the operating status of the RAC, T 1 and T 0 The operating status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C ware the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC ,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

[0065] Preferably, the total load change determination unit performs demand response based on the load side demand response model, the functional penetration rate division result and the air conditioning participation peak load regulation model to determine the total load change of the load side demand response, including:

[0066]

[0067] in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR 1 , BR 2 , B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

[0068] The present invention provides a method and system for load-side demand response of a low-carbon city integrated energy system, including: establishing a load-side demand response model of a low-carbon city integrated energy system; dividing the air-conditioning function penetration rate on the load side of the low-carbon city integrated energy system, and determining the function penetration rate division result; establishing a model for air-conditioning participation in peak-shaving of a low-carbon city integrated energy system; performing demand response based on the load-side demand response model, the function penetration rate division result and the air-conditioning participation in peak-shaving model, and determining the total load change of the load-side demand response. The present invention divides the load demand response into price-sensitive demand response and energy conversion demand response according to the energy consumption characteristics of the load side of the low-carbon city integrated energy system; proposes that the air-conditioning on the load side of the low-carbon city integrated energy system participates in the peak-shaving of the power grid, and provides a quantitative calculation method for the air-conditioning function penetration rate, which can provide a reference for the low-carbon operation optimization strategy of the low-carbon city integrated energy system. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] A more complete understanding of exemplary embodiments of the present invention may be obtained by referring to the following drawings:

[0070] Figure 1 It is a flow chart of a load-side demand response method 100 of a low-carbon urban integrated energy system according to an embodiment of the present invention;

[0071] Figure 2 It is a schematic diagram of a load-side demand response process of a low-carbon urban integrated energy system taking into account the participation of electric energy flow in grid peak regulation according to an embodiment of the present invention;

[0072] Figure 3 A topological structure diagram of a low-carbon urban integrated energy system according to an embodiment of the present invention;

[0073] Figure 4 is a schematic diagram of an active output curve after the air conditioner participates in peak shaving and valley filling according to an embodiment of the present invention;

[0074] Figure 5 Schematic diagram of the structure of a load-side demand response system 500 of a low-carbon urban integrated energy system according to an embodiment of the present invention. DETAILED DESCRIPTION

[0075] Now, exemplary embodiments of the present invention are described with reference to the accompanying drawings. However, the present invention can be implemented in many different forms and is not limited to the embodiments described herein. These embodiments are provided to disclose the present invention in detail and completely and to fully convey the scope of the present invention to those skilled in the art. The terms used in the exemplary embodiments shown in the accompanying drawings are not intended to limit the present invention. In the accompanying drawings, the same units / elements are marked with the same reference numerals.

[0076] Unless otherwise specified, the terms (including technical terms) used herein have the commonly understood meanings to those skilled in the art. In addition, it is understood that the terms defined in commonly used dictionaries should be understood to have the same meanings as those in the context of the relevant fields, and should not be understood as idealized or overly formal meanings.

[0077] Figure 1 FIG. 1 is a flow chart of a method 100 for responding to a load side demand of a low-carbon urban integrated energy system according to an embodiment of the present invention. Figure 1 As shown, the load-side demand response of the low-carbon city integrated energy system provided by the present invention participates in the peak-shaving model to perform demand response, and determines the total load change in the load-side demand response. The present invention divides the load demand response into price-sensitive demand response and energy conversion demand response according to the energy consumption characteristics of the load side of the low-carbon city integrated energy system; it proposes that the load-side air conditioner of the low-carbon city integrated energy system participates in the peak-shaving of the power grid, and provides a quantitative calculation method for the air-conditioning function penetration rate, which can provide a reference for the low-carbon operation optimization strategy of the low-carbon city integrated energy system. The load-side demand response method 100 of the low-carbon city integrated energy system provided by the embodiment of the present invention starts from step 101. In step 101, a load-side demand response model of the low-carbon city integrated energy system is established.

[0078] Preferably, the load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model;

[0079] The load reduction model includes:

[0080]

[0081] The adjustable load model comprises:

[0082]

[0083] The energy conversion demand response sub-model includes:

[0084]

[0085] in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL (t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

[0086] Combination Figure 2 As shown, in the present invention, it is first necessary to input the demand response classification result and each demand response parameter and the initial state quantity of the load.

[0087] Among them, according to the load response characteristics of the power grid, the demand response of the demand side of the low-carbon city integrated energy system is divided into price-sensitive demand response and energy conversion demand response. The demand response parameters include curtailable load demand response parameters, adjustable load demand response parameters and energy conversion load demand response parameters. Among them, curtailable load demand response parameters and load initial state quantities include: initial load at time t The change in electricity price at time k after DR Δψ k , initial electricity price at time k The flexibility coefficient ε of the load at time t to the electricity price at time k t,k , the load price demand flexibility matrix can be reduced CL (t,k), user satisfaction coefficient ζ Cus Among them, the adjustable load demand response parameters and load initial state quantities include: adjustable load price demand flexibility matrix Ξ AL (t, k), initial adjustable load at time t Among them, the energy conversion load demand response parameters and load initial state quantities include: energy conversion coefficient ζ; the unit calorific value of electric energy and thermal energy is H ε , H h ; The energy utilization rates of electrical energy and thermal energy are η ε and η h .

[0088] In the present invention, a load-side demand response model for a low-carbon urban integrated energy system is established.

[0089] According to the characteristics of grid load response, the demand side demand response of low-carbon urban integrated energy system is divided into price-sensitive demand response and energy conversion demand response:

[0090] The loads participating in price-sensitive demand response are divided into curtailable loads and adjustable loads. The curtailable load model uses the price demand flexibility matrix to represent the curtailable load demand response characteristics:

[0091] The element ε in the tth row and kth column of the flexibility matrix Ξ(t,k) t,k It represents the flexibility coefficient of load at time t to electricity price at time k:

[0092]

[0093] In the formula, is the load change at time t after DR, is the initial load at time t, Δψ k is the change in electricity price at time k after DR, is the initial electricity price at time k.

[0094] The change in load that can be reduced after the load demand response is:

[0095]

[0096] In the formula, is the initial load reduction at time t, CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient.

[0097] The adjustable load model uses the price demand flexibility matrix to represent the adjustable load demand response characteristics. The change in adjustable load at time t after demand response for:

[0098]

[0099] In the formula, is the initial adjustable load at time t; AL (t,k) is the adjustable load price demand flexibility matrix.

[0100] The energy conversion load demand response model is:

[0101]

[0102] In the formula, are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

[0103] In step 102, the air conditioning function penetration rate is divided on the load side of the low-carbon city comprehensive energy system to determine the function penetration rate division result.

[0104] Preferably, the step of dividing the air conditioning function penetration rate on the load side of the low-carbon city comprehensive energy system and determining the function penetration rate division result comprises:

[0105] Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system

[0106] Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ;

[0107] use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS :

[0108] use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

[0109] In the present invention, the air conditioning function penetration rate on the load side of the low-carbon city comprehensive energy system is divided according to the air conditioning load utilization rate on the previous day. Specifically, the function penetration rate division steps are as follows:

[0110] (1) Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system

[0111] (2) Collect the latest grid-connected rated power value P of the air-conditioning load of the low-carbon urban integrated energy system AC,Rated ;

[0112] (3) Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS , and satisfy:

[0113]

[0114] (4) Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR , and satisfy:

[0115]

[0116] In step 103, a peak load regulation model of air conditioning in a low-carbon urban integrated energy system is established.

[0117] Preferably, the establishment of a low-carbon city comprehensive energy system air conditioning participation peak load regulation model includes:

[0118] Establish the air conditioning thermodynamic model, including:

[0119]

[0120] Build an indoor air temperature model, including:

[0121]

[0122] Among them, the polynomial f 1 、f 2 、f 3 and f 4 satisfy:

[0123]

[0124] The power model of air conditioner participating in power grid peak load regulation at time t is established as follows:

[0125]

[0126] Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W.V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C.W is the water vapor temperature coefficient, is the operating status of the RAC, T 1 and T 0 The operating status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,qz The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C w are the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

[0127] In the present invention, a peak load regulation model of air conditioners in a low-carbon city comprehensive energy system is established. Specifically, it includes:

[0128] Establish the thermodynamic model of air conditioning:

[0129]

[0130] Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W.V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity. 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 They are constants with values ​​of 1.07, 0.2, -0.65, -2.7, 6.105, 17.27 and 237.7 respectively.

[0131] The indoor air temperature model is established as:

[0132]

[0133] Among them, Q p is the rated cooling capacity of the coolant group, γ C.W is the water vapor temperature coefficient, is the operating status of the RAC, T 1 and T 0 The operating status and Initial room temperature.

[0134] Among them, the polynomial f 1 、f 2 、f 3 and f 4 satisfy:

[0135]

[0136] Among them, σ a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z are the heat dissipation per unit area of ​​equipment and lighting respectively. x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the aggregation coefficient. m w , C w are the mass and specific heat of chilled water respectively.

[0137] Among them, the cooling power when the air conditioner reaches steady state is:

[0138]

[0139] Among them, R a is the equivalent thermal resistance of RAC.

[0140] The power model of air conditioner participating in power grid peak load regulation at time t is established as follows:

[0141]

[0142] Among them, P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

[0143] In step 104, based on the load-side demand response model, the functional penetration rate division result and the air-conditioning participation peak load regulation model, demand response is performed to determine the total load change of the load-side demand response.

[0144] Preferably, the step of performing demand response based on the load-side demand response model, the functional penetration rate division result and the air-conditioning peak load regulation model to determine the total load change of the load-side demand response includes:

[0145]

[0146] in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR 1 , BR 2 , B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

[0147] Combination Figure 2 As shown, according to the established model and the input response parameters and initial load state quantity, calculation is performed, and the total load change of the load side demand response is output.

[0148] Among them, the total load change of the load side demand response at time t is:

[0149]

[0150] in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR 1 , BR 2 , B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

[0151] In the present invention, a typical low-carbon urban integrated energy system is established based on the method of the present invention, and the topology is as follows: Figure 3 As shown. Among them, electric energy and gas energy are supplied by the upper power grid and gas grid respectively. Gas purchased from the upper gas grid is used to supply combined heat and power (CHP) and gas boiler (GB), and the remaining electric energy can be sold to the upper power grid; the energy coupling equipment includes CHP, heat pump (HP) and GB, which can realize the bidirectional flow of electric and thermal energy; CHP is composed of gas turbine (GT), waste heat boiler (WHB) and low-temperature waste heat power generation device based on organic Rankine cycle (ORC), and the operation mode is thermal-electric decoupling, which can adapt to different operating conditions of the system; since gas boilers need power support during startup and operation, HP and gas boilers can absorb wind power and bear part of the heat load. Among them, wind power is transmitted from the large-scale new energy base of Shagohuang to the low-carbon urban integrated energy system coupled with the power of the receiving grid through direct current transmission. Smart buildings on the load side are gradually participating in the low-carbon optimization operation of low-carbon integrated energy systems, mainly represented by electric vehicles and cold storage air conditioners. Introducing demand response can smooth out load curve fluctuations, achieve interactive coupling of electricity and heat, shaving peaks and filling valleys, and reduce operating costs. The active output curve of the air conditioner after participating in peak shaving and filling valleys is as follows: Figure 4 shown.

[0152] The method of the present invention proposes a load-side demand response model for a low-carbon city integrated energy system that takes into account the participation of electric energy flow in grid peak regulation; proposes a method for calculating the air-conditioning function penetration rate applicable to a low-carbon city integrated energy system; proposes an air-conditioning load peak regulation model applicable to a low-carbon city integrated energy system, solves the technical problem of lack of low-carbon operation scheduling of modern low-carbon city integrated energy systems, and can provide a reference for low-carbon operation optimization strategies of low-carbon city integrated energy systems.

[0153] Figure 5 FIG. 5 is a schematic diagram of the structure of a load-side demand response system 500 of a low-carbon urban integrated energy system according to an embodiment of the present invention. Figure 5 As shown, the load-side demand response system 500 of the low-carbon urban integrated energy system provided in the embodiment of the present invention includes: a demand response model establishment unit 501, a functional penetration rate division unit 502, a peak-shaving model establishment unit 503 and a total load change determination unit 504.

[0154] Preferably, the demand response model establishing unit 501 is used to establish a load-side demand response model for a low-carbon urban integrated energy system.

[0155] Preferably, the load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model;

[0156] The load reduction model includes:

[0157]

[0158] The adjustable load model comprises:

[0159]

[0160] The energy conversion demand response sub-model includes:

[0161]

[0162] in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL(t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

[0163] Preferably, the functional penetration rate division unit 502 is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city comprehensive energy system and determine the functional penetration rate division result.

[0164] Preferably, the functional penetration rate division unit 502 is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city integrated energy system, and determine the functional penetration rate division result, including:

[0165] Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system

[0166] Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ;

[0167] use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS :

[0168] use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

[0169] Preferably, the peak load regulation model establishing unit 503 is used to establish a peak load regulation model in which air conditioning participates in the low-carbon city integrated energy system.

[0170] Preferably, the peak load model establishing unit 503 establishes a peak load model for air conditioners in a low-carbon city integrated energy system, including:

[0171] Establish the air conditioning thermodynamic model, including:

[0172]

[0173] Build an indoor air temperature model, including:

[0174]

[0175] Among them, the polynomial f 1 、f 2 、f 3 and f4 satisfy:

[0176]

[0177] The power model of air conditioner participating in power grid peak load regulation at time t is established as follows:

[0178]

[0179] Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W.V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C.W is the water vapor temperature coefficient, is the operating status of the RAC, T 1 and T 0 The operating status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C w are the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

[0180] Preferably, the total load change determination unit 504 is used to perform demand response based on the load-side demand response model, the functional penetration rate division result and the air-conditioning participation peak load regulation model, and determine the total load change of the load-side demand response.

[0181] Preferably, the total load change determination unit 504 performs demand response based on the load side demand response model, the function penetration rate division result and the air conditioning participation peak load regulation model to determine the total load change of the load side demand response, including:

[0182]

[0183] in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR 1 , BR 2 , B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

[0184] The load-side demand response system 500 of the low-carbon urban integrated energy system of the embodiment of the present invention corresponds to the load-side demand response method 100 of the low-carbon urban integrated energy system of another embodiment of the present invention, which will not be repeated here.

[0185] The present invention has been described with reference to a few embodiments. However, it is known to those skilled in the art that other embodiments than the one disclosed above are equally within the scope of the present invention.

[0186] Generally, all terms used in the present invention are interpreted according to their ordinary meaning in the technical field, unless otherwise explicitly defined therein. All references to "a / said / the [device, component, etc.]" are open to interpretation as at least one instance of the device, component, etc., unless otherwise explicitly stated. The steps of any method disclosed herein do not necessarily have to be performed in the exact order disclosed, unless explicitly stated.

[0187] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0188] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0189] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0190] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0191] Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents, and any modifications or equivalent replacements that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the present invention.

Claims

1. A load-side demand response method for a low-carbon urban integrated energy system, characterized in that: The method comprises: Establish a load-side demand response model for low-carbon urban integrated energy systems; Carry out air conditioning function penetration rate division on the load side of the low-carbon city comprehensive energy system and determine the function penetration rate division results; Establish a peak load regulation model for air conditioning in low-carbon urban integrated energy systems; Based on the load-side demand response model, the functional penetration rate division result and the air-conditioning participation peak load regulation model, demand response is performed to determine the total load change of the load-side demand response.

2. The method according to claim 1, characterized in that The load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model; The load reduction model includes: The adjustable load model comprises: The energy conversion demand response sub-model includes: in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL (t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

3. The method according to claim 1, characterized in that The above-mentioned classification of air conditioning function penetration rate on the load side of the low-carbon city comprehensive energy system and determination of the function penetration rate classification results include: Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ; use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS : use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

4. The method according to claim 1, characterized in that The establishment of a low-carbon city comprehensive energy system air conditioning participation peak load regulation model includes: Establish the air conditioning thermodynamic model, including: Build an indoor air temperature model, including: Among them, polynomials f1, f2, f3 and f4 satisfy: The power model of air conditioner participating in power grid peak load regulation at time t is established as follows: Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W.V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C·W is the water vapor temperature coefficient, is the running status of RAC, T1 and T0 are the running status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C w are the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and I RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

5. The method according to claim 1, characterized in that The step of performing demand response based on the load-side demand response model, the functional penetration rate division result and the air-conditioning participation peak load regulation model, and determining the total load change amount of the load-side demand response includes: in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR1, BR2, B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

6. A load-side demand response system for a low-carbon urban integrated energy system, characterized in that: The system comprises: Demand response model building unit, used to build a load-side demand response model for a low-carbon city integrated energy system; The functional penetration rate division unit is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city comprehensive energy system and determine the functional penetration rate division result; The peak load regulation model establishment unit is used to establish a peak load regulation model for air conditioners in a low-carbon city comprehensive energy system; The total load change determination unit is used to perform demand response based on the load side demand response model, the functional penetration rate division result and the air conditioning participation peak load regulation model, and determine the total load change of the load side demand response.

7. The system according to claim 6, characterized in that The load-side demand response model of the low-carbon city integrated energy system includes: a price-sensitive demand response sub-model and an energy conversion demand response sub-model, and the price-sensitive demand response sub-model includes: a curtailable load model and an adjustable load model; The load reduction model includes: The adjustable load model comprises: The energy conversion demand response sub-model includes: in, is the initial load reduction at time t, is the initial load reduction at time t; CL (t,k) is the flexible matrix of load price demand that can be reduced, ψ k is the electricity price at time k, ζ Cus is the user satisfaction coefficient, is the initial electricity price at time k; is the initial adjustable load at time t, AL (t,k) is the adjustable load price demand flexibility matrix, is the initial adjustable load at time t; are the replaceable electrical load and the corresponding thermal load respectively; ζ is the energy conversion coefficient; H ε , H h are the unit calorific values ​​of electrical energy and thermal energy respectively; η ε , η h are the energy utilization rates of electrical energy and thermal energy, respectively.

8. The system according to claim 6, characterized in that The functional penetration rate division unit is used to divide the air conditioning functional penetration rate on the load side of the low-carbon city comprehensive energy system and determine the functional penetration rate division result, including: Collect the maximum power of building cooling / heating in low-carbon urban integrated energy system Collect the latest grid-connected rated power value P of air conditioning load in low-carbon urban integrated energy system AC,Rated ; use Calculate the peak load utilization rate η of the comprehensive energy system in low-carbon cities AC,PS : use Calculate the demand response utilization rate η of the low-carbon city comprehensive energy system AC,DR .

9. The system according to claim 6, characterized in that The peak load model establishment unit establishes a peak load model in which air conditioners in a low-carbon urban integrated energy system participate in the peak load, including: Establish the air conditioning thermodynamic model, including: Build an indoor air temperature model, including: Among them, polynomials f1, f2, f3 and f4 satisfy: The power model of air conditioner participating in power grid peak load regulation at time t is established as follows: Among them, Θ H is the indoor body temperature, T in (t) is the indoor air temperature at time t, P W·V (t) is the water vapor pressure, S(t) is the wind speed at time t, σ air is the relative air humidity; Δ 11 , Δ 12 , Δ 13 , Δ 14 , Δ 21 , Δ 22 and Δ 23 are all constants; Q p is the rated cooling capacity of the coolant group, γ C·W is the water vapor temperature coefficient, is the running status of RAC, T1 and T0 are the running status and The initial room temperature at a is the relative humidity of air, V B is the building volume, κ S is the heat storage coefficient of the building’s interior wall, S w_in is the interior wall area of ​​the building, S top and S w are the thermal conductivity of the building roof and wall respectively. wc , T wj are the outlet and inlet temperatures of cooling water, T out (t) is the outdoor temperature at time t, κ e , κ z , κ x are equipment cooling load, lighting cooling load and human visual cooling load coefficient, q e ,q z The heat dissipation per unit area of ​​equipment and lighting is q x ,q n are the sensible heat and latent heat dissipation of the human body, Q L (t) is the indoor cooling load at time t, Q r (t) is the cooling load of personnel at time t, S C is the cooling area, and are the passenger flow number and passenger flow number boundary, ε Clu is the clustering coefficient; m w , C w are the mass and specific heat capacity of chilled water respectively; Q RAC R is the cooling power when the air conditioner reaches steady state; a is the equivalent thermal resistance of RAC; P RAC (t) is the peak shaving potential at time t, N RAC is the number of refrigeration and air conditioning systems in IES, Q RAC,i is the cooling power of the i-th refrigeration and air conditioning system, Δ RAC1 is the main factor of the electric power of the refrigeration and air conditioning system, Δ RAC2 and l RAC They are the main coefficient and constant term of the cooling power of the refrigeration and air-conditioning systems respectively.

10. The system according to claim 6, characterized in that The total load change determination unit performs demand response based on the load side demand response model, the function penetration rate division result and the air conditioning participation peak load regulation model to determine the total load change of the load side demand response, including: in, is the total load change of the load side demand response at time t; is the load change that can be reduced after the load demand response at time t; is the change in load demand response that can be adjusted at time t; are the replaceable electrical load and the corresponding thermal load respectively; P RAC (t) is the peak shaving potential of air conditioning at time t; B CL , B AL , BR1, BR2, B RAC They are respectively the reducible load demand response logic control factor, the adjustable load demand response logic control factor, the replaceable electric load logic control factor, the replaceable thermal load logic control factor and the RAC peak load regulation logic control factor.

Images