A method for tracking energy flow distribution of an integrated energy system

Abstract

The application discloses a kind of integrated energy system flow distribution analytical tracking methods, comprising: determining the network topology information of target integrated energy system, multi-energy pipeline parameters, multi-energy load position, multi-energy end position, energy station equipment parameters and operation mode;Through integrated energy system table system, non-equilibrium node load data is obtained, and integrated energy system equivalent flow mechanism model is established;According to direct flow calculation method, the flow distribution of integrated energy system is calculated, and integrated energy system road is established;On road, integrated energy system network lossless processing is carried out;The relationship between source end-node of integrated energy system is analyzed, the relationship between total flow of node i, flow out of node i and load flow on node i is constructed, and node-load flow distribution matrix D is obtained;The flow relationship between source end-node-load of integrated energy system is directly derived, and source end-load, source end-branch flow correlation matrix and distribution relationship are obtained.

Description

Technical Field

[0001] This invention belongs to the field of energy quality analysis technology for integrated energy systems, and relates to the fields of integrated energy systems, effective energy utilization of integrated energy systems, power systems, thermal systems, and gas systems. It particularly relates to an integrated energy system. Stream allocation parsing and tracing method. Background Technology

[0002] Currently, a unified national energy market has not yet been fully established, and numerous bottlenecks and pain points remain. Regional and industry-specific market segmentation remains a significant issue. The construction of a unified energy market is increasingly highlighting the commodity attributes of energy, necessitating the improvement of a fairer and more reasonable energy pricing mechanism. Furthermore, the advancement of "dual-carbon" energy conservation has led to the large-scale integration of various renewable energy sources, but the proportion and pathways for their integration remain unclear.

[0003] In past research, traditional power flow tracking algorithms have provided strong support for electricity market pricing and green energy consumption. However, with the continuous evolution of the energy internet and the development of multi-energy complementarity technologies, integrated energy systems are gradually becoming the main form of energy carrier for human society, placing higher demands on the pricing of multiple energy sources and the coordinated consumption of various renewable energy sources. Due to factors such as unreasonable cost-sharing mechanisms and conflicts of interest between energy sectors, the electricity, natural gas, and heat markets are basically priced separately, hindering the coordinated development of energy interconnection from an economic perspective. At the same time, both non-renewable and renewable energy sources vary in quality. A large amount of green electricity, green hydrogen, geothermal energy, and other renewable energy sources are being developed and utilized, with green electricity, mainly from photovoltaics and wind turbines, accounting for the largest proportion. How to promote the coordinated consumption of multiple renewable energy sources is a key issue in the process of advancing dual-carbon goals. Leveraging integrated energy systems... The flow allocation analysis and tracing method integrates the energy system's source-grid-load-storage components. The flow distribution relationship is presented intuitively to energy suppliers, energy users, and decision-making bodies, providing important technical support for the high-quality development of my country's energy system and the establishment of a unified energy market.

[0004] In summary, existing power flow tracing calculation methods have limitations in integrated energy systems involving the coupling of multiple energy sources such as electricity, gas, and heat. As a common attribute of multiple energy forms, it can uniformly characterize the effective energy of a comprehensive energy system and clarify the flow traces of effective energy in the network within a comprehensive energy system. Summary of the Invention

[0005] The purpose of this invention is to overcome the limitations of existing technologies that fail to clearly trace the flow of effective energy in a network within an integrated energy system, based on equivalent... The flow mechanism model establishes a comprehensive energy system A flow tracing model is proposed for an integrated energy system. The flow allocation parsing tracing method of this invention can be used to calculate the effective energy of an integrated energy system. The movement traces within the network inform energy pricing within an integrated energy system, constructing a system based on... It provides a basis for the integrated energy market of trading partners, clarifies the pathways and proportions for renewable energy consumption, and provides technical support for promoting renewable energy consumption.

[0006] The objective of this invention is achieved through the following technical solution:

[0007] An integrated energy system Stream allocation parsing and tracing methods include:

[0008] S1. Determine the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system;

[0009] S2. Through integrated energy systems The table system obtains the load of unbalanced nodes. Data, establishing an equivalent integrated energy system Flow mechanism model;

[0010] (201) The integrated energy system The metering system consists of devices installed at key nodes in an integrated energy system, capable of directly measuring the values ​​at each key node. Parameters, the Parameters include but are not limited to Momentum, Flow, load Source end

[0011] (202) The equivalent of the integrated energy system Flow mechanism models, including power system equivalents Flow mechanism model, natural gas system equivalent Flow mechanism model, thermodynamic system equivalent Flow mechanism model and energy station model:

[0012] S3. Based on direct Stream computing methods for calculating integrated energy systems Flow distribution and the establishment of an integrated energy system. Road; the integrated energy system A road is defined as consisting of branch roads, impedance, source, Each load Composed of components Flow loop;

[0013] S4. Integrated Energy System Road to integrate energy system network Damage reduction processing is performed to ensure the integrity of the line at both the beginning and end. Flow consistency;

[0014] S5. Construct the total number of nodes in the integrated energy system Flow vector, injection of all nodes Flow vector and source end on all nodes The relationship between flow vectors is established, forming a node-load relationship. Flow assignment matrix D;

[0015] S6. Derive the source-load and source-branch relationships of the integrated energy system. Stream allocation relationship.

[0016] Furthermore, in step S1, the network topology information of the integrated energy system includes the topology of the power system, natural gas system, hydrogen energy system, heating system, and cooling system; the multi-energy pipeline parameters of the integrated energy system include the line resistance, reactance, length, and model of the power system; the pipeline length, diameter, and model of the natural gas system; and the length, diameter, thermal conductivity, and roughness of the heating and cooling pipelines; the multi-energy load locations of the integrated energy system include the positions of the power load, natural gas load, hydrogen energy load, heating load, and cooling load in the topology; the multi-energy terminal locations of the integrated energy system include the positions of the power source, natural gas gate station, hydrogen energy gate station, heating source, and cooling source in the topology; the energy station equipment parameters and operating modes of the integrated energy system include the configuration, model, and output ratio of the energy conversion units inside the energy station; and the power source includes power plants and distributed new energy sources.

[0017] Furthermore, in step (202),

[0018] The power system equivalent In the flow mechanism model, power lines Flow is represented as nodes Pressure and the flow rate of medium in power system branches The product form:

[0019] In the formula: For power nodes Voltage, equal to node voltage

[0020] power lines Loss is represented as Pressure difference and flow rate of medium in power system branches The product form:

[0021]

[0022] In the formula: and These are the nodes at both ends of the line. The voltages are equal to the voltages at both ends of the node. and

[0023] Power node injection Flow e e,q Represented as nodes The product of voltage and node injected current:

[0024]

[0025] In the formula: Assuming the conjugate of the injected current into the node, and selecting the direction of current injection from the power system into the node as positive, then the corresponding e for the power source, electrical load, and connection node... e,q These represent negative values, positive values, and 0, respectively.

[0026] n e,l Power lines Flow vector e e n e,l Power lines Loss column vector Δe e n e,n Dimensional power node injection Flow vector e e,q They are represented as follows:

[0027]

[0028]

[0029]

[0030] In the formula: A e,- and A e These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the power system, both of order n. e,n ×n e,l ; and n e,n Power Node Press column vector and n e,n Conjugate column vector of injected current at a power node; For n e,l A conjugate column vector of current in a power line;

[0031] The natural gas system equivalent In the flow mechanism model, natural gas pipelines Flow is represented as nodes Pressure and natural gas system medium flow rate m g The product form:

[0032] e g =U g m g (37)

[0033] In the formula: U g For natural gas nodes Pressure, equal to node potential p g ;

[0034] Natural gas pipeline Loss is represented as Pressure differential and natural gas system medium flow rate m g The product form:

[0035] Δe g =(U g1 -U g2 )m g (38)

[0036] In the formula: U g1 and U g2 These are the nodes at both ends of the pipeline. The pressure is equal to the pressure at both ends of the node. potential p g1 and p g2 ;

[0037] Natural gas node injection Flow e g,q Represented as nodes The product of pressure and nodal injection gas flow rate:

[0038] e g,q =U g m g,q (39)

[0039] Where: m g,q To determine the gas flow rate for the node, the positive direction is chosen as the gas flow from the natural gas system into the node. At this point, the gas source, gas load, and the corresponding e of the connected node are... g,q These represent negative values, positive values, and 0, respectively.

[0040] n g,p Wei Natural Gas Pipeline Flow vector e g n g,p Wei Natural Gas Pipeline Loss column vector Δe g ng,n Natural gas node injection Flow vector e g,q They are represented as follows:

[0041]

[0042]

[0043] e g,q =diag(U g )m g,q (42)

[0044] In the formula: A g,- and A g These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the natural gas system, both of order n. g,n ×n g,p ;U g and m g,q n g,n Wei Natural Gas Node Press column vector and n g,n Column vector of natural gas node injection flow rate;

[0045] The thermal system is equivalent to In the flow mechanism model, the thermal branch Flow is represented as nodes Pressure and corresponding thermodynamic system branch medium flow rate m h The product of:

[0046] e h =U h m h (43)

[0047] Thermal branch Loss is represented by the nodes at both ends of the branch. Pressure difference and corresponding flow rate m of the medium in the thermal system branch h The product of:

[0048] Δe h =(U h1 -U h2 )m h (44)

[0049] In the formula: U h1 and U h2 The nodes at both ends of the thermal branch are respectively Pressure;

[0050] Thermal node injection Flow e h,q Represented as nodes The product of pressure and mass flow rate of injected water at the node:

[0051] e h,q =(p s -p r )m h,q =U h m h,q (45)

[0052] Where: m h,q The mass flow rate of water injected into the node, when the node is a heat load node, is m. h,q For heat load water flow rate m h,L When this node is a heat source node, m h,q The heat source water flow rate m h,S When the node is a connection node, m h,q =0;

[0053] n h,p Weiheli branch Flow vector e h n h,p Weiheli branch Loss column vector Δe h n h,n Thermal node injection Flow vector e h,q They are represented as follows:

[0054]

[0055]

[0056] e h,q =diag(U h )m h,q (48)

[0057] In the formula: A h,- and A h These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the water supply network, respectively; U h For n h,n Thermal nodes Compression column vector; m h,q m is the column vector of water flow rates at thermal nodes. h This corresponds to the flow rate vector of the pipeline.

[0058] The energy station model is based on nodes. Damage indicates energy station Flow model, this part Loss equals the electrical energy consumed or supplied by the energy station in the power system. e EH,e Fuel consumed or supplied in natural gas systems eEH,g Heat consumed or supplied in a thermal system e EH,h The sum, i.e., the nodal injection of the equivalent electricity, natural gas, and heat nodes of the energy station. Sum of flows:

[0059]

[0060] Furthermore, in step S4, the integrated energy system network is without... Damage treatment, The loss is borne by the nodes at both ends of the pipeline, as shown below:

[0061]

[0062] In the formula: e Node,noloss For none Damaged nodes Stream vector; eN is none Nodes before degradation Flow vector; A is the node-branch correlation matrix of the integrated energy system; Δe loss For pipelines Loss vector;

[0063] none Source node of integrated energy system after damage treatment Flow vector and load end node The flow vectors are represented as follows:

