A method and apparatus for optimizing the simulation operation of combined heat and power production

Abstract

This invention relates to the field of new energy power generation technology, specifically providing a method and apparatus for optimizing the simulation operation of combined heat and power (CHP) production. The method includes: solving a pre-constructed CHP simulation operation optimization model to obtain target values ​​for CHP simulation operation parameters; and using the target values ​​of the CHP simulation operation parameters as the optimization result. The pre-constructed CHP simulation operation optimization model includes an objective function with the goal of maximizing new energy power generation and its corresponding constraints. The technical solution provided by this invention can provide technical support for the assessment and calculation of the new energy absorption capacity of provincial power grids containing heating systems.

Description

Technical Field

[0001] This invention relates to the field of new energy power generation technology, specifically to a method and apparatus for optimizing the simulation operation of combined heat and power production. Background Technology

[0002] With the increasing proportion of renewable energy, the risk of renewable energy curtailment on the power grid is constantly intensifying. To meet winter heating demands, thermal power units need to maintain high output levels and a large number of units in operation during the heating season, thus compressing the space for renewable energy generation. As a result, renewable energy curtailment during the heating season has become the main cause of renewable energy curtailment in the "Three Norths" region (Northeast, North, and Northwest China). To improve the utilization rate of renewable energy during the heating season, fully utilizing the thermal inertia and heat storage capacity of the heating system, and reducing the operating capacity of thermal power units through combined heat and power optimization, are important technical approaches to improve the utilization rate of renewable energy during the heating season.

[0003] Conducting annual 8760-hour renewable energy production simulation calculations based on time-series production simulation technology is an important technical means to assess the renewable energy absorption capacity of the power grid, rationally optimize the medium- and long-term operation modes of the power grid and power sources, and promote renewable energy absorption. Current renewable energy absorption capacity assessment technologies do not fully consider the role of combined heat and power (CHP) optimization, making it difficult to accurately quantify and assess the renewable energy absorption capacity during the heating season. The heating system has a complex structure, including three main parts: heat source (heating units), heating network, and heat users. According to the principles of thermodynamics and energy conservation, the heat balance equation can be characterized by a series of complex nonlinear, time-delayed dynamic differential equations. Due to the large number of heating units and heat users in the power grid and the complex structure of the heating network, the complex dynamic differential equation model cannot meet the needs of rapid calculation for annual 8760-hour production simulation. Therefore, it is urgent to propose a CHP joint operation optimization modeling method suitable for annual time-series production simulation calculations to meet the needs of renewable energy absorption capacity assessment calculations. Summary of the Invention

[0004] To overcome the above-mentioned shortcomings, this invention proposes a method and apparatus for optimizing the simulation operation of combined heat and power production.

[0005] Firstly, a method for optimizing the simulation operation of combined heat and power (CHP) production is provided, the method comprising:

[0006] Solve the pre-built simulation operation optimization model of combined heat and power production to obtain the target values ​​of the simulation operation parameters of combined heat and power production.

[0007] The target values ​​of the simulated operation parameters of the combined heat and power production are used as the optimization results of the combined heat and power production simulation operation.

[0008] The pre-constructed combined heat and power (CHP) production simulation and optimization model includes an objective function aimed at maximizing new energy power generation and its corresponding constraints.

[0009] Preferably, the simulated operation parameters of the combined heat and power production include at least one of the following: power generation of extraction steam heating units, heating power of extraction steam heating units, number of operating extraction steam heating units, power generation of back-pressure heating units, heating power of back-pressure heating units, number of operating back-pressure heating units, power generation of new energy sources, heating power input to the heating network, heating power output from the heating network, heat exchange between indoor air and building walls, heat exchange between exterior walls and outdoor air, indoor temperature of heat users, and indoor wall temperature of heat users.

[0010] Preferably, the objective function is calculated as follows:

[0011]

[0012] In the above formula, obj is the objective value of the objective function, Δt is the unit time interval, and p w (t) represents the power generation of new energy sources during time period t, and T represents the control period.

[0013] Preferably, the constraints include power balance constraints, the mathematical model of which is as follows:

[0014]

[0015] In the above formula, Let be the power generation of the i-th type extraction steam heating unit during time period t. p represents the power generation of the j-th type back-pressure heating unit during time period t. d (t) represents the electrical load during time period t, I represents the number of types of extraction steam heating units, and J represents the number of types of back-pressure heating units.

