A method and system for analyzing the joint state of multiple energy flows in an integrated energy system

Abstract

The present invention belongs to the technical field of integrated energy system planning and provides a method and system for analyzing the combined state of multiple energy flows in an integrated energy system. The method comprises obtaining a heat transfer equivalent circuit between two nodes of a thermal system during steady-state operation based on a preset thermoelectric equivalent rule, thereby constructing a computational model of the thermal system; treating the power grid and heat network in the integrated energy system as a whole, combining the computational models of the power system and the thermal system, and constructing a unified thermoelectric model of the integrated energy system; and calculating the state quantities and node powers of the integrated energy system based on the unified thermoelectric model of the integrated energy system and initial state quantity data of the integrated energy system for use in integrated energy system planning.

Description

Technical Field

[0001] The present invention belongs to the technical field of integrated energy system planning, and in particular relates to a method and system for analyzing the joint state of multiple energy flows in an integrated energy system. Background Art

[0002] The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.

[0003] To achieve the planning, design, and operational optimization of integrated energy systems, it is crucial to conduct joint analysis and calculation of multiple energy flows. While the modeling of power networks is currently well established, with mature flow calculation methods, the heating network has yet to develop a unified calculation method. Existing technologies employ a power-hydraulic-thermal solution approach, decomposing and solving the system's hydraulic, thermal, and power flow equations or combining them into a comprehensive system of equations.

[0004] The inventors found that due to the heterogeneity of thermoelectric energy flows, the analysis and calculation methods of the thermoelectric coupling network cannot be unified, resulting in a large number of state variables in the entire system, high dimensionality and strong nonlinearity of the equation group, high calculation complexity, a large amount of calculation for each iteration, and difficulty in solving, which reduces the efficiency of the joint state analysis of multiple energy flows in the integrated energy system. Summary of the Invention

[0005] In order to solve the technical problems existing in the above-mentioned background technology, the present invention provides a method and system for joint state analysis of multi-energy flows in an integrated energy system. Based on the unified thermal and electric model of the integrated energy system, the method and system realize the unified analysis of heterogeneous multi-energy flows, which can reduce the computational complexity of the joint analysis of the power grid and the heat network while meeting the calculation accuracy, and lay the foundation for research on the planning, operation, collaborative optimization and scheduling of the integrated energy system.

[0006] In order to achieve the above object, the present invention adopts the following technical solutions:

[0007] A first aspect of the present invention provides a method for analyzing the joint state of multiple energy flows in an integrated energy system.

[0008] A method for analyzing the joint state of multiple energy flows in an integrated energy system, comprising:

[0009] Based on the preset thermoelectric equivalent rules, the heat transfer equivalent circuit between two nodes of the thermal system in steady-state operation is obtained, and the calculation model of the thermal system is constructed;

[0010] Considering the power grid and heat network in the integrated energy system as a whole, the calculation models of the power system and the thermal system are combined to construct a unified thermal and power model of the integrated energy system;

[0011] Based on the unified thermal power model of the integrated energy system and the initial data of the state quantity of the integrated energy system, the state quantity and node power of the integrated energy system are calculated for the planning of the integrated energy system.

[0012] As an implementation method, the preset thermoelectric equivalent rule is:

[0013] The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters;

[0014] The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters;

[0015] The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters;

[0016] The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters.

[0017] As an implementation method, in the preset thermal-electric equivalent rule, the heating network nodes are divided into heat load nodes, heat source nodes and intermediate nodes.

[0018] As an implementation method, for a heat load node, the heat load power and the load outlet water temperature are known;

[0019] For the heat source node, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged;

[0020] For the middle node, there is no temperature difference on both sides of the node and the thermal power is 0.

[0021] As an implementation method, the Newton-Raphson method is used to calculate the state quantity and node power of the integrated energy system.

[0022] The second aspect of the present invention provides a multi-energy flow joint state analysis system for an integrated energy system.

[0023] A multi-energy flow joint state analysis system for an integrated energy system, comprising:

[0024] The calculation model construction module of the thermal system is used to obtain the heat transfer equivalent circuit between two nodes of the thermal system during steady-state operation based on the preset thermoelectric equivalent rules, and to construct the calculation model of the thermal system;

[0025] A unified thermal power model construction module, which treats the power grid and thermal network in the integrated energy system as a whole, combines the computational models of the power system and the thermal system, and constructs a unified thermal power model of the integrated energy system;

[0026] The state analysis and calculation module is used to calculate the state quantity and node power of the integrated energy system based on the unified thermal power model of the integrated energy system and the initial state quantity data of the integrated energy system for use in the planning of the integrated energy system.

[0027] As an implementation method, the Newton-Raphson method is used to calculate the state quantity and node power of the integrated energy system.

[0028] As an implementation method, the preset thermoelectric equivalent rule is:

[0029] The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters;

[0030] The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters;

[0031] The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters;

[0032] The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters.

[0033] As an implementation method, in the preset thermal-electric equivalent rule, the heating network nodes are divided into heat load nodes, heat source nodes and intermediate nodes.

[0034] As an implementation method, for a heat load node, the heat load power and the load outlet water temperature are known;

[0035] For the heat source node, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged;

[0036] For the middle node, there is no temperature difference on both sides of the node and the thermal power is 0.

[0037] A third aspect of the present invention provides a computer-readable storage medium.

[0038] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps in the above-mentioned method for analyzing the joint state of multiple energy flows in an integrated energy system.

[0039] A fourth aspect of the present invention provides an electronic device.

[0040] An electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the above-described method for analyzing the joint state of multiple energy flows in an integrated energy system are implemented.

[0041] Compared with the prior art, the present invention has the following beneficial effects:

[0042] The present invention constructs a calculation model of the thermal system based on preset thermoelectric equivalence rules, and then combines the calculation models of the power system and the thermal system to construct a unified thermoelectric model of the integrated energy system, establishing connections between different energy fields and realizing unified analysis of heterogeneous multi-energy flows. Compared with traditional analysis methods, this method can reduce the computational complexity of the joint analysis of the power grid and the heat network while meeting the calculation accuracy, laying the foundation for research on the planning and operation, collaborative optimization and scheduling of the integrated energy system.

[0043] Advantages of additional aspects of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] The accompanying drawings, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0045] Figure 1 This is a flow chart of a method for analyzing the joint state of multiple energy flows in an integrated energy system according to an embodiment of the present invention;

[0046] Figure 2 is a heat transfer equivalent circuit between two nodes of a thermal system during steady-state operation according to an embodiment of the present invention;

[0047] Figure 3 This is a multi-energy flow joint analysis and calculation process based on the Newton-Raphson method in an embodiment of the present invention;

[0048] Figure 4 It is a structural diagram of a multi-energy flow joint state analysis system for an integrated energy system according to an embodiment of the present invention. DETAILED DESCRIPTION

[0049] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0050] It should be noted that the following detailed descriptions are illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which the present invention belongs.

[0051] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.

[0052] Example 1

[0053] Reference Figure 1 This embodiment provides a method for analyzing the joint state of multiple energy flows in an integrated energy system, which includes:

[0054] S101: Based on a preset thermoelectric equivalent rule, a heat transfer equivalent circuit between two nodes of a thermal system during steady-state operation is obtained, and a calculation model of the thermal system is constructed.

[0055] In the specific implementation process, the preset thermoelectric equivalent rule is:

[0056] The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters;

[0057] The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters;

[0058] The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters;

[0059] The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters.

[0060] The preset thermoelectric equivalence rule, shown in Table 1, rationally demonstrates the equivalent relationships between the basic thermoelectric state parameters and the equivalent relationships between thermoelectric power and energy parameters. This thermoelectric equivalence rule is more reasonable and comprehensive than the traditional thermoelectric equivalence rule, facilitating energy coordination analysis between thermoelectric and thermal power.

[0061] Table 1 Circuit comparison of heating network

[0062]