[0064]

[0065] In the formula: e sourse For source node Stream vector; e load For load end nodes Stream vector;

[0066] none Integrated energy system pipelines after damage treatment The flow vector is represented as:

[0067]

[0068] In the formula, e noloss For none Damaged pipeline Stream vector; e is zero Pipeline before damage Stream vector.

[0069] Furthermore, step S5 is detailed as follows:

[0070] Ensure node outflow flowing Balance, total outflow from nodes The vector relationship of the flow is represented as:

[0071] e Node -A2e noloss =e load (53)

[0072] In the formula: e Node For none Damaged nodes Stream vector; e noloss For none Damaged pipeline Flow vector; A2 is the node-branch incidence matrix of the outflow node; e load For load end nodes The flow vector can be zero;

[0073] Expanding equation (23), we get:

[0074]

[0075] Right now,

[0076] De Node =e load (55)

[0077] In the formula, D represents the node-load ratio. Stream allocation matrix.

[0078] Furthermore, step S6 is detailed as follows:

[0079] Integrated energy system source end The flow column vector is written as:

[0080] e source =diag(e source ) -1 [1;1;…;1] (56)

[0081] Obtain the source of the integrated energy system Stream vectors and nodes The relationship between the flow vectors is as follows:

[0082] e source =diag(e source ) -1 diag(e Node ) -1 e Node (57)

[0083] The node-load defined in step S5 The flow allocation matrix contains nodes. Flow and load Relationship e Node =(D) -1 e load e load For load end nodes The flow vector then yields the source end. and load The relationship between them is:

[0084] e source =diag(e source ) -1 diag(e Node ) -1 (D) -1 e load (58)

[0085] Therefore, the source-load of the integrated energy system The correlation matrix of a flow is represented as:

[0086] R S-L =diag(e source ) -1 diag(e Node ) -1 (D) -1 (59)

[0087] For branch roads If the flow allocation coefficient and the outflow node allocation coefficient are the same, then the source-branch correlation matrix is:

[0088] R S-P =-R S-L A1 (60)

[0089] In the formula: A1 is the node-branch correlation matrix of the injected node; based on the source-load... The correlation matrix R of the flow S-L Source-branch correlation matrix R S-P It can calculate the source-load and source-branch parameters of the integrated energy system without requiring any allocation principles. Flow allocation ratio.

[0090] This invention also provides an integrated energy system. The stream allocation parsing and tracing apparatus includes:

[0091] The data storage unit is used to store the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system.

[0092] Integrated energy system equivalent The flow mechanism model generation unit is used to generate models based on integrated energy systems. The table system obtains the load of unbalanced nodes. Data generates integrated energy system equivalent Flow mechanism model, equivalent of integrated energy system Flow mechanism models include power system equivalents Flow mechanism model, natural gas system equivalent Flow mechanism model, thermodynamic system equivalent Flow mechanism model and energy station model;

[0093] Calculation unit, used for direct Stream computing methods for calculating integrated energy systems Flow distribution and the establishment of an integrated energy system. road;

[0094] Processing unit for integrated energy system Road to integrate energy system network Damage reduction processing is performed to ensure the integrity of the line at both the beginning and end. Flow consistency;

[0095] Solver element, used to solve based on the total number of nodes Flow vector, outflow node Flow vectors and load on nodes The relationship between flow vectors forms the node-load relationship. The flow distribution matrix D ultimately leads to the derivation of the source-load and source-branch relationships of the integrated energy system. Stream allocation relationship.

[0096] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the program, implements the integrated energy system. The steps of the stream allocation parsing tracing method.

[0097] The present invention also provides a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when executed by a processor, implements the integrated energy system. The steps of the stream allocation parsing tracing method.

[0098] Compared with the prior art, the beneficial effects of the technical solution of the present invention are:

[0099] 1. Compared with traditional power flow tracing methods, this invention takes into account the quantitative and qualitative differences in different forms of energy within a comprehensive energy system, based on... As a common attribute of various energy forms, it can uniformly characterize the effective energy of a comprehensive energy system.

[0100] 2. This invention is based on equivalence A comprehensive energy system based on flow mechanism model. The flow distribution analysis and tracing method can clearly identify the flow traces of effective energy in the network of an integrated energy system, and can calculate the effective energy of each source end of the integrated energy system to branches and loads. The proportion of contributions allocated.

[0101] 3. This invention has broad application prospects. It comprehensively considers both the "quantity" and "quality" of energy to establish an equivalent integrated energy system. A flow mechanism model is proposed for a comprehensive energy system. The flow assignment analytical tracing model can solve for known non-equilibrium nodes in practical engineering. Integrated energy system under certain conditions Energy flow distribution relationships; with increasingly close energy interconnection, the energy market is placing higher demands on energy pricing mechanisms, taking into account the premise of meeting energy demand. Possessing commodity attributes, decision-makers As the trading partner, utilizing integrated energy systems The flow allocation analysis and tracking method establishes a reasonable, fair, and effective cost allocation method, breaks down industry barriers, and serves as a technical support for setting energy prices.

[0102] 4. Just as non-renewable energy sources vary in energy quality, so too do renewable energy sources. Currently, a large amount of renewable energy, such as green electricity, green hydrogen, and geothermal energy, is being developed and utilized. Among these, green electricity, primarily from photovoltaic and wind turbines, accounts for the largest proportion and also boasts the highest energy quality, making it ideal for integrated energy systems. The flow distribution analytical tracing method calculates the effective energy distribution relationship of the source-grid-load in the integrated energy system. It can clearly determine the distribution ratio and distribution path of different quality energy supplied by renewable energy to each link of the integrated energy system, so as to guide the coordinated consumption of multiple renewable energy sources. It also gives a "green label" to a large number of renewable energy sources such as green electricity, green hydrogen, and geothermal energy.