[0016] Preferably, the constraints include heating unit heat balance constraints, and the mathematical model for the heating unit heat balance constraints is as follows:

[0017]

[0018] In the above formula, Let be the heating power of the i-th type extraction steam heating unit during time period t. h represents the heating power of the j-th type back-pressure heating unit during time period t. d (t) represents the heating power input to the heating network during time period t.

[0019] Furthermore, the constraints include heating network heat balance constraints, the mathematical model of which is as follows:

[0020]

[0021] In the above formula, h d (t-ΔT) represents the heating power input to the heating network during the time period t-ΔT, h o (t) represents the heat output of the heating network to the user side during time period t, η represents the heat loss of the heating network, and ΔT represents the heat delay of the heating network.

[0022] Furthermore, the formula for calculating the heat loss of the heating network is as follows:

[0023]

[0024] In the above formula, L is the pipe length, k, k1, and k2 are the first, second, and third coefficients, respectively, and t e The external temperature of the pipeline is denoted as u, and the flow velocity in the pipeline is denoted as u.

[0025] The calculation formulas for the first coefficient, the second coefficient, and the third coefficient are as follows:

[0026]

[0027] In the above formula, α1 is the convective heat transfer coefficient between the water flow and the pipe wall, L1 is the inner circumference of the pipe wall, F is the cross-sectional area of ​​the flow channel, c is the specific heat of water, ρ is the density of water, and F p Let c be the cross-sectional area of ​​the pipe. p ρ is the specific heat of the pipe wall. p α is the density of the pipe wall, α2 is the heat transfer coefficient of the pipe insulation layer, and L2 is the outer circumference of the pipe wall.

[0028] Furthermore, the formula for calculating the thermal delay of the heating network is as follows:

[0029]

[0030] In the above formula, L is the length of the pipe, r is the inner diameter of the pipe, and G... t This represents the water flow rate per unit time in the pipeline.

[0031] Furthermore, the constraints include heat user heat supply balance constraints, the mathematical model of which is as follows:

[0032]

[0033] In the above formula, H w (t) represents the heat exchanged between indoor air and building walls during time period t, H. l (t) represents the heat exchanged between the exterior wall and the outdoor air during time interval t. a (t) and T a (t+1) represents the indoor temperature at time t and time t+1, respectively. o(t) represents the outdoor temperature during time period t, T w (t) and T w (t+1) represents the wall temperature at time t and time t+1, respectively, and s represents the total heating area of ​​the building. and These represent the upper and lower limits of indoor temperature, ΔT. a K1 represents the upper limit of indoor temperature change in adjacent time periods, K2 represents the indoor heat transfer coefficient, K3 represents the wall heat storage coefficient, and K4 represents the outdoor heat transfer coefficient.

[0034] Furthermore, the calculation formulas for the indoor heat transfer coefficient K1, wall heat storage coefficient K2, outdoor heat transfer coefficient K3, and wall heat dissipation coefficient K4 are as follows:

[0035]

[0036] In the above formula, h is the height of each floor of the building, and c a It is the specific heat capacity of air, ρ a Let be the density of air, b be the ratio of the building's external surface area to its volume, α be the ratio of the building's interior wall area to its exterior wall area, x be the heat transfer coefficient between the indoor air and the walls, Δt be the time interval, and c be the density of air. w It is the specific heat capacity of the wall, ρ w For the density of the wall, δ w This refers to the thickness of the wall.

[0037] Furthermore, the constraints include thermoelectric coupling constraints for the extraction steam heating unit, and the mathematical model for these constraints is as follows:

[0038]

[0039] In the above formula, Let P be the number of type i extraction steam heating units in operation during time period t. A P represents the power generation at vertex A in the operating range of the extraction steam heating unit. B P represents the power generation at vertex B in the operating range of the extraction steam heating unit. C P represents the power generation at vertex C in the operating range of the extraction steam heating unit. D h represents the power generation at vertex D in the operating range of the extraction steam heating unit. A h represents the heating power at vertex A in the operating range of the extraction steam heating unit. B h represents the heating power at vertex B in the operating range of the extraction steam heating unit. C h represents the heating power at vertex C in the operating range of the extraction steam heating unit. D This represents the heating power at vertex D in the operating range of the extraction steam heating unit.

[0040] Furthermore, the mathematical model for the thermoelectric coupling constraint of the back-pressure heating unit is as follows:

[0041]