[0063] In Table 1, heat capacity flow CP = c p m. In the formula, c p is the specific heat capacity of water, in J / (kg·K); m is the mass flow rate of water in the pipe, in kg / s.

[0064] Based on the above thermoelectric analogy rules, the heat transfer equivalent circuit between two nodes of the thermal system during steady-state operation is as follows: Figure 2 shown.

[0065] T i is the node temperature at the head end of the heat network pipeline, T j is the temperature of the terminal node of the heat network pipeline, CP i is the heat capacity flow of the heat network pipeline head node, CP j is the heat capacity flow at the end node of the heat network pipeline.

[0066] According to the temperature drop equation of the heat transfer pipe and the thermoelectric analogy rules in Table 1, the node admittance matrix of the thermal system can be derived: for:

[0067]

[0068] Where, is the equivalent conductance (S) and susceptance (S) of the ij branch, T i 、CP i is the temperature (K) and heat capacity flow (J / sK) of the starting node of the heat network pipeline, T a is the ambient temperature (K), λ is the total heat transfer coefficient per unit length of the pipeline (W·m -1 K -1 ), L ij is the pipe length of the ij branch (m).

[0069] Based on the above equivalent rules, heating network nodes can be divided into heat load nodes, heat source nodes and intermediate nodes.

[0070] For the heat load node, the heat load power and load outlet water temperature are usually known in the actual thermal system.

[0071] For heat source nodes, the regional heating network cannot meet the heating needs of all heat users by only adjusting the outlet water temperature of the heat source. Therefore, the quantity regulation method is adopted, that is, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged.

[0072] For an unloaded intermediate node, there is no temperature difference on both sides of the node and the thermal power is 0.

[0073] The thermal system nodes are classified according to the known variables, as shown in Table 2.

[0074] Table 2 Node classification of thermal system

[0075] Node Classification Given amount Amount to be demanded Heat load node Thermal load power, load return water temperature Load water supply temperature, heat capacity flow rate Heat source node Heat source water outlet temperature Thermal output power, heat capacity flow rate Intermediate Node Thermal power Supply / return water temperature, heat capacity flow rate

[0076] Assume that the thermal system has n nodes, including m1 heat load nodes, m2 heat source nodes, and (n-m1-m2) intermediate nodes.

[0077] Each heating network node considers two variables: temperature T and heat capacity flow CP, so two equations are required. Analogous to the traditional DC power grid flow calculation model, combined with the thermoelectric analogy rules shown in Table 1, the thermal network calculation model can be derived as follows:

[0078]

[0079] The thermal network model is a universal steady-state calculation model for each heating network node. Compared to AC power networks, thermal networks do not include reactive power and phase angle, making their network calculation model simpler.

[0080] The thermal calculation model based on Newton's method can be obtained from formula (2):

[0081]

[0082] Where, ΔF h is the imbalance vector of the thermodynamic system; ΔCP i is the heat capacity flow imbalance at node i; Δφ i is the thermal power imbalance of node i; is the given amount of heat capacity flow at node i; is the thermal power given by node i.

[0083] The state variable x of the thermodynamic system h for:

[0084]

[0085] State correction value Δx of the thermal system h for:

[0086]

[0087] Jacobian matrix of the heat network:

[0088]

[0089] Where I is the identity matrix.

[0090] The thermal system analysis and calculation process is as follows:

[0091] a. Input the original data of the thermal system, including the node number of each pipe section and the parameters of the heat transfer equivalent circuit, the node number of each heat load, its return water temperature and the thermal power used, the node number of each heat source, its outlet water temperature, the node number of each intermediate node and its thermal power, etc.

[0092] b. Set the initial value of the node temperature to be calculated and initial value of heat capacity flow That is, setting the initial value of the state variables of the thermal system The permissible value ε of the given state variable error.

[0093] c. Calculate the nodal admittance matrix of the thermal system by formula (1)

[0094] d. Calculate the thermal power of the heating network node using formula (2)

[0095] e. Calculate the system imbalance using formula (3)

[0096] f. Calculate the Jacobian matrix from formula (6)

[0097] g. Calculate the state correction:

[0098] h. Update state:

[0099] i. Determine whether If not, return to step c and continue iterating until the accuracy requirement is met; if satisfied, end the iteration and Calculate the thermal power φ of the heating network node i , and finally output the calculation results.

[0100] In a combined heat and power energy system, a heat pump is a coupling device between the power system and the thermal system. It can transfer heat energy from a low-temperature heat source to a high-temperature heat source by consuming electricity to achieve heating. Its mathematical model is:

[0101] φ hp =η COP P hp (7)

[0102] Where, P hp is the input electrical power of the heat pump, φ hp is the heat output power of the heat pump, η COP is the heat energy conversion efficiency of the heat pump.

[0103] S102: Consider the power grid and heat network in the integrated energy system as a whole, combine the calculation models of the power system and the thermal system, and construct a unified thermal and electrical model of the integrated energy system.

[0104] Using a simultaneous solution approach, the power and heat networks are considered as a whole, constructing an unbalanced set of equations for the integrated energy network. The Jacobian matrix of the thermal-electric coupling is then derived and solved using the Newton-Raphson method. This simultaneous solution has a high model dimension, a long single iteration time, and a small number of iterations, but requires high initial values ​​and good convergence.

[0105] The equations of the power system and the thermal system are combined for unified iteration. The iterative formula of the Newton-Raphson method is as follows:

[0106]

[0107] Where k is the number of iterations;

[0108] x is the state variable of the system:

[0109]

[0110] θ i is the voltage phase angle of node i (deg), V i is the voltage amplitude at node i (kV).

[0111] Δx is the state correction of the system:

[0112]

[0113] ΔF is the system imbalance vector:

[0114]

[0115] ΔP i , ΔQ i is the active unbalance (MW) and reactive unbalance (MVar) of node i; is the active power injection (MW) and reactive power injection (MVar) given by node i; V i 、V j is the voltage amplitude of node i and node j (kV); G ij 、B ij is the equivalent conductance (S) and susceptance (S) of the ij branch; θ ij is the phase angle difference (deg) between the voltage at the beginning and end of the ij branch.

[0116] The Jacobian matrix J is derived from the imbalance matrix ΔF.

[0117]

[0118] The diagonal block elements J of the Jacobian matrix J e and J h Respectively represent the relationship between the power flow of the electrical and thermal subsystems and their own state quantities. The calculation formula is the same as when each subsystem is calculated separately. eh Indicates the impact of the power grid on the heat network through the thermoelectric coupling device, the matrix J he Indicates the impact of the thermal network on the power grid through thermoelectric coupling equipment.

[0119] The submatrix J representing the power system e for:

[0120]

[0121] For the off-diagonal elements of each block matrix (j≠i), we have:

[0122]

[0123] For the diagonal elements (j=i) of each block matrix, we have:

[0124]

[0125] The coupling device of the thermoelectric coupling link is a heat pump, and its mathematical model, equation (7), can be obtained:

[0126] Thp CP hp =η COP P hp (16)

[0127] Since the district heating network adopts the quantity regulation method, the outlet water temperature of the heat pump is T hp Keep unchanged, only its heat capacity flow rate CP hp Adjustment is made, so the electric-to-thermal sub-matrix J eh for:

[0128]

[0129] In the grid-connected mode, any power shortage in the grid is provided by the grid connection point, so the power balance node has no coupling with the thermal power. The partial derivative of the thermal power deviation with respect to the power variable is zero, so the thermal-to-electric submatrix J he for:

[0130]

[0131] Submatrix J representing the thermal system h for:

[0132]

[0133] S103: Based on the unified thermal power model of the integrated energy system and the initial state quantity data of the integrated energy system, the state quantity and node power of the integrated energy system are calculated for use in the planning of the integrated energy system.

[0134] During the specific implementation process, the Newton-Raphson method is used to calculate the state quantity and node power of the integrated energy system.

[0135] The calculation process of multi-energy flow joint analysis based on Newton-Raphson method is as follows: Figure 3 As shown:

[0136] 1) Input the raw data of the thermal electron system.

[0137] 2) Set the initial value of the system state variable x (0) ; Given the allowable value ε of the state variable error.

[0138] 3) Form the node admittance matrix of the power system and thermal system.

[0139] 4) Apply the initial value of the state variable x (0) Calculate the active and reactive power of electrical systems and the thermal power of thermal systems.

[0140] 5) Calculate the system imbalance ΔF using formula (11): (k) .

[0141] 6) Calculate the overall Jacobian matrix J by formula (12) (k) .

[0142] 7) Calculate the state correction value Δx by formula (8): (k+1) .

[0143] 8) Update state x (k+1) .

[0144] 9) Determine whether ||Δx is satisfied (k+1) ||<ε, if not satisfied, return to step 3) and continue iterating until the accuracy requirement is met; if satisfied, end the iteration and use x (k+1) Calculate the reactive power of the grid PV nodes, the active and reactive power of the grid balancing nodes, and the thermal power of the heating network nodes, and finally output the calculation results.

[0145] For operating integrated energy systems, analytical calculations can determine whether the grid bus voltage, branch current, and power exceed limits, and whether the heating network node temperature, flow rate, and thermal power are within specified ranges. If these do not meet requirements, timely measures should be taken to adjust operational methods. For integrated energy systems under planning, analytical calculations can provide a basis for selecting power supply plans and electrical equipment for the power system, and heating plans and thermal equipment for the thermal system. This analytical calculation method can also provide raw data for integrated energy system fault analysis and operational optimization.

[0146] This example studies a unified calculation method for thermal and power systems to facilitate the connection between different energy sectors and achieve a unified analysis of heterogeneous multi-energy flows. Compared with traditional analysis methods, this method reduces the computational complexity of the joint analysis of power and heat networks while maintaining computational accuracy, laying the foundation for research on the planning, operation, and coordinated optimization of integrated energy systems.

[0147] Example 2

[0148] Reference Figure 4 This embodiment provides a multi-energy flow joint state analysis system for an integrated energy system, which specifically includes the following modules:

[0149] (1) A computational model construction module for the thermal system, which is used to obtain the heat transfer equivalent circuit between two nodes of the thermal system during steady-state operation based on the preset thermoelectric equivalent rules, and to construct a computational model of the thermal system.

[0150] In the specific implementation process, the preset thermoelectric equivalent rule is:

[0151] The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters;

[0152] The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters;

[0153] The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters;

[0154] The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters.

[0155] Among them, in the preset thermal-electric equivalent rule, the heating network nodes are divided into heat load nodes, heat source nodes and intermediate nodes.

[0156] For the heat load node, the heat load power and load outlet water temperature are known;

[0157] For the heat source node, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged;

[0158] For the middle node, there is no temperature difference on both sides of the node and the thermal power is 0.

[0159] (2) A unified thermal power model construction module, which is used to treat the power grid and thermal network in the integrated energy system as a whole, combine the calculation models of the power system and the thermal system, and construct a unified thermal power model of the integrated energy system.

[0160] (3) A state analysis and calculation module, which is used to calculate the state quantity and node power of the integrated energy system based on the unified thermal power model of the integrated energy system and the initial state quantity data of the integrated energy system for use in the planning of the integrated energy system.

[0161] It should be noted here that the various modules in this embodiment correspond one-to-one to the various steps in Example 1, and the specific implementation processes are the same, which will not be repeated here.

[0162] Example 3

[0163] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the steps in the method for analyzing the joint state of multiple energy flows in an integrated energy system as described above are implemented.

[0164] Example 4

[0165] This embodiment provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for analyzing the joint state of multiple energy flows in an integrated energy system as described above are implemented.

[0166] The present invention is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products of the embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams and the combination of processes and / or blocks in the flowcharts 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, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0167] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection of the present invention.

Claims

