CN114004047B - Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system - Google Patents
Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system Download PDFInfo
- Publication number
- CN114004047B CN114004047B CN202111148312.3A CN202111148312A CN114004047B CN 114004047 B CN114004047 B CN 114004047B CN 202111148312 A CN202111148312 A CN 202111148312A CN 114004047 B CN114004047 B CN 114004047B
- Authority
- CN
- China
- Prior art keywords
- natural gas
- node
- time
- operation model
- branch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005485 electric heating Methods 0.000 title claims abstract description 24
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 320
- 239000003345 natural gas Substances 0.000 claims abstract description 160
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 239000007789 gas Substances 0.000 claims abstract description 43
- 238000009826 distribution Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 74
- 230000006835 compression Effects 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 14
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 230000017105 transposition Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000006698 induction Effects 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract description 3
- 230000010354 integration Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Fluid Mechanics (AREA)
- Operations Research (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a matrix operation model modeling method of an electric heating gas interconnection multi-energy system, which is used for acquiring data information of the electric heating gas interconnection multi-energy system; establishing a matrixing operation model of the power system, wherein the matrixing operation model comprises a direct current power flow matrixing operation model and an alternating current power flow matrixing operation model; establishing a matrixing operation model of a thermodynamic system, wherein the matrixing operation model comprises a hydraulic model matrixing operation model, a node pressure and pressure drop matrixing operation model and a flow and temperature matrixing operation model; establishing a matrixing operation model of the natural gas system, wherein the matrixing operation model comprises a flow balance matrixing operation model and a pressure distribution matrixing operation model; and acquiring data information output by the matrixing operation models of all the electric power systems, the thermodynamic systems and the natural gas systems. The model constructed by the invention has strong universality and wide application range, provides theoretical guidance and reference for rapid topology structure analysis, rapid data arrangement and induction of an electric heating gas interconnection multi-energy system, and saves the on-site time of engineering service and manpower and material resources.
Description
Technical Field
The invention relates to a matrix operation model modeling method of an electrothermal gas interconnection multi-energy system, and belongs to the technical field of electrothermal gas interconnection multi-energy systems.
Background
In recent years, the development of energy systems presents diversified, intelligent and informationized trends, and the combined coordinated operation of multiple or multiple energy conversion coupling devices provides rich possibility for the operation management, control, adjustment and the like of an electric-thermal gas interconnection multi-energy system. The energy utilization is developed towards the directions of multi-energy coordination, multi-energy complementary utilization, deep cascade utilization of energy coordination and multi-source coupling integration. In order to further improve the comprehensive energy efficiency level, cultivate the new mode amateur state of energy supply and marketing, reduce the operation and maintenance economic cost and realize the win-win situation of multiple aspects, strengthen the coordinated interaction and integration of multiple energy sources, networks, charges and storage depths, construct an interconnected multi-energy system which is one of important technical means for realizing multiple aspects of interconnection and intercommunication, complementary interaction of energy flows and open interconnection of multiple types of energy sources of a novel energy system in the future, and develop the basic premise of constructing the interconnected multi-energy system that the electric heating gas interconnected multi-energy system has strong universality and wide application range.
In the aspects of energy production, transmission, conversion, distribution, storage, utilization, market transaction and the like, the whole interconnected multi-energy system needs to be analyzed and researched by comprehensively considering strategies and methods of coupling interconnection, unified integration, collaborative scheduling and optimal management, the coupling interconnection and integration among multiple energy sources are enhanced, and the coordination complementation and collaborative optimization of multiple energy sources are promoted to become the necessary trend of the development of a novel energy system in the future. It should be noted that the multi-heterogeneous energy interconnection system with the multi-energy coordination and complementation characteristic becomes an effective scheme for solving the problems of high operation economic cost, pollution discharge amplification, low comprehensive energy efficiency level and the like of a single energy system. Particularly, the electric, thermal and gas interconnection energy system coordinates and complements a plurality of energy sources, is beneficial to realizing economic and low-carbon targets, is beneficial to realizing bidirectional complementary interaction and collaborative optimization of a plurality of energy sources, and is one of effective means for improving the flexibility of system operation scheduling and resource coordination configuration.
However, the type of the energy conversion device in the electric heating gas interconnection multi-energy system is complex and various, the coupling interconnection mode of the heterogeneous energy subsystem is complex, and the like, so that the general modeling is complex. Therefore, the universal matrixing modeling method for the broadcast and television hot gas interconnection multi-energy system with strong universality and application range is very necessary and urgent in the energy subsystem level.
Disclosure of Invention
The purpose is as follows: in order to overcome the defects in the prior art, the invention provides a modeling method for a matrixing operation model of an electric heating gas interconnection multi-energy system, which aims to realize rapid topological structure analysis, rapid data arrangement and induction, establish a matrixing operation model with strong universality and wide application range on a system level, further achieve the purposes of saving engineering service field time and manpower and material resources and providing theoretical guidance and reference for modeling and operation of the electric heating gas interconnection multi-energy system.
The technical scheme is as follows: in order to solve the technical problems, the invention adopts the following technical scheme:
a matrix operation model modeling method of an electric heating gas interconnection multi-energy system comprises the following steps:
(1) The method comprises the steps of obtaining data information of an electric heating gas interconnection multi-energy system, wherein the data information comprises data information of coupling interconnection architecture of the electric heating gas interconnection multi-energy system, energy conversion coupling equipment, upper and lower output power limits of a power supply, a heat source and a natural gas source, a power system topological structure, a thermodynamic system topological structure, a natural gas system topological structure, power system line parameters, thermodynamic system pipeline branch parameters and natural gas system pipeline branch parameters.
(2) And establishing a matrixing operation model of the power system, wherein the matrixing operation model comprises a direct current load flow matrixing operation model and an alternating current load flow matrixing operation model.
(3) And establishing a matrixing operation model of the thermodynamic system, wherein the matrixing operation model comprises a hydraulic model matrixing operation model, a node pressure and pressure drop matrixing operation model and a flow and temperature matrixing operation model.
(4) And establishing a matrixing operation model of the natural gas system, wherein the matrixing operation model comprises a flow balance matrixing operation model and a pressure distribution matrixing operation model.
(5) Acquiring data information output by matrix operation models of all power systems, thermodynamic systems and natural gas systems, wherein the data information comprises power supply output, node electric power, node reactive power, node voltage phase angle and node voltage data information of the power systems; the heat source output, the node mass flow, the node thermal power, the node pressure and the pipeline pressure loss data information of the thermodynamic system; the gas source output, node pressure, pipeline natural gas flow, node natural gas flow, compressor consumption flow and compressor compression ratio data information of the natural gas system.
Further, the step (2) of establishing a matrixing operation model of the power system includes:
the DC power flow matrixing operation model has the following specific expression:
Wherein: to inject a power vector at time t node, where element/> Here, theAnd/>Generator output and load at time t node i,/>, respectivelyAs susceptance matrix, element B ij=-1/xij,/>, of the susceptance matrixB ij、xij is susceptance and reactance of a power network line or branch ij respectively, ij is a line or branch, and the numbers of the initial end nodes and the final end nodes are i and j respectively; /(I)Is the voltage phase angle vector at the node at time t; subscripts i and j are the numbers of different nodes of the power network.
The alternating current power flow matrixing operation model has the following specific expression:
Wherein: p t、Qt is the active power vector and the reactive power vector of the node in the power system at time t respectively; u t is the node voltage vector in the power system at time tV; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; * And (5) performing conjugate operation on the orientation quantity.
Further, the node number of the thermodynamic system is from 1 to node based on the graph theory idea, namely the total number of nodes of the thermodynamic system is node; the pipes or branches of the thermodynamic system are numbered from 1 to b, i.e. the total number of branches of the thermodynamic system is b.
Further, the step (3) of establishing a matrixing operation model of the thermodynamic system includes:
the hydraulic model matrixing operation model has the following specific expression:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) The mass flow vector of the thermal working medium of the heating power pipeline; 0 is zero matrix vector; /(I)A basic circuit matrix associated with a closed circuit of the thermodynamic system; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working fluid head vector raised for the pressure circulation pump at time t.
The node pressure and pressure drop matrixing operation model has the following specific expression:
Wherein:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) Is a node pressure vector in the thermodynamic system at time t; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working medium pressure head vector raised by the pressure circulation pump at time t; the superscript b is the total number of branches of the thermodynamic system; superscript' is vector transposition operation; Is the pressure at node in the thermodynamic system at time tdata; /(I) Is the pressure loss of the thermodynamic pipe b in the thermodynamic system at time t; /(I)The hot working medium pressure head is raised by the pressure circulating pump in the thermodynamic system at time t.
The flow and temperature matrixing operation model has the following specific expression:
Wherein: The input and output thermal power vectors of the thermal pipeline in the thermodynamic system at time t are respectively; /(I) Is a node thermal load vector in the thermodynamic system at time t; c p is the specific heat capacity of the hot working medium; /(I)The method comprises the steps of heating a working medium mass flow vector for a pipeline terminal node of a thermodynamic system at time t; /(I)The temperature vectors of the inlet port and the outlet port at the junction of the thermodynamic pipeline and the thermodynamic pipeline at time t are respectively; /(I)Respectively a node-mass flow inflow pipeline starting point correlation matrix and a node-mass flow outflow pipeline ending point correlation matrix in the thermodynamic system at time t.
Further, the natural gas system comprises a medium-pressure and high-pressure transmission natural gas system, the air flow in the natural gas pipeline is closely related to the pressure intensity of nodes at two sides of the pipeline and the physical condition of pipeline transmission, and the mathematical relationship is as follows:
wherein: subscripts i, j and ij are respectively the head end node, the tail end node and the pipeline number of the natural gas pipeline; Natural gas pressure at natural gas system nodes i, j, respectively, at time t; /(I) Is a natural gas system natural gas flow direction variable at time t; /(I)Is a characteristic constant of the pipeline; /(I)Is the air flow in the duct at time t; k ij is equivalent characteristic physical parameter of the pipeline branch, and is also called natural gas pipeline branch impedance; /(I)The equivalent pressure square difference of the head and the tail of the branch is obtained.
Based on the theory of the graph, the nodes in the topological relation of the medium and high pressure natural gas system are numbered, the number of the nodes is from 1 to node, namely the total number of the nodes is node; the natural gas pipeline or the topological branch is numbered, and the branch number is from 1 to b, namely the total number of the branches is b. Then, the compression ratio of the compressors of each branch may be respectively corresponding to the branch numbersTo/>The flow direction variable of the natural gas system pipeline can respectively correspond to the serial number of the branchTo/>The equivalent characteristic parameters of the branches can respectively correspond to the branch numbers K 1 to K b; the equivalent pressure squared differences at the head and the tail of the branch can respectively correspond to the serial numbers of the branch as/>To/>The injected natural gas flow of each node can respectively correspond to the node number of/>To/>The natural gas flow of each branch can be respectively corresponding to the branch number of/>To/>The natural gas flow consumed by the compressors on each branch can be respectively corresponding to the branch number/>To the point of
Further, the step (4) of establishing a matrixing operation model of the natural gas system includes: the flow balance matrixing operation model and the pressure distribution matrixing operation model comprise:
the flow balance matrixing operation model has the following specific expression:
Wherein:
Wherein: The node-branch complete incidence matrix of the natural gas system; f t is the natural gas flow vector through the pipeline in the natural gas network at time t; /(I) The method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; /(I)Is the natural gas flow vector input into the bypass compressor in the natural gas network at time t; q t is the node injection flow vector in the natural gas system at time t; the superscript b is the total number of the branches of the natural gas pipeline; the subscript node is the total number of topological nodes of the natural gas system; superscript' is vector transposition operation; /(I)Is the natural gas flow of branch b in the natural gas system at time t; /(I)Is the natural gas flow consumed by the compressor at time t branch b; /(I)For the injection of natural gas flow at node t.
The pressure distribution matrixing operation model has the following specific expression:
Wherein:
Wherein: the method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; is a diagonal matrix in the natural gas network related to the compression ratio of the bypass compressors at time t; /(I) The node-branch complete endpoint correlation matrix of the natural gas system; /(I)Square vector of node pressure in natural gas system at time t; pi t is the equivalent pressure square difference vector of the first end and the second end of the branch in the natural gas system at time t; superscript' is vector transposition operation; /(I)Is the compressor compression ratio of branch b at time t; /(I)Is the square of the pressure at node in the natural gas system at time t; /(I)Is the natural gas flow direction of branch b in the natural gas system at time t; k b is the physical parameter of the pipeline branch b; /(I)Is the natural gas flow through line b at time t.
Further, the natural gas system further comprises:
① When no compressor exists in the ith branch, the compression ratio is set at the moment And set the natural gas flow consumed by the compressor/>In particular, for electric compressors, always set/>
② The gas flow direction of a natural gas system is usually determined, and particularly in a high-pressure natural gas transmission system, the gas flow direction of various regulating valves and practical engineering applications cannot be easily changed, so that the natural gas topological relation is usually determined, namely, the flow direction variable of a pipeline of the natural gas system is usually a known quantity.
③ When the ith node is an air source nodeThe value is positive, which is the time of the gas load node/>The value is negative and is the intermediate node/>The value is 0.
④ The matrix equation of the medium-high pressure transmission natural gas system is also applicable to the low-pressure gas distribution system, and at the moment
The beneficial effects are that: compared with the prior art, the modeling method for the matrix operation model of the electric heating gas interconnection multi-energy system provided by the invention comprehensively considers the electric heating gas heterogeneous energy subsystem at the same time, and can provide theoretical guidance for modeling and operation of the multi-energy system; the invention provides a matrixing operation model of an electric power system, a thermodynamic system and a natural gas system, which has strong universality and wide application range; the invention can effectively provide theoretical guidance and reference for rapid topological structure analysis, rapid data arrangement and induction of the electric heating gas interconnection multi-energy system on the system level, and save the engineering service field time and manpower and material resources.
Drawings
FIG. 1 is a flow chart of a method implementation of the present invention.
Detailed Description
The invention will be further described with reference to specific examples.
Example 1:
a matrix operation model modeling method of an electric heating gas interconnection multi-energy system, as shown in figure 1, comprises the following steps:
(1) Acquiring data information of electric-thermal interconnection multi-energy system
The method comprises the steps of obtaining data information of an electric heating gas interconnection multi-energy system, wherein the data information comprises data such as coupling interconnection architecture of the electric heating gas interconnection multi-energy system, energy conversion coupling equipment, upper and lower output power limits of a power supply, a heat source and a natural gas source, a power system topological structure, a thermodynamic system topological structure, a natural gas system topological structure, power system line parameters, thermodynamic system pipeline branch parameters, natural gas system pipeline branch parameters and the like.
(2) Establishing a matrixing operation model of a power system
The method comprises the steps of establishing a direct current power flow equation for a power system, wherein the specific expression of a matrix form is as follows:
Wherein: to inject a power vector at time t node, where element/> Here/>And/>Respectively generating output and load of a generator at a time t node i, wherein subscripts i and j are the numbers of different nodes of the power network; Is susceptance matrix, wherein the element B ij=-1/xij,/> B. x is susceptance and reactance of a power network line or branch respectively, ij is the line or branch, and the numbers of the head and tail end nodes are i and j respectively; /(I)Is the voltage phase angle vector at time t.
The mathematical model described by the classical alternating current power flow model in the power system in a matrix form is specifically as follows:
Wherein: p t、Qt is the active power vector and the reactive power vector of the node in the power system at time t respectively; u t is the node voltage vector in the power system at time tV; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; * And (5) performing conjugate operation on the orientation quantity.
(3) Establishing a matrixing operation model of a thermodynamic system
Aiming at a thermodynamic system, based on the graph theory, the node numbers of the thermodynamic system are from 1 to node, namely the total number of nodes of the thermodynamic system is node; the pipes or branches of the thermodynamic system are numbered from 1 to b, i.e. the total number of branches of the thermodynamic system is b.
The specific expression of the hydraulic model matrixing operation of the thermodynamic system is as follows:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) The mass flow vector of the thermal working medium of the heating power pipeline; 0 is zero matrix vector; /(I)A basic circuit matrix associated with a closed circuit of the thermodynamic system; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working fluid head vector raised for the pressure circulation pump at time t.
When the node pressure and pressure drop in the thermodynamic system are subjected to matrix operation, the specific expression is as follows:
Wherein:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) Is a node pressure vector in the thermodynamic system at time t; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working medium pressure head vector raised by the pressure circulation pump at time t; the superscript b is the total number of branches of the thermodynamic system; /(I)Is the pressure at node in the thermodynamic system at time tdata; /(I)Is the pressure loss of the thermodynamic pipe b in the thermodynamic system at time t; /(I)The hot working medium pressure head is raised by the pressure circulating pump in the thermodynamic system at time t.
The matrixing operation model for expressing the flow and the temperature of the thermodynamic system by using the matrix equation set is specifically as follows:
Wherein: The input and output thermal power vectors of the thermal pipeline in the thermodynamic system at time t are respectively; /(I) Is a node thermal load vector in the thermodynamic system at time t; c p is the specific heat capacity of the hot working medium; /(I)The method comprises the steps of heating a working medium mass flow vector for a pipeline terminal node of a thermodynamic system at time t; /(I)The temperature vectors of the inlet port and the outlet port at the junction of the thermodynamic pipeline and the thermodynamic pipeline at time t are respectively; /(I)Respectively a node-mass flow inflow pipeline starting point correlation matrix and a node-mass flow outflow pipeline ending point correlation matrix in the thermodynamic system at time t.
(4) Establishing a matrixing operation model of a natural gas system
For a medium-high pressure natural gas transmission system, the air flow in a natural gas pipeline is closely related to the pressure intensity of nodes at two sides of the pipeline and the physical condition of pipeline transmission, and the mathematical relationship is as follows:
wherein: subscripts i, j and ij are respectively the head end node, the tail end node and the pipeline number of the natural gas pipeline; Natural gas pressure at natural gas system nodes i, j, respectively, at time t; /(I) Is a natural gas system natural gas flow direction variable at time t; /(I)Is a characteristic constant of the pipeline; /(I)Is the air flow in the duct at time t; k ij is equivalent characteristic physical parameter of the pipeline branch, and is also called natural gas pipeline branch impedance; /(I)The equivalent pressure square difference of the head and the tail of the branch is obtained.
Based on the theory of the graph, analyzing the topological relation of the medium-high pressure transmission natural gas system, numbering the topological nodes of the natural gas system, wherein the number of the nodes is from 1 to node, namely the total number of the nodes is node; the natural gas pipeline or the topological branch is numbered, and the branch number is from 1 to b, namely the total number of the branches is b. Then, the compression ratio of the compressors of each branch may be respectively corresponding to the branch numbersTo/>The flow direction variable of the natural gas system pipeline can respectively correspond to the serial number of the branchTo/>The equivalent characteristic parameters of the branches can respectively correspond to the branch numbers K 1 to K b; the equivalent pressure squared differences at the head and the tail of the branch can respectively correspond to the serial numbers of the branch as/>To/>The injected natural gas flow of each node can respectively correspond to the node number of/>To/>The natural gas flow of each branch can be respectively corresponding to the branch number of/>To/>The natural gas flow consumed by the compressors on each branch can be respectively corresponding to the branch number/>To/>
The natural gas system flow balance matrixing equation expression is as follows:
Wherein:
Wherein: The node-branch complete incidence matrix of the natural gas system; f t is the natural gas flow vector through the pipeline in the natural gas network at time t; /(I) The method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; /(I)Is the natural gas flow vector input into the bypass compressor in the natural gas network at time t; q t is the node injection flow vector in the natural gas system at time t; the superscript b is the total number of the branches of the natural gas pipeline; the subscript node is the total number of topological nodes of the natural gas system; superscript' is vector transposition operation; /(I)Is the natural gas flow of branch b in the natural gas system at time t; /(I)Is the natural gas flow consumed by the compressor at time t branch b; /(I)For the injection of natural gas flow at node t.
The pressure distribution matrixing equation expression of the natural gas system is as follows:
Wherein:
Wherein: the method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; is a diagonal matrix in the natural gas network related to the compression ratio of the bypass compressors at time t; /(I) The node-branch complete endpoint correlation matrix of the natural gas system; /(I)Square vector of node pressure in natural gas system at time t; pi t is the equivalent pressure square difference vector of the first end and the second end of the branch in the natural gas system at time t; superscript' is vector transposition operation; /(I)Is the compressor compression ratio of branch b at time t; /(I)Is the square of the pressure at node in the natural gas system at time t; /(I)Is the natural gas flow direction of branch b in the natural gas system at time t; k b is the physical parameter of the pipeline branch b; /(I)Is the natural gas flow through line b at time t.
The following description is given to a matrixing model of a natural gas system:
① When no compressor exists in the ith branch, the compression ratio is set at the moment And set the natural gas flow consumed by the compressor/>In particular, for electric compressors, always set/>
② The gas flow direction of a natural gas system is usually determined, and particularly in a high-pressure natural gas transmission system, the gas flow direction of various regulating valves and practical engineering applications cannot be easily changed, so that the natural gas topological relation is usually determined, namely, the flow direction variable of a pipeline of the natural gas system is usually a known quantity.
③ When the ith node is an air source nodeThe value is positive, which is the time of the gas load node/>The value is negative and is the intermediate node/>The value is 0.
④ The matrix equation of the medium-high pressure transmission natural gas system is also applicable to the low-pressure gas distribution system, and at the moment
(5) Outputting matrixed model data information
Outputting matrix model data information, including data information such as power supply output, node electric power, node reactive power, node voltage phase angle, node voltage and the like of the power system; the heat source output, the node mass flow, the node thermal power, the node pressure intensity, the pipeline pressure intensity loss and other data information of the thermodynamic system; and the data information such as gas source output, node pressure, pipeline natural gas flow, node natural gas flow, compressor consumption flow, compressor compression ratio and the like of the natural gas system.
The foregoing is only a preferred embodiment of the invention, it being noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present invention, and such modifications and adaptations are intended to be comprehended within the scope of the invention.
Claims (4)
1. A matrix operation model modeling method of an electric heating gas interconnection multi-energy system is characterized by comprising the following steps of: the method comprises the following steps:
acquiring data information of an electric-thermal interconnection multi-energy system;
according to the data information, a matrixing operation model of the power system is established, wherein the matrixing operation model comprises a direct current power flow matrixing operation model and an alternating current power flow matrixing operation model;
According to the data information, a matrixing operation model of the thermodynamic system is established, wherein the matrixing operation model comprises a hydraulic model matrixing operation model, a node pressure and pressure drop matrixing operation model and a flow and temperature matrixing operation model;
According to the data information, a matrixing operation model of the natural gas system is established, wherein the matrixing operation model comprises a flow balance matrixing operation model and a pressure distribution matrixing operation model;
the DC power flow matrixing operation model has the following specific expression:
Wherein: to inject a power vector at time t node, where element/> Here/>AndGenerator output and load at time t node i,/>, respectivelyAs susceptance matrix, element B ij=-1/xij,/>, of the susceptance matrixB ij、xij is susceptance and reactance of a power network line or branch ij respectively, ij is a line or branch, and the numbers of the initial end nodes and the final end nodes are i and j respectively; /(I)Is the voltage phase angle vector at the node at time t; subscripts i and j are the numbers of different nodes of the power network;
the alternating current power flow matrixing operation model has the following specific expression:
Wherein: p t、Qt is the active power vector and the reactive power vector of the node in the power system at time t respectively; u t is the node voltage vector in the power system at time tV; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; * Performing conjugate operation on the orientation quantity;
The node numbers of the thermodynamic system are from 1 to node based on the graph theory idea, namely the total number of the nodes of the thermodynamic system is node; the number of the pipelines or the branches of the thermodynamic system is from 1 to b, namely the total number of the branches of the thermodynamic system is b;
the hydraulic model matrixing operation model has the following specific expression:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) The mass flow vector of the thermal working medium of the heating power pipeline; 0 is zero matrix vector; /(I)A basic circuit matrix associated with a closed circuit of the thermodynamic system; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working medium pressure head vector raised by the pressure circulation pump at time t;
The node pressure and pressure drop matrixing operation model has the following specific expression:
Wherein:
Wherein: is a complete correlation matrix associated with thermodynamic system nodes and thermodynamic pipes; /(I) Is a node pressure vector in the thermodynamic system at time t; /(I)Is a thermodynamic pipeline pressure loss vector at time t; /(I)A hot working medium pressure head vector raised by the pressure circulation pump at time t; the superscript b is the total number of branches of the thermodynamic system; superscript' is vector transposition operation; /(I)Is the pressure at node in the thermodynamic system at time tdata; /(I)Is the pressure loss of the thermodynamic pipe b in the thermodynamic system at time t; /(I)A hot working medium pressure head elevated by a pressure circulation pump in the thermodynamic system at time t;
the flow and temperature matrixing operation model has the following specific expression:
Wherein: The input and output thermal power vectors of the thermal pipeline in the thermodynamic system at time t are respectively; /(I) Is a node thermal load vector in the thermodynamic system at time t; c p is the specific heat capacity of the hot working medium; /(I)The method comprises the steps of heating a working medium mass flow vector for a pipeline terminal node of a thermodynamic system at time t; /(I)The temperature vectors of the inlet port and the outlet port at the junction of the thermodynamic pipeline and the thermodynamic pipeline at time t are respectively; /(I)Respectively a node-mass flow inflow pipeline starting point correlation matrix and a node-mass flow outflow pipeline ending point correlation matrix in the thermodynamic system at time t;
The natural gas system comprises a medium-pressure and high-pressure natural gas transmission system, the air flow in the natural gas pipeline is closely related to the pressure intensity of nodes at two sides of the pipeline and the physical condition of pipeline transmission, and the mathematical relationship is as follows:
wherein: subscripts i, j and ij are respectively the head end node, the tail end node and the pipeline number of the natural gas pipeline; Natural gas pressure at natural gas system nodes i, j, respectively, at time t; /(I) Is a natural gas system natural gas flow direction variable at time t; kappa ij is the pipeline characteristic constant; /(I)Is the air flow in the duct at time t; k ij is equivalent characteristic physical parameter of the pipeline branch, and is also called natural gas pipeline branch impedance; /(I)The equivalent pressure square difference of the head and the tail of the branch is obtained;
Based on the theory of the graph, the nodes in the topological relation of the medium and high pressure natural gas system are numbered, the number of the nodes is from 1 to node, namely the total number of the nodes is node; numbering natural gas pipelines or topological branches, wherein the number of branches is from 1 to b, namely the total number of branches is b; then, the compression ratio of the compressors of each branch may be respectively corresponding to the branch numbers To/>The flow direction variable of the natural gas system pipeline can respectively correspond to the serial number of the branchTo/>The equivalent characteristic parameters of the branches can respectively correspond to the branch numbers K 1 to K b; the equivalent pressure squared differences at the head and the tail of the branch can respectively correspond to the serial numbers of the branch as/>To/>The injected natural gas flow of each node can respectively correspond to the node number of/>To/>The natural gas flow of each branch can respectively correspond to the branch numbers f 1 t to/>The natural gas flow consumed by the compressors on each branch can be respectively corresponding to the branch number/>To/>
The flow balance matrixing operation model has the following specific expression:
Wherein:
Wherein: The node-branch complete incidence matrix of the natural gas system; f t is the natural gas flow vector through the pipeline in the natural gas network at time t; /(I) The method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; Is the natural gas flow vector input into the bypass compressor in the natural gas network at time t; q t is the node injection flow vector in the natural gas system at time t; b is the total number of branches of the natural gas pipeline; the subscript node is the total number of topological nodes of the natural gas system; superscript' is vector transposition operation; /(I) Is the natural gas flow of branch b in the natural gas system at time t; /(I)Is the natural gas flow consumed by the compressor at time t branch b; /(I)For the injection of natural gas flow at node t;
the pressure distribution matrixing operation model has the following specific expression:
Wherein:
Wherein: the method comprises the steps of (1) associating a matrix for a complete starting point in a natural gas network, wherein the complete starting point is related to a network node and a pipeline branch; /(I) Is a diagonal matrix in the natural gas network related to the compression ratio of the bypass compressors at time t; /(I)The node-branch complete endpoint correlation matrix of the natural gas system; /(I)Square vector of node pressure in natural gas system at time t; pi t is the equivalent pressure square difference vector of the first end and the second end of the branch in the natural gas system at time t; superscript' is vector transposition operation; /(I)Is the compressor compression ratio of branch b at time t; /(I)Is the square of the pressure at node in the natural gas system at time t; /(I)Is the natural gas flow direction of branch b in the natural gas system at time t; k b is the physical parameter of the pipeline branch b; /(I)Is the natural gas flow through line b at time t.
2. The modeling method for the matrixing operation model of the electric heating gas interconnection multi-energy system according to claim 1, which is characterized by comprising the following steps: the electric heating gas interconnection multi-energy system data information comprises an electric heating gas interconnection multi-energy system coupling interconnection framework, energy conversion coupling equipment, a power supply, upper and lower limits of output power of a heat source and a natural gas source, an electric power system topological structure, a thermodynamic system topological structure, a natural gas system topological structure, electric power system line parameters, thermodynamic system pipeline branch parameters and natural gas system pipeline branch parameter data information.
3. The modeling method for the matrixing operation model of the electric heating gas interconnection multi-energy system according to claim 1, which is characterized by comprising the following steps: further comprises:
Acquiring data information output by matrix operation models of all power systems, thermodynamic systems and natural gas systems, wherein the data information comprises power supply output, node electric power, node reactive power, node voltage phase angle and node voltage data information of the power systems; the heat source output, the node mass flow, the node thermal power, the node pressure and the pipeline pressure loss data information of the thermodynamic system; the gas source output, node pressure, pipeline natural gas flow, node natural gas flow, compressor consumption flow and compressor compression ratio data information of the natural gas system.
4. The modeling method for the matrixing operation model of the electric heating gas interconnection multi-energy system according to claim 1, which is characterized by comprising the following steps: further comprises:
When no compressor exists in the ith branch, the compression ratio is set at the moment And the natural gas flow consumed by the compressor is setFor the electric compressor, always set/>
The gas flow direction of the natural gas system is determined, and in a high-pressure natural gas transmission system, the gas flow direction of various regulating valves and practical engineering application cannot be easily changed, so that the natural gas topological relation is determined, namely the flow direction variable of a pipeline of the natural gas system is a known quantity;
when the ith node is an air source node The value is positive, which is the time of the gas load node/>The value is negative and is the intermediate node/>The value is 0;
The matrix equation of the medium-high pressure transmission natural gas system is also applicable to the low-pressure gas distribution system, and at the moment
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111148312.3A CN114004047B (en) | 2021-09-29 | 2021-09-29 | Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111148312.3A CN114004047B (en) | 2021-09-29 | 2021-09-29 | Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114004047A CN114004047A (en) | 2022-02-01 |
CN114004047B true CN114004047B (en) | 2024-05-28 |
Family
ID=79921977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111148312.3A Active CN114004047B (en) | 2021-09-29 | 2021-09-29 | Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114004047B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115562029B (en) * | 2022-10-17 | 2023-06-16 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200662A1 (en) * | 2018-04-20 | 2019-10-24 | 东北大学 | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
CN111563315A (en) * | 2020-04-08 | 2020-08-21 | 重庆大学 | Topological analysis-based steady-state energy flow calculation method for electricity-gas comprehensive energy system |
CN112016033A (en) * | 2020-08-07 | 2020-12-01 | 大连理工大学 | Electric-heat-gas comprehensive energy system load flow calculation method based on forward-backward substitution method |
CN112994020A (en) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | Multi-energy system load flow decomposition calculation method |
CN113487171A (en) * | 2021-07-02 | 2021-10-08 | 安徽帅尔信息科技有限公司 | Business resource allocation method based on improved slime mold optimization algorithm |
-
2021
- 2021-09-29 CN CN202111148312.3A patent/CN114004047B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200662A1 (en) * | 2018-04-20 | 2019-10-24 | 东北大学 | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
CN111563315A (en) * | 2020-04-08 | 2020-08-21 | 重庆大学 | Topological analysis-based steady-state energy flow calculation method for electricity-gas comprehensive energy system |
CN112016033A (en) * | 2020-08-07 | 2020-12-01 | 大连理工大学 | Electric-heat-gas comprehensive energy system load flow calculation method based on forward-backward substitution method |
CN112994020A (en) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | Multi-energy system load flow decomposition calculation method |
CN113487171A (en) * | 2021-07-02 | 2021-10-08 | 安徽帅尔信息科技有限公司 | Business resource allocation method based on improved slime mold optimization algorithm |
Non-Patent Citations (4)
Title |
---|
Zhong yongjie,et al..Hierarchical Multi-Objective Fuzzy Collaborative Optimization of Integrated Energy System under Off-Design Performance.Energies.2019,第1-27页. * |
含电—热—气—冷子系统的区域综合能源系统多场景优化调度;钟永洁等;电力系统自动化;20190625;第43卷(第12期);第76-84页 * |
多能源网络的广义电路分析理论――(一)支路模型;杨经纬;张宁;康重庆;;电力系统自动化;20200510(第09期);全文 * |
电力-天然气综合能源系统建模与仿真;常晓勇;王玉婷;陈福锋;李玉平;徐程骥;;浙江电力;20180725(第07期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114004047A (en) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109978625B (en) | Comprehensive energy system multi-target operation optimization method considering electric and thermal gas network | |
CN110210747B (en) | Flexible scheduling method for electric-heating-gas interconnection energy system | |
CN108960503B (en) | Multi-scene optimization analysis method of comprehensive energy system based on interior point method | |
CN106056251B (en) | A kind of Optimization Scheduling of electric-thermal coupling multipotency streaming system | |
CN106096777B (en) | A kind of Optimization Scheduling of electric-gas coupling multipotency streaming system | |
CN111815042A (en) | Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model | |
CN107800158B (en) | A kind of economy of taking into account couples multipotency streaming system Optimization Scheduling with the electric-thermal of efficiency | |
CN109347107A (en) | One kind be incorporated into the power networks electric heating interconnection integrated energy system tidal current computing method | |
CN111899120A (en) | Method for establishing comprehensive energy planning and operation combined optimization model | |
CN107910871A (en) | A kind of energy internet is provided multiple forms of energy to complement each other control method | |
CN110502791B (en) | Comprehensive energy system steady-state modeling method based on energy concentrator | |
CN114004047B (en) | Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system | |
CN112035984A (en) | Electric-gas-storage area comprehensive energy system collaborative planning method | |
CN110968827A (en) | Optimal configuration method for multi-region comprehensive energy system | |
CN108717598A (en) | Index calculating method temporarily drops in the network voltage of electric-gas interconnection integrated energy system | |
CN111062598A (en) | Distributed optimal scheduling method and system for comprehensive energy system | |
CN112330493B (en) | Energy system modeling and comprehensive analysis method, device and storage medium | |
CN109066695A (en) | Two-stage linearization electrical optimal energy flow calculation method | |
CN111310310A (en) | Thermodynamic system static power flow fast decoupling calculation method for quantity adjustment | |
CN114462163B (en) | exergy hub model considering energy quality | |
Li et al. | Exergy Hub: A novel energy hub model considering energy quality | |
CN115545340A (en) | Comprehensive energy station-network collaborative planning method considering double uncertainties | |
Hu et al. | Multi-Energy Flow Calculation Model and Solution Method for Regional Integrated Energy System | |
Ayele et al. | Pseudo-dynamic simulation on a district energy system made of coupling technologies | |
Ma et al. | Optimal Energy Flow Calculation of Electricity-heat Integrated Energy System with Heat Pump |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |