A comprehensive demand response method for thermal and electrical coupled parks

Abstract

The present invention discloses a comprehensive demand response method for a thermal-electric coupled park. First, according to different energy consumption types, industrial users participating in the demand response are divided into traditional single type, dual-energy coupled type and multi-energy flexible type. By analyzing the response characteristics of each type of user, a demand response model is established; secondly, various types of equipment owned by the industrial park operator are modeled; finally, according to the interaction process between the industrial park operator and the superior power grid, a comprehensive demand response model of the thermal-electric coupled park operator that responds to the peak-shaving demand of the superior power grid is established. The objective functions are the response rate of the park operator to the peak-shaving demand of the superior power grid, the renewable energy consumption rate and the maximum total profit of the park operator. Considering multiple constraints, by coordinating the output of each equipment in the park and the comprehensive demand response of various industrial users, the optimal scheduling plan and the response strategy of each industrial user that can fully respond to the peak-shaving demand of the superior power grid are obtained, thereby improving the renewable energy consumption level in the park.

Description

Technical Field

[0001] The present invention relates to the field of new energy technology, and in particular to a comprehensive demand response method for a thermal-electrical coupling park. Background Art

[0002] Against the backdrop of the development of the Energy Internet, energy coupling technology has further improved, and a large number of energy coupling devices have emerged, making it possible to achieve multi-energy coupling coordination within a given area. Industrial parks, with their large load volumes and diverse energy consumption, are a key scenario for multi-energy integration systems. As the physical carrier of the Energy Internet, integrated energy systems are a new type of multi-energy integration system that achieves complementary power generation from multiple sources, including electricity, gas, and heat, and optimizes the source-grid-load-storage link. They are also the primary form of future energy system development. The evolution of industrial parks, traditionally powered by electricity, towards integrated energy systems with complementary power sources, can improve energy efficiency, achieve green, low-carbon, and sustainable industrial development while ensuring economic returns.

[0003] The development of industrial parks into integrated energy systems is driving the diversification of user energy demands, and the gradual evolution of traditional electricity demand response (TDR) towards integrated demand response (IDR). At the same time, users possess a more diverse range of flexible loads, allowing them to actively respond to signals from higher-level authorities by adjusting their own resources, participating in regulating the supply and demand balance of various energy sources within the system and promoting two-way interaction between supply and demand.

[0004] Therefore, further developing the potential of comprehensive demand response in industrial parks, categorizing and analyzing their response characteristics, and establishing refined models can not only fully tap the user response potential, but also improve the comprehensive utilization rate of equipment, thereby delaying or reducing equipment investment and construction costs, which is of great engineering significance.

[0005] Currently, research on demand response in industrial parks has yielded some results, but most studies focus on the characteristics of electricity demand response under different electricity pricing mechanisms. For integrated demand response with vertical energy substitution characteristics, some scholars have established simple models, primarily based on energy price signals or incentive mechanisms. Others have considered the response speeds of different types of integrated demand response and divided them into different response time scales. However, none of these studies consider the coupling and conversion of multiple energy sources such as electricity, heat, and gas during industrial production, nor the different energy consumption characteristics of various industrial users within the park. They fail to fully tap into the users' inherent regulation potential and lack consideration for their energy quality requirements. Consequently, existing integrated demand response models are not well adapted to the energy supply structure of industrial parks and lack operability in real-world scenarios. Urgent questions remain: how to fine-tune the integrated demand response resources of various industrial users within industrial parks to further tap into their response potential, and how to formulate integrated demand response strategies for industrial parks. Summary of the Invention

[0006] The purpose of the present invention is to provide a comprehensive demand response method for a thermal-electric coupling park in response to the problems in the prior art that there is little research on comprehensive demand response for industrial parks, industrial parks lack unified energy optimization and scheduling management, and their demand response potential is difficult to be further explored and utilized. The present invention starts from the perspective of thermal-electric coupling park operators, considers the coupling and conversion of multiple energy sources such as electricity, heat, and gas in the industrial production process, as well as the different energy consumption characteristics of various industrial users in the park, classifies industrial users in the park, and establishes a refined model of comprehensive demand response resources for various types of industrial users. On this basis, with the response rate of the park operator to the peak-shaving demand of the upper power grid, the renewable energy consumption rate and the total revenue of the park operator as optimization targets, a comprehensive demand response model for the thermal-electric coupling park is constructed to maximize the response rate of the park operator to the peak-shaving demand of the upper power grid and the renewable energy consumption level in the park.

[0007] In order to achieve the above object, the technical solution adopted by the present invention is:

[0008] A comprehensive demand response method for a thermal-electrical coupled park comprises the following steps:

[0009] Step 1: Based on different energy consumption types, industrial users are divided into traditional single type, dual-energy coupling type, and multi-energy flexible type, and demand response models are established for each type of industrial user;

[0010] Establish energy conversion models for various energy production equipment within the thermal-electrical coupled park;

[0011] Step 2: Establish a comprehensive demand response model for the thermal-electric coupled park that responds to the peak-shaving demand of the upper power grid, with the park operator's response rate ζ to the upper power grid's peak-shaving demand, the renewable energy consumption rate ψ, and the park operator's total revenue as the objective function;

[0012] maxΓ=φ1ζ+φ2ψ+φ3(R P -C)

[0013] Γ represents the target value, R P represents the revenue from energy sales, C represents the total cost, φ1, φ2, and φ3 represent three priority factors, and φ1>φ2>φ3;

[0014] Energy balance constraints, equipment output constraints, tie line constraints, schedulable production task constraints, and IDR response constraints for various industrial users are used as constraints.

[0015] Step 3: Solve the comprehensive demand response model of the thermal-electrical coupled park to obtain a comprehensive demand response plan for the thermal-electrical coupled park.

[0016] Preferably, in step 1, the demand response model of a traditional single-type industrial user is:

[0017]

[0018]

[0019]

[0020]

[0021]

[0022] Where n represents the schedulable production task index, and there are N production tasks in total; t n It represents the starting period of the schedulable production task n in the actual production process; and is a binary variable, In the actual production process, the schedulable production task n starts at time period t. Indicates that the schedulable production task n in the original production plan is in operation during time period t. In the actual production process, the schedulable production task n is in operation during time period t; T represents the scheduling period, T n represents the start period of the schedulable production task n in the original production plan; T n，s represents the continuous running time of the schedulable production task n; T n，max and T n，min They represent the upper and lower limits of the start-up period of the schedulable production task n; Le，t1 represents the electric load demand in time period t after only the production task can be dispatched, L e，n represents the power demand of the schedulable production task n; L e，t0 Represents the electricity load demand in time period t in the original production plan.

[0023] Preferably, in step 1, the demand response model of the response model of the electric-thermal coupling type industrial user is:

[0024] The electric-thermal coupling load model is:

[0025]

[0026] Where, Indicates the change in electrical load due to the change in thermal load during period t; represents the change in heat load during period t; η eh represents the conversion factor of thermoelectric units; ξ represents the thermoelectric coupling coefficient;

[0027] Heat load demand during non-demand response period:

[0028]

[0029]