[0042] In the above formula, t represents the number of back-pressure heating units of type j operating during time period t, and h represents the height of each floor of the building.

[0043] Secondly, a combined heat and power (CHP) production simulation and optimization device is provided, the CHP production simulation and optimization device comprising:

[0044] The acquisition module is used to solve the pre-built optimization model for combined heat and power production simulation operation and obtain the target values ​​of the parameters for combined heat and power production simulation operation.

[0045] The simulation result output module is used to take the target values ​​of the combined heat and power production simulation operation parameters as the optimization results of the combined heat and power production simulation operation.

[0046] The pre-constructed combined heat and power (CHP) production simulation and optimization model includes an objective function aimed at maximizing new energy power generation and its corresponding constraints.

[0047] Thirdly, a storage medium is provided, the storage medium including a stored program, wherein, when the program is running, the device where the storage medium is located is controlled to execute the combined heat and power production simulation operation optimization method.

[0048] Fourthly, a processor is provided for running a program, wherein the program executes the combined heat and power production simulation operation optimization method during operation.

[0049] The above-described technical solutions of the present invention have at least one or more of the following beneficial effects:

[0050] This invention provides a method and apparatus for optimizing the simulation operation of combined heat and power (CHP) production, comprising: solving a pre-constructed CHP simulation operation optimization model to obtain target values ​​of CHP simulation operation parameters; and using the target values ​​of the CHP simulation operation parameters as the optimization result of the CHP simulation operation. The pre-constructed CHP simulation operation optimization model includes an objective function with the goal of maximizing renewable energy power generation and its corresponding constraints. The CHP simulation operation optimization model proposed in this invention simplifies the description of the operation process of complex heating systems. The model is a mixed-integer linear programming model, which can meet the practical production simulation calculation requirements of provincial-scale power grids. Attached Figure Description

[0051] Figure 1 This is a schematic diagram of the main steps of the combined heat and power production simulation operation optimization method according to an embodiment of the present invention;

[0052] Figure 2 This is a schematic diagram of the operating range of the extraction steam heating unit according to an embodiment of the present invention;

[0053] Figure 3 This is a schematic diagram of the thermoelectric coupling operation range of the S extraction steam heating units according to an embodiment of the present invention;

[0054] Figure 4 This is a schematic diagram of the thermoelectric coupling operation zone of the back-pressure heating unit according to an embodiment of the present invention;

[0055] Figure 5 This is a schematic diagram of the thermoelectric coupling operation range of the S-type back-pressure heating unit according to an embodiment of the present invention;

[0056] Figure 6 This is a main structural block diagram of the combined heat and power production simulation operation optimization device according to an embodiment of the present invention. Detailed Implementation

[0057] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0058] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0059] This invention proposes a method for simulating and optimizing the operation of combined heat and power (CHP) production. Heating systems involve complex operating characteristics of heat sources, heat networks, and heat users, exhibiting large-scale, nonlinear, and time-delayed characteristics, making modeling challenging. This invention extracts key features of the heating system, such as thermoelectric coupling characteristics, heat storage, heat dissipation, and thermal inertia, and utilizes techniques like equivalent aggregation to establish a CHP simulation optimization model that meets the needs of time-series production simulation. The proposed CHP simulation optimization model is a mixed-integer linear optimization model, capable of meeting the application requirements for rapid calculation in production simulation.

[0060] See appendix Figure 1 , Figure 1 This is a schematic flowchart illustrating the main steps of a combined heat and power (CHP) production simulation and optimization method according to an embodiment of the present invention. Figure 1 As shown, the combined heat and power production simulation operation optimization method in this embodiment of the invention mainly includes the following steps:

[0061] Step S101: Solve the pre-constructed combined heat and power production simulation operation optimization model to obtain the target values ​​of the combined heat and power production simulation operation parameters;

[0062] Step S102: Use the target values ​​of the combined heat and power production simulation operation parameters as the optimization results of the combined heat and power production simulation operation;

[0063] The pre-constructed combined heat and power (CHP) production simulation and optimization model includes an objective function aimed at maximizing new energy power generation and its corresponding constraints.

[0064] In this embodiment, the simulated operation parameters of the combined heat and power production include at least one of the following: power generation of extraction steam heating units, heating power of extraction steam heating units, number of operating extraction steam heating units, power generation of back-pressure heating units, heating power of back-pressure heating units, number of operating back-pressure heating units, power generation of new energy sources, heating power input to the heating network, heating power output from the heating network, heat exchange between indoor air and building walls, heat exchange between exterior walls and outdoor air, indoor temperature of heat users, and indoor wall temperature of heat users.