[0103] 5. With the advancement of the Energy Internet and the increasing demands of users for energy quality, higher requirements are being placed on the integrated energy supply side for the efficient maintenance of basic data for integrated energy system state estimation, as well as the online security analysis and multi-energy operation scheduling capabilities of integrated energy systems. The Energy Internet has proposed an integrated energy system energy management system (IES-EMS). In the future, IES-EMS will further consider improvements in energy quality or the addition of new energy sources. The stream analysis module will result in a large number of Data monitoring equipment and the complex defects it generates Data emerges, bad Data detection and identification will become a crucial function in integrated energy system state estimation, serving as a vital guarantee for online safety analysis and operational control of integrated energy systems. To ensure the accuracy of subsequent integrated energy system effective state estimation, [the following methods can be used]. Flow tracing calculations identify defective measurement data, through Flow tracing computation can filter and correct for defective measurements, reducing their impact on topology error identification and improving its accuracy. Attached Figure Description

[0104] Figure 1 For integrated energy systems Flowchart of the flow allocation parsing tracing calculation method.

[0105] Figure 2 A comprehensive energy system topology diagram for an example.

[0106] Figure 3 Integrated Energy System as an Example Road and Flow distribution diagram.

[0107] Figure 4 For example, the integrated energy system has no After damage Road and Flow distribution diagram. Detailed Implementation

[0108] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only for explaining the present invention and are not intended to limit the present invention.

[0109] This embodiment provides an integrated energy system. The flow allocation parsing tracing method, its calculation method is as follows: Figure 1 As shown.

[0110] S1. Determine the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system. The network topology information of the integrated energy system includes the topology structure of multiple energy systems such as the power system, natural gas (hydrogen) system, and heat (cold) system. The multi-energy pipeline parameters of the integrated energy system include the line resistance, reactance, length, and model of the power system; the pipeline length, diameter, and model of the natural gas system; and the length, diameter, thermal conductivity, and roughness of the heat (cold) pipeline. The multi-energy load locations of the integrated energy system include the positions of the power load, natural gas (hydrogen) load, and heat (cold) load in the topology. The multi-energy terminal locations of the integrated energy system include the positions of the power source (power plant, distributed new energy), natural gas (hydrogen) gate station, and heat (cold) source in the topology. The energy station equipment parameters and operating modes of the integrated energy system include the configuration, model, and output ratio of the energy conversion units within the energy station.

[0111] S2. Through integrated energy systems The table system obtains the load of unbalanced nodes. Data, establishing an equivalent integrated energy system Flow mechanism model.

[0112] (202) The integrated energy system The metering system consists of devices installed at key nodes in an integrated energy system, capable of directly measuring the values ​​at each key node. Parameters (including but not limited to) Momentum, Flow, load Source end );

[0113] (202) The equivalent of the integrated energy system Flow mechanism models, including power system equivalents Flow mechanism model, natural gas system equivalent Flow mechanism model, thermodynamic system equivalent Flow mechanism model and energy station model:

[0114] (2021) Power System Equivalent In the flow mechanism model, power lines A flow can be represented as a node Pressure and medium flow rate The product form of (power system branch line currents):

[0115]

[0116] In the formula: For power nodes Voltage, equal to node voltage

[0117] power lines Loss can be expressed as Pressure difference and medium flow rate The product form of (power system branch line currents):

[0118]

[0119] In the formula: and These are the nodes at both ends of the line. The voltages are equal to the voltages at both ends of the node. and

[0120] Power node injection Flow e e,q Can be represented as nodes The product of voltage and node injected current:

[0121]

[0122] In the formula: Assuming the conjugate of the injected current into the node, and selecting the direction of current injection from the power system into the node as positive, then the corresponding e for the power source, electrical load, and connection node... e,q The values ​​are negative, positive, and 0, respectively.

[0123] n e,l Power lines Flow vector e e n e,l Power lines Loss column vector Δe e n e,n Dimensional power node injection Flow vector e e,q They can be represented as:

[0124]

[0125]

[0126]

[0127] In the formula: A e,- and A e These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the power system, both of order n. e,n ×n e,l ; and n e,n Power Node Press column vector and ne,n Conjugate column vector of injected current at a power node; For n e,l A conjugate column vector of current in a power line.

[0128] (2022) Natural Gas System Equivalent In the flow mechanism model, natural gas pipelines A flow can be represented as a node Pressure and medium flow rate m g The product form of (natural gas system branch flow rates):

[0129] e g =U g m g *67)

[0130] In the formula: U g For natural gas nodes Compression, equal to node potential p g .

[0131] Natural gas pipeline Loss can be expressed as Pressure difference and medium flow rate m g The product form of (natural gas system branch flow rates):

[0132] Δe g =(U g1 -U g2 )m g (68)

[0133] In the formula: U g1 and U g2 These are the nodes at both ends of the pipeline. The pressure is equal to the pressure at both ends of the node. potential p g1 and p g2 .

[0134] Natural gas node injection Flow e g,q Can be represented as nodes The product of pressure and nodal injection gas flow rate:

[0135] e g,q =U g m g,q (69)

[0136] Where: m g,q Similar to the reference direction of water flow in a thermal system, this paper selects the direction of gas flow from the natural gas system into the node as positive, and considers the gas source, gas load, and e corresponding to the connecting node as the positive direction. g,q The values ​​are negative, positive, and 0, respectively.

[0137] n g,p Wei Natural Gas Pipeline Flow vector e g n g,p Wei Natural Gas Pipeline Loss column vector Δe g n g,n Natural gas node injection Flow vector e g,q They can be represented as:

[0138]

[0139]

[0140] e g,q =diag(U g )m g,q (72)

[0141] In the formula: A g,- and A g These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the natural gas system, both of order n. g,n ×n g,p U g and m g,q n g,n Wei Natural Gas Node Press column vector and n g,n Column vector of natural gas node injection flow rate. (2023) Equivalent thermal system In the flow mechanism model, the thermal branch A flow can be represented as a node Pressure and corresponding medium flow rate m h The product of (the flow rates of the branch water in the thermal system):

[0142] e h =U h m h (73)

[0143] Thermal branch Loss can be represented as the nodes at both ends of a branch. Pressure difference and corresponding medium flow rate m h The product of (the flow rates of the branch water in the thermal system):

[0144] Δe h =(U h1 -U h2 )m h (74)

[0145] In the formula: U h1 and Uh2 The nodes at both ends of the thermal branch are respectively Pressure.

[0146] Thermal node injection Flow e h,q Can be represented as nodes The product of pressure and mass flow rate of injected water at the node:

[0147] e h,q =(p s -p r )m h,q =U h m h,q (75)

[0148] Where: m h,q The mass flow rate of water injected into the node. When the node is a heat load node, m h,q For heat load water flow rate m h,L When this node is a heat source node, m h,q The flow rate of the heat source water is m h,S When the node is a connection node, m h,q It is 0.

[0149] n h,p Weiheli branch Flow vector e h n h,p Weiheli branch Loss column vector Δe h n h,n Thermal node injection Flow vector e h,q They can be represented as:

[0150]

[0151]

[0152] e h,q =diag(U h )m h,q (78)

[0153] In the formula: A h,- and A h These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the water supply network, respectively; U h For n h,n Thermal nodes Compression column vector; m h,q m is the column vector of water flow rates at thermal nodes. h This is the corresponding pipe flow rate vector.

[0154] (2024) Energy Station Model, based on nodes Damage indicates energy station Flow model, this part Loss equals the electrical energy consumed or supplied by the energy station in the power system. e EH,e Fuel consumed or supplied in natural gas systems e EH,g Heat consumed or supplied in a thermal system e EH,h The sum, i.e., the nodal injection of the equivalent electricity, natural gas, and heat nodes of the energy station. Sum of flows:

[0155]

[0156] S3. Based on direct Stream computing methods for calculating integrated energy systems Flow distribution and the establishment of an integrated energy system. Road; the integrated energy system A road, analogous to the analysis approach of an electrical system based on electrical "circuits," is defined as a branch circuit. impedance, source, Load, etc. Composed of components Flow loop.

[0157] S4. Road to integrate energy system network Damage reduction processing is performed to ensure the integrity of the line at both the beginning and end. Flow consistency;

[0158] Integrated energy system network Damage treatment, The loss is borne by the nodes at both ends of the pipeline, which can be represented as:

[0159]

[0160] In the formula: e Node,noloss For none Damaged nodes Stream vector; e N For none Nodes before degradation Flow vector; A is the node-branch correlation matrix of the integrated energy system; Δe loss For pipelines Loss vector.

[0161] none Source node of integrated energy system after damage treatment Flow vector and load end node The flow vectors are represented as follows:

[0162]

[0163] In the formula: e sourse For source node Stream vector; e load For load end nodes Stream vector.

[0164] none Integrated energy system pipelines after damage treatment The flow vector is represented as:

[0165]

[0166] In the formula, e noloss For none Damaged pipeline Stream vector; e is zero Pipeline before damage Stream vector.

[0167] S5. Total number of nodes constructed Flow vector, outflow node Flow vectors and load on nodes The relationship between flow vectors is used to form node-load relationships. Flow allocation matrix D; includes: guaranteeing node outflow flowing Balance, total outflow from nodes The vector relationship of a flow can be represented as:

[0168] e Node -A2e noloss =e load (83)

[0169] In the formula: e Node For none Damaged nodes Flow vector; A2 is the node-branch incidence matrix of the outflow node; e load It can be zero.

[0170] Expanding equation (53), we get:

[0171]

[0172] Right now,

[0173] De Node =e load (85)

[0174] In the formula, D represents the node-load ratio. Stream allocation matrix.

[0175] S6. Derive the source-load and source-branch relationships of the integrated energy system. Flow distribution relationship. Specifically: source end of integrated energy system. The flow vector can be written as:

[0176] e source =diag(e source ) -1 [1;1;…;1] (86)

[0177] Furthermore, the source end of the integrated energy system is obtained. Stream vectors and nodes The relationship between the flow vectors is as follows:

[0178] e source =diag(e source ) -1 diag(e Node ) -1 e Node (87)

[0179] The node-load defined in step S5 The flow allocation matrix contains nodes. Flow and load Relationship e Node =(D) -1 e load Then the source end can be obtained. and load The relationship between them is:

[0180] e source =diag(e source ) -1 diag(e Node ) -1 (D) -1 e load (88)

[0181] Therefore, the source-load of the integrated energy system The correlation matrix of a flow can be represented as:

[0182] R S-L =diag(e source ) -1 diag(e Node ) -1 (D) -1 (89)

[0183] For branch roads If the flow allocation coefficient and the outflow node allocation coefficient are the same, then the source-branch correlation matrix is:

[0184] R S-P =-R S-L A1(90)

[0185] In the formula: A1 is the node-branch association matrix of the injected node.

[0186] According to source-load The correlation matrix R of the flow S-L Source-branch correlation matrix R S-P This allows for the calculation of the source-load and source-branch parameters of the integrated energy system without requiring any allocation principles. Flow allocation ratio.

[0187] Specifically, the following examples illustrate the aforementioned integrated energy system. The method for stream allocation parsing and tracing is further explained below:

[0188] The integrated energy system proposed in this embodiment Measurement methods and An example of the table system applied to an integrated energy system with 6 power nodes, 8 thermal nodes, and 5 natural gas nodes is presented to illustrate the integrated energy system. Measurement and verification of the effectiveness of the present invention.

[0189] To unify the electricity-gas-heat system of the integrated energy system Flow tracing is used to perform an equivalent transformation on a typical integrated energy system example. The integrated energy system topology used in this embodiment is as follows: Figure 2 As shown, the system includes a 6-node power system, a 5-node natural gas system, and an 8-node heating system. The power system is based on modified IEEE 33-node distribution network parameters, with a source voltage of 12.66 kV. The medium-pressure natural gas system is based on modified 5-node natural gas system parameters, with a gas source pressure of 5 bar, a theoretical natural gas combustion temperature of 1973℃, and a calorific value of 45.75 MJ / m³. The heating system is based on modified 8-node heating system parameters, with dual heat sources supplying heat at 100℃ for both sources and a load outlet temperature of 50℃. The example includes two energy stations: nodes E6, G9, and H12 are coupled to a combined heat and power (CHP) unit in energy station 1, while nodes E5 and H13 are coupled to an electric boiler (EB) in energy station 2. The gas-to-electricity and gas-to-heat conversion efficiencies of the CHP in energy station 1 are 0.3 and 0.4, respectively, while the electricity-to-heat conversion efficiency of the EB in energy station 2 is 0.95. The system operates using a "heat-driven power generation" method. The pipeline parameters for the power, natural gas, and heating systems are shown in Tables 1-3 of Appendix A, respectively. Electrical load... Natural gas load Heat load The data is shown in Appendix Table 4-6, as... Input data for the stream direct computation method.

[0190] Table 1 Power Line Parameters

[0191]

[0192] Table 2 Parameters of Thermal Pipelines

[0193]

[0194] Table 3 Natural Gas Pipeline Parameters

[0195]

[0196] Table 4 Electrical Load

[0197]

[0198] Table 5 Heat Load

[0199]

[0200]

[0201] Table 6 Natural Gas Load

[0202]

[0203] The integrated energy system proposed in this invention The stream allocation parsing and tracing method first performs the following steps on the implementation example: Stream computing yields a comprehensive energy system. Road and Flow distribution as Figure 3 As shown, Energy Station 1 and Energy Station 2 The losses are 687.931kW and 694.164kW respectively, with the numbers in parentheses representing the power supplied from the source. or pipeline The numbers near the boxes indicate the pipeline. Loss, the number pointed to by the arrow represents the load consumption.

[0204] Secondly, none Integrated energy system after damage Road and Flow distribution as Figure 4 As shown, the energy station's The loss is considered as the load on nodes E5 and E6, and on all pipelines. All losses are borne by adjacent nodes in the thermal system. Damage is considered a load After no The overall integrated energy system and its individual nodes after damage treatment The flow remains balanced and can adapt to the principle of proportional distribution. Stream tracing can also be used for verification. The correctness of flow tracing.

[0205] according to Figure 1 The aforementioned calculation steps yield the integrated energy system. The calculation results of the flow allocation parsing tracing method are shown in Tables 7 and 8. It can be seen that all source ends The sum of the distribution factors for any load is always equal to 1. For example, for node H19: 0.371 + 0.629 = 1. This indicates that load node H19 receives its load from source nodes E1 and G7. The flow rates were 37.1% and 62.9%, respectively. These results verify the effectiveness and correctness of the present invention.

[0206] Table 7 Source of the stream allocation parsing tracing method -load Allocation coefficient

[0207]

[0208]

[0209] Table 8 Source-branch of the flow distribution parser Flow distribution (kW)

[0210]

[0211]

[0212] Preferably, this embodiment also provides an integrated energy system. The stream allocation parsing and tracing apparatus includes:

[0213] The data storage unit is used to store the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system.

[0214] Integrated energy system equivalent The flow mechanism model generation unit is used to generate models based on integrated energy systems. The table system obtains the load of unbalanced nodes. Data generates integrated energy system equivalent Flow mechanism model, equivalent of integrated energy system Flow mechanism models include power system equivalents Flow mechanism model, natural gas system equivalent Flow mechanism model, thermodynamic system equivalent Flow mechanism model and energy station model;

[0215] Calculation unit, used for direct Stream computing methods for calculating integrated energy systems Flow distribution and the establishment of an integrated energy system. road;

[0216] Processing unit for integrated energy system Road to integrate energy system network Damage reduction processing is performed to ensure the integrity of the line at both the beginning and end. Flow consistency;

[0217] Solver element, used to solve based on the total number of nodes Flow vector, outflow node Flow vectors and load on nodes The relationship between flow vectors forms the node-load relationship. The flow distribution matrix D ultimately leads to the derivation of the source-load and source-branch relationships of the integrated energy system. Stream allocation relationship.

[0218] Preferably, embodiments of this application also provide an integrated energy system capable of implementing the above embodiments. A specific implementation of an electronic device for all steps in the stream allocation parsing and tracing method, the electronic device specifically including the following:

[0219] Processor, memory, communications interface, and bus;

[0220] The processor, memory, and communication interface communicate with each other via a bus; the communication interface is used to realize information transmission between server-side devices, metering devices, and user-side devices.

[0221] The processor is used to call computer programs stored in memory, and when the processor executes the computer programs, it implements the integrated energy system in the above embodiments. All steps in the stream allocation parsing tracing method.

[0222] Preferably, embodiments of this application also provide an integrated energy system capable of implementing the above embodiments. A computer-readable storage medium containing all steps of the flow allocation parsing and tracing method, wherein a computer program is stored on the computer-readable storage medium, and when executed by a processor, the computer program implements all steps of the integrated energy system flow allocation parsing and tracing method in the above embodiments.