[0030] Where, represents the medium-grade heat consumed by the steam-driven equipment during period t; Indicates the low-grade heat recovery efficiency of steam-driven equipment; Q h，t represents the total heat load demand of electric heating industrial users during period t; Respectively represent the medium-grade heat load demand and low-grade heat load demand for direct utilization; η d ,η low They represent the conversion efficiency of heat supply to direct steam driven loads and low-grade thermal power equipment respectively;

[0031] Heat load demand during demand response period:

[0032]

[0033] Where Q hDR，t represents the total heat load demand of electric heating type industrial users in the post-response period t; They represent the response changes of the medium-grade thermal power directly utilized and the changes of the recovered and converted low-grade thermal power respectively;

[0034] The electric load demand during the demand response period is:

[0035]

[0036] Where, L e，t2 represents the electric load demand in time period t after only the electric-thermal coupling type user responds, L e，t0 Represents the electricity load demand in time period t in the original production plan.

[0037] Preferably, in step 1, the demand response model of the response model of the electrically coupled industrial user is:

[0038] The electrical coupling load model is:

[0039] ΔL etr，t =-η ge ΔL gtr，t

[0040] Where ΔL etr，t , Δ Lgtr，t The increment of electricity and gas load caused by electricity / gas substitution during actual operation; η ge Indicates the conversion coefficient of equal calorific value of electricity and natural gas;

[0041] The electrical load substitution rate curve expression is as follows:

[0042]

[0043] Where λ ge represents the load substitution rate, Δp ge represents the electricity price difference, a ge 、b ge and λ ge，max They represent the electrical substitution dead zone threshold, saturation zone threshold and maximum electrical load substitution rate respectively;

[0044] The fitted electric load and gas load in period t after taking into account the response of alternative loads are expressed as:

[0045] L e，t3 =L e，t0 -λ ge L e，t0 +ΔL etr，t

[0046] L g，t1 =L g，t0 +λ ge L e，t0 / η ge +ΔL gtr，t

[0047] Where, L e，t3 represents the electric load demand in time period t after only the electrically coupled users respond, L e，t0 represents the electricity load demand in period t in the original production plan; L g，t1 represents the gas load demand in period t after only the electrical coupling type users respond, L g，t0Represents the gas load demand in time period t in the original production plan.

[0048] Preferably, in step 1, the demand response model of the multi-energy flexible industrial user's response model is:

[0049] ΔL e，t =ΔL g，t L NG η e / Δt

[0050] ΔQ h，t =ΔL e，t / ξη eh

[0051] L e，t4 =L e，t0 -ΔL e，t

[0052] Q h，t =Q h，t0 -ΔQ h，t

[0053] L g，t2 =L g，t0 +ΔL g，t

[0054] Where ΔL e，t Indicates the change in electricity purchased from the park; ΔL g，t Indicates the change in natural gas required by industrial users to meet the change in electricity; L NG Indicates the lower calorific value of natural gas; η e represents the efficiency of natural gas power generation; Δt represents the time step; ΔQ h，t represents the change in heat load; ξ represents the thermoelectric coupling coefficient; η eh Indicates the conversion factor of thermoelectric units; L e，t4 represents the electric load demand in time period t after only the multi-energy flexible users respond; L e，t0 represents the electricity load demand in period t in the original production plan; Q h，t represents the heat load demand in time period t after the response; Q h，t0 represents the heat load demand in period t of the original production plan; L g，t2 represents the gas load demand in time period t after only the multi-energy flexible users respond; L g，t0 Represents the gas load demand in time period t in the original production plan.

[0055] Preferably, in step 2, the response rate ζ of the park operator to the peak-shaving demand of the upper-level power grid is:

[0056]

[0057] Where: ΔPref，t P represents the peak load demand issued by the upper power grid during period t, that is, the power purchased by the park operator needs to be reduced; e，t 、P eDR，t The electricity purchase volume before and after the IDR of the park operator respectively;

[0058] The renewable energy consumption rate ψ is:

[0059]

[0060] Where: P WT，t and P PV，t Respectively represent the wind turbine and photovoltaic output power during period t; P Ab，t Indicates the amount of wind and solar power abandoned during period t.

[0061] Preferably, in step 2, the energy balance constraints include electric power balance, thermal power balance and gas power balance, which can be expressed as:

[0062] P MTe，t +P PV，t +P wT，t +P be，t +(P ES，dis，t -P ES，ch，t )-P P2Ge，t -P EB，t =L e，t

[0063] Q MTh，t +Q GB，t +Q EB，t +(Q HS，dis，t -Q HS，ch，t )=L h，t

[0064] P bg，t +P P2Gg，t -P MTg，t -P GBg，t =L g，t

[0065] Where, L e，t , L h，t and L g，t are the electricity, heat and gas load demands during period t respectively; P MTe，t represents the output power of the gas turbine during period t; P PV，t represents the photovoltaic output power during period t; P WT，t represents the wind power output power during period t, P be，t P represents the amount of electricity purchased from the upper power grid during period t; ES，ch，t and P ES，dis，t They represent the charging and discharging power of the energy storage during period t; P P2Ge，t and P EB，tThey represent the power consumption of the power-to-gas equipment and the electric boiler during period t; Q MTh，t , Q GB，t , Q EB，t They represent the heat output of gas turbine, gas boiler and electric boiler during period t respectively; g HS，ch，t and Q HS，dis，t P represents the thermal energy storage charging and discharging power during period t; bg，t represents the gas purchase volume during period t, P P2Gg，t P represents the output of the power-to-gas equipment during period t; MTg，t and P GBg，t are the gas consumption of gas turbine and gas boiler during period t respectively.

[0066] Preferably, in step 2, the energy storage device constraint condition is:

[0067] E m，t =(1-μ m )E m，t-1 +(P m，cb，t η m，ch -P m，dis，t / η m，dis )Δt

[0068]

[0069]

[0070] E m，mi n≤E m，t ≤E m，max

[0071] 0≤I m，ch，t +I m，dis，t ≤1

[0072] E m，0 =E m，T

[0073] Where, E m，t-1 and E m，t is the energy storage capacity of the mth type energy storage device during periods t-1 and t; η m，ch and η m，dis represents the charging and discharging efficiency of the mth type of energy storage device; μ m represents the energy loss rate of the mth type energy storage device; P m，ch，t and P m，dis，t Respectively represent the charging and discharging power of the energy storage device; I m，ch，t and I m，dis，t It is a binary variable, and its value is 1, indicating that the m-th type of energy storage device is charging or discharging during period t. m，ch，t and I m，dis，t Cannot be 1 at the same time; and is the maximum charge and discharge energy of the mth type energy storage device; E m，max and E m，min Respectively represent the upper and lower limits of the storage capacity of the mth type of energy storage device; E m，0 、E m，T They represent the stored energy of the mth type of energy storage device at the initial moment and after one cycle respectively.

[0074] Preferably, in step 2, the schedulable production task constraints are:

[0075]

[0076]

[0077]

[0078]

[0079]

[0080]

[0081] Where, T back Indicates the period during which the schedulable production task n2 lags behind the schedulable production task n1.

