CN114004047A - Electric heating gas interconnection multi-energy system matrix operation model modeling method - Google Patents
Electric heating gas interconnection multi-energy system matrix operation model modeling method Download PDFInfo
- Publication number
- CN114004047A CN114004047A CN202111148312.3A CN202111148312A CN114004047A CN 114004047 A CN114004047 A CN 114004047A CN 202111148312 A CN202111148312 A CN 202111148312A CN 114004047 A CN114004047 A CN 114004047A
- Authority
- CN
- China
- Prior art keywords
- natural gas
- node
- operation model
- time
- pressure
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005485 electric heating Methods 0.000 title claims description 6
- 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
- 239000007789 gas Substances 0.000 claims abstract description 49
- 238000009826 distribution Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 74
- 230000005540 biological transmission Effects 0.000 claims description 15
- 230000006835 compression Effects 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 15
- 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
- 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 3
- 230000000295 complement effect Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000004927 fusion 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
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
Images
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 modeling method of a matrix operation model 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 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; 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 system, the thermodynamic system and the natural gas system. The model constructed by the invention has strong universality and wide application range, provides theoretical guidance and reference for quick topological structure analysis, quick data arrangement and induction of the electric-heating-gas interconnected multi-energy system, and saves the time, manpower and material resources of an engineering service site.
Description
Technical Field
The invention relates to a modeling method of a matrix operation model of an electric-heat-gas interconnection multi-energy system, and belongs to the technical field of electric-heat interconnection multi-energy systems.
Background
In recent years, the development of energy systems shows diversified, intelligent and informatization trends, and the combined coordinated operation of multiple or multiple energy conversion coupling devices provides abundant possibilities for operation management and control, regulation and the like of the electric-thermal-interconnection multi-energy system. The energy utilization develops towards the direction of multi-energy coordination, multi-energy complementary utilization, energy coordination deep cascade utilization and multi-source coupling integration. In order to further improve comprehensive energy efficiency grade, cultivate new energy supply and marketing mode state, reduce operation and maintenance operation economic cost, realize win-win situation, strengthen source, network, load and storage depth coordination interaction and fusion of various energy sources, the construction of the interconnected multi-energy system is one of important technical means for realizing multi-aspect interconnection and intercommunication, energy flow complementation and mutual aid and multi-type energy open interconnection of a novel energy system in the future, and the development of the modeling with strong universality and wide application range of the electric-thermal-electric-interconnection multi-energy system is a basic premise for constructing the interconnected multi-energy system.
In multiple aspects of energy production, transmission, conversion, allocation, storage, utilization, market trading and the like, the whole interconnected multi-energy system needs to be analyzed and researched by comprehensively considering the strategies and methods of coupling interconnection, unified integration, cooperative scheduling and optimized management, coupling interconnection and integration among multiple energy sources are enhanced, the coordination and complementation of multiple energy sources are promoted, and the cooperative optimization becomes the inevitable trend of the development of a novel energy system in the future. It should be noted that the heterogeneous energy interconnection system with the characteristics of multi-energy coordination and complementation becomes an effective scheme for solving the problems of high economic cost, pollution discharge and amplification, low comprehensive energy efficiency level and the like of single energy system operation. Particularly, the electricity, heat and gas interconnected energy system coordinates and complements multiple energy sources, is favorable for realizing economic and low-carbon targets, is favorable for realizing bidirectional complementary interaction and collaborative optimization of the multiple energy sources, and is one of effective means for improving the flexibility of system operation scheduling and resource coordination configuration.
However, the general modeling is complicated due to the fact that the types of the energy conversion equipment in the electric-heating-gas interconnection multi-energy system are complex and various, the coupling interconnection mode of the heterogeneous energy subsystem is complex, and the like. Therefore, it is necessary and urgent to provide a universal matrix modeling method for a multi-energy system with strong universality and wide application range 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 of a matrix operation model of an electric-heat-gas interconnected multi-energy system, which aims to realize quick topological structure analysis, quick data arrangement and induction on a system level and establish the matrix operation model with strong universality and wide application range, thereby achieving the purposes of saving engineering service site time, manpower and material resources and providing theoretical guidance and reference for modeling and operation of the electric-heat-gas interconnected multi-energy system.
The technical scheme is as follows: in order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a modeling method for a matrix operation model 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-heat-gas interconnection multi-energy system, wherein the data information comprises a coupling interconnection framework of the electric-heat-gas interconnection multi-energy system, energy conversion coupling equipment, power supplies, upper and lower limits of output power of 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 circuit parameters, thermodynamic system pipeline branch parameters and natural gas system pipeline branch parameter data information.
(2) And 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.
(3) 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 a matrixing operation model of all electric power systems, thermodynamic systems and natural gas systems, wherein the data information comprises power supply output of the electric power systems, node electric power, node reactive power, node voltage phase angle and node voltage data information; data information of heat source output, node mass flow, node thermal power, node pressure and pipeline pressure loss of the thermodynamic system; the natural gas system comprises the data information of gas source output, node pressure, pipeline natural gas flow, node natural gas flow, compressor consumption flow and compressor compression ratio.
Further, the step (2) of establishing a matrixing operation model of the power system comprises:
the direct current power flow matrixing operation model has the following specific expression:
in the formula:
for injecting a power vector at a time t node, where the elements
Here, the
And
respectively generator output and load at time t node i,
being a susceptance matrix, element B of the susceptance matrixij＝-1/xij、
Bij、xijRespectively the susceptance and reactance of a power network circuit or a branch circuit ij, wherein ij is the circuit or the branch circuit, and the serial numbers of the nodes at the head end and the tail end are i and j respectively;
is the voltage phase angle vector at time t node; subscripts i and j are numbers of different nodes of the power network.
The alternating current power flow matrixing operation model has the following specific expression:
in the formula: pt、QtRespectively are active power vectors and reactive power vectors of nodes in the power system at time t; u shapetIs a node voltage vector in the power system at time t; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; is an orientation quantity conjugate operation.
Further, the thermodynamic system is based on the idea of graph theory, the node numbers of the thermodynamic system are from 1 to the nodes, that is, the total number of the nodes of the thermodynamic system is the nodes; the number of the pipes or branches of the thermodynamic system is 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 comprises:
the hydraulic model matrixing operation model has the following specific expression:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a thermal medium mass flow vector of a thermal pipeline; 0 is a zero matrix vector;
a basic loop matrix related to a thermodynamic system closed loop;
is the thermodynamic pipeline pressure loss vector at time t;
is the vector of the head of the hot working medium raised by the pressure circulation pump at time t.
The node pressure and pressure drop matrixing operation model has the following specific expression:
wherein:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a node pressure vector in the thermodynamic system at time t;
is the thermodynamic pipeline pressure loss vector at time t;
is the hot working medium pressure head vector raised by the pressure circulating pump at the time t; the superscript b is the total number of branches of the thermodynamic system; superscript' is a vector transposition operation;
is the pressure at the node in the thermodynamic system at time t;
is the pressure loss in thermodynamic pipe b at time t;
is the head of the hot working medium 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:
in the formula:
respectively inputting and outputting thermal power vectors of a thermal pipeline in a thermal system at time t;
is a node thermal load vector in the thermodynamic system at time t; cpThe specific heat capacity of the hot working medium;
the thermal medium mass flow vector is the thermal medium mass flow vector of the pipeline terminal node of the thermodynamic system at time t;
respectively are temperature vectors of an inlet port and an outlet port at a junction of the thermal pipeline and the thermal pipeline at time t;
respectively is a node-mass flow inflow pipeline starting point incidence matrix and a node-mass flow outflow pipeline end point incidence matrix in the thermodynamic system at time t.
Further, the natural gas system comprises a medium-pressure and high-pressure natural gas transmission system, the gas flow in the natural gas pipeline is closely related to the pressure intensity of nodes on two sides of the pipeline and the physical conditions of pipeline transmission, and the mathematical relationship is as follows:
in the formula: subscripts i, j and ij are a natural gas pipeline head end node, a pipeline tail end node and a pipeline number respectively;
respectively the natural gas pressure at the natural gas system nodes i and j at the time t;
is the natural gas flow direction variable of the natural gas system at time t;
is a pipeline characteristic constant;
is the amount of air flow in the duct at time t; kijIs a pipeline branch equivalent characteristic physical parameter, also called natural gas pipeline branch impedance;
is the square difference of the equivalent pressure at the head end and the tail end of the branch.
Numbering nodes in the topological relation of the medium-pressure and high-pressure natural gas systems based on the idea of graph theory, wherein the number of the nodes is from 1 to the node, namely the total number of the nodes is the node; numbering natural gas pipelines or topological branches, and the number of the branchesFrom 1 to b, i.e. the total number of branches is b. Then, the compression ratio of the compressor of each branch can be respectively numbered as
To
The natural gas system pipeline flow direction variables can respectively correspond to branch numbers of
To
The branch equivalent characteristic parameters can respectively correspond to branch numbers K1To Kb(ii) a The equivalent pressure square differences of the head and the tail ends of the branch can respectively correspond to the serial numbers of the branches as
To
The injected natural gas flow rate of each node can be respectively numbered as
To
The natural gas flow of each branch can be respectively numbered as
To
The natural gas flow consumed by the compressor on each branch can be respectively numbered as
To
Further, the step (4) of establishing a matrixing operation model of the natural gas system comprises the following steps: 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:
in the formula:
a natural gas system node-branch complete incidence matrix; f. oftIs the natural gas flow vector through the pipeline in the natural gas network at time t;
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
the natural gas flow vector input into the branch compressor in the natural gas network at time t; q. q.stInjecting flow vectors for nodes in the natural gas system at time t; the superscript b is the total number of the natural gas pipeline branches; subscript node is the total number of topological nodes of the natural gas system; superscript' is a vector transposition operation;
is the natural gas flow rate in branch b in the natural gas system at time t;
is the natural gas flow consumed by the compressor on time t leg b;
the flow of natural gas is injected at time tnode.
The pressure distribution matrixing operation model has the following specific expression:
wherein:
in the formula:
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
is a diagonal matrix relating to the compression ratio of the branch compressors in the natural gas network at time t;
a natural gas system node-branch complete terminal incidence matrix;
is the node pressure flat direction quantity in the natural gas system at time t; II typetThe equivalent pressure square error vector of the head end and the tail end of a branch in the natural gas system at the time t is shown; superscript' is a vector transposition operation;
the compressor compression ratio for leg b at time t;
is the square of the pressure at the node in the natural gas system at time t;
is the natural gas flow direction of branch b in the natural gas system at time t; kbIs the physical parameter of the pipeline branch b;
is the flow of natural gas through line b at time t.
Further, the natural gas system further comprises:
when there is no compressor in the ith branch, setting compression ratio
And setting the natural gas flow consumed by the compressor
In particular, for the electric compressor, provision is always made for
Secondly, the gas flow direction of the natural gas system is usually determined, and particularly in a high-pressure natural gas transmission system, the gas flow direction cannot be easily changed in various regulating valves and practical engineering application, 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 node
When the value is positive, it is the gas load node
When the value is negative, as the intermediate node
The value is 0.
The matrixing equation of the medium-high pressure natural gas transmission system is also suitable for the low-pressure gas distribution system, and at the moment
Has the advantages that: compared with the prior art, the electric heating gas interconnection multi-energy system matrixing operation model modeling method comprehensively considers the electric heating gas heterogeneous energy subsystem at the same time, and can provide theoretical guidance for modeling and running 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 interconnected multi-energy system on the system level, and saves the time of engineering service site and manpower and material resources.
Drawings
FIG. 1 is a flow chart of a method of the present invention.
Detailed Description
The present invention will be further described with reference to the following examples.
Example 1:
a modeling method of a matrix operation model of an electric-heating-gas interconnection multi-energy system is shown in figure 1 and comprises the following steps:
(1) acquiring data information of electric heating and gas interconnection multi-energy system
The method comprises the steps of obtaining data information of the electric-heat-gas interconnection multi-energy system, wherein the data information comprises data such as a coupling interconnection framework of the electric-heat-gas interconnection multi-energy system, energy conversion coupling equipment, output power upper and lower 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 circuit parameters, thermodynamic system pipeline branch parameters, natural gas system pipeline branch parameters and the like.
(2) Establishing a matrixing operation model of an electric power system
Establishing a direct current power flow equation for the power system, wherein a specific expression of a matrix form is as follows:
in the formula:
for injecting a power vector at a time t node, where the elements
Here, the
And
the output power and the load of the generator at a time t node i are respectively shown, and subscripts i and j are numbers of different nodes of the power network;
is a susceptance matrix in which the element Bij＝-1/xij、
B. x is susceptance and reactance of the power network circuit or branch circuit respectively, ij is the circuit or branch circuit, and the serial numbers of the nodes at the head end and the tail end of the circuit or branch circuit are i and j respectively;
is a vector of voltage phase angles at node t at time t.
The classical alternating current power flow model in the power system is described by a matrix form, and the mathematical model is concretely as follows:
in the formula: pt、QtRespectively are active power vectors and reactive power vectors of nodes in the power system at time t; u shapetIs a node voltage vector in the power system at time t; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; is an orientation quantity conjugate operation.
(3) Establishing a thermodynamic system matrixing operation model
Aiming at a thermodynamic system, based on the idea of graph theory, the node number of the thermodynamic system is from 1 to node, namely the total number of the nodes of the thermodynamic system is the node; the number of the pipes or branches of the thermodynamic system is 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:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a thermal medium mass flow vector of a thermal pipeline; 0 is a zero matrix vector;
a basic loop matrix related to a thermodynamic system closed loop;
is the thermodynamic pipeline pressure loss vector at time t;
is the vector of the head of the hot working medium raised by the pressure circulation pump at time t.
When the node pressure and the pressure drop in the thermodynamic system are subjected to matrixing operation, the specific expression is as follows:
wherein:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a node pressure vector in the thermodynamic system at time t;
is the thermodynamic pipeline pressure loss vector at time t;
is the hot working medium pressure head vector raised by the pressure circulating pump at the time t; the superscript b is the total number of branches of the thermodynamic system;
is the pressure at the node in the thermodynamic system at time t;
is the pressure loss in thermodynamic pipe b at time t;
is the head of the hot working medium raised by the pressure circulating pump in the thermodynamic system at time t.
The matrix operation model for expressing the flow and the temperature of the thermodynamic system by using a matrix equation set is concretely as follows:
in the formula:
respectively inputting and outputting thermal power vectors of a thermal pipeline in a thermal system at time t;
is a node thermal load vector in the thermodynamic system at time t; cpThe specific heat capacity of the hot working medium;
the thermal medium mass flow vector is the thermal medium mass flow vector of the pipeline terminal node of the thermodynamic system at time t;
respectively are temperature vectors of an inlet port and an outlet port at a junction of the thermal pipeline and the thermal pipeline at time t;
respectively is a node-mass flow inflow pipeline starting point incidence matrix and a node-mass flow outflow pipeline end point incidence matrix in the thermodynamic system at time t.
(4) Establishing a natural gas system matrixing operation model
For a medium-high pressure natural gas transmission system, the gas flow in a natural gas pipeline is closely related to the pressure of nodes on two sides of the pipeline and the physical conditions of pipeline transmission, and the mathematical relationship is as follows:
in the formula: subscripts i, j and ij are a natural gas pipeline head end node, a pipeline tail end node and a pipeline number respectively;
respectively the natural gas pressure at the natural gas system nodes i and j at the time t;
is the natural gas flow direction variable of the natural gas system at time t;
is a pipeline characteristic constant;
is the amount of air flow in the duct at time t; kijIs a pipeline branch equivalent characteristic physical parameter, also called natural gas pipeline branch impedance;
is the square difference of the equivalent pressure at the head end and the tail end of the branch.
Analyzing the topological relation of the medium and high-voltage natural gas transmission system based on the idea of graph theory, numbering topological nodes of the natural gas transmission system, and counting the number of the nodes from 1 to the node, namely counting the total number of the nodes as the node; numbering the natural gas pipelines or topological branches, wherein the number of the branches is from 1 to b, namely the total number of the branches is b. Then, the compression ratio of the compressor of each branch can be respectively numbered as
To
The natural gas system pipeline flow direction variables can respectively correspond to branch numbers of
To
The branch equivalent characteristic parameters can respectively correspond to branch numbers K1To Kb(ii) a The equivalent pressure square differences of the head and the tail ends of the branch can respectively correspond to the serial numbers of the branches as
To
The injected natural gas flow rate of each node can be respectively numbered as
To
The natural gas flow of each branch can be respectively numbered as
To
The natural gas flow consumed by the compressor on each branch can be respectively numbered as
To
The expression of the flow balance matrixing equation of the natural gas system is as follows:
wherein:
in the formula:
a natural gas system node-branch complete incidence matrix; f. oftFor flowing through pipes in natural gas networks at times tA natural gas flow vector;
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
the natural gas flow vector input into the branch compressor in the natural gas network at time t; q. q.stInjecting flow vectors for nodes in the natural gas system at time t; the superscript b is the total number of the natural gas pipeline branches; subscript node is the total number of topological nodes of the natural gas system; superscript' is a vector transposition operation;
is the natural gas flow rate in branch b in the natural gas system at time t;
is the natural gas flow consumed by the compressor on time t leg b;
the flow of natural gas is injected at time tnode.
The expression of the pressure distribution matrixing equation of the natural gas system is as follows:
wherein:
in the formula:
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
is a diagonal matrix relating to the compression ratio of the branch compressors in the natural gas network at time t;
a natural gas system node-branch complete terminal incidence matrix;
is the node pressure flat direction quantity in the natural gas system at time t; II typetThe equivalent pressure square error vector of the head end and the tail end of a branch in the natural gas system at the time t is shown; superscript' is a vector transposition operation;
the compressor compression ratio for leg b at time t;
is the square of the pressure at the node in the natural gas system at time t;
is the natural gas flow direction of branch b in the natural gas system at time t; kbIs the physical parameter of the pipeline branch b;
is the flow of natural gas through line b at time t.
The matrixing model for the natural gas system is explained as follows:
when there is no compressor in the ith branch, setting compression ratio
And setting the natural gas flow consumed by the compressor
In particular, for the electric compressor, provision is always made for
Secondly, the gas flow direction of the natural gas system is usually determined, and particularly in a high-pressure natural gas transmission system, the gas flow direction cannot be easily changed in various regulating valves and practical engineering application, 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 node
When the value is positive, it is the gas load node
When the value is negative, as the intermediate node
The value is 0.
The matrixing equation of the medium-high pressure natural gas transmission system is also suitable for the low-pressure gas distribution system, and at the moment
(5) Outputting matrixed model data information
Outputting matrixing model data information including data information such as power supply output of an electric power system, node electric power, node reactive power, node voltage phase angle, node voltage and the like; data information of heat source output, node mass flow, node thermal power, node pressure, pipeline pressure loss and the like of the thermodynamic system; the natural gas system comprises 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.
The above description is only of the preferred embodiments of the present invention, and it should be 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 invention and these are intended to be within the scope of the invention.
Claims (9)
1. A modeling method for a matrix operation model 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 heating and gas interconnection multi-energy system;
establishing a matrixing operation model of the power system according to the data information, 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 the thermodynamic system according to the data information, 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;
and establishing a matrixing operation model of the natural gas system according to the data information, wherein the matrixing operation model comprises a flow balance matrixing operation model and a pressure distribution matrixing operation model.
2. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model of claim 1 is characterized in that: the data information of the electric-heat-gas interconnection multi-energy system comprises a coupling interconnection framework of the electric-heat-gas interconnection multi-energy system, energy conversion coupling equipment, power supplies, upper and lower limits of output power of 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 circuit parameters, thermodynamic system pipeline branch parameters and natural gas system pipeline branch parameter data information.
3. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model of claim 1 is characterized in that: further comprising:
acquiring data information output by a matrixing operation model of all electric power systems, thermodynamic systems and natural gas systems, wherein the data information comprises power supply output of the electric power systems, node electric power, node reactive power, node voltage phase angle and node voltage data information; data information of heat source output, node mass flow, node thermal power, node pressure and pipeline pressure loss of the thermodynamic system; the natural gas system comprises the data information of gas source output, node pressure, pipeline natural gas flow, node natural gas flow, compressor consumption flow and compressor compression ratio.
4. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model according to claim 1 or 3, characterized in that: the specific expression of the direct current power flow matrixing operation model is as follows:
in the formula:
for injecting a power vector at a time t node, where the elements
Here, the
And
respectively generator output and load at time t node i,
being a susceptance matrix, element B of the susceptance matrixij＝-1/xij、
Bij、xijRespectively the susceptance and reactance of a power network line or a branch ij, wherein ij is the line or the branch, and the serial numbers of the nodes at the head end and the tail end are respectivelyi、j；
Is the voltage phase angle vector at time t node; subscripts i and j are numbers of different nodes of the power network;
the specific expression of the alternating current power flow matrixing operation model is as follows:
in the formula: pt、QtRespectively are active power vectors and reactive power vectors of nodes in the power system at time t; u shapetIs a node voltage vector in the power system at time t; y is a node admittance matrix; re and Im are respectively the real part and the imaginary part of the vector; is an orientation quantity conjugate operation.
5. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model according to claim 1 or 3, characterized in that: the thermodynamic system is based on the idea of graph theory, the node number of the thermodynamic system is from 1 to the node, namely the total number of the nodes of the thermodynamic system is the node; the number of the pipes or branches of the thermodynamic system is from 1 to b, i.e. the total number of branches of the thermodynamic system is b.
6. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model according to claim 5 is characterized in that: the hydraulic model matrixing operation model has the following specific expression:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a thermal medium mass flow vector of a thermal pipeline; 0 is a zero matrix vector;
a basic loop matrix related to a thermodynamic system closed loop;
is the thermodynamic pipeline pressure loss vector at time t;
is the hot working medium pressure head vector raised by the pressure circulating pump at the time t;
the node pressure and pressure drop matrixing operation model has the following specific expression:
wherein:
in the formula:
is a complete correlation matrix related to the thermodynamic system nodes and the thermodynamic pipelines;
is a node pressure vector in the thermodynamic system at time t;
is the thermodynamic pipeline pressure loss vector at time t;
is the hot working medium pressure head vector raised by the pressure circulating pump at the time t; the superscript b is the total number of branches of the thermodynamic system; superscript' is a vector transposition operation;
is the pressure at the node in the thermodynamic system at time t;
is the pressure loss in thermodynamic pipe b at time t;
is the head of the hot working medium 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:
in the formula:
respectively inputting and outputting thermal power vectors of a thermal pipeline in a thermal system at time t;
is a node thermal load vector in the thermodynamic system at time t; cpThe specific heat capacity of the hot working medium;
the thermal medium mass flow vector is the thermal medium mass flow vector of the pipeline terminal node of the thermodynamic system at time t;
respectively at time tmoriThe temperature vectors of an inlet port and an outlet port at the junction of the pipeline and the heat distribution pipeline;
respectively is a node-mass flow inflow pipeline starting point incidence matrix and a node-mass flow outflow pipeline end point incidence matrix in the thermodynamic system at time t.
7. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model according to claim 1 or 3, characterized in that: the natural gas system comprises a medium-pressure and high-pressure natural gas transmission system, wherein the gas flow in the natural gas pipeline is closely related to the pressure intensity of nodes on two sides of the pipeline and the physical conditions of pipeline transmission, and the mathematical relationship is as follows:
in the formula: subscripts i, j and ij are a natural gas pipeline head end node, a pipeline tail end node and a pipeline number respectively;
respectively the natural gas pressure at the natural gas system nodes i and j at the time t;
is the natural gas flow direction variable of the natural gas system at time t;
is a pipeline characteristic constant;
is the amount of air flow in the duct at time t; kijIs a pipeline branch equivalent characteristic physical parameter, also called natural gas pipeline branch impedance;
the square difference of the equivalent pressure at the head end and the tail end of the branch is taken as the pressure difference of the equivalent pressure at the head end and the tail end of the branch;
numbering nodes in the topological relation of the medium-pressure and high-pressure natural gas systems based on the idea of graph theory, wherein the number of the nodes is from 1 to the node, namely the total number of the nodes is the node; numbering natural gas pipelines or topological branches, wherein the number of the branches is from 1 to b, namely the total number of the branches is b; then, the compression ratio of the compressor of each branch can be respectively numbered as
To
The natural gas system pipeline flow direction variables can respectively correspond to branch numbers of
To
The branch equivalent characteristic parameters can respectively correspond to branch numbers K1To Kb(ii) a The equivalent pressure square differences of the head and the tail ends of the branch can respectively correspond to the serial numbers of the branches as
To
The injected natural gas flow rate of each node can be respectively numbered as
To
The natural gas flow of each branch can be respectively numbered as
To
The natural gas flow consumed by the compressor on each branch can be respectively numbered as
To
8. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model of claim 7 is characterized in that: the flow balance matrixing operation model has the following specific expression:
wherein:
in the formula:
a natural gas system node-branch complete incidence matrix; f. oftIs the natural gas flow vector through the pipeline in the natural gas network at time t;
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
at time t daysThe natural gas flow vector input into the branch compressor in the natural gas network; q. q.stInjecting flow vectors for nodes in the natural gas system at time t; the superscript b is the total number of the natural gas pipeline branches; subscript node is the total number of topological nodes of the natural gas system; superscript' is a vector transposition operation;
is the natural gas flow rate in branch b in the natural gas system at time t;
is the natural gas flow consumed by the compressor on time t leg b;
is the flow of injected natural gas at time tnode;
the pressure distribution matrixing operation model has the following specific expression:
wherein:
in the formula:
a complete starting point incidence matrix related to network nodes and pipeline branches in the natural gas network;
is a diagonal matrix relating to the compression ratio of the branch compressors in the natural gas network at time t;
a natural gas system node-branch complete terminal incidence matrix;
is the node pressure flat direction quantity in the natural gas system at time t; II typetThe equivalent pressure square error vector of the head end and the tail end of a branch in the natural gas system at the time t is shown; superscript' is a vector transposition operation;
the compressor compression ratio for leg b at time t;
is the square of the pressure at the node in the natural gas system at time t;
is the natural gas flow direction of branch b in the natural gas system at time t; kbIs the physical parameter of the pipeline branch b;
is the flow of natural gas through line b at time t.
9. The modeling method of the electric-thermal-gas interconnection multi-energy system matrixing operation model of claim 7 is characterized in that: further comprising:
when no compressor is arranged in the ith branch, the compression ratio is set at the moment
And setting the natural gas flow consumed by the compressor
In particular, for the electric compressor, provision is always made for
The gas flow direction of a natural gas system is usually determined, particularly in a high-pressure natural gas transmission system, the gas flow direction cannot be easily changed in various regulating valves and practical engineering application, 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 node
When the value is positive, it is the gas load node
When the value is negative, as the intermediate node
The value is 0;
the matrixing equation of the medium-high pressure natural gas transmission system is also suitable for 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 true CN114004047A (en) | 2022-02-01 |
| CN114004047B 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | 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, 2 March 2019 (2019-03-02), pages 1 - 27 * |
| 常晓勇;王玉婷;陈福锋;李玉平;徐程骥;: "电力-天然气综合能源系统建模与仿真", 浙江电力, no. 07, 25 July 2018 (2018-07-25) * |
| 杨经纬;张宁;康重庆;: "多能源网络的广义电路分析理论――(一)支路模型", 电力系统自动化, no. 09, 10 May 2020 (2020-05-10) * |
| 钟永洁等: "含电—热—气—冷子系统的区域综合能源系统多场景优化调度", 电力系统自动化, vol. 43, no. 12, 25 June 2019 (2019-06-25), pages 76 - 84 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114004047B (en) | 2024-05-28 |
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 | |
| CN111815042A (en) | Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model | |
| CN111899120A (en) | Method for establishing comprehensive energy planning and operation combined optimization model | |
| CN109347107A (en) | One kind be incorporated into the power networks electric heating interconnection integrated energy system tidal current computing method | |
| CN107947245B (en) | Equivalent optimal power flow model construction method considering natural gas system constraint | |
| CN107910871A (en) | A kind of energy internet is provided multiple forms of energy to complement each other control method | |
| Zhang et al. | Optimization of distributed integrated multi-energy system considering industrial process based on energy hub | |
| CN109102104A (en) | A kind of multi-energy system synergism low-carbon planing method based on carbon emission flow | |
| CN112035984A (en) | Electric-gas-storage area comprehensive energy system collaborative planning method | |
| CN110968827A (en) | Optimal configuration method for multi-region comprehensive energy system | |
| CN114004047B (en) | Modeling method for matrixing operation model of electric heating gas interconnection multi-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 | |
| CN111310310B (en) | Thermodynamic system static state trend rapid decoupling calculation method for quantity adjustment | |
| CN112214883A (en) | Micro-energy-grid comprehensive energy system operation optimization method and system | |
| CN114462163B (en) | exergy hub model considering energy quality | |
| CN115545340A (en) | Comprehensive energy station-network collaborative planning method considering double uncertainties | |
| CN109002963A (en) | Electric heating energy market equilibrium method for solving and device | |
| Li et al. | Exergy Hub: A novel energy hub model considering energy quality | |
| CN115034453A (en) | Flexible operation method and computing platform of comprehensive energy system considering carbon trading market | |
| CN114004075A (en) | Multi-objective integrated design and optimal control method and system for comprehensive energy system | |
| Hu et al. | Multi-Energy Flow Calculation Model and Solution Method for Regional Integrated Energy System | |
| CN114066056A (en) | Optimal scheduling method and system considering flexibility of thermoelectric cooperative comprehensive energy system | |
| Ayele et al. | Pseudo-dynamic simulation on a district energy system made of coupling technologies |
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 |