[0065] The constraints include at least one of the following: power balance constraints, heating unit heat supply balance constraints, heating network heat supply balance constraints, heat user heat supply balance constraints, extraction-type heating unit thermoelectric coupling constraints, and back-pressure heating unit thermoelectric coupling constraints.

[0066] Furthermore, the constraints may also include: system backup demand constraints, number of operating heating units constraints, output range constraints of heating units constraints, output ramp-up constraints of heating units constraints, and output range constraints of new energy sources.

[0067] In one implementation, the objective function is calculated as follows:

[0068]

[0069] In the above formula, obj is the objective value of the objective function, Δt is the unit time interval, and p w (t) represents the power generation of new energy sources during time period t, and T represents the control period.

[0070] In one implementation, the constraint conditions include a power balance constraint, the mathematical model of which is as follows:

[0071]

[0072] In the above formula, Let be the power generation of the i-th type extraction steam heating unit during time period t. p represents the power generation of the j-th type back-pressure heating unit during time period t. d (t) represents the electrical load during time period t, I represents the number of types of extraction steam heating units, and J represents the number of types of back-pressure heating units.

[0073] In one embodiment, the constraints include heating unit heat balance constraints, the mathematical model of which is as follows:

[0074]

[0075] In the above formula, Let be the heating power of the i-th type extraction steam heating unit during time period t. h represents the heating power of the j-th type back-pressure heating unit during time period t. d (t) represents the heating power input to the heating network during time period t.

[0076] In one embodiment, the constraints include heating network heat balance constraints, the mathematical model of which is as follows:

[0077]

[0078] In the above formula, h d (t-ΔT) represents the heating power input to the heating network during the time period t-ΔT, h o (t) represents the heat output of the heating network to the user side during time period t, η represents the heat loss of the heating network, and ΔT represents the heat delay of the heating network.

[0079] A heating network pipeline can be represented as a heating pipeline with a multi-layered insulation structure. Based on thermodynamic principles, the formula for calculating the heat loss of a heating network is as follows:

[0080]

[0081] In the above formula, L is the pipe length, k, k1, and k2 are the first, second, and third coefficients, respectively, and t e The external temperature of the pipeline is denoted as u, and the flow velocity in the pipeline is denoted as u.

[0082] The calculation formulas for the first coefficient, the second coefficient, and the third coefficient are as follows:

[0083]

[0084] In the above formula, α1 is the convective heat transfer coefficient between the water flow and the pipe wall, L1 is the inner circumference of the pipe wall, F is the cross-sectional area of ​​the flow channel, c is the specific heat of water, ρ is the density of water, and F p Let c be the cross-sectional area of ​​the pipe. p ρ is the specific heat of the pipe wall. pα is the density of the pipe wall, α2 is the heat transfer coefficient of the pipe insulation layer, and L2 is the outer circumference of the pipe wall.

[0085] The flow time of the liquid in the heating pipeline is approximately equal to the heat network delay time. Specifically, the calculation formula for the heat delay of the heating network is as follows:

[0086]

[0087] In the above formula, L is the length of the pipe, r is the inner diameter of the pipe, and G... t This represents the water flow rate per unit time in the pipeline.

[0088] In one implementation, the constraints include heat user heat supply balance constraints, the mathematical model of which is as follows:

[0089]

[0090] In the above formula, H w (t) represents the heat exchanged between indoor air and building walls during time period t, H. l (t) represents the heat exchanged between the exterior wall and the outdoor air during time interval t. a (t) and T a (t+1) represents the indoor temperature at time t and time t+1, respectively. o (t) represents the outdoor temperature during time period t, T w (t) and T w (t+1) represents the wall temperature at time t and time t+1, respectively, and s represents the total heating area of ​​the building. and These represent the upper and lower limits of indoor temperature, ΔT. a K1 represents the upper limit of indoor temperature change in adjacent time periods, K2 represents the indoor heat transfer coefficient, K3 represents the wall heat storage coefficient, and K4 represents the outdoor heat transfer coefficient.

[0091] The calculation formulas for the indoor heat transfer coefficient K1, wall heat storage coefficient K2, outdoor heat transfer coefficient K3, and wall heat dissipation coefficient K4 are as follows:

[0092]

[0093] In the above formula, h is the height of each floor of the building, and c a It is the specific heat capacity of air, ρ a Let be the density of air, b be the ratio of the building's external surface area to its volume, α be the ratio of the building's interior wall area to its exterior wall area, x be the heat transfer coefficient between the indoor air and the walls, Δt be the time interval, and c be the density of air. w It is the specific heat capacity of the wall, ρ w For the density of the wall, δ w This refers to the thickness of the wall.

[0094] Formulas 1 through 4 in the heat user heat supply balance constraint describe the thermal and heat dissipation characteristics of the heat user's building, reflecting the heat exchange process between the air inside the building and the outdoor air, including two aspects: heat exchange between the air inside the building and the walls; and heat exchange between the walls and the outdoor air. Formulas 5 and 6 in the heat user heat supply balance constraint reflect the heat user's thermal comfort needs.

[0095] In one embodiment, the constraints include thermoelectric coupling constraints for extraction-type heating units, and the mathematical model for these constraints is as follows:

[0096]

[0097] In the above formula, Let P be the number of type i extraction steam heating units in operation during time period t. A P represents the power generation at vertex A in the operating range of the extraction steam heating unit. B P represents the power generation at vertex B in the operating range of the extraction steam heating unit. C P represents the power generation at vertex C in the operating range of the extraction steam heating unit. D h represents the power generation at vertex D in the operating range of the extraction steam heating unit. A h represents the heating power at vertex A in the operating range of the extraction steam heating unit. B h represents the heating power at vertex B in the operating range of the extraction steam heating unit. C h represents the heating power at vertex C in the operating range of the extraction steam heating unit. D This represents the heating power at vertex D in the operating range of the extraction steam heating unit.