[0082] Preferably, in step 3, the standard form of the solution model is:

[0083]

[0084] In the formula, stg(x)=0 is the constraint condition, x i is a continuous variable, x i The minimum and maximum values of x are min and x max , x j It is an integer 0-1 variable.

[0085] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0086] 1. The present invention takes into account the differences in energy consumption types and demand response forms of industrial users, and divides industrial users into three categories: traditional single type, dual-energy coupling type, and multi-energy flexible type. It also establishes refined demand response models for their comprehensive demand response resources, which can further tap the response potential of industrial users and improve the energy flexibility of each industrial user.

[0087] 2. The present invention constructs a comprehensive demand response strategy generation model for thermal-electric coupled park operators with the goal of maximizing the park operator's response rate to the peak-shaving demand of the superior power grid, the renewable energy consumption rate and the total profit of the park operator, and proposes a refined response strategy for thermal-electric coupled park operators, which enables the park operator to fully respond to the peak-shaving demand of the superior power grid while improving the renewable energy consumption level within the park. BRIEF DESCRIPTION OF THE DRAWINGS

[0088] Figure 1 It is the relationship curve between the electrical load substitution rate and the electrical price difference.

[0089] Figure 2 This is the structural diagram of the thermal and electrical coupling park.

[0090] Figure 3 It is a flow chart of the interaction between the superior power grid, park operators and users.

[0091] Figure 4 It is the electricity load demand curve of various industrial users in the park on a typical day.

[0092] Figure 5 It is a heat load demand curve of various industrial users in the park on a typical day.

[0093] Figure 6 It is a gas load demand curve of various industrial users in the park on a typical day.

[0094] Figure 7 This is a typical solar photovoltaic and wind turbine output curve in the park.

[0095] Figure 8 This is the result diagram of the park's original planned electric power scheduling.

[0096] Figure 9 This is the result diagram of the park's original planned gas power scheduling.

[0097] Figure 10 This is the result diagram of the park's original planned thermal power scheduling.

[0098] Figure 11 This is a diagram of the operation status of the park's originally planned schedulable production tasks.

[0099] Figure 12 This is the electric power dispatch result diagram after the park responds.

[0100] Figure 13 This is the gas power dispatch result diagram after the park responds.

[0101] Figure 14 This is the thermal power dispatch result diagram after the park responds.

[0102] Figure 15It is a diagram of the operation status of schedulable production tasks after the park responds.

[0103] Figure 16 This is a diagram of the electricity load response of various industrial users in the park.

[0104] Figure 17 It is the heat load response of various industrial users in the park.

[0105] Figure 18 It is the gas load response of various industrial users in the park. DETAILED DESCRIPTION

[0106] The present invention will be described in detail below with reference to the accompanying drawings.

[0107] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0108] Example 1

[0109] A comprehensive demand response method for thermal-electrically coupled industrial parks first categorizes participating industrial users into traditional single-energy, dual-energy coupled, and multi-energy flexible types based on their energy usage. Dual-energy coupled includes both thermal-electric coupled and electrical-electric coupled types. By analyzing the response characteristics of each type of user, a corresponding demand response model is established.

[0110] Secondly, various types of equipment owned by industrial park operators are modeled, including distributed energy such as wind turbines and photovoltaic cells, energy conversion devices such as gas turbines and electric boilers, and energy storage devices.

[0111] Finally, based on the interaction process between industrial park operators and the upper-level power grid, a comprehensive demand response model for thermal and electrical coupled industrial park operators was established to respond to the upper-level power grid's peak-shaving needs. Taking the park operator's response rate to the upper-level power grid's peak-shaving needs, the renewable energy consumption rate, and the maximization of the park operator's total revenue as the objective functions, and considering constraints such as energy balance, equipment output, upper and lower limits on tie lines, and dispatchable production tasks, by coordinating the output of various equipment within the park and the comprehensive demand responses of various industrial users, an optimal dispatch plan and response strategies for each industrial user were derived that fully responded to the upper-level power grid's peak-shaving needs, while also improving the park's renewable energy consumption level.

[0112] The specific steps include:

[0113] 1 Industrial User Classification and Response Model

[0114] Industrial park users may have different energy coupling devices based on production needs. These users have a certain ability to adjust their own energy needs and can serve as important demand response resources for participating in higher-level peak load regulation. Based on their energy usage, industrial users are categorized as traditional single-energy users, dual-energy coupled users, and multi-energy flexible users.

[0115] 1.1 Traditional single response model

[0116] Traditional single-use industrial users have high electricity consumption but low heat and gas load demands. They can only respond to electricity as a single energy source and typically adjust and shift production tasks according to production plans, production processes, and the specific sequence between different industrial production tasks to participate in demand-side management to maximize profits. Considering the coupling between industrial production tasks and the startup and operation status of production tasks, the following model can be established to describe schedulable production tasks:

[0117]

[0118]

[0119]

[0120]

[0121]

[0122] Where n represents the schedulable production task index, and there are N production tasks in total; t n It represents the starting period of the schedulable production task n in the actual production process; and is a binary variable, In the actual production process, the schedulable production task n starts at time period t. Indicates that the schedulable production task n in the original production plan is in operation during time period t. In the actual production process, the schedulable production task n is in operation during time period t; T represents the scheduling period, T n represents the start period of the schedulable production task n in the original production plan; T n，s represents the continuous running time of the schedulable production task n; I n，max and I n，min They represent the upper and lower limits of the start-up period of the schedulable production task n; L e，t1 represents the electric load demand in time period t after only the production task can be dispatched, L e，n represents the power demand of the schedulable production task n; L e，t0 Represents the electricity load demand in time period t in the original production plan.

[0123] 1.2 Dual-energy coupling response model

[0124] Dual-energy coupled industrial users are primarily divided into two types: electric-thermal coupled loads and electrical-coupled loads. These loads can achieve multi-energy complementarity based on their own energy-coupled production equipment, such as electric boilers and power-to-gas (P2G) equipment. Adjustments in demand for one energy source can affect the supply balance of the other. Therefore, by balancing the demands of both energy sources, energy shortages can be alleviated and peak-shaving and valley-filling can be achieved.

[0125] (1) Electrothermal coupling type

[0126] Industrial users with coupled electricity and heat have production equipment that exhibits thermal-electric coupling. Therefore, these users can adjust their heat load to alter their electricity demand. When the upper-level power grid demands peak load reduction, two approaches can be used to increase gas turbine power generation and reduce the park's power purchases from the grid: 1) reducing user electricity purchase demand; 2) increasing user heat demand, leading to increased power generation from the park's gas turbines. The first approach involves industrial users leveraging the plant's multi-energy flow coupling mechanism to adjust the energy usage strategies of each device based on peak load reduction instructions issued by the park operator, thereby reducing power purchases during peak load reduction periods and aligning them closely with the peak load reduction target. The second approach involves the park operator using incentives and subsidies to encourage industrial users to increase their heat load demand, thereby increasing gas turbine power generation while also increasing heat output. These two peak load reduction strategies can be implemented simultaneously to maximize response potential.