[0223] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. In particular, hardware + program embodiments are relatively simple in description because they are fundamentally similar to method embodiments; relevant parts can be referred to the descriptions in the method embodiments.

[0224] The foregoing has described specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims may be performed in a different order than that shown in the embodiments and may still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require the specific or sequential order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are possible or may be advantageous.

[0225] While this application provides method operation steps as shown in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only execution order. In actual device or client product execution, the method can be executed sequentially as shown in the embodiments or drawings, or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0226] 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.

[0227] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment 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.

[0228] Unless otherwise specified, the model numbers of the various devices in this embodiment of the invention are not limited, and any device that can perform the above functions is acceptable.

[0229] Those skilled in the art will understand that the accompanying drawings are merely schematic diagrams of a preferred embodiment, and the sequence numbers of the above embodiments of the present invention are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0230] This invention is not limited to the embodiments described above. The above description of specific embodiments is intended to illustrate and explain the technical solutions of this invention. The specific embodiments described above are merely illustrative and not restrictive. Without departing from the spirit and scope of the claims, those skilled in the art can make many specific modifications based on the teachings of this invention, and these modifications all fall within the scope of protection of this invention.

Claims

1. A method for analyzing and tracing the distribution of energy in a comprehensive energy system, characterized in that, include: S1. Determine the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system; S2. Obtain load data of unbalanced nodes through the integrated energy system metering system, and establish an equivalent flow mechanism model of the integrated energy system; (201) The integrated energy system metering system is a device installed at each key node in the integrated energy system, which can directly measure the energy parameters of each key node. The energy parameters include, but are not limited to, energy potential, energy current, load energy, and source energy. (202) The equivalent slack flow mechanism model of the integrated energy system includes the equivalent slack flow mechanism model of the power system, the equivalent slack flow mechanism model of the natural gas system, the equivalent slack flow mechanism model of the thermal system, and the energy station model: S3. Calculate the current distribution of the integrated energy system according to the direct current calculation method, and establish the current circuit of the integrated energy system; the current circuit of the integrated energy system is defined as a current loop composed of branches, current impedance, current source, and current load components; S4. Perform lossless processing on the integrated energy system network to ensure that the current at the beginning and end of the line is consistent; S5. Construct the relationship between the total current vector of all nodes in the integrated energy system, the injected current vector of all nodes, and the source current vector of all nodes, and form a node-load current distribution matrix. ; S6. Derive the current distribution relationships between the source end and load, and between the source end and branch in the integrated energy system; Specifically as follows: The source-end flow column vector of the integrated energy system can be written as: (26) The relationship between the source flow column vector and the node flow column vector of the integrated energy system is as follows: (27) The node-load flow allocation matrix defined in step S5 establishes a relationship between node flow and load flow. , Let the load-end node flow vector be , then the relationship between the source node flow and the load flow is: (28) Therefore, the correlation matrix of the source-load flow of the integrated energy system is expressed as: (29) If the distribution coefficients of the branch flow and the outflow nodes are the same, then the source-branch correlation matrix is: (30) In the formula: The node-branch association matrix of the injected node; Based on the correlation matrix of source-load current Source-branch correlation matrix It can calculate the current distribution ratio between the source end and the load, and between the source end and the branch in the integrated energy system without any allocation principle.

2. The integrated energy system flow distribution analysis and tracking method according to claim 1, characterized in that, In step S1, the network topology information of the integrated energy system includes the topology of the power system, natural gas system, hydrogen energy system, heating system, and cooling system; the multi-energy pipeline parameters of the integrated energy system include the line resistance, reactance, length, and model of the power system; the pipeline length, diameter, and model of the natural gas system; and the length, diameter, thermal conductivity, and roughness of the heating and cooling pipelines; the multi-energy load locations of the integrated energy system include the positions of the power load, natural gas load, hydrogen energy load, heating load, and cooling load in the topology; the multi-energy terminal locations of the integrated energy system include the positions of the power source, natural gas gate station, hydrogen energy gate station, heating source, and cooling source in the topology; the energy station equipment parameters and operating modes of the integrated energy system include the configuration, model, and output ratio of the energy conversion units inside the energy station; the power source includes power plants and distributed new energy sources.

