CN113343531A - Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference - Google Patents

Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference Download PDF

Info

Publication number
CN113343531A
CN113343531A CN202110684263.9A CN202110684263A CN113343531A CN 113343531 A CN113343531 A CN 113343531A CN 202110684263 A CN202110684263 A CN 202110684263A CN 113343531 A CN113343531 A CN 113343531A
Authority
CN
China
Prior art keywords
mass flow
node
pressure
equal
natural gas
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
Application number
CN202110684263.9A
Other languages
Chinese (zh)
Other versions
CN113343531B (en
Inventor
王程
吴科宏
毕天姝
王敬尧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202110684263.9A priority Critical patent/CN113343531B/en
Publication of CN113343531A publication Critical patent/CN113343531A/en
Application granted granted Critical
Publication of CN113343531B publication Critical patent/CN113343531B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an electricity-gas comprehensive energy system dynamic energy flow obtaining method based on explicit difference, which comprises the steps of firstly obtaining the linear relations between the mass flow and pressure of any node in different areas in a natural gas pipeline difference grid, boundary conditions and initial conditions according to a difference equation of the mass flow and the pressure; constructing a two-port model of mass flow and pressure of any node on the natural gas pipeline differential grid with respect to boundary conditions and initial conditions according to the obtained linear relation; and embedding the two-port model into a dynamic energy flow calculation model of the electric-gas comprehensive energy system to obtain the dynamic energy flow of the electric-gas comprehensive energy system. The method converts the partial differential equation of the natural gas pipeline flow into a time domain linear algebraic model, is convenient to embed a time domain optimization model, can greatly reduce the state variable scale and the model calculation cost compared with the traditional finite difference method, and improves the calculation efficiency.

Description

Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference
Technical Field
The invention relates to the technical field of power systems, in particular to a method for acquiring dynamic energy flow of an electricity-gas comprehensive energy system based on explicit difference.
Background
The electricity-gas comprehensive energy system is a novel energy system integrating generation, transmission, distribution and consumption of heterogeneous energy such as electric energy and natural gas, can effectively improve the comprehensive energy efficiency and the wind-light power generation consumption level of a user side terminal, is an important component of a clean, low-carbon, safe and efficient energy system in China, and is a beneficial way for realizing the carbon peak reaching, carbon neutralization and medium-long term targets in China. The calculation of the electricity-gas coupling energy flow is the basis for realizing the efficient operation of the regional electricity-gas comprehensive energy system, the active power flow distribution rule of the power system can be approximately depicted by a direct current power flow model with low calculation cost and high accuracy, the energy flow distribution rule of the natural gas system is described by a partial differential equation, and the calculation and analysis cost is high.
At present, the existing research adopts an implicit finite difference method to convert a partial differential equation into a linear algebraic equation to realize the high-precision calculation of the electricity-gas energy flow, but because a large amount of auxiliary variables and constraints are introduced at the same time, the calculation cost is higher; meanwhile, a plurality of research works try to construct a natural gas dynamic model with high precision and low complexity in a frequency domain and a complex frequency domain, but an equivalent conversion condition of a time domain-frequency domain/a complex frequency domain needs to be additionally added during calculation of the electric-gas coupling energy flow, so that the scale of the electric-gas coupling energy flow calculation model is increased, and the calculation efficiency is reduced.
Disclosure of Invention
The invention aims to provide a method for acquiring dynamic energy flow of an electricity-gas integrated energy system based on explicit difference.
The purpose of the invention is realized by the following technical scheme:
an explicit differential-based method for acquiring dynamic energy flow of an electricity-gas integrated energy system, the method comprising:
step 1, obtaining linear relations between mass flow and pressure of any node in different areas in a natural gas pipeline differential grid, boundary conditions and initial conditions according to a differential equation of the mass flow and the pressure;
step 2, constructing a two-port model of mass flow and pressure of any node on the natural gas pipeline differential grid with respect to boundary conditions and initial conditions according to the linear relation obtained in the step 1;
and 3, embedding the two-port model into a dynamic energy flow calculation model of the electric-gas comprehensive energy system to obtain the dynamic energy flow of the electric-gas comprehensive energy system.
According to the technical scheme provided by the invention, the partial differential equation of the natural gas pipeline flow is converted into the time domain linear algebraic model by the method, so that the time domain optimization model is conveniently embedded, and compared with the traditional finite difference method, the method can greatly reduce the state variable scale and the model calculation cost and improve the calculation efficiency.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for acquiring a dynamic energy flow of an electricity-gas integrated energy system based on an explicit difference according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a dynamic three-layer explicit differential grid of a natural gas pipeline according to an embodiment of the present invention;
fig. 3 is a schematic diagram of region division of the differential mesh according to the embodiment of the present invention;
FIG. 4 is a schematic diagram of the differential type of the regions 1 and 2 according to the embodiment of the present invention;
FIG. 5 is a schematic diagram of the differential form of the regions 3 and 4 according to the embodiment of the present invention;
FIG. 6 is a schematic view of an exemplary electro-pneumatic energy system topology according to the present invention;
FIG. 7 is a schematic view of a single pipe according to an exemplary embodiment of the present invention;
FIG. 8 is a graph showing a comparison of single pipe results for two processes in accordance with an exemplary embodiment of the present invention;
FIG. 9 is a schematic view of topology and pipeline parameters of an exemplary natural gas testing system in accordance with the present invention;
FIG. 10 is a schematic view of a dynamic curve of gas pressures at nodes and mass flows at ends of pipes in a natural gas system according to an embodiment of the present invention;
FIG. 11 is a graph illustrating a comparison of the load of the power system and the output of the generator according to an exemplary embodiment of the present invention;
FIG. 12 is a graph illustrating the comparison of gas pressure and gas turbine consumption for a natural gas system according to an exemplary embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention are clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the present invention will be further described in detail with reference to the accompanying drawings, and as shown in fig. 1, a flow chart of a method for acquiring a dynamic power flow of an electricity-gas integrated energy system based on an explicit difference provided by the embodiment of the present invention is schematically shown, where the method includes:
step 1, obtaining linear relations between mass flow and pressure of any node in different areas in a natural gas pipeline differential grid, boundary conditions and initial conditions according to a differential equation of the mass flow and the pressure;
in this step, first, the second-order hyperbolic partial differential equation of the mass flow is expressed as:
Figure BDA0003123813320000031
the second order hyperbolic partial differential equation for pressure is expressed as:
Figure BDA0003123813320000032
wherein M is natural gas mass flow, pi is natural gas pressure, c is gas sound velocity, lambda is friction coefficient, D is pipeline diameter, x is space variable, t is time variable,
Figure BDA0003123813320000033
in order to be the average flow rate,
Figure BDA0003123813320000034
is the sign of the partial derivative;
then, the solution is regionalized into discrete points according to the natural gas pipeline differential grid, as shown in fig. 2, which is a schematic diagram of a dynamic three-layer explicit differential grid of the natural gas pipeline according to an embodiment of the present invention,
Figure BDA0003123813320000035
representing the mass flow M at node (j, k)(x＝j,t＝k)；
Discretizing the differential equation into a differential equation, approximating a second order partial derivative using a second order center difference quotient, expressed as:
Figure BDA0003123813320000036
Figure BDA0003123813320000037
wherein tau is a time infinitesimal; h is a spatial infinitesimal; u (x)j,tk) Representing the pressure or mass flow at the point with the space serial number j (j is more than or equal to 0 and less than or equal to N) and the time serial number k (T is more than or equal to 0 and less than or equal to T) in the differential grid;
the first order partial derivative of approximate time using the forward difference quotient is expressed as:
Figure BDA0003123813320000038
substituting equations (3), (4) and (5) into the second-order hyperbolic partial differential equation (1) of the mass flow rate, the algebraic form of the difference with respect to the mass flow rate is obtained as:
Figure BDA0003123813320000041
reusing spatial direction transformations
Figure BDA0003123813320000042
Combined vertical type (5) and (6) are obtained
Figure BDA0003123813320000043
In the formula, K1、K2、K3、K4Is a difference coefficient whose value is:
Figure BDA0003123813320000044
by giving boundary conditions and initial conditions, solving a differential expression (8) of a second-order hyperbolic partial differential equation about mass flow, and for a section of natural gas pipeline, assuming that the head-end air pressure and the tail-end mass flow are known, the following steps are provided:
Figure BDA0003123813320000045
in the formula (I), the compound is shown in the specification,
Figure BDA0003123813320000046
indicating the gas pressure at the inlet (x ═ N) of the duct(x＝N,t)；Πin(t) is a function of the duct inlet pressure over time;
Figure BDA0003123813320000047
denotes the mass flow M at the outlet of the pipe (x ═ 0)(x＝0,t)；Mout(t) is the initial value of the mass flow at the outlet of the pipeline;
assuming that the system is in a steady state at the initial moment, the mass flow is satisfied
Figure BDA0003123813320000048
Considering that the state quantity of equation (8) contains only mass flow, the pressure boundary condition of equation (9) needs to be converted into a mass flow boundary condition, for which the first-order partial derivative of the forward difference quotient approximation space is used:
Figure BDA0003123813320000049
the united type (1), (5) and (11) equivalently converts the pressure boundary condition (9) into a mass flow boundary condition:
Figure BDA00031238133200000410
in the formula (I), the compound is shown in the specification,
Figure BDA00031238133200000411
indicating the gas pressure pi at the node (j, k)(x＝j,t＝k)；
The boundary conditions of the three-layer explicit difference equation for the mass flow are thus obtained, namely equations (9), (12); and initial conditions, i.e. formulae (9), (10);
in a similar way, firstly, according to a second-order hyperbolic partial differential equation (2) of the pressure, a differential algebraic form of the pressure is obtained in a joint type (3), (4) and (5):
Figure BDA0003123813320000051
secondly, the initial conditions of the formula (13) are the air pressure of each node at the initial time and the air pressure change rate at the initial time;
the joint type (5) and (11) obtain the air pressure of each node at the initial moment as follows:
Figure BDA0003123813320000052
assuming that the system is in a steady state at the initial time, the air pressure change rate at the initial time is:
Figure BDA0003123813320000053
finally, considering that the state quantity of equation (13) contains only pressure, it is necessary to convert the mass flow boundary into a pressure boundary, which is satisfied at the pipe outlet according to equation (9):
Figure BDA0003123813320000054
boundary conditions for three-layer explicit differential equations for pressure are thus obtained, namely equations (9), (16); and the initial conditions, i.e., equations (14), (15).
Step 2, constructing a two-port model of mass flow and pressure of any node on the natural gas pipeline differential grid with respect to boundary conditions and initial conditions according to the linear relation obtained in the step 1;
in this step, firstly, according to the difference in distance between the point to be solved and the boundary point and the initial point on the natural gas pipeline differential grid, the area to be solved is divided into four areas, as shown in fig. 3, which is a schematic diagram of area division of the differential grid according to the embodiment of the present invention, and the respective areas are:
1) region 1(k ═ 1, j < N); 2) region 2(k ═ 1, j ═ N); 3) region 3 (k is more than or equal to 2 and less than or equal to T, j is more than or equal to 1 and less than N); 4) region 4(2 ≦ k ≦ T, j ═ N);
firstly, a two-port model of mass flow of any node on a natural gas pipeline differential grid with respect to boundary conditions and initial conditions is constructed, specifically:
(1) region 1: when k is 1 and j is less than N, as shown in fig. 4, it is a schematic diagram of a differential form of the regions 1 and 2 according to the embodiment of the present invention, and the differential form of the mass flow rate is as shown in fig. 4(a), specifically:
Figure BDA0003123813320000055
assuming that the mass flow of each node of the pipeline at the initial moment is equal to the mass flow at the tail end of the pipeline, when k is less than 2,
Figure BDA0003123813320000056
and then, satisfy:
Figure BDA0003123813320000057
bringing formula (18) into formula (17) yields:
Figure BDA0003123813320000061
in the formula (I), the compound is shown in the specification,
Figure BDA0003123813320000062
to represent
Figure BDA0003123813320000063
And
Figure BDA0003123813320000064
the linear correlation coefficient of (a) has a value of:
Figure BDA0003123813320000065
(2) region 2: when k is 1 and j is N, the differential form of the mass flow rate is as shown in fig. 4(b), specifically:
Figure BDA0003123813320000066
in the formula:
Figure BDA0003123813320000067
when the region 1 and the region 2 are combined into the region a and the equations (19) and (20) are observed, the mass flow satisfies the following conditions when k is 1 and j is 1 ≦ N (i.e., the region 1+ the region 2):
Figure BDA0003123813320000068
the matrix form of the formula (21) is
Figure BDA0003123813320000069
In the formula:
Figure BDA00031238133200000610
Figure BDA00031238133200000611
(3) region 3: when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than N, as shown in FIG. 5, the difference form of the areas 3 and 4 in the embodiment of the present invention is schematically shown, and the difference form of the mass flow is shown in FIG. 5(a), and satisfies the following conditions:
Figure BDA00031238133200000612
(4) region 4: when k is 2 ≦ T, and j is N, the differential form of the mass flow is as shown in fig. 5(b), and satisfies:
Figure BDA00031238133200000613
when the region 3 and the region 4 are merged into the region B, and the formula (23) and the formula (24) are observed, when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than or equal to N (namely the region 3+ the region 4), the mass flow satisfies:
Figure BDA00031238133200000614
so far, a two-port form of mass flow rates of the area A and the area B is obtained, the area A and the area B are combined, and when k is more than or equal to 1 and less than or equal to T and j is more than or equal to 1 and less than or equal to N, the mass flow rates satisfy:
Figure BDA0003123813320000071
in the formula, Aj,BjDifference coefficient matrices, m, of (T × T) respectivelyjThe mass flow for a node with spatial index j is a (T × 1) matrix:
Figure BDA0003123813320000072
equation (26) is a two-port model of the mass flow of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition;
then, a two-port model of the pressure of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition is constructed, specifically:
(1) region 1: when k is 1 and j < N, the pressure difference satisfies the following conditions:
Figure BDA0003123813320000073
as can be seen from equation (27) and the air pressure at each node at the initial time of synchronization, when k is 0 and j is equal to or less than 1 and equal to or less than N, the air pressure satisfies:
Figure BDA0003123813320000074
in the formula
Figure BDA0003123813320000075
The rate of change of the gas pressure at the initial time and the boundary conditions at the outlet of the pipeline are brought into (27) to obtain:
Figure BDA0003123813320000076
the formula (29) is arranged as
Figure BDA0003123813320000077
Wherein:
Figure BDA0003123813320000081
the matrix form of equation (30) is:
Figure BDA0003123813320000082
in the formula:
Figure BDA0003123813320000083
Figure BDA0003123813320000084
(2) region 2: when k is 1 and j is N, the pressure satisfies:
Figure BDA0003123813320000085
in the formula:
Figure BDA0003123813320000086
(3) region 3: when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than N, the pressure satisfies the following conditions:
Figure BDA0003123813320000087
(4) region 4: when k is more than or equal to 2 and less than or equal to T, j is equal to N, the pressure satisfies the following conditions:
Figure BDA0003123813320000088
combining 4 areas, and when k is more than or equal to 1 and less than or equal to T and j is more than or equal to 1 and less than or equal to N, the pressure is satisfied as follows:
Figure BDA0003123813320000089
in the formula, Cj,DjDifference coefficient matrices, n, of (T x T) respectivelyjThe pressure at a node with spatial index j is a (T1) matrix:
Figure BDA00031238133200000810
thus, a two-port model of the pressure of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition is obtained, namely equation (35).
In the specific implementation, because the state quantities of the head end and the tail end of the natural gas pipeline are only related to the boundary conditions and the initial conditions, the time domain two-port model does not contain state variables of other space micro elements, and compared with a traditional finite difference method, the state variable scale and the model calculation cost are greatly reduced.
In addition, the two-port model of the mass flow and the pressure of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition can be further represented by combining:
Figure BDA0003123813320000091
when the initial conditions are head end air pressure and tail end mass flow, the head end mass flow and the tail end air pressure of the pipeline meet the following conditions:
Figure BDA0003123813320000092
thereby obtaining a time domain two-port model (37) of the gas pressure and the mass flow at the head end and the tail end of the natural gas pipeline.
And 3, embedding the two-port model into a dynamic energy flow calculation model of the electric-gas comprehensive energy system to obtain the dynamic energy flow of the electric-gas comprehensive energy system.
In this step, the operation cost of the regional electricity-gas integrated energy system includes the traditional unit power generation cost and the natural gas source gas supply cost, and the objective function of the dynamic energy flow is the minimum operation cost, namely:
Figure BDA0003123813320000093
in the formula (I), the compound is shown in the specification,
Figure BDA0003123813320000094
calculating a set of time periods for the powerflow;
Figure BDA0003123813320000095
is a traditional unit set; c. CgcIs the g thcCost coefficient of the conventional unit;
Figure BDA0003123813320000096
is the g thcTraditional unitAt a time period teInternal output force;
Figure BDA0003123813320000097
calculating the number of time segments for the natural gas energy flow;
Figure BDA0003123813320000098
is a medium-pressure air source set; c. CwThe cost coefficient of the w gas source;
Figure BDA0003123813320000099
for the w gas source during the time period tgMass flow rate inside;
for example, as shown in fig. 6, a topological schematic diagram of an electric-gas integrated energy system according to an example of the present invention is composed of a 13-node natural gas system and a 12-node power system, and includes three power generators (one conventional thermal power generating unit and two gas power generating units, which are respectively connected to nodes 8 and 9 of the natural gas system); 10 electrical loads; three conventional gas loads; two gas turbine gas loads.
The energy flow model of the natural gas system satisfies the following constraints:
(1) pipeline air pressure-mass flow constraint:
according to the time domain two-port model of the gas pressure and the mass flow at the head end and the tail end of the natural gas pipeline, the gas pressure and the mass flow at the head end and the tail end of the p-th pipeline meet the following conditions:
Figure BDA0003123813320000101
in the formula (I), the compound is shown in the specification,
Figure BDA0003123813320000102
respectively representing mass flow and air pressure matrixes at the head end of the p-th pipeline in all time periods;
Figure BDA0003123813320000103
respectively representing the mass flow and the air pressure matrix of the p-th pipeline tail end in all time periods;
(2) air pressure and mass flow restraint of an air source:
Figure BDA0003123813320000104
(3) node mass flow balance constraint:
for natural gas node ngAnd the sum of the mass flow of all pipelines connected with the node and the net injection mass flow of the node is zero:
Figure BDA0003123813320000105
in the formula, ngIs a natural gas node;
Figure BDA0003123813320000106
for accessing natural gas node ngThe gas source set of (2);
Figure BDA0003123813320000107
is the end and natural gas node ngA set of connected pipes;
Figure BDA0003123813320000108
is a head end and natural gas node ngA set of connected pipes; d is an element of Dg(ng) To be located at natural gas node ngA load set of (a);
(4) and (3) node air pressure constraint:
Figure BDA0003123813320000109
in the formula (I), the compound is shown in the specification,Π、
Figure BDA00031238133200001010
the lower limit and the upper limit of the natural gas node air pressure are set;
the power flow model of the power system satisfies the following constraints:
because the power network of the regional electricity-gas comprehensive energy system belongs to a power distribution network, the topology of the regional electricity-gas comprehensive energy system is generally radial, and the power Flow of the power system can be approximated by combining a Branch Flow power Flow model and a second-order conical convex relaxation technology;
(1) node active and reactive power balance constraint:
Figure BDA00031238133200001011
Figure BDA00031238133200001012
in the formula, neIs a power system node;
Figure BDA00031238133200001013
for node n of the power systemeA collection of legacy units;
Figure BDA00031238133200001014
for node n of the power systemeA set of gas turbine units;
Figure BDA00031238133200001015
for node n of the power systemeA set of electrical loads;
Figure BDA0003123813320000111
for the line head section and the node n of the power systemeA set of connected power lines of (c);
Figure BDA0003123813320000112
for line end and power system node neA set of connected power lines of (c);
Figure BDA0003123813320000113
the active power and the reactive power of the traditional unit are obtained;
Figure BDA0003123813320000114
the active power and the reactive power of the thermal power generating unit are obtained;
Figure BDA0003123813320000115
active and reactive loads for the power system; plt、QltThe active power and the reactive power of the line are obtained; r isl、xlThe resistance and reactance of the circuit are shown; i isltIs the square of the branch current;
Figure BDA0003123813320000116
is the ground conductance of the node;
Figure BDA0003123813320000117
is the square of the node voltage;
(2) the line voltage drop equation and the line current equation of the second-order cone-convex relaxation are as follows:
Figure BDA0003123813320000118
Figure BDA0003123813320000119
in the formula, vl-tIs the line end voltage; v. ofl+tIs the line first section voltage;
(3) and (3) generator capacity constraint:
the active power and reactive power upper and lower limits of the generator are constrained as follows:
Figure BDA00031238133200001110
in the formula (I), the compound is shown in the specification,
Figure BDA00031238133200001111
for node n of the power systemeA set of generator sets;P、
Figure BDA00031238133200001112
the upper limit and the lower limit of active power of the generator are set;Q、
Figure BDA00031238133200001113
the upper limit and the lower limit of the reactive power of the generator are set;
the constraints after the electrical-to-electrical coupling are:
(1) the active power and gas consumption of gas power generation are restricted:
Figure BDA00031238133200001114
in the formula, eta is energy conversion efficiency;
(2) the time coupling constraint of the power energy flow calculation variable and the natural gas energy flow calculation variable is as follows:
Figure BDA00031238133200001115
in the formula (I), the compound is shown in the specification,
Figure BDA00031238133200001116
is an upper rounding function; and N is the ratio of the time resolution of the electric and gas decision variables.
The following describes the advantages of the above dynamic power flow acquisition method with specific examples:
1) the results of the calculations of the Implicit finite Difference Method (IDM) and the Method described in this example (Explicit Difference model, Explicit Difference Method, EDM) were compared for a single natural gas pipeline. Fig. 7 is a schematic diagram of a single pipeline according to an exemplary embodiment of the present invention, where the length of the natural gas pipeline is 7783 m, the diameter of the pipeline is 0.5 m, the head-end pressure is maintained at a constant pressure of 0.4MPa, and the end load is changed at intervals of 150 seconds. The length of the defined temporal and spatial bins is 0.5 seconds and 181 meters, respectively.
As shown in fig. 8, which is a comparison diagram of results of a single pipeline in two methods according to the embodiment of the present invention, when the mass flow at the outlet of the pipeline has a step change, the mass flow at the inlet obtained by the two calculation methods does not change immediately, and the change has an obvious "time lag" phenomenon. In the aspect of calculation time, a computer configured as an Intel i 5-84002.8 GHz, 16GB memory is used, the calculation time of the implicit differential model is 27.36 seconds, while the calculation time of the time domain two-port model proposed herein is 1.07 seconds, and the speed is increased by more than 25 times, which embodies the advantage of the method in the aspect of calculation efficiency.
2) The validity of the natural gas pipeline time domain two-port model constructed in the text is verified by adopting a 4-node natural gas system example, and meanwhile, the calculation time is consumed by comparing with an implicit differential model, as shown in fig. 9, a schematic diagram of the topology and pipeline parameters of the natural gas testing system of the example is shown, wherein a node 1 is a gas source (GW1), the gas pressure is constant at 0.4MPa, and nodes 3 and 4 are gas loads (GL1 and GL 2).
Fig. 10 is a schematic diagram showing dynamic curves of pressure at each node and mass flow at the head and tail ends of each pipeline in the natural gas system according to the example of the present invention, where the boundary conditions are Π 1, Md1, and Md2, which respectively represent pressure at node 1, load Md1, and Md2 in the natural gas system. II 2-EDM, II 3-EDM and II 4-EDM respectively represent the air pressure of the nodes 2, 3 and 4 calculated according to the EDM; min1-EDM, Min2-EDM, Min3-EDM represent the mass flow into pipes 1#, 2#, 3# calculated according to EDM, respectively. By taking IDM as comparison, the maximum deviation of the node air pressure and the pipeline mass flow obtained by EDM solution with the same differential precision is 0.29 percent and 0.47 percent respectively. It should be noted that the two calculation results are not identical because the difference formats used by the two methods are different, and in order to eliminate the influence of the boundary condition values on the calculation speed, several groups of loads are randomly generated, and the two methods are respectively used for calculating the natural gas system load flow, and the lengths of the limited time and space elements are 1 second and 362 meters respectively. The average time consumption of the method is 4.21 seconds, and the average time consumption of the implicit differential model is 131.38 seconds, which represents the great advantage of the technology in the aspect of solving efficiency.
3) Dynamic energy flow calculation example for regional electricity-gas comprehensive energy system
The topology of the regional electricity-gas comprehensive energy system is shown in fig. 6, and the regional electricity-gas comprehensive energy system is composed of a 13-node natural gas system and a 12-node electric power system and comprises three generators; a traditional thermal power generating unit. The two gas turbine units are respectively connected with natural gas system nodes 8 and 9; a gas well; the system comprises ten electric loads; three conventional gas loads; two gas turbine gas loads. The time interval length of the power system model is selected to be 15 minutes, energy flow calculation results of the regional electricity-gas comprehensive energy system within 6 hours are considered, the time interval and the space interval are 24 time intervals, the length of a limited time infinitesimal and the length of a limited space infinitesimal in the natural gas system model are respectively 9 seconds and 3303 meters, 40000 variables are obtained, 56450 constraint conditions are obtained, and the model solving time is 8.66 seconds.
Fig. 11 is a graph showing the comparison between the load of the power system and the output of the generator according to the example of the present invention, fig. 12 is a graph showing the comparison between the gas pressure of the natural gas system and the gas consumption of the gas turbine according to the example of the present invention, and fig. 11 and 12 show that: the method provided by the embodiment of the invention can rapidly depict the dynamic energy flow of the electricity-gas integrated energy system, and the dynamic model can reflect the dynamic process of the air pressure in the natural gas system and can better reflect the real condition of the operation of the combined system.
It is noted that those skilled in the art will recognize that embodiments of the present invention are not described in detail herein.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (5)

1. An explicit differential-based method for acquiring dynamic power flow of an electric-gas integrated energy system, the method comprising:
step 1, obtaining linear relations between mass flow and pressure of any node in different areas in a natural gas pipeline differential grid, boundary conditions and initial conditions according to a differential equation of the mass flow and the pressure;
step 2, constructing a two-port model of mass flow and pressure of any node on the natural gas pipeline differential grid with respect to boundary conditions and initial conditions according to the linear relation obtained in the step 1;
and 3, embedding the two-port model into a dynamic energy flow calculation model of the electric-gas comprehensive energy system to obtain the dynamic energy flow of the electric-gas comprehensive energy system.
2. The method for acquiring the dynamic power flow of the electricity-gas integrated energy system based on the explicit difference as claimed in claim 1, wherein the process of the step 1 is specifically as follows:
first, the second-order hyperbolic partial differential equation of mass flow is expressed as:
Figure FDA0003123813310000011
the second order hyperbolic partial differential equation for pressure is expressed as:
Figure FDA0003123813310000012
wherein M is natural gas mass flow, pi is natural gas pressure, c is gas sound velocity, lambda is friction coefficient, D is pipeline diameter, x is space variable, t is time variable,
Figure FDA0003123813310000013
in order to be the average flow rate,
Figure FDA0003123813310000014
is the sign of the partial derivative;
and then regionalizing the solution into discrete points according to the natural gas pipeline differential grid, wherein:
Figure FDA0003123813310000015
representing the mass flow M at node (j, k)(x＝j,t＝k)；
Discretizing the differential equation into a differential equation, approximating a second order partial derivative using a second order center difference quotient, expressed as:
Figure FDA0003123813310000016
Figure FDA0003123813310000017
wherein tau is a time infinitesimal; h is a spatial infinitesimal; u (x)j,tk) Representing the pressure or mass flow at the point with the space serial number j (j is more than or equal to 0 and less than or equal to N) and the time serial number k (T is more than or equal to 0 and less than or equal to T) in the differential grid;
the first order partial derivative of approximate time using the forward difference quotient is expressed as:
Figure FDA0003123813310000021
substituting equations (3), (4) and (5) into the second-order hyperbolic partial differential equation (1) of the mass flow rate, the algebraic form of the difference with respect to the mass flow rate is obtained as:
Figure FDA0003123813310000022
reusing spatial direction transformations
Figure FDA0003123813310000023
Combined vertical type (5) and (6) are obtained
Figure FDA0003123813310000024
In the formula, K1、K2、K3、K4Is a difference coefficient whose value is:
Figure FDA0003123813310000025
by giving boundary conditions and initial conditions, solving a differential expression (8) of a second-order hyperbolic partial differential equation about mass flow, and for a section of natural gas pipeline, assuming that the head-end air pressure and the tail-end mass flow are known, the following steps are provided:
Figure FDA0003123813310000026
in the formula (I), the compound is shown in the specification,
Figure FDA0003123813310000027
indicating the gas pressure at the inlet (x ═ N) of the duct(x＝N,t)；Πin(t) is a function of the duct inlet pressure over time;
Figure FDA0003123813310000028
denotes the mass flow M at the outlet of the pipe (x ═ 0)(x＝0,t)；Mout(t) is the initial value of the mass flow at the outlet of the pipeline;
assuming that the system is in a steady state at the initial moment, the mass flow is satisfied
Figure FDA0003123813310000029
Considering that the state quantity of equation (8) contains only mass flow, the pressure boundary condition of equation (9) needs to be converted into a mass flow boundary condition, for which the first-order partial derivative of the forward difference quotient approximation space is used:
Figure FDA00031238133100000210
the united type (1), (5) and (11) equivalently converts the pressure boundary condition (9) into a mass flow boundary condition:
Figure FDA00031238133100000211
in the formula (I), the compound is shown in the specification,
Figure FDA00031238133100000212
indicating the gas pressure pi at the node (j, k)(x＝j,t＝k)；
The boundary conditions of the three-layer explicit difference equation for the mass flow are thus obtained, namely equations (9), (12); and initial conditions, i.e. formulae (9), (10);
in a similar way, firstly, according to a second-order hyperbolic partial differential equation (2) of the pressure, a differential algebraic form of the pressure is obtained in a joint type (3), (4) and (5):
Figure FDA0003123813310000031
secondly, the initial conditions of the formula (13) are the air pressure of each node at the initial time and the air pressure change rate at the initial time;
the joint type (5) and (11) obtain the air pressure of each node at the initial moment as follows:
Figure FDA0003123813310000032
assuming that the system is in a steady state at the initial time, the air pressure change rate at the initial time is:
Figure FDA0003123813310000033
finally, considering that the state quantity of equation (13) contains only pressure, it is necessary to convert the mass flow boundary into a pressure boundary, which is satisfied at the pipe outlet according to equation (9):
Figure FDA0003123813310000034
boundary conditions for three-layer explicit differential equations for pressure are thus obtained, namely equations (9), (16); and the initial conditions, i.e., equations (14), (15).
3. The method for acquiring the dynamic power flow of the electricity-gas integrated energy system based on the explicit difference as claimed in claim 1, wherein in step 2, the region to be solved is divided into four regions according to the difference between the point to be solved and the boundary point and the initial point on the natural gas pipeline differential grid, and the four regions are respectively:
1) region 1(k ═ 1, j < N); 2) region 2(k ═ 1, j ═ N); 3) region 3 (k is more than or equal to 2 and less than or equal to T, j is more than or equal to 1 and less than N); 4) region 4(2 ≦ k ≦ T, j ═ N);
firstly, a two-port model of mass flow of any node on a natural gas pipeline differential grid with respect to boundary conditions and initial conditions is constructed, specifically:
(1) region 1: when k is 1 and j is less than N, the difference form of the mass flow is as follows:
Figure FDA0003123813310000035
assuming that the mass flow of each node of the pipeline at the initial moment is equal to the mass flow at the tail end of the pipeline, when
Figure FDA0003123813310000036
And then, satisfy:
Figure FDA0003123813310000041
bringing formula (18) into formula (17) yields:
Figure FDA0003123813310000042
in the formula (I), the compound is shown in the specification,
Figure FDA0003123813310000043
to represent
Figure FDA0003123813310000044
And
Figure FDA0003123813310000045
Figure FDA0003123813310000046
the linear correlation coefficient of (a) has a value of:
Figure FDA0003123813310000047
(2) region 2: when k is 1 and j is N, the difference form of the mass flow is specifically:
Figure FDA0003123813310000048
in the formula:
Figure FDA0003123813310000049
when the region 1 and the region 2 are combined into the region a and the equations (19) and (20) are observed, when k is 1 and j is not less than 1 and not more than N, that is, the region 1+ the region 2, the mass flow satisfies:
Figure FDA00031238133100000410
the matrix form of the formula (21) is
Figure FDA00031238133100000411
In the formula:
Figure FDA00031238133100000412
Figure FDA00031238133100000413
(3) region 3: when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than N, the difference form of the mass flow satisfies the following conditions:
Figure FDA00031238133100000414
(4) region 4: when k is more than or equal to 2 and less than or equal to T, j is equal to N, the difference form of the mass flow satisfies:
Figure FDA00031238133100000415
combining the area 3 and the area 4 into an area B, and observing the formula (23) and the formula (24), when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than or equal to N, namely the area 3+ the area 4, the mass flow satisfies:
Figure FDA00031238133100000416
so far, a two-port form of mass flow rates of the area A and the area B is obtained, the area A and the area B are combined, and when k is more than or equal to 1 and less than or equal to T and j is more than or equal to 1 and less than or equal to N, the mass flow rates satisfy:
Figure FDA0003123813310000051
in the formula, Aj,BjDifference coefficient matrices, m, of (T × T) respectivelyjThe mass flow for a node with spatial index j is a (T × 1) matrix:
Figure FDA0003123813310000052
equation (26) is a two-port model of the mass flow of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition;
then, a two-port model of the pressure of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition is constructed, specifically:
(1) region 1: when k is 1 and j < N, the pressure difference satisfies the following conditions:
Figure FDA0003123813310000053
as can be seen from equation (27) and the air pressure at each node at the initial time of synchronization, when k is 0 and j is equal to or less than 1 and equal to or less than N, the air pressure satisfies:
Figure FDA0003123813310000054
in the formula
Figure FDA0003123813310000055
The rate of change of the gas pressure at the initial time and the boundary conditions at the outlet of the pipeline are brought into (27) to obtain:
Figure FDA0003123813310000056
the formula (29) is arranged as
Figure FDA0003123813310000057
Wherein:
Figure FDA0003123813310000061
the matrix form of equation (30) is:
Figure FDA0003123813310000062
in the formula:
Figure FDA0003123813310000063
Figure FDA0003123813310000064
(2) region 2: when k is 1 and j is N, the pressure satisfies:
Figure FDA0003123813310000065
in the formula:
Figure FDA0003123813310000066
(3) region 3: when k is more than or equal to 2 and less than or equal to T and j is more than or equal to 1 and less than N, the pressure satisfies the following conditions:
Figure FDA0003123813310000067
(4) region 4: when k is more than or equal to 2 and less than or equal to T, j is equal to N, the pressure satisfies the following conditions:
Figure FDA0003123813310000068
combining 4 areas, and when k is more than or equal to 1 and less than or equal to T and j is more than or equal to 1 and less than or equal to N, the pressure is satisfied as follows:
Figure FDA0003123813310000069
in the formula, Cj,DjDifference coefficient matrices, n, of (T x T) respectivelyjThe pressure at a node with spatial index j is a (T1) matrix:
Figure FDA00031238133100000610
thus, a two-port model of the pressure of any node on the natural gas pipeline differential grid with respect to the boundary condition and the initial condition is obtained, namely equation (35).
4. The method for acquiring dynamic power flow of an electric-gas integrated energy system based on explicit difference as claimed in claim 3, wherein in step 2, the two-port model of mass flow and pressure of any node on the natural gas pipeline difference grid with respect to boundary condition and initial condition is further represented by:
Figure FDA0003123813310000071
when the initial conditions are head end air pressure and tail end mass flow, the head end mass flow and the tail end air pressure of the pipeline meet the following conditions:
Figure FDA0003123813310000072
thereby obtaining a time domain two-port model (37) of the gas pressure and the mass flow at the head end and the tail end of the natural gas pipeline.
5. The method for acquiring the dynamic energy flow of the electric-gas integrated energy system based on the explicit difference as claimed in claim 4, wherein in step 3, the operation cost of the regional electric-gas integrated energy system includes the conventional unit power generation cost and the natural gas source gas supply cost, and the objective function of the dynamic energy flow is to minimize the operation cost, that is:
Figure FDA0003123813310000073
in the formula (I), the compound is shown in the specification,
Figure FDA0003123813310000074
calculating a set of time periods for the powerflow;
Figure FDA0003123813310000075
is a traditional unit set; c. CgcIs the g thcCost coefficient of the conventional unit;
Figure FDA00031238133100000713
is the g thcThe conventional unit is in time period teInternal output force;
Figure FDA0003123813310000076
calculating the number of time segments for the natural gas energy flow;
Figure FDA0003123813310000077
is a medium-pressure air source set; c. CwThe cost coefficient of the w gas source;
Figure FDA0003123813310000078
for the w gas source during the time period tgMass flow rate inside;
the energy flow model of the natural gas system satisfies the following constraints:
(1) pipeline air pressure-mass flow constraint:
according to the time domain two-port model of the gas pressure and the mass flow at the head end and the tail end of the natural gas pipeline, the gas pressure and the mass flow at the head end and the tail end of the p-th pipeline meet the following conditions:
Figure FDA0003123813310000079
in the formula (I), the compound is shown in the specification,
Figure FDA00031238133100000710
respectively representing mass flow and air pressure matrixes at the head end of the p-th pipeline in all time periods;
Figure FDA00031238133100000711
respectively representing the mass flow and the air pressure matrix of the p-th pipeline tail end in all time periods;
(2) air pressure and mass flow restraint of an air source:
Figure FDA00031238133100000712
(3) node mass flow balance constraint:
for natural gas node ngAnd the sum of the mass flow of all pipelines connected with the node and the net injection mass flow of the node is zero:
Figure FDA0003123813310000081
in the formula, ngIs a natural gas node;
Figure FDA0003123813310000082
for accessing natural gas node ngThe gas source set of (2);
Figure FDA0003123813310000083
is the end and natural gas node ngA set of connected pipes;
Figure FDA0003123813310000084
is a head end and natural gas node ngAre connected with each otherA set of pipes of (a); d is an element of Dg(ng) To be located at natural gas node ngA load set of (a);
(4) and (3) node air pressure constraint:
Figure FDA0003123813310000085
in the formula (I), the compound is shown in the specification,Π、
Figure FDA00031238133100000820
the lower limit and the upper limit of the natural gas node air pressure are set;
the power flow model of the power system satisfies the following constraints:
(1) node active and reactive power balance constraint:
Figure FDA0003123813310000087
Figure FDA0003123813310000088
in the formula, neIs a power system node;
Figure FDA0003123813310000089
for node n of the power systemeA collection of legacy units;
Figure FDA00031238133100000810
for node n of the power systemeA set of gas turbine units;
Figure FDA00031238133100000811
for node n of the power systemeA set of electrical loads;
Figure FDA00031238133100000812
as a lineFirst section and electric power system node neA set of connected power lines of (c);
Figure FDA00031238133100000813
for line end and power system node neA set of connected power lines of (c);
Figure FDA00031238133100000814
the active power and the reactive power of the traditional unit are obtained;
Figure FDA00031238133100000815
the active power and the reactive power of the thermal power generating unit are obtained;
Figure FDA00031238133100000816
active and reactive loads for the power system; plt、QltThe active power and the reactive power of the line are obtained; r isl、xlThe resistance and reactance of the circuit are shown; i isltIs the square of the branch current;
Figure FDA00031238133100000817
is the ground conductance of the node;
Figure FDA00031238133100000818
is the square of the node voltage;
(2) the line voltage drop equation and the line current equation of the second-order cone-convex relaxation are as follows:
Figure FDA00031238133100000819
Figure FDA0003123813310000091
in the formula (I), the compound is shown in the specification,
Figure FDA00031238133100000910
is the line end voltage;
Figure FDA0003123813310000099
is the line first section voltage;
(3) and (3) generator capacity constraint:
the active power and reactive power upper and lower limits of the generator are constrained as follows:
Figure FDA0003123813310000092
in the formula (I), the compound is shown in the specification,
Figure FDA0003123813310000093
for node n of the power systemeA set of generator sets;
Figure FDA0003123813310000094
the upper limit and the lower limit of active power of the generator are set;
Figure FDA0003123813310000095
the upper limit and the lower limit of the reactive power of the generator are set;
the constraints after the electrical-to-electrical coupling are:
(1) the active power and gas consumption of gas power generation are restricted:
Figure FDA0003123813310000096
in the formula, eta is energy conversion efficiency;
(2) the time coupling constraint of the power energy flow calculation variable and the natural gas energy flow calculation variable is as follows:
Figure FDA0003123813310000097
in the formula (I), the compound is shown in the specification,
Figure FDA0003123813310000098
is an upper rounding function; and N is the ratio of the time resolution of the electric and gas decision variables.
CN202110684263.9A 2021-06-21 2021-06-21 Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference Active CN113343531B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110684263.9A CN113343531B (en) 2021-06-21 2021-06-21 Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110684263.9A CN113343531B (en) 2021-06-21 2021-06-21 Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference

Publications (2)

Publication Number Publication Date
CN113343531A true CN113343531A (en) 2021-09-03
CN113343531B CN113343531B (en) 2023-04-18

Family

ID=77477957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110684263.9A Active CN113343531B (en) 2021-06-21 2021-06-21 Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference

Country Status (1)

Country Link
CN (1) CN113343531B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114065668A (en) * 2021-11-25 2022-02-18 重庆大学 Graph theory-based quantitative calculation method for water distribution system flow and water head pressure along line
CN114285037A (en) * 2022-03-07 2022-04-05 天津大学 Method for determining control parameter stability region of regional electricity-gas integrated energy system
CN114841075A (en) * 2022-05-19 2022-08-02 华北电力大学 Optimal energy flow calculation modeling method for electricity-gas comprehensive energy system
CN115169950A (en) * 2022-07-26 2022-10-11 山东大学 Distributed cooperation method and system for electric-gas system based on multi-parameter planning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150347651A1 (en) * 2014-06-02 2015-12-03 University Of Windsor System and Method for Determining Heat and Fluid Flow in or Around Objects
CN107291990A (en) * 2017-05-24 2017-10-24 河海大学 Energy stream emulation mode based on electrical interconnection integrated energy system transient Model
CN109740242A (en) * 2018-12-29 2019-05-10 重庆大学 Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method
CN111695269A (en) * 2020-06-23 2020-09-22 中国电力科学研究院有限公司 Multi-time-interval electricity-gas comprehensive energy system state estimation method, system and device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150347651A1 (en) * 2014-06-02 2015-12-03 University Of Windsor System and Method for Determining Heat and Fluid Flow in or Around Objects
CN107291990A (en) * 2017-05-24 2017-10-24 河海大学 Energy stream emulation mode based on electrical interconnection integrated energy system transient Model
CN109740242A (en) * 2018-12-29 2019-05-10 重庆大学 Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method
CN111695269A (en) * 2020-06-23 2020-09-22 中国电力科学研究院有限公司 Multi-time-interval electricity-gas comprehensive energy system state estimation method, system and device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵霞等: "电-气综合能源系统能流计算的改进方法", 《电工技术学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114065668A (en) * 2021-11-25 2022-02-18 重庆大学 Graph theory-based quantitative calculation method for water distribution system flow and water head pressure along line
CN114065668B (en) * 2021-11-25 2024-04-05 重庆大学 Quantitative calculation method for flow and water head pressure along water distribution system based on graph theory
CN114285037A (en) * 2022-03-07 2022-04-05 天津大学 Method for determining control parameter stability region of regional electricity-gas integrated energy system
CN114841075A (en) * 2022-05-19 2022-08-02 华北电力大学 Optimal energy flow calculation modeling method for electricity-gas comprehensive energy system
CN115169950A (en) * 2022-07-26 2022-10-11 山东大学 Distributed cooperation method and system for electric-gas system based on multi-parameter planning
CN115169950B (en) * 2022-07-26 2023-03-24 山东大学 Distributed cooperation method and system for electric-gas system based on multi-parameter planning