[0127] The medium- and low-temperature exhaust heat generated by industrial park operators' gas turbine power generation is recycled through waste heat boilers and sold to industrial users within the park as medium-grade hot steam (170°C-550°C), serving as the primary heat source for their industrial production processes. There are two main ways for general industrial users to utilize this heat: 1) directly using the medium-grade heat purchased from the park in their industrial production processes; 2) using steam-driven equipment to recycle this medium-grade hot steam into low-grade hot steam (ambient temperature -170°C), which is then used as a heat source for factory refrigeration equipment or to generate hot water.

[0128] Therefore, the following model can be established for the electrothermal coupling load:

[0129]

[0130] Where, Indicates the change in electrical load due to the change in thermal load during period t; represents the change in heat load during period t; η ehIt represents the conversion coefficient of thermoelectric units and can be taken as 0.7; ξ represents the thermoelectric coupling coefficient, which characterizes the coupling strength of the thermoelectric load. A positive value indicates positive thermoelectric coupling, that is, the electric and thermal loads change in the same direction; a negative value indicates negative thermoelectric coupling, that is, the electric and thermal loads change in opposite directions.

[0131] Modeling is performed for the electric and heating load demands of this type of industrial users during non-demand response periods and demand response periods.

[0132] 1) Heat load demand during non-demand response period:

[0133]

[0134]

[0135] Where, represents the medium-grade heat consumed by the steam-driven equipment during period t; Indicates the low-grade heat recovery efficiency of steam-driven equipment; Q h，t represents the total heat load demand of electric heating industrial users during period t; Respectively represent the medium-grade heat load demand and low-grade heat load demand for direct utilization; η d ,η low They represent the conversion efficiency of heat supply to direct steam driven loads and low-grade thermal power equipment respectively;

[0136] 2) Heat load demand during demand response period:

[0137]

[0138] Where Q hDR，t represents the total heat load demand of electric heating type industrial users in the post-response period t; They represent the response changes of the medium-grade thermal power directly utilized and the changes of the recovered and converted low-grade thermal power respectively;

[0139] Therefore, according to formula (6), the electric load demand during the demand response period can be obtained as:

[0140]

[0141] Where, L e，t2 represents the electric load demand in time period t after only the electric-thermal coupling type user responds, L e，t0 Represents the electricity load demand in time period t in the original production plan.

[0142] (2) Electrical coupling type

[0143] Electrically coupled industrial users own electrically coupled equipment, such as air-conditioning equipment based on electric / gas hybrid refrigeration, P2G, etc. Park operators can sign contracts with such industrial users to encourage their electrically coupled replaceable loads to use natural gas to meet energy demand during peak electricity consumption periods, and to use electricity instead of gas during off-peak periods. Therefore, industrial users can obtain certain compensation benefits without affecting their energy comfort.

[0144] ΔL etr，t =-η ge ΔL gtr，t (11)

[0145] Where ΔL etr，t , ΔL gtr，t The increment of electricity and gas load caused by electricity / gas substitution during actual operation; η ge Indicates the conversion coefficient of equal calorific value between electricity and natural gas.

[0146] At the same time, industrial users can independently coordinate their own energy-consuming equipment based on the electricity price difference to reduce their energy costs. Based on the principles of consumer psychology, the load substitution rate is defined to characterize the response degree of industrial users. Its response model can be approximately fitted into a piecewise linear function, and finally the electrical load of each period can be fitted according to the load substitution rate of different periods. The uncertainty of the replaceable load response based on the principles of consumer psychology is as follows: Figure 1 As shown, the horizontal axis represents the electricity price difference Δp ge , the vertical axis represents the load substitution rate λ ge The load replacement rate is defined as the ratio of the amount of electricity replaced by the user's natural gas during that period to the amount of electricity consumed during that period.

[0147] Depend on Figure 1 The electrical load substitution rate curve expression is as follows:

[0148]

[0149] Where λ ge represents the load substitution rate, Δp ge represents the electricity price difference, a ge 、b ge and λ ge，max They represent the electrical substitution dead zone threshold, saturation zone threshold and maximum electrical load substitution rate respectively;

[0150] The fitted electric load and gas load in period t after taking into account the response of alternative loads are expressed as:

[0151] L e，t3 =L e，t0 -λ ge L e，t0 +ΔL etr，t (13)

[0152] L g，t1 =L g，t0 +λ ge L e，t0 / η ge +ΔL gtr，t (14)

[0153] Where, L e，t3 represents the electric load demand in time period t after only the electrically coupled users respond, L e，t0 represents the electricity load demand in period t in the original production plan; L g，t1 represents the gas load demand in period t after only the electrical coupling type users respond, L g，t0 Represents the gas load demand in time period t in the original production plan.

[0154] 1.3 Multifunctional and flexible response model

[0155] Multi-energy flexible industrial users have multi-energy coupling conversion equipment such as combined heat and power (CHP) units and gas turbines. Their energy consumption types are flexible and they can respond to electricity, heat and gas energy.

[0156] ΔL e，t =ΔL g，t L NG η e / Δt (15)

[0157] ΔQ h，t =ΔL e，t / ξη eh (16)

[0158] L e，t4 =L e，t0 -ΔL e，t (17)

[0159] Q h，t =Q h，t0 -ΔQ h，t (18)

[0160] L g，t2 =L g，t0 +ΔL g，t (19)

[0161] Where ΔL e，t Indicates the change in electricity purchased from the park; ΔL g，t Indicates the change in natural gas required by industrial users to meet the change in electricity; L NG Indicates the lower calorific value of natural gas, generally 9.7kWh / m 3 ;η e represents the efficiency of natural gas power generation; Δt represents the time step; ΔQh，t represents the change in heat load; ξ represents the thermoelectric coupling coefficient; η eh Indicates the conversion factor of thermoelectric units; L e，t4 represents the electric load demand in time period t after only the multi-energy flexible users respond; L e，t0 represents the electricity load demand in period t in the original production plan; Q h，t represents the heat load demand in time period t after the response; Q h，t0 represents the heat load demand in period t of the original production plan; L g，t2 represents the gas load demand in time period t after only the multi-energy flexible users respond; L g，t0 Represents the gas load demand in time period t in the original production plan.

[0162] 2 Thermal and electrical coupling park equipment model

[0163] This technology uses a typical thermal-electric coupled industrial park as a scenario to model the energy production equipment, energy conversion equipment, and energy storage equipment owned by the industrial park operator. The source side includes renewable energy equipment such as wind turbines and photovoltaic cells, and the energy conversion equipment includes gas turbines, gas boilers, electric boilers, P2G, etc., and is equipped with power and heat storage devices. The structural diagram of the thermal-electric coupled park is shown below. Figure 2 shown.

[0164] 2.1 Wind turbine

[0165] The output power P of the wind turbine WT，t The approximate relationship between it and wind speed v can be expressed by the following function:

[0166]

[0167] Where, P r is the rated power of the wind turbine; v represents the actual wind speed; v ci 、v r 、v co They are the cut-in wind speed, rated wind speed and cut-out wind speed of the wind turbine respectively.

[0168] 2.2 Photovoltaic cells