[0098] Among them, the operating range of extraction steam heating units is as follows Figure 2 As shown, the horizontal axis represents the unit's heating power, and the vertical axis represents the unit's power generation power. AB c represents the slope of line AB. BC c represents the slope of line BC. DC This represents the slope of line DC. When the unit is in the start-up state, the operating range of the extraction steam heating unit is the area enclosed by ABCD, which is... Figure 2 It can be seen that, while maintaining a certain heating power, the power generation can be adjusted within a certain range. For example, when the heating power is equal to h, the corresponding minimum power generation is P. E The maximum power generation capacity is P F That is, the generating capacity of the unit can be within the range [P] E P F[Internal variation. When the unit is in a shutdown state, the unit's operating range is limited to zero, that is, both heating power and power generation power are 0.]

[0099] Regarding the thermoelectric coupling constraint of the extraction steam heating unit in the mathematical model, assuming a total of N extraction steam heating units, S units are in the on-state during time period t, and N / S units are in the off-state. The operating states of the S on-state units are aggregated, and... Figure 3 For the combined thermoelectric coupling operating range of S extraction steam heating units, compared to the operating range of a single unit, the coordinates of the vertices of the ABCD interval become (p A S,h A S), (P B S,h B S), (P C S,h C S), (P D S,h D S).

[0100] Furthermore, the mathematical model for the thermoelectric coupling constraint of the back-pressure heating unit is as follows:

[0101]

[0102] In the above formula, t represents the number of back-pressure heating units of type j operating during time period t, and h represents the height of each floor of the building.

[0103] Figure 4 A schematic diagram of the thermoelectric coupling operation range of a back-pressure heating unit is shown. As can be seen from the diagram, the power generation of the back-pressure heating unit is limited by the heating power; thermoelectricity and power are mutually constrained, and the power generation and heating power have a linear relationship, meaning the heating power ranges from 0 to h. B The range of power generation variation is P A To P B For back-pressure heating units of the same capacity and type, an equivalent aggregated operation model is established. Assuming that S units are in operation during time period t, the operating status of these S units is aggregated, and... Figure 5 A schematic diagram of the thermoelectric coupling operating range of the S-type back-pressure heating units is shown.

[0104] Furthermore, the combined heat and power production simulation operation optimization model provided in this embodiment of the invention can be solved by calling the mathematical programming solver Cplex software.

[0105] Based on the same inventive concept, the present invention also provides a combined heat and power production simulation operation optimization device, such as... Figure 6 As shown, the combined heat and power production simulation operation optimization device includes:

[0106] The acquisition module is used to solve the pre-built optimization model for combined heat and power production simulation operation and obtain the target values ​​of the parameters for combined heat and power production simulation operation.

[0107] The simulation result output module is used to take the target values ​​of the combined heat and power production simulation operation parameters as the optimization results of the combined heat and power production simulation operation.

[0108] The pre-constructed combined heat and power (CHP) production simulation and optimization model includes an objective function aimed at maximizing new energy power generation and its corresponding constraints.

[0109] Preferably, the simulated operation parameters of the combined heat and power production include at least one of the following: power generation of extraction steam heating units, heating power of extraction steam heating units, number of operating extraction steam heating units, power generation of back-pressure heating units, heating power of back-pressure heating units, number of operating back-pressure heating units, power generation of new energy sources, heating power input to the heating network, heating power output from the heating network, heat exchange between indoor air and building walls, heat exchange between exterior walls and outdoor air, indoor temperature of heat users, and indoor wall temperature of heat users.

[0110] Furthermore, the constraints include at least one of the following: power balance constraints, heating unit heat supply balance constraints, heating network heat supply balance constraints, heat user heat supply balance constraints, extraction-type heating unit thermoelectric coupling constraints, and back-pressure heating unit thermoelectric coupling constraints.

[0111] Furthermore, the constraints include at least one of the following: system backup demand constraints, number of operating heating units constraints, output range constraints of heating units, output ramp-up constraints of heating units, and output range constraints of new energy sources.

[0112] Furthermore, the objective function is calculated as follows:

[0113]

[0114] In the above formula, obj is the objective value of the objective function, Δt is the unit time interval, and p w (t) represents the power generation of new energy sources during time period t, and T represents the control period.

[0115] Furthermore, the constraints include power balance constraints, the mathematical model of which is as follows:

[0116]

[0117] In the above formula, Let be the power generation of the i-th type extraction steam heating unit during time period t. p represents the power generation of the j-th type back-pressure heating unit during time period t. d (t) represents the electrical load during time period t, I represents the number of types of extraction steam heating units, and J represents the number of types of back-pressure heating units.

[0118] Furthermore, the constraints include heating unit heat balance constraints, the mathematical model of which is as follows:

[0119]

[0120] In the above formula, Let be the heating power of the i-th type extraction steam heating unit during time period t. h represents the heating power of the j-th type back-pressure heating unit during time period t. d (t) represents the heating power input to the heating network during time period t.

[0121] Furthermore, the constraints include heating network heat balance constraints, the mathematical model of which is as follows:

[0122]

[0123] In the above formula, h d (t-ΔT) represents the heating power input to the heating network during the time period t-ΔT, h o (t) represents the heat output of the heating network to the user side during time period t, η represents the heat loss of the heating network, and ΔT represents the heat delay of the heating network.

[0124] Furthermore, the formula for calculating the heat loss of the heating network is as follows:

[0125]

[0126] In the above formula, L is the pipe length, k, k1, and k2 are the first, second, and third coefficients, respectively, and t e The external temperature of the pipeline is denoted as u, and the flow velocity in the pipeline is denoted as u.

[0127] The calculation formulas for the first coefficient, the second coefficient, and the third coefficient are as follows:

[0128]

[0129] In the above formula, α1 is the convective heat transfer coefficient between the water flow and the pipe wall, L1 is the inner circumference of the pipe wall, F is the cross-sectional area of ​​the flow channel, c is the specific heat of water, ρ is the density of water, and F p Let c be the cross-sectional area of ​​the pipe. p ρ is the specific heat of the pipe wall. p α is the density of the pipe wall, α2 is the heat transfer coefficient of the pipe insulation layer, and L2 is the outer circumference of the pipe wall.

[0130] Furthermore, the formula for calculating the thermal delay of the heating network is as follows:

[0131]

[0132] In the above formula, L is the length of the pipe, r is the inner diameter of the pipe, and G... t This represents the water flow rate per unit time in the pipeline.

[0133] Furthermore, the constraints include heat user heat supply balance constraints, the mathematical model of which is as follows:

[0134]

[0135] In the above formula, H w (t) represents the heat exchanged between indoor air and building walls during time period t, H. l (t) represents the heat exchanged between the exterior wall and the outdoor air during time interval t. a (t) and T a (t+1) represents the indoor temperature at time t and time t+1, respectively. o (t) represents the outdoor temperature during time period t, T w (t) and T w (t+1) represents the wall temperature at time t and time t+1, respectively, and s represents the total heating area of ​​the building. and These represent the upper and lower limits of indoor temperature, ΔT. a K1 represents the upper limit of indoor temperature change in adjacent time periods, K2 represents the indoor heat transfer coefficient, K3 represents the wall heat storage coefficient, and K4 represents the outdoor heat transfer coefficient.

[0136] Furthermore, the calculation formulas for the indoor heat transfer coefficient K1, wall heat storage coefficient K2, outdoor heat transfer coefficient K3, and wall heat dissipation coefficient K4 are as follows:

[0137]

[0138] In the above formula, h is the height of each floor of the building, and c a It is the specific heat capacity of air, ρ a Let be the density of air, b be the ratio of the building's external surface area to its volume, α be the ratio of the building's interior wall area to its exterior wall area, x be the heat transfer coefficient between the indoor air and the walls, Δt be the time interval, and c be the density of air. w It is the specific heat capacity of the wall, ρ w For the density of the wall, δ w This refers to the thickness of the wall.

[0139] Furthermore, the constraints include thermoelectric coupling constraints for the extraction steam heating unit, and the mathematical model for these constraints is as follows:

[0140]

[0141] In the above formula, Let P be the number of type i extraction steam heating units in operation during time period t. A P represents the power generation at vertex A in the operating range of the extraction steam heating unit. B P represents the power generation at vertex B in the operating range of the extraction steam heating unit. C p represents the power generation at vertex C in the operating range of the extraction steam heating unit. D h represents the power generation at vertex D in the operating range of the extraction steam heating unit. A h represents the heating power at vertex A in the operating range of the extraction steam heating unit. B h represents the heating power at vertex B in the operating range of the extraction steam heating unit. C h represents the heating power at vertex C in the operating range of the extraction steam heating unit. D This represents the heating power at vertex D in the operating range of the extraction steam heating unit.

[0142] Furthermore, the mathematical model for the thermoelectric coupling constraint of the back-pressure heating unit is as follows:

[0143]

[0144] In the above formula, t represents the number of back-pressure heating units of type j operating during time period t, and h represents the height of each floor of the building.

[0145] Furthermore, the present invention provides a storage medium comprising a stored program, wherein, when the program is executed, the device containing the storage medium is controlled to perform the combined heat and power production simulation operation optimization method.

[0146] Furthermore, the present invention provides a processor for running a program, wherein the program executes the aforementioned cogeneration simulation operation optimization method during operation.

[0147] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0148] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0149] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0150] These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable apparatus for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0151] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A method for optimizing the simulation operation of combined heat and power production, characterized in that, The method includes: Solve the pre-built simulation operation optimization model of combined heat and power production to obtain the target values ​​of the simulation operation parameters of combined heat and power production. The target values ​​of the simulated operation parameters of the combined heat and power production are used as the optimization results of the combined heat and power production simulation operation. The pre-constructed combined heat and power (CHP) production simulation and optimization model includes: an objective function and constraints aimed at maximizing new energy power generation; The objective function is calculated as follows: In the above formula, The objective value of the objective function. For a unit of time period, for t The renewable energy power generation capacity during a given time period, where T is the control period; The constraints include the heating unit's heat supply balance constraint, and the mathematical model for the heating unit's heat supply balance constraint is as follows: In the above formula, For the first i Extraction-type heating units in t Heating power during the period For the first j Back-pressure type heating unit in t Heating power during the period for t The heating power input to the time-period heating network; The constraints include thermoelectric coupling constraints for extraction-type heating units, and the mathematical model for these constraints is as follows: In the above formula, For the first Extraction-type heating units in Number of units in operation during a time period This refers to the power generation at vertex A in the operating range of the extraction steam heating unit. This refers to the power generation at vertex B in the operating range of the extraction steam heating unit. This refers to the power generation at vertex C in the operating range of the extraction steam heating unit. This refers to the power generation at vertex D in the operating range of the extraction steam heating unit. This refers to the heating power supplied at vertex A in the operating range of the extraction steam heating unit. This refers to the heating power at vertex B in the operating range of the extraction steam heating unit. This refers to the heating power at vertex C in the operating range of the extraction steam heating unit. The heating power at vertex D in the operating range of the extraction steam heating unit; The mathematical model for the thermoelectric coupling constraint of the back-pressure heating unit is as follows: In the above formula, For the first Back-pressure type heating unit in Number of machines started during a given time period The height of each floor of the building.