1. A method for analyzing the joint state of multiple energy flows in an integrated energy system, characterized in that: include: Based on the preset thermoelectric equivalent rules, the heat transfer equivalent circuit between two nodes of the thermal system during steady-state operation is obtained, and the calculation model of the thermal system is constructed. The thermal network calculation model is: ; is the node temperature at the head end of the heat network pipeline, is the temperature of the terminal node of the heat network pipeline, is the heat capacity flow of the node at the head end of the heat network pipeline; is the nodal admittance matrix of the thermal system; The power grid and heat network in the integrated energy system are considered as a whole, and the computational models of the power system and thermal system are combined to construct a unified thermal and electrical model of the integrated energy system. A simultaneous solution method is used to consider the power grid and heat network as a whole, construct the unbalanced equations of the integrated energy network, and then derive the Jacobian matrix of the thermal and electrical coupling, which is solved using the Newton-Raphson method. Based on the unified thermal power model of the integrated energy system and the initial state data of the integrated energy system, the state quantities and node powers of the integrated energy system are calculated for use in the planning of the integrated energy system. The preset thermoelectric equivalent rule is: The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters; The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters; The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters; The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters.

2. The method for analyzing the multi-energy flow state of an integrated energy system according to claim 1, characterized in that: In the preset thermal-electric equivalent rule, the heating network nodes are divided into heat load nodes, heat source nodes and intermediate nodes.

3. The method for analyzing the multi-energy flow state of an integrated energy system according to claim 2, characterized in that: For the heat load node, the heat load power and load outlet water temperature are known; For the heat source node, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged; For the middle node, there is no temperature difference on both sides of the node and the thermal power is 0.

4. The method for analyzing the multi-energy flow state of an integrated energy system according to claim 1, wherein: The Newton-Raphson method is used to calculate the state quantity and node power of the integrated energy system.

5. A multi-energy flow joint state analysis system for an integrated energy system, characterized in that: include: The calculation model construction module of the thermal system is used to obtain the heat transfer equivalent circuit between two nodes of the thermal system during steady-state operation based on the preset thermoelectric equivalent rules, and construct the calculation model of the thermal system; the thermal network calculation model is: ; is the node temperature at the head end of the heat network pipeline, is the temperature of the terminal node of the heat network pipeline, is the heat capacity flow of the node at the head end of the heat network pipeline; is the nodal admittance matrix of the thermal system; A unified thermal power model construction module treats the power grid and thermal network in an integrated energy system as a whole, combines the computational models of the power system and thermal system, and constructs a unified thermal power model of the integrated energy system. This module uses a simultaneous solution method to treat the power grid and thermal network as a whole, constructs the unbalanced equations of the integrated energy network, and then derives the Jacobian matrix of the thermal power coupling, which is solved using the Newton-Raphson method. A state analysis and calculation module is used to calculate the state quantities and node powers of the integrated energy system based on the unified thermal power model of the integrated energy system and the initial state quantity data of the integrated energy system for use in the planning of the integrated energy system; The preset thermoelectric equivalent rule is: The voltage in the circuit parameters is analogous to the temperature in the thermal circuit parameters; The current in the circuit parameters is analogous to the heat capacity flow in the thermal circuit parameters; The electrical power in the circuit parameters is analogous to the thermal power in the thermal circuit parameters; The electrical energy in the circuit parameters is analogous to the thermal energy in the thermal circuit parameters; Or the Newton-Raphson method can be used to calculate the state quantity and node power of the integrated energy system.

6. The integrated energy system multi-energy flow joint state analysis system according to claim 5, characterized in that: In the preset thermal-electric equivalent rule, the heating network nodes are divided into heat load nodes, heat source nodes and intermediate nodes; For the heat load node, the heat load power and load outlet water temperature are known; For the heat source node, the mass flow rate of the pipeline is adjusted, and the outlet water temperature of the heat source node remains unchanged; For the middle node, there is no temperature difference on both sides of the node and the thermal power is 0.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps in the method for analyzing the joint state of multiple energy flows in an integrated energy system as described in any one of claims 1 to 4 are implemented.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the method for analyzing the joint state of multiple energy flows in an integrated energy system as described in any one of claims 1 to 4 are implemented.

Images