[0169] The actual output power of a photovoltaic cell, which is related to the light intensity and ambient temperature, can be obtained based on test parameters such as its maximum output power under standard test conditions:

[0170]

[0171] Where, P PV，t Indicates the actual output power of the photovoltaic cell during the t period during operation; f PVIt represents the output power derating factor, which is equal to the ratio of the actual output power of the photovoltaic cell during operation to the output power of the photovoltaic cell under rated conditions. It can generally be taken as f PV =0.95; P STC Indicates the maximum output power of photovoltaic cells under standard test conditions; G t Indicates the actual light intensity during period t; G STC Indicates the light intensity received by photovoltaic cells under standard test conditions, which can be G STC =1kW / m 2 ; k represents the temperature coefficient; T t Indicates the surface temperature of the photovoltaic cell during the t period during the actual operation; T STC Indicates the temperature of the photovoltaic cell under standard test environment, T STC =25℃.

[0172] 2.3 Gas Turbine

[0173] The typical physical model of a gas turbine can be expressed as:

[0174]

[0175]

[0176] Q MTh，t =η r C MTh Q MT，t (6)

[0177] Where, P MT，t P represents the output power of the gas turbine during period t; MTg，t represents the natural gas consumption of the gas turbine during period t; η MTe Indicates the power generation efficiency of the gas turbine; Q MT，t Indicates the waste heat discharged by the gas turbine; η L Indicates the heat loss coefficient; Q MTh，t Indicates the heating power output of the waste heat recovery unit; η r Indicates the flue gas recovery rate; C MTh Indicates the heating coefficient of the waste heat recovery unit.

[0178] 2.4 Gas Boiler

[0179] A gas boiler uses the input natural gas to generate heat energy at a certain efficiency. The mathematical model of a gas boiler is:

[0180] Q GB，t =η GB ·P GBg，t (7)

[0181] Where Q GB，tand P GBg，t They represent the heat output and natural gas input of the gas boiler during period t; η GB Indicates the efficiency of the gas boiler.

[0182] 2.5 Electric Boiler

[0183] The electric boiler inputs electricity and generates cooling energy. The mathematical model of the electric refrigerator can be expressed as follows:

[0184] Q EB，t =η EB ·P EB，t (8)

[0185] Where Q EB，t and P EB，t They represent the output thermal power and input electrical power of the electric boiler during period t respectively; η EB Indicates the heating efficiency of electric boiler.

[0186] 2.6 Power-to-gas equipment

[0187] The power-to-gas equipment model is as follows:

[0188] P P2Gg，t =η P2G ·P P2Ge，t (9)

[0189] Where, P P2Gg，t and P P2Ge，t They represent the P2G output gas power and input electric power during period t respectively; η P2G Indicates the P2G conversion efficiency.

[0190] 2.7 Energy Storage Devices

[0191] Energy storage devices have certain common characteristics. Their energy storage and release characteristics can be characterized by three parameters: device storage capacity, self-consumption rate, and storage and release efficiency. The energy storage state change of energy storage devices is expressed by the following formula:

[0192]

[0193] Where m represents the type of energy storage device, m=ES, HS represents electrical energy storage and thermal energy storage respectively; SOC m，t and SOC m，t-1 Respectively represent the energy storage status of the mth type energy storage device in the tth period and the t-1th period during operation; δ m represents the storage loss rate of the mth type of energy storage; and They represent the energy storage and discharge efficiency of the mth type of energy storage device respectively; and Binary variables representing the operating status of energy storage devices, P m，ch，t and Pm，dis，t They represent the charging and discharging power of the energy storage device respectively.

[0194] 3. Response strategies for thermal and electrical coupled park operators

[0195] The park operator is responsible for the interaction and coordination of the entire park's optimized scheduling process. When the upper power grid issues a peak-shaving instruction, the park operator decomposes the peak-shaving demand according to the response contract signed in advance with the industrial user, sends incentive signals to various users, and obtains the optimal scheduling plan and response strategy by coordinating the output of its own equipment and guiding various demand response resources. The interaction process between the upper power grid, park operators and users is as follows: Figure 3 shown.

[0196] 3.1 Objective Function

[0197] The thermal-electric coupled park operator response strategy aims to coordinate the output of various equipment within the park and the response plans of various industrial users, allowing the park operator to respond as fully as possible to the peak-shaving demand of the upper power grid while improving the park's renewable energy absorption capacity and the park operator's total revenue. Therefore, this method selects the park operator's response rate to the upper power grid's peak-shaving demand, the renewable energy absorption rate, and the park operator's total revenue as optimization objectives. Among them, if t0-t1 is the peak-shaving period, the park operator's response rate to the upper power grid's peak-shaving demand is:

[0198]

[0199] Where: Δ Pef，t P represents the peak load demand issued by the upper power grid during period t, that is, the power purchased by the park operator needs to be reduced; e，t 、P eDR，t The electricity purchase volume before and after the IDR of the park operator respectively;

[0200] The renewable energy consumption rate ψ is:

[0201]

[0202] Where: P WT，t and P PV，t Respectively represent the wind turbine and photovoltaic output power during period t; P Ab，t Indicates the amount of wind and solar power abandoned during period t.

[0203] The total revenue of the park operator is “energy sales revenue R P -Total cost C", where total cost C mainly includes the energy purchase cost C P , operation and maintenance costs C OM , carbon emission cost C CE and industrial user IDR compensation cost C IDR , as shown in equations (31)-(36).

[0204]

[0205] C=C P +C OM +C CE +C IDR (12)

[0206]

[0207]

[0208]

[0209]

[0210] Where, T represents the total operating hours, which can be 24 hours; r E，t 、r G，t and r H，t P is the electricity, gas and heat prices sold by operators to users during period t; be，t 、P bg，t They represent the amount of electricity and gas purchased by the park operator from the upper energy grid during period t; P se，t 、P sg，t 、P sh，t They represent the electricity, gas and heat power sold by the park operator to users during period t; c E，t 、c G，t They represent the electricity purchase price and gas purchase price in period t respectively; c om，j P represents the unit operation and maintenance cost of the jth device; j，t represents the output of the jth device in time period t; x represents the processing cost per unit mass of CO2; μ e and μ g Represent the equivalent carbon emission coefficients of electricity purchase and gas purchase respectively; γ sf , γ hd , γ etr and γ mf They respectively represent the cost incurred by shifting the dispatchable production task for one period of time, the cost of the additional unit heat load demand of the thermoelectric coupling user, the unit energy cost of the electrically replaceable load, and the unit energy conversion cost of the multi-energy flexible user.

[0211] Considering that the main purpose of this method is for the thermal-electric coupling park operator to respond to the peak-shaving demand of the upper power grid as much as possible and to improve the level of renewable energy consumption, the three optimization objectives mentioned above are prioritized. Among them, improving the response rate of the park operator to the peak-shaving demand of the upper power grid is the first-level optimization goal; improving the renewable energy consumption rate is the second-level optimization goal; and improving the total revenue of the park operator is the third-level optimization goal. Therefore, three priority factors are set: φ1 = 10000000, φ2 = 1000000, and φ3 = 1. Since the response rate and consumption rate are small, the priority factor has a large value, and the objective function can be expressed as:

[0212] maxΓ=φ1ζ+φ2ψ+φ3(R P -C) (37)

[0213] Where Γ represents the target value.

[0214] 3.2 Constraints

[0215] (1) Energy balance constraints

[0216] Including electrical power balance, thermal power balance and gas power balance.

[0217] P MTe，t +P PV，t +P WT，t +P be，t +(P ES，dis，t -P ES，ch，t )-P P2Ge，t -P EB，t =L e，t (38)

[0218] Q MTh，t +Q GB，t +Q EB，t +(Q HS，dis，t -Q HS，ch，t )=L h，t (39)

[0219] P bg，t +R P2Gg，t -P MTg，t -P GBg，t =L g，t (40)

[0220] Where, L e，t , L h，t and L g，t are the electricity, heat and gas load demands during period t; P MTe，t represents the output power of the gas turbine during period t; P PV，t represents the photovoltaic output power during period t; P WT，t represents the wind power output power during period t, Pbe，t P represents the amount of electricity purchased from the upper power grid during period t; ES，ch，t and P ES，dis，t They represent the charging and discharging power of the energy storage during period t; P P2Ge，t and P EB，t They represent the power consumption of the power-to-gas equipment and the electric boiler during period t; Q MTh，t , Q GB，t , Q EB，t They represent the heat output of gas turbine, gas boiler and electric boiler during period t respectively; Q HS，ch，t and Q HS，dis，t P represents the thermal energy storage charging and discharging power during period t; bg，t represents the gas purchase volume during period t, P P2Gg，t P represents the output of the power-to-gas equipment during period t; MTg，t and P GBg，t are the gas consumption of gas turbine and gas boiler during period t respectively.

[0221] (2) Equipment output constraints

[0222] 1) Energy storage device constraints:

[0223] E m，t =(1-μ m )E m，t-1 +(P m，ch，t η m，ch -P m，dis，t / η m，dis )Δt (41)

[0224]

[0225]

[0226] E m，min ≤E m，t ≤E m，max (44)

[0227] 0≤I m，ch，t +I m，dis，t ≤1 (45)

[0228] E m，0 =E m，T (46)

[0229] Where, E m，t-1 and E m，t is the energy storage capacity of the mth type energy storage device during periods t-1 and t; η m，ch and η m，dis represents the charging and discharging efficiency of the mth type of energy storage device; μ m represents the energy loss rate of the mth type energy storage device; P m，ch，t and P m，dis，tRespectively represent the charging and discharging power of the energy storage device; I m，ch，t and I m，dis，t It is a binary variable, and its value is 1, indicating that the m-th type of energy storage device is charging or discharging during period t. m，ch，t and I m，dis，t Cannot be 1 at the same time; and is the maximum charge and discharge energy of the mth type energy storage device; E m，max and E m，min Respectively represent the upper and lower limits of the storage capacity of the mth type of energy storage device; E m，0 、E m，T They represent the stored energy of the mth type of energy storage device at the initial moment and after one cycle respectively.

[0230] 2) Output constraints of other equipment:

[0231] P j，min ≤P j，t ≤P j，max (17)

[0232] Where, P j，t P represents the output of the j-th type of equipment in period t; j，max and P j，min They represent the upper and lower limits of the output of the jth type of equipment respectively.

[0233] (3) Tie line constraints

[0234] Park operators are restricted in purchasing energy from the power grid and gas network:

[0235] P be，min ≤P be，t ≤P be，max (18)

[0236] P bg，min ≤P bg，t ≤P bg，max (19)

[0237] Where, P be，max and P be，min Respectively represent the upper and lower limits of electricity purchase; P bg，max and P bg，min Respectively represent the upper and lower limits of gas purchase volume.

[0238] (4) Schedulable production task constraints

[0239] There may be production timing issues between different production tasks of industrial users, resulting in a time coupling relationship. Therefore, the production timing constraints for schedulable production tasks are established as follows:

[0240]

[0241]

[0242]

[0243]

[0244]

[0245]

[0246] Where, T back Indicates the time period that the schedulable production task n2 lags behind the schedulable production task n1. Formula (50) indicates that the schedulable production task n2 needs to be started after the schedulable production task n1 is started, and formulas (51)-(52) indicate that the schedulable production task n2 needs to be started after the schedulable production task n1 is completed. back The schedulable production task n2 and the schedulable production task n1 are started at the same time. Formula (53) means that the schedulable production task n2 and the schedulable production task n1 are started at the same time. Formula (54) means that the schedulable production task n1 and the schedulable production task n2 are run at the same time. Formula (55) means that the schedulable production task n1 and the schedulable production task n2 cannot be run at the same time.

[0247] (5) Constraints on IDR response volume for various industrial users

[0248] Electrically coupled industrial users:

[0249]

[0250]

[0251] Where, Indicates the upper and lower limits of the load change that can be replaced by electricity and gas during period t.

[0252] Electric and thermal coupling industrial users:

[0253]

[0254]

[0255]

[0256] Where, Indicates the maximum and minimum values converted into low-grade thermal power; Indicates the maximum and minimum values of the medium-grade thermal power that can be directly utilized; Indicates the maximum and minimum values that low-grade thermal power can respond to.

[0257] Multi-energy coupling industrial users:

[0258]

[0259]

[0260]

[0261] Where, Indicates the maximum and minimum values of the change in the amount of electricity purchased by users from the park; They represent the maximum and minimum changes in the amount of natural gas required caused by changes in the amount of electricity purchased by users; They respectively represent the maximum and minimum changes in the amount of waste heat generated by natural gas power generation.

[0262] 4 Solution method

[0263] The model built by this technology is essentially a mixed integer linear programming model, which can be solved by calling the YALMIP toolbox and CPLEX solver in the commercial software MATLAB. The standard form of the solution model is shown in formula (64):

[0264]

[0265] In the formula, maxΓ is the objective function. The model built by this technology takes the maximum total revenue of the park operator as the objective function; stg(x) = 0 is the constraint condition, including energy balance constraint, equipment output constraint, industrial user IDR response constraint and other constraints; x is the optimization variable, including the continuous variable x i and integer 0-1 variable x j , x i The minimum and maximum values of x are min and x max .

[0266] Example 2

[0267] In order to verify the effectiveness of the refined response strategy generation method for thermal-electrical coupled park operators proposed in the present invention, an industrial park is used as an example for analysis.

[0268] 1. Scenario Description

[0269] Taking an industrial park as an example, the electrical and thermal load demand curves of various industrial users in the park on a typical day, as well as the photovoltaic and wind turbine output curves are as follows: Figure 4-Figure 7Tables 1 and 2 show the equipment parameters and time-of-use electricity prices within the industrial park, respectively. Table 3 shows the dispatchable production task parameters for a traditional single-use industrial user. The upper-level power grid issues peak-shaving instructions: between 7:00 PM and 8:00 PM, the peak-shaving demands are 8,000 kW and 8,500 kW, respectively. A traditional single-use industrial user has 10 dispatchable production tasks, with a production task scheduling fee of 100 yuan per hour. These dispatchable production tasks must meet the following six constraints: each task must be completed within 24 hours; n2 must start after n1 or simultaneously; n3 and n4 must start simultaneously; n6 must start at least one hour after n5 has completed; n7 and n8 must not operate simultaneously; and n9 and n10 must operate simultaneously.