2. The method as described in claim 1, characterized in that, The simulated operating parameters for combined heat and power production include at least one of the following: power generation of extraction steam heating units, heating power of extraction steam heating units, number of operating extraction steam heating units, power generation of back-pressure heating units, heating power of back-pressure heating units, number of operating back-pressure heating units, power generation of new energy sources, heating power input to the heating network, heating power output from the heating network, heat exchange between indoor air and building walls, heat exchange between exterior walls and outdoor air, indoor temperature of heat users, and indoor wall temperature of heat users.

3. The method as described in claim 1, characterized in that, The constraints include power balance constraints, the mathematical model of which is as follows: In the above formula, For the first i Extraction-type heating units in t Power generation during the period For the first j Back-pressure type heating unit in t Power generation during the period for t Electricity load during a given time period I This refers to the number of types of extraction-type heating units. J This refers to the number of back-pressure heating unit types.

4. The method as described in claim 1, characterized in that, The constraints include heating network heat balance constraints, and the mathematical model for these constraints is as follows: In the above formula, for t-△T The heating power input to the time-period heating network. For heating pipe network in The heat output to the user side during the time period For heat loss in the heating network, This refers to the thermal delay of the heating network.

5. The method as described in claim 4, characterized in that, The formula for calculating the heat loss of the heating network is as follows: In the above formula, For the length of the pipe, , , They are the first coefficient, the second coefficient, and the third coefficient, respectively. The external ambient temperature of the pipeline. The flow velocity in the pipe; The calculation formulas for the first coefficient, the second coefficient, and the third coefficient are as follows: In the above formula, The convective heat transfer coefficient between the water flow and the pipe wall. The inner circumference of the pipe wall. The cross-sectional area of ​​the flow channel. For the specific heat of water, The density of water, The cross-sectional area of ​​the pipe. For the specific heat of the pipe wall, The density of the pipe wall, The heat transfer coefficient of the pipe insulation layer, It is the outer circumference of the pipe wall.