3. The integrated energy system flow distribution analysis and tracking method according to claim 1, characterized in that, In step (202), In the equivalent backflow mechanism model of the power system, the backflow of power lines is expressed as the node backflow and the flow rate of the medium in the power system branches. The product form: (1) In the formula: The voltage at the power node is equal to the node voltage. ; Power line losses are expressed as voltage difference and the flow rate of the medium in the power system branches. The product form: (2) In the formula: and These are the voltages at both ends of the line, respectively, and are equal to the voltages at both ends of the line. and ; Power node injection current Expressed as the product of node voltage and node injected current: (3) In the formula: The conjugate of the injected current to the node is chosen, with the positive direction being the injection of current from the power system into the node. At this point, the power source, electrical load, and the corresponding connection node... These represent negative values, positive values, and 0, respectively. Dimensional power line flow vector , Power line loss column vector , Dimensional power node injection flow column vector They are represented as follows: (4) (5) (6) In the formula: and These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the power system, both of order [number missing]. ; and They are respectively dimensional power node pressure column vector sum Conjugate column vector of injected current at a power node; for A conjugate column vector of current in a power line; In the equivalent backflow mechanism model of the natural gas system, the backflow of natural gas pipelines is represented by the node backflow pressure and the flow rate of the natural gas system medium. The product form: (7) In the formula: The nodal pressure of natural gas is equal to the nodal potential. ; The loss in a natural gas pipeline is expressed as the pressure difference and the flow rate of the natural gas system medium. The product form: (8) In the formula: and The pressure at each of the two ends of the pipeline is equal to the pressure at each of the two ends of the pipeline. and ; Natural gas node injection flow Expressed as the product of nodal pressure and nodal injection gas rate: (9) In the formula: To determine the gas flow rate for the node, the positive direction is chosen as the gas flow from the natural gas system into the node. At this point, the gas source, gas load, and the corresponding connection node... These represent negative values, positive values, and 0, respectively. Dimensional natural gas pipeline flow vector , Natural gas pipeline loss column vector , Dimensional natural gas node injection flow column vector They are represented as follows: (10) (11) (12) In the formula: and These are the outflow node-branch correlation matrix and the node-branch correlation matrix of the natural gas system, both of order [number missing]. ; and They are respectively Dimensional natural gas node pressure column vector sum Column vector of natural gas node injection flow rate; In the equivalent swirling mechanism model of the thermal system, the swirling flow in the thermal branch is expressed as the node swirling pressure and the corresponding flow rate of the medium in the thermal system branch. The product of: (13) The thermal branch loss is expressed as the pressure difference between the nodes at both ends of the branch and the corresponding flow rate of the medium in the thermal system branch. The product of: (14) In the formula: and These are the node pressures at both ends of the thermal branch; Thermal node injection flow Represented as nodal pressure and nodal injection water mass flow rate The product of: (15) In the formula: The mass flow rate of water injected into the node is determined when the node is a heat load node. For heat load water flow rate When the node is a heat source node, For heat source water flow rate When the node is a connected node, =0; Thermal branch flow column vector , Thermal branch loss column vector , Thermal node injection flow column vector They are represented as follows: (16) (17) (18) In the formula: and These are the outflow node-branch correlation matrix and the water supply network node-branch correlation matrix, respectively. for Dimensional thermal node pressure column vector; This is a column vector of water flow rates at thermal nodes; This corresponds to the flow rate vector of the pipeline. The energy station model uses node losses to represent the energy station current model. This portion of the loss is equal to the electrical energy consumed or supplied by the energy station in the power system. Fuel consumed or supplied in the natural gas system The heat consumed or supplied in the thermal system The sum, i.e., the sum of the nodal injection currents of the equivalent electricity, natural gas, and heat nodes of the energy station: (19)。 4. The integrated energy system flow distribution analysis and tracing method according to claim 1, characterized in that, In step S4, the integrated energy system network lossless processing distributes the losses to the nodes at both ends of the pipeline, as follows: (20) In the formula: This is the lossless node flow vector; The node flow vector before lossless transformation; The node-branch correlation matrix of the integrated energy system; This is the pipeline loss vector; The lossless energy flow vectors of the source-end nodes and load-end nodes of the integrated energy system are expressed as follows: (21) In the formula: The source node's flow vector; The flow vector of the load-end node; The lossless processing of the integrated energy system pipeline flow vector is expressed as: (22) In the formula, The pipeline flow vector after lossless transformation; This is the pipeline flow vector before lossless transformation.

5. The integrated energy system flow distribution analysis and tracking method according to claim 1, characterized in that, Step S5 is as follows: To ensure the balance of outflow from nodes, the vector relationship of the total outflow from nodes can be expressed as: (23) In the formula: This is the lossless node flow vector; The pipeline flow vector after lossless transformation; The node-branch association matrix for outflow nodes; This is the flow vector at the load end node, and can be zero. Expanding equation (23), we get: (24) Right now, (25) In the formula, Assign a node-load flow distribution matrix.

6. A comprehensive energy system flow distribution analysis and tracking device, characterized in that, include: The data storage unit is used to store the network topology information, multi-energy pipeline parameters, multi-energy load locations, multi-energy terminal locations, energy station equipment parameters, and operating modes of the target integrated energy system. The integrated energy system equivalent current mechanism model generation unit is used to generate an integrated energy system equivalent current mechanism model based on the load data of unbalanced nodes obtained from the integrated energy system metering system. The integrated energy system equivalent current mechanism model includes an equivalent current mechanism model for a power system, an equivalent current mechanism model for a natural gas system, an equivalent current mechanism model for a heating system, and an energy station model. The calculation unit is used to calculate the current distribution of the integrated energy system according to the direct current calculation method, and to establish the current path of the integrated energy system; The processing unit is used to perform lossless processing of the integrated energy system circuit to ensure that the current at the beginning and end of the line is consistent; The solving unit is used to form a node-load flow distribution matrix based on the relationship between the total flow vector of the node, the flow vector flowing out of the node, and the load flow vector on the node. Ultimately, the current distribution relationships between the source and load, and between the source and branch of the integrated energy system are derived. Specifically as follows: The source-end flow column vector of the integrated energy system can be written as: (26) The relationship between the source flow column vector and the node flow column vector of the integrated energy system is as follows: (27) The node-load flow allocation matrix defined in step S5 establishes a relationship between node flow and load flow. , Let the load-end node flow vector be , then the relationship between the source node flow and the load flow is: (28) Therefore, the correlation matrix of the source-load flow of the integrated energy system is expressed as: (29) If the distribution coefficients of the branch flow and the outflow nodes are the same, then the source-branch correlation matrix is: (30) In the formula: The node-branch correlation matrix for the injected node; based on the source-load flow correlation matrix. Source-branch correlation matrix It can calculate the current distribution ratio between the source end and the load, and between the source end and the branch in the integrated energy system without any allocation principle.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps of the integrated energy system flow allocation analysis and tracing method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the integrated energy system flow allocation analysis and tracing method according to any one of claims 1 to 5.

Images