[0270] Table 1 Parameters of equipment in the industrial park

[0271]

[0272] Table 2 Time-of-use electricity prices

[0273]

[0274] Table 3 Scheduling production task parameters

[0275]

[0276] 2 Optimization results and analysis

[0277] 2.1 Optimization operation results

[0278] When the upper power grid does not issue peak shaving requirements, the original planned scheduling strategy of the park operator and the original planned operation of the dispatchable production tasks of traditional single industrial users are as follows: Figures 8-11 shown.

[0279] When the upper power grid issues peak load demand, the operator of the thermal and electrical coupled park will adopt a refined and comprehensive demand response strategy such as Figures 12-18 As shown, Figure 12-14 Provide the output of each device for the park operator, Figures 15-18 The response of various industrial users to the peak-shaving needs of park operators.

[0280] Depend on Figure 12-14As can be seen, during peak load shaving periods, the park operator coordinates the demand response resources of various industrial users and adjusts the output of various equipment to meet 100% of the upstream grid's peak load shaving needs, while also fully absorbing renewable energy. During off-peak hours (1:00-6:00 and midnight), wind turbine output is high. Due to the thermoelectric coupling characteristics of gas turbines, heat generated during these periods not only meets the thermal load but also generates electricity. During these periods, the electrical load and electricity required for electric boilers are primarily met by wind turbines and gas turbines. Because upstream grid electricity prices are lower during these periods, the park operator purchases electricity for P2G gas production and energy storage charging. During the normal hours (7:00-10:00), photovoltaic output gradually increases, coinciding with a gradual increase in electrical load demand and a decrease in thermal load demand. During these periods, the electrical load is primarily met by photovoltaics, wind turbines, and gas turbines, while thermal load demand is primarily met by the gas turbine alone. During normal hours, from 6:00 PM to 5:00 PM, the thermal load is met by heat released from thermal energy storage and the gas turbine. During this period, when grid electricity prices are low, the park operator purchases electricity from the upstream grid to meet the load and charge the energy storage. During peak hours, from 11:00 AM to 2:00 PM and from 7:00 PM to 9:00 PM, the energy storage is discharged while the gas turbine charges the thermal energy storage after meeting the thermal load. This increases the gas turbine output and reduces the amount of electricity purchased from the grid, allowing the park operator to meet the peak load regulation needs of the upstream grid. Natural gas required for gas loads, gas turbines, and gas turbine boilers is primarily met through purchased gas. P2G technology partially meets gas load demand during low electricity demand periods.

[0281] When the park operator adjusts the output of each device to respond to the peak shaving demand of the upper power grid, the response of various industrial users is required to be as follows: Figures 15-18 As shown. Figure 11 and Figure 15 It can be seen that the schedulable production tasks n2, n5, n6, n8, n9 and n10 of traditional single industrial users have been transferred. Figure 8 The original plan for electric power scheduling had wind and solar curtailment during the 1:00 period. The park operator shifted the start-up period of production tasks n2, n5, and n8 to the 1:00 period. This was mainly because the original planned start-up period of production tasks n2, n5, and n8 was close to the 1:00 period and the duration was short. Therefore, they were shifted to the 1:00 period, thereby improving the renewable energy consumption level in the park. At the same time, the production tasks n6, n9, and n10, which were originally planned to be in the peak period, were shifted to other periods. Among them, production task n6 was shifted to 21:00-24:00 due to the minimum start-up time limit, and n9 and n10 were shifted to another peak period of 11:00-14:00, in order to maximize the response rate to the peak-shaving demand of the upper power grid and ensure the economic efficiency of the park operation. Figure 16-Figure 18It can be seen that during peak-shaving periods, industrial users with electricity-heat coupling simultaneously increase their demand for medium- and low-grade heat loads, increasing gas turbine output and reducing the park operator's electricity purchases. Electrically coupled industrial users, based on electricity price signals, independently substitute electricity to meet their energy needs. Simultaneously, the park operator implements incentive-based electricity substitution for these industrial users, substituting electricity for gas during off-peak periods to reduce wind and solar curtailment and increasing gas substitution for electricity during peak periods to meet the response requirements of the upper-level grid. Multi-energy flexible industrial users increase their demand for heat loads during peak-shaving periods, while simultaneously increasing their demand for gas to reduce their demand for electricity, thereby meeting the response requirements of the upper-level grid.

[0282] 2.2 Comparative Analysis

[0283] In order to verify the effectiveness and economic superiority of the model proposed in this invention, three modes are set for simulation comparison and analysis:

[0284] Mode 1: No industrial user classification is performed, and only the electricity demand response within the industrial park is considered;

[0285] Mode 2: No industrial user classification is performed, and comprehensive demand response for electricity, heating, and gas within the industrial park is considered;

[0286] Model 3: Classify industrial users and consider the comprehensive demand response of electricity, heat and gas of various industrial users in the park, which is the complete model proposed in this invention.

[0287] Table 4 shows the response and operating costs of the park operator under the three models. It can be seen that under Model 1, the park operator can only meet 72.7% of the upstream grid's peak-shaving demand, and its renewable energy absorption rate is 98.8%. Model 2, which considers the park's comprehensive electricity and heat demand response, enables the park operator to fully respond to the upstream grid's peak-shaving demand, increasing the renewable energy absorption rate to 99.2%. The park operator's total daily operating revenue increases by 5.26% compared to Model 1. Model 3, which classifies each industrial user within the park and implements a refined response, achieves a 100% response to the upstream grid's peak-shaving demand while increasing its total daily operating revenue by 5.51% compared to Model 2 and 11.1% compared to Model 1. Furthermore, the park's renewable energy absorption rate under Model 3 increases to 100%. This shows that refined response based on the classification of industrial users within the park enables the park operator to fully respond to the upstream grid's peak-shaving demand, improves the park's renewable energy absorption level, and enhances the economic efficiency of system operation.

[0288] Table 4 Total revenue and response of park operators under different modes

[0289]

[0290]

[0291] 3 Conclusion

[0292] From the perspective of thermal-electrical coupled park operators, this paper categorizes industrial users within the park into three types: traditional single-energy type, dual-energy coupled type, and multi-energy flexible type. Response models are then established for each type. Based on these models, a comprehensive demand response strategy model for park operators is constructed, with the goals of maximizing the park operator's response rate to the peak-shaving demand of the upper power grid, the renewable energy absorption rate, and the park operator's total revenue. Ultimately, a refined response strategy for thermal-electrical coupled park operators is derived, leading to the following conclusions:

[0293] (1) Compared with the traditional demand response in which only electricity is involved, considering the comprehensive demand response of heat and electricity in the park can better meet the peak load demand of the upper power grid and improve the renewable energy consumption level of the park.

[0294] (2) Refined response to the classification of industrial users within the park can enable park operators to fully respond to the peak-shaving needs of the upper power grid, while significantly improving the level of renewable energy consumption and the economic efficiency of system operation.