Also Published As

Publication number Publication date
CN113343531B (en) 2023-04-18

Similar Documents

Publication Publication Date Title
CN113343531B (en) Method for acquiring dynamic energy flow of electricity-gas integrated energy system based on explicit difference
CN107291990B (en) Energy flow simulation method based on transient model of electricity-gas interconnection comprehensive energy system
CN102012956B (en) A kind of wind energy turbine set equivalence method considering the random fluctuation of wind energy turbine set input wind speed and direction
CN106096269B (en) The Interval Power Flow calculation method of natural gas grid in a kind of electric-gas coupled system
CN110797896A (en) Hybrid energy storage capacity optimal configuration method in micro-grid wind and solar energy storage system
CN109636009B (en) Method and system for establishing neural network model for determining line loss of power grid
CN110429648B (en) Small interference stability margin probability evaluation method considering wind speed random fluctuation
CN105930980B (en) A kind of electricity turns the integrated energy system linear multi likelihood energy stream method of gas
CN110866213A (en) Multi-network steady-state energy flow analysis method and device for electricity-gas integrated energy system
CN109697308A (en) A kind of natural gas transmission systematic steady state modeling method considering pipe network time-delay characteristics
CN105162116B (en) A kind of section economic load dispatching Nonlinear Dual optimization method of the solution containing wind-powered electricity generation
CN113627021A (en) Electrical interconnection system optimal energy flow calculation method based on sequence convex programming
CN114298455A (en) Comprehensive energy hybrid modeling method based on GAN technology
CN111563315B (en) Topology analysis-based steady-state energy flow calculation method for electric-gas comprehensive energy system
CN104993503A (en) Island microgrid frequency control method
CN111797537A (en) Reliability evaluation method and device for power distribution system containing wind power
Sedighnejad et al. Performance evaluation of a hybrid wind-diesel-compressed air energy storage system
CN108964143B (en) Natural gas network static equivalent model of electricity-gas comprehensive energy system
CN113065981B (en) Source-storage-network coordination planning method considering energy storage comprehensive income
CN114139349A (en) Combined optimization method of day-ahead units of gas-electricity combined system considering pipeline climbing
CN114841075B (en) Optimal energy flow calculation modeling method of electricity-gas integrated energy system
Tang et al. Reliability evaluation of EGIES with highly nonlinear modeling of gas system components
CN109980649A (en) It is a kind of meter and multiple stable point the probability load flow calculation method based on saddle point approximation method
CN110707705B (en) Power flow sequence analysis model of electric-gas coupling comprehensive energy system
CN117094613B (en) Model construction method and device applied to comprehensive energy system

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