6. The method as described in claim 4, characterized in that, The formula for calculating the heat delay of the heating network is as follows: In the above formula, The length of the pipe, r The inner diameter of the pipe. This represents the water flow rate per unit time in the pipeline.

7. The method as described in claim 1, characterized in that, The constraints include heat user heat supply balance constraints, and the mathematical model for these constraints is as follows: In the above formula, For indoor air and building walls in Heat exchange during a period of time For the exterior walls and outdoor air in Heat exchange during a period of time and They are respectively Time period and Indoor temperature during the period for Outdoor temperature during the time period, and They are respectively Time period and Wall temperature over time The total heating area of ​​the building. and These are the upper and lower limits of indoor temperature, respectively. This represents the upper limit of indoor temperature variation in adjacent time periods. The indoor heat transfer coefficient, The heat storage coefficient of the wall. The outdoor heat transfer coefficient, This is the heat dissipation coefficient of the wall.

8. The method as described in claim 7, characterized in that, The indoor heat transfer coefficient Wall heat storage coefficient Outdoor heat transfer coefficient Wall heat dissipation coefficient The calculation formula is as follows: In the above formula, The height of each floor of the building, It is the specific heat capacity of air. For the density of air, It is the ratio of a building's external surface area to its volume. This represents the ratio of the area of ​​the building's interior walls to the area of ​​its exterior walls. The heat transfer coefficient between indoor air and walls. For time intervals, It is the specific heat capacity of the wall. For the density of the wall, This refers to the thickness of the wall.

9. An apparatus based on the combined heat and power production simulation operation optimization method according to any one of claims 1-8, characterized in that, The device includes: The acquisition module is used to solve the pre-built optimization model for combined heat and power production simulation operation and obtain the target values ​​of the parameters for combined heat and power production simulation operation. The simulation result output module is used to take the target values ​​of the combined heat and power production simulation operation parameters as the optimization results of the combined heat and power production simulation operation. The pre-constructed combined heat and power (CHP) production simulation and optimization model includes an objective function aimed at maximizing new energy power generation and its corresponding constraints.

10. A storage medium, characterized in that, The storage medium includes a stored program, wherein, when the program is executed, it controls the device containing the storage medium to perform the combined heat and power production simulation operation optimization method according to any one of claims 1 to 8.

11. A processor, characterized in that, The processor is used to run a program, wherein the program executes the cogeneration simulation operation optimization method according to any one of claims 1 to 8.

Images