[0295] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A comprehensive demand response method for a thermal-electrical coupled park, characterized in that: The following steps are involved: Step 1: Based on different energy consumption types, industrial users are divided into traditional single type, dual-energy coupling type, and multi-energy flexible type, and demand response models are established for each type of industrial user; Establish energy conversion models for various energy production equipment within the thermal-electrical coupled park; Step 2: Establish a comprehensive demand response model for the thermal and electrical coupled park that responds to the peak-shaving demand of the upper power grid, and use the response rate of the park operator to the peak-shaving demand of the upper power grid. , renewable energy consumption rate and the total revenue of the park operator as the objective function; represents the target value, Indicates the revenue from energy sales, represents the total cost, 、 、 represents three priority factors, and > > ; Energy balance constraints, equipment output constraints, tie line constraints, schedulable production task constraints, and IDR response constraints for various industrial users are used as constraints. Step 3: Solve the comprehensive demand response model of the thermal-electrical coupled park to obtain a comprehensive demand response plan for the thermal-electrical coupled park; In step 1, the demand response model of traditional single industrial users is: Where, n Represents the schedulable production task index, there are N production tasks; Represents schedulable production tasks n The start-up period in the actual production process; 、 、 and is a binary variable, , Indicates the schedulable production tasks in the actual production process n In the period t start up, Indicates the schedulable production tasks in the original production plan n In the period t In running state, Indicates the schedulable production tasks in the actual production process n In the period t In operation; T represents the scheduling period, Represents schedulable production tasks n The start-up period in the original production plan; Represents schedulable production tasks n The duration of operation; and Represent the schedulable production tasks n Upper and lower limits of the start period; Indicates that only production tasks can be scheduled after the response period t The electrical load demand, Represents schedulable production tasks n electricity demand; Indicates the period in the original production plan t Electric load demand; The demand response model of the response model of the electric-thermal coupling industrial user is: The electric-thermal coupling load model is: Where, express t The change in electrical load due to changes in heat load during a period of time; express t Heat load variation during the period; Indicates the conversion factor of thermoelectric units; represents the thermoelectric coupling coefficient; Heat load demand during non-demand response period: Where, Indicates time period t Medium-grade heat consumed by steam-driven equipment; Indicates the low-grade heat recovery efficiency of steam-driven equipment; Indicates time period t Total heat load demand of electric heating industrial users; 、 They represent the medium-grade heat load demand and low-grade heat load demand for direct utilization respectively; 、 They represent the conversion efficiency of heat supply to direct steam driven loads and low-grade thermal power equipment respectively; Heat load demand during demand response period: Where, Indicates the post-response period t Total heat load demand of electric heating industrial users; 、 They represent the response changes of the medium-grade thermal power directly utilized and the changes of the recovered and converted low-grade thermal power respectively; The electric load demand during the demand response period is: Where, Indicates the time period after the response of the electrothermal coupling type user t The electrical load demand, Indicates the period in the original production plan t Electric load demand; The demand response model of the response model of electrically coupled industrial users is: The electrical coupling load model is: Where, 、 The increase in electricity and gas loads caused by electricity / gas substitution during actual operation; Indicates the conversion coefficient of equal calorific value of electricity and natural gas; The electrical load substitution rate curve expression is as follows: Where, represents the load substitution rate, represents the electrical price difference, 、 and They represent the electrical substitution dead zone threshold, saturation zone threshold and maximum electrical load substitution rate respectively; After accounting for alternative load response t The time period fitting electric load and gas load are expressed as: Where, Indicates the period after the response of the electrically coupled user only t The electrical load demand, Indicates the period in the original production plan t Electric load demand; Indicates the period after the response of the electrically coupled user only t Gas load demand, Indicates the period in the original production plan t Gas load demand; The demand response model of the multi-energy flexible industrial user response model is: Where, Indicates the change in electricity purchased from the park; Indicates the change in natural gas required by industrial users to meet the change in electricity; Indicates the lower calorific value of natural gas; represents the efficiency of natural gas power generation; represents the time step; Indicates the change in heat load; represents the thermoelectric coupling coefficient; Indicates the conversion factor of thermoelectric units; Indicates the period after which only multi-functional flexible users respond t Electric load demand; Indicates the period in the original production plan t Electric load demand; Indicates the time period after the response t heat load requirements; Indicates the period in the original production plan t heat load requirements; Indicates the period after which only multi-functional flexible users respond t Gas load demand; Indicates the period in the original production plan t Gas load demand.

2. The comprehensive demand response method for a thermal-electrical coupled park according to claim 1, characterized in that: In step 2, the response rate of the park operator to the peak load demand of the upper power grid for: Where: Indicates that the upper power grid is t Peak load demand during the time period, i.e. the amount of power purchased by the park operator that needs to be reduced; 、 The electricity purchase volume before and after the IDR of the park operator respectively; Renewable energy consumption rate for: Where: and Respectively t Wind turbine and photovoltaic output power during the time period; express t The amount of wind and solar power abandoned during the period.

3. The comprehensive demand response method for a thermal-electrical coupled park according to claim 1, characterized in that: In step 2, the energy balance constraints include electric power balance, thermal power balance, and gas power balance, which can be expressed as: Where, 、 and They are t Electricity, heat and gas load demands during the time period; express t The gas turbine output electrical power during the period; express t Photovoltaic output power during the period; express t Wind power output power during the period, express t Purchase electricity from the higher-level power grid during the time period; and Respectively t The charging and discharging power of the time-slot electric energy storage; and Respectively t Power consumption of power-to-gas equipment and electric boilers during time periods; 、 、 Represents gas turbine, gas boiler and electric boiler respectively t Heat production during the period; and express t Thermal energy storage charging and discharging power during the time period; express t Gas purchase volume for the period, express t Output of power-to-gas equipment during time period; and They are t Gas consumption of gas turbines and gas boilers during each period.

4. The comprehensive demand response method for a thermal-electrical coupled park according to claim 1, characterized in that: In step 2, the energy storage device constraint conditions are: Where, and For the m Energy storage device t -1 and t Energy storage during the time period; and Indicates the m The charging and discharging efficiency of similar energy storage devices; Indicates the m Energy loss rate of energy storage devices; and Respectively represent the charging and discharging power of the energy storage device; and It is a binary variable, and its value is 1. t Time period m Energy storage devices charge or discharge energy, and Cannot be 1 at the same time; and For the m Maximum charge and discharge energy of energy storage devices; and Respectively represent m The upper and lower limits of the storage capacity of energy storage devices; 、 Respectively represent m The energy storage of a similar energy storage device at the initial moment and after one cycle.

5. The comprehensive demand response method for a thermal-electrical coupled park according to claim 1, characterized in that: In step 2, the schedulable production task constraints are: Where, Represents schedulable production tasks n 2. Lagging behind schedulable production tasks n 1 period.

6. The comprehensive demand response method for a thermal-electrical coupled park according to claim 1, characterized in that: In step 3, the standard form of the solution model is: Where, is the constraint condition, is a continuous variable, The minimum and maximum values of and , It is an integer 0-1 variable.

Images