CN113486532A - Dynamic safety control method for electric heating comprehensive energy system - Google Patents
Dynamic safety control method for electric heating comprehensive energy system Download PDFInfo
- Publication number
- CN113486532A CN113486532A CN202110828652.4A CN202110828652A CN113486532A CN 113486532 A CN113486532 A CN 113486532A CN 202110828652 A CN202110828652 A CN 202110828652A CN 113486532 A CN113486532 A CN 113486532A
- Authority
- CN
- China
- Prior art keywords
- node
- formula
- temperature
- transmission
- power
- 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
- 238000005485 electric heating Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000005540 biological transmission Effects 0.000 claims abstract description 126
- 239000011159 matrix material Substances 0.000 claims abstract description 75
- 230000035945 sensitivity Effects 0.000 claims abstract description 58
- 230000008878 coupling Effects 0.000 claims abstract description 30
- 238000010168 coupling process Methods 0.000 claims abstract description 30
- 238000005859 coupling reaction Methods 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 238000013507 mapping Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 41
- 238000002347 injection Methods 0.000 claims description 29
- 239000007924 injection Substances 0.000 claims description 29
- 239000013598 vector Substances 0.000 claims description 28
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 12
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 239000008400 supply water Substances 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 2
- 238000011217 control strategy Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000003068 static effect Effects 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/20—Design optimisation, verification or simulation
-
- 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
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- 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
Abstract
The invention discloses a dynamic safety control method of an electric heating comprehensive energy system, which comprises the following steps: 10) deducing a transmission matrix of boundary conditions and initial conditions based on a heat supply network dynamic model of a central implicit differential format, and establishing a pipeline equivalent model; 20) based on a pipeline equivalent model, combining a network topological equation, and establishing a source-load equation of heat supply network dynamics; 30) based on a power grid load flow model and a heat supply network dynamic source-load equation, a sensitivity calculation formula of an electric power network and a heat supply network is established, a coupling equipment model is combined, a global mapping relation of an electric heating comprehensive energy system is established, and a dynamic safety control method is established.
Description
Technical Field
The invention relates to the field of energy system modeling and operation analysis, in particular to a dynamic safety control method for an electric heating comprehensive energy system.
Background
The increasing environmental problems and energy crisis have promoted the revolution of energy technology, and the electric-thermal comprehensive energy system is a main form of energy revolution, is helpful to improve the complementarity of electric energy and heat energy utilization, and is widely applied to practical engineering. However, due to the difference in the multi-time scale characteristics of the power system and the thermodynamic system, the model of the comprehensive electric heating energy system is essentially a complex set of high-dimensional constants and partial differential equations, and the solution of the complex set of constants and partial differential equations is very complex. In addition, since electric energy propagates in the power system at the speed of light, it tends to stabilize in milliseconds; the heat energy is transmitted in the thermodynamic system through working medium flow, and can be stabilized only after a control instruction is issued, so that the control time scale in the electric power system is generally short, and the control time scale of the thermodynamic system is generally long.
In the existing research, a safety control strategy is generally formulated according to an energy flow calculation result at a control moment, but because the control time interval of a thermodynamic system is long, the dynamic process evolution between the control time intervals is often neglected, so that the formulation of the control strategy is relatively one-sided. In addition, the existing research generally adopts a steady-state model for modeling analysis, neglects the dynamic and multi-time scale characteristics in the electric heating comprehensive energy system, and has larger difference between the obtained result and the real situation. In addition, due to the coupling of heating power and water power, the existing research generally analyzes the static safety problem of the electric heating comprehensive energy system through an optimization method or an alternate iteration solving method, lacks an accurate and quantitative safety control strategy, and adjusts the state overrun which may occur in the operation process, thereby ensuring the operation of the system.
Disclosure of Invention
In order to solve the defects mentioned in the background technology, the invention aims to provide a dynamic safety control method of an electric heating comprehensive energy system, which deduces the transmission coefficients of the initial condition and the boundary condition in a differential format heat supply network dynamic model, establishes a pipeline equivalent model of state quantity distribution, constructs a source-charge equation of heat supply network dynamics by combining a topological equation, reveals a source-charge mapping relation, and provides a dynamic safety control strategy through a global sensitivity calculation formula of the electric heating comprehensive energy system. Provides theoretical guidance for the safe operation of the electric heating comprehensive energy system.
The purpose of the invention can be realized by the following technical scheme:
a dynamic safety control method for an electric heating comprehensive energy system comprises the following steps:
step 10) deducing transmission coefficients of boundary conditions and initial conditions based on a heat supply network dynamic model of a central implicit differential format, and establishing a pipeline equivalent model;
step 20) establishing a dynamic source-load equation of the heat supply network based on the pipeline equivalent model and in combination with a network topology equation;
and step 30) establishing a sensitivity calculation formula of the electric power and thermal power network based on the power flow model and the dynamic source-load equation of the heat supply network, establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model, and establishing a dynamic safety control method.
Further, the step 10) is specifically as follows:
step 101) according to the central implicit difference format, establishing a heat supply network dynamic model as follows:
in the formula, i and j represent a space step and a time step of the temperature distribution of the pipeline, respectively, T is the temperature of the pipeline with the ambient temperature as a reference value,
represents the temperature, μ, of the segment at time j at the ith of the pipe1-3The method comprises the following steps of constructing a constant coefficient term for expressing temperature transmission characteristics, specifically:
wherein Δ x and Δ t are the space and time step of the difference; v is the flow velocity of the pipeline, and r is the thermal resistance coefficient of the pipeline;
step 102) expressing the initial conditions, the boundary conditions and the temperature of the end of the pipeline in a vector form, and introducing a permutation and combination calculation formula
Defining reduced coefficient terms simultaneously
Initial condition T0Boundary condition T0And pipe end temperature TMThe vector is shown in equation (5):
in the formula, M and N are the space and time segment number of the pipeline temperature respectively;
permutation and combination calculation formula
Expressed as:
step 103) arbitrary boundary conditions
And status point
There are two types of transmission paths:
(1)
first experiences a spatial transmission loss to
Then experiences M-1 section space transmission loss and k-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-1 section space transmission loss and k-j-1 section time transmission loss to
αkjThe transmission relationship of the boundary condition at the moment j to the branch terminal temperature at the moment k is shown as follows:
in the formula, the 1 st term and the 2 nd term on the right of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively,
and
denotes alphakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 104) any initial conditions
And status point
There are two types of transmission paths:
(1)
after a period of transmission loss to
Then experiences the M-j section space transmission loss and the N-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-j-1 space transmission loss and k-1 time transmission loss to
βkjAnd (3) representing the transmission relation of the j-th section initial condition to the branch terminal temperature at the k moment, and represented as:
in the formula, the 1 st and 2 nd terms on the right of the equal sign correspond to the transmission coefficients, ψ, of the path (1) and the path (2), respectivelykj1And psikj2Is represented by betakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 105) combining the transmission matrices of the initial conditions and the boundary conditions, the pipe end temperature at a single moment is expressed as:
expanding the formula (11) from a single moment k to a full period, wherein the temperature vector at the tail end of the pipeline is expressed as a linear combination of an initial condition and a boundary condition, and is expressed as a formula (12):
TM＝αT0+βT0 (12)
where α and β are the transmission matrices for the boundary and initial conditions, respectively.
Further, the step 20) is specifically as follows:
step 201) constructing a multi-type network topological equation comprising a pipeline initial end temperature-node temperature incidence matrix ApnPipeline flow-inflow node incidence matrix AinPipeline flow-outflow node incidence matrix Aout(ii) a The block matrix is an N-dimensional square matrix, NbAnd NhThe number of the pipelines and the number of the nodes of the heat supply network are respectively, and each topological equation is specifically expressed as follows:
(1) pipeline initial end temperature-node temperature correlation matrix Apn：ApnContaining Nb×NhA square matrix, if the initial end temperature of the pipeline i is the same as the node j, ApnThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Ain：AinContaining Nh×NbA square matrix, if the flow flows into the node i from the pipeline j, AinThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(3) pipe flow-inflow node incidence matrix Aout：AoutContaining Nh×NbA square matrix, if the flow flows into the pipeline j from the node i, AoutThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
step 202) combining a pipeline equivalent model to deduce a dynamic source-load equation of the heat supply network; based on pipeline initial end temperature vector T0And node temperature TnBy means of the correlation matrix ApnThe relationship, expressed as:
T0＝ApnTn (13)
the tube temperature mixing equation is expressed as:
in the formula, Tn,iDenotes the temperature, T, of node iM,bDenotes the end temperature, m, of the pipe bbDenotes the mass flow of the pipe b, Es,iRepresenting a set of pipes with node i as the head node, Ee,iRepresenting a set of pipes, V, with pipe i as the last nodehThe node set node outflow flow of the heat supply network comprises two parts of node outflow to load and node outflow to pipeline, and the formula (14) is popularized to the system by using a topological equation and expressed as follows:
AoutmT0+dTn＝AinmTM (15)
in the formula, m is a flow vector in the dynamic model, and d is a node injection flow vector;
bringing formula (12) and formula (13) into formula (15) in this order, yields:
MoutTn＝JTn+B (16)
in the formula, MoutThe method comprises the following steps of A, taking an outflow flow matrix of nodes in a dynamic model, J taking a transmission matrix of different node temperatures in the dynamic model, and B taking a temperature component of an initial condition action; mouJ and B are respectively:
Mout＝diag(AoutmApn+d),J＝AinmαApn,B＝AinmβΤ0 (17)
Vh,srand Vh,nsRespectively, the source node and the non-source node of the thermodynamic system, then equation (17) is expanded as:
in the formula (I), the compound is shown in the specification,
and
outgoing traffic matrices, T, being a set of source nodes and a set of non-source nodes, respectivelyn,srAnd Tn,nsTemperature vectors for source and non-source nodes, J11、J12、J21And J22Respectively, a block matrix, B, in a coefficient matrix JsrAnd BnsTemperature components acting on the source node and the non-source node for initial conditions;
according to a Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network is obtained and expressed as:
further, the step 30) is specifically as follows:
step 301), establishing a grid sensitivity calculation formula, including node voltage amplitude, phase angle sensitivity about node injection active power and reactive power, and line transmission power sensitivity to node injection active power and reactive power, as shown in formula (20) and formula (21), respectively:
in the formula, U and theta are the node voltage amplitude and phase angle respectively, P and Q are the node injected active power and reactive power respectively, and PlAnd QlActive and reactive power, S, respectively, for the transmission of network branchesUPSensitivity, S, representing the magnitude of the node voltage with respect to the active power injected into the nodeθPSensitivity, S, representing the phase angle of the node voltage with respect to the active power injected into the nodeUQIndicating the sensitivity of the node voltage amplitude with respect to the node injected reactive power, SθQRepresenting the sensitivity of the node voltage phase angle to the node injected reactive power;
step 302) establishing a heat supply network sensitivity calculation formula according to a source-load equation in the heat supply network dynamic model, wherein the heat supply network sensitivity calculation formula comprises the sensitivity of the water supply temperature at a source node relative to the water supply temperature at a non-source node, the sensitivity of the water supply temperature at a node relative to the return water temperature at the node, and the sensitivity of the water supply temperature at the node relative to the thermal power of the node, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the formula (I), the compound is shown in the specification,
indicating the temperature of the water supply at node i at time k,
represents the return water temperature of the node i at the moment k,
represents the thermal power of node i at time k, diInjection flow for node i, CρIs specific heat capacity of working medium, Vh,ldFor a set of load nodes in a heat network, K1,ji,nkIs K1Elements of the jth row and the ith column in the block matrix of the jth row and the ith column;
step 303) establishing a global mapping relation of the electric heating comprehensive energy system by combining a coupling equipment model, wherein the coupling equipment model is expressed as:
φCHP＝Ch-cuPCHP (25)
φHP＝ckPHP (26)
in the formula, phiCHPAnd PCHPRespectively representing the thermal power and the electric power output by the cogeneration plant, ChTo the maximum output thermal power of the cogeneration plant, cuIs the heat-to-power ratio of a cogeneration plantHPAnd PHPRespectively representing the thermal power output and the electric power consumed by the thermoelectric conversion device, ckIs the thermoelectric ratio of the thermoelectric conversion device;
according to the formula (22) and the formula (24), obtaining the sensitivity of the thermal power of the source node in the heat supply network relative to the thermal power of the load node; and (3) obtaining the sensitivities of the injection power of the node in the power grid, the voltage amplitude and the branch transmission power relative to the thermal power of the load node in the heat supply network by combining a coupling equipment model, wherein the sensitivities are respectively expressed as a formula (28) and a formula (29):
in the formula (I), the compound is shown in the specification,
representing the injected active power of node i at time k,
representing the transmission power, V, of the branch at time kcpRepresenting a coupling node set in the electric heating integrated energy system, eta represents the conversion coefficient of the coupling equipment, and eta is-1/c for the electric heating cogeneration equipmentu(ii) a For a thermoelectric conversion device, η is 1/ck；
Step 304) according to the coupling relation of the electric heating system, a dynamic safety control method is established, which comprises the steps of adjusting the node injection power to control the node voltage amplitude and the branch transmission active power, adjusting the source node water supply temperature to control the load node water supply temperature, adjusting the source node water supply temperature to control the node voltage amplitude and the branch transmission active power, and respectively entering (30) to (33):
ΔP＝(U-Ulim+dU)/SUP (30)
wherein, Δ P represents the control quantity of the node injection power,
control quantity, U, representing supply water temperature at node i at time klimAnd Pl,limRepresenting node voltage amplitude threshold and branch transmission active power threshold, dU and dPlIs the margin of the node voltage amplitude.
The invention has the beneficial effects that:
the method establishes the analytical relationship between the dynamic characteristic of the heat supply network state quantity and the initial condition and the boundary condition, deduces a source-load direct mapping function, quantitatively describes the influence of the dynamic characteristic of the heat supply network on the operation of the power grid, is favorable for accurately formulating the dynamic safety control strategy of the electric heating comprehensive energy system considering the dynamic characteristic of the heat supply network, and ensures the stable operation of the system.
Drawings
The invention will be further described with reference to the accompanying drawings.
FIG. 1 is a schematic flow chart of a dynamic safety control method of an electric heating comprehensive energy system according to an embodiment of the invention;
FIG. 2 is a block diagram of a thermodynamic system employed in an embodiment of the present invention;
FIG. 3 is a comparison of the pipe end temperatures before and after control in an embodiment of the present invention;
fig. 4 is a comparison of active power injected into nodes before and after control in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be 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 of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1: a dynamic safety control method for an electric heating comprehensive energy system comprises the following steps:
step 10) deducing transmission coefficients of boundary conditions and initial conditions based on a heat supply network dynamic model of a central implicit differential format, and establishing a pipeline equivalent model;
step 20) establishing a dynamic source-load equation of the heat supply network based on the pipeline equivalent model and in combination with a network topology equation;
and step 30) establishing a sensitivity calculation formula of the electric power and thermal power network based on the power flow model and the dynamic source-load equation of the heat supply network, establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model, and establishing a dynamic safety control strategy.
The step 10) specifically comprises the following steps:
step 101) according to the central implicit difference format, establishing a heat supply network dynamic model as follows:
in the formula, i and j represent a space step and a time step of the temperature distribution of the pipeline, respectively, T is the temperature of the pipeline with the ambient temperature as a reference value,
represents the temperature, μ, of the segment at time j at the ith of the pipe1-3The method comprises the following steps of constructing a constant coefficient term for expressing temperature transmission characteristics, wherein the constant coefficient term specifically comprises the following steps:
wherein Δ x and Δ t are the space and time step of the difference; v is the flow velocity of the pipeline, and r is the thermal resistance coefficient of the pipeline;
step 102) representing the initial conditions, the boundary conditions and the pipeline end temperature into a vector form, and introducing a permutation and combination calculation formula
Defining reduced coefficient terms simultaneously
Initial condition T0Boundary condition T0And pipe end temperature TMThe vector is shown in equation (5). In the formula, M and N are the space and time segment numbers of the pipeline temperature respectively.
Permutation and combination calculation formula
Can be expressed as:
step 103) arbitrary boundary conditions
And status point
There are two types of transmission paths: (1)
first experiences a spatial transmission loss to
Then experiences M-1 section space transmission loss and k-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Experiences the spatial transmission loss of the M-1 sectionAnd k-j-1 time transmission loss to
Definition of alphakjThe transmission relationship of the boundary condition at the time j to the branch end temperature at the time k can be represented as follows:
in the formula, the 1 st term and the 2 nd term on the right of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively,
and
denotes alphakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 104) any initial conditions
And status point
There are two types of transmission paths: (1)
after a period of transmission loss to
Then experiences the M-j section space transmission loss and the N-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-j-1 space transmission loss and k-1 time transmission loss to
Definition of betakjThe transmission relation of the j-th section initial condition to the branch end temperature at the k moment can be represented as follows:
in the formula, the 1 st and 2 nd terms on the right of the equal sign correspond to the transmission coefficients, ψ, of the path (1) and the path (2), respectivelykj1And psikj2Is represented by betakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 105) combining the transmission matrices of the initial conditions and the boundary conditions, the pipe end temperature at a single moment can be expressed as:
expanding the formula (11) from a single time k to a full time, the pipeline end temperature vector can be expressed as a linear combination of an initial condition and a boundary condition, and can be expressed as a formula (12), wherein alpha and beta are transmission matrixes of the boundary and the initial condition respectively.
TM＝αT0+βT0 (12)
The step 20) comprises the following steps:
step 201) constructing a multi-type network topological equation comprising a pipeline initial end temperature-node temperature incidence matrix ApnPipeline flow-inflow node incidence matrix AinPipeline flow-outflow node incidence matrix Aout. Defining a blocking matrix as an N-dimensional square matrix, NbAnd NhThe number of pipes and the number of nodes of the heat supply network are respectively, and each topological equation can be specifically expressed as follows:
(1) pipeline initial end temperature-node temperature correlation matrix Apn：ApnContaining Nb×NhA square matrix, if the initial end temperature of the pipeline i is the same as the node j, ApnThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Ain：AinContaining Nh×NbA square matrix, if the flow flows from the pipeline j to the node i, AinThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Aout：AoutContaining Nh×NbA square matrix, if the flow flows into the pipeline j from the node i, AoutThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
step 202) combining the pipeline equivalent model to deduce a dynamic source-load equation of the heat supply network. Based on pipeline initial end temperature vector T0And node temperature TnCan pass through the incidence matrix ApnThe relationship, expressed as:
T0＝ApnTn (13)
the pipe temperature mixing equation can be expressed as:
in the formula, Tn,iDenotes the temperature, T, of node iM,bIndicating the end temperature of the pipe b，mbDenotes the mass flow of the pipe b, Es,iRepresenting a set of pipes with node i as the head node, Ee,iRepresenting a set of pipes, V, with pipe i as the last nodehIs a node assembly of a heat supply network. Considering that node outgoing traffic includes two parts, node outgoing to load and node outgoing to pipe, equation (14) can be generalized to the system using topological equations, which are expressed as:
AoutmT0+dTn＝AinmTM (15)
in the formula, m is a flow vector in the dynamic model, and d is a node injection flow vector. By bringing formula (12) and formula (13) into formula (15) in this order, it is possible to obtain:
MoutTn＝JTn+B (16)
in the formula, MoutThe matrix is an outflow flow matrix of nodes in the dynamic model, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by an initial condition. MouJ and B are respectively:
Mout＝diag(AoutmApn+d),J＝AinmαApn,B＝AinmβΤ0 (17)
definition Vh,srAnd Vh,nsRespectively, the source node and the non-source node of the thermodynamic system, equation (17) can be expanded as follows:
in the formula (I), the compound is shown in the specification,
and
outgoing traffic matrices, T, being a set of source nodes and a set of non-source nodes, respectivelyn,srAnd Tn,nsTemperature vectors for source and non-source nodes, J11、J12、J21And J22Respectively, a block matrix, B, in a coefficient matrix JsrAnd BnsThe temperature components acting on the source node and the non-source node for the initial condition.
According to the Gaussian elimination method, a source-load equation in the dynamic model of the heat supply network can be obtained and expressed as follows:
the step 30) includes:
step 301), establishing a grid sensitivity calculation formula, including node voltage amplitude, phase angle sensitivity about node injection active power and reactive power, and line transmission power sensitivity to node injection active power and reactive power, as shown in formula (20) and formula (21), respectively:
in the formula, U and theta are the node voltage amplitude and phase angle respectively, P and Q are the node injected active power and reactive power respectively, and PlAnd QlActive and reactive power, S, respectively, for the transmission of network branchesUPSensitivity, S, representing the magnitude of the node voltage with respect to the active power injected into the nodeθPSensitivity, S, representing the phase angle of the node voltage with respect to the active power injected into the nodeUQRepresenting the sensitivity of the node voltage amplitude to the node injected reactive power; sθQRepresenting the sensitivity of the node voltage phase angle to the node injected reactive power.
Step 302) establishing a heat supply network sensitivity calculation formula according to a source-load equation in the heat supply network dynamic model, wherein the heat supply network sensitivity calculation formula comprises the sensitivity of the water supply temperature at a source node relative to the water supply temperature at a non-source node, the sensitivity of the water supply temperature at a node relative to the return water temperature at the node, and the sensitivity of the water supply temperature at the node relative to the thermal power of the node, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the formula (I), the compound is shown in the specification,
indicating the temperature of the water supply at node i at time k,
represents the return water temperature of the node i at the moment k,
represents the thermal power of node i at time k, diInjection flow for node i, CρIs specific heat capacity of working medium, Vh,ldFor a set of load nodes in a heat network, K1,ji,nkIs K1The jth row and the ith column of the block matrix.
Step 303) establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model. The coupling device model may be expressed as:
φCHP＝Ch-cuPCHP (25)
φHP＝ckPHP (26)
in the formula, phiCHPAnd PCHPRespectively representing the thermal power and the electric power output by the cogeneration plant, ChTo the maximum output thermal power of the cogeneration plant, cuIs the heat-to-power ratio of a cogeneration plantHPAnd PHPRespectively representing the thermal power output and the electric power consumed by the thermoelectric conversion device, ckIs the thermoelectric ratio of the thermoelectric conversion device.
According to the formula (22) and the formula (24), the sensitivity of the thermal power of the source node in the heat supply network relative to the thermal power of the load node can be obtained. By combining the coupling equipment model, the sensitivities of the node injection power, the voltage amplitude and the branch transmission power in the power grid to the thermal power of the load node in the heat supply network can be obtained, and the sensitivities are respectively shown as a formula (28) and a formula (29).
In the formula (I), the compound is shown in the specification,
representing the injected active power of node i at time k,
representing the transmission power, V, of the branch at time kcpRepresenting a coupling node set in the electric heating integrated energy system, eta represents the conversion coefficient of the coupling equipment, and eta is-1/c for the electric heating cogeneration equipmentu(ii) a For a thermoelectric conversion device, η is 1/ck。
Step 304) constructing a dynamic safety control strategy according to the coupling relation of the electric heating system, wherein the dynamic safety control strategy comprises adjusting the voltage amplitude of a node injection power control node and the active power transmitted by a branch, adjusting the water supply temperature of a source node to control the water supply temperature of a load node, and adjusting the voltage amplitude of the node water supply temperature control node and the active power transmitted by the branch, which are respectively shown in formulas (30) to (33).
ΔP＝(U-Ulim+dU)/SUP (30)
Wherein, Δ P represents the control quantity of the node injection power,
control quantity, U, representing supply water temperature at node i at time klimAnd Pl,limRepresenting node voltage amplitude threshold and branch transmission active power threshold, dU and dPlIs the margin of the node voltage amplitude.
The application example is as follows:
the system shown in fig. 2 is taken as an example for explanation.
As shown in fig. 1, an embodiment of the present invention provides a dynamic safety control method for an electric heating integrated energy system, including the following steps:
step 10) deducing transmission coefficients of boundary conditions and initial conditions based on a heat supply network dynamic model of a central implicit differential format, and establishing a pipeline equivalent model;
step 20) establishing a dynamic source-load equation of the heat supply network based on the pipeline equivalent model and in combination with a network topology equation;
and step 30) establishing a sensitivity calculation formula of the electric power and thermal power network based on the power flow model and the dynamic source-load equation of the heat supply network, establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model, and establishing a dynamic safety control strategy.
In the above embodiment, the step 10) specifically includes:
step 101) according to the central implicit difference format, establishing a heat supply network dynamic model as follows:
in the formula, i and j represent a space step and a time step of the temperature distribution of the pipeline, respectively, T is the temperature of the pipeline with the ambient temperature as a reference value,
represents the temperature, μ, of the segment at time j at the ith of the pipe1-3The method comprises the following steps of constructing a constant coefficient term for expressing temperature transmission characteristics, wherein the constant coefficient term specifically comprises the following steps:
wherein Δ x and Δ t are the space and time step of the difference; v is the flow velocity of the pipeline, and r is the thermal resistance coefficient of the pipeline;
step 102) representing the initial conditions, the boundary conditions and the pipeline end temperature into a vector form, and introducing a permutation and combination calculation formula
Defining reduced coefficient terms simultaneously
Initial condition T0Boundary condition T0And pipe end temperature TMThe vector is shown in equation (5). In the formula, M and N are the space and time segment numbers of the pipeline temperature respectively.
Permutation and combination calculation formula
Can be expressed as:
step 103) arbitrary boundary conditions
And status point
There are two types of transmission paths: (1)
first experiences a spatial transmission loss to
Then experiences M-1 section space transmission loss and k-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-1 section space transmission loss and k-j-1 section time transmission loss to
Definition ofαkjThe transmission relationship of the boundary condition at the time j to the branch end temperature at the time k can be represented as follows:
in the formula, the 1 st term and the 2 nd term on the right of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively,
and
denotes alphakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 104) any initial conditions
And status point
There are two types of transmission paths: (1)
after a period of transmission loss to
Then experiences the M-j section space transmission loss and the N-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-j-1 space transmission loss and k-1 time transmission loss to
Definition of betakjThe transmission relation of the j-th section initial condition to the branch end temperature at the k moment can be represented as follows:
in the formula, the 1 st and 2 nd terms on the right of the equal sign correspond to the transmission coefficients, ψ, of the path (1) and the path (2), respectivelykj1And psikj2Is represented by betakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 105) combining the transmission matrices of the initial conditions and the boundary conditions, the pipe end temperature at a single moment can be expressed as:
expanding the formula (11) from a single time k to a full time, the pipeline end temperature vector can be expressed as a linear combination of an initial condition and a boundary condition, and can be expressed as a formula (12), wherein alpha and beta are transmission matrixes of the boundary and the initial condition respectively.
TM＝αT0+βT0 (12)
In the above embodiment, the step 20) specifically includes:
step 201) constructing a multi-type network topological equation comprising pipesTemperature-node temperature correlation matrix A at the beginning of the trackpnPipeline flow-inflow node incidence matrix AinPipeline flow-outflow node incidence matrix Aout. Defining a blocking matrix as an N-dimensional square matrix, NbAnd NhThe number of pipes and the number of nodes of the heat supply network are respectively, and each topological equation can be specifically expressed as follows:
(1) pipeline initial end temperature-node temperature correlation matrix Apn：ApnContaining Nb×NhA square matrix, if the initial end temperature of the pipeline i is the same as the node j, ApnThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Ain：AinContaining Nh×NbA square matrix, if the flow flows from the pipeline j to the node i, AinThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Aout：AoutContaining Nh×NbA square matrix, if the flow flows into the pipeline j from the node i, AoutThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
step 202) combining the pipeline equivalent model to deduce a dynamic source-load equation of the heat supply network. Based on pipeline initial end temperature vector T0And node temperature TnCan pass through the incidence matrix ApnThe relationship, expressed as:
T0＝ApnTn (13)
the pipe temperature mixing equation can be expressed as:
in the formula, Tn,iDenotes the temperature, T, of node iM,bDenotes the end temperature, m, of the pipe bbDenotes the mass flow of the pipe b, Es,iRepresenting a set of pipes with node i as the head node, Ee,iIs shown as a tubePipeline set with channel i as end node, VhIs a node assembly of a heat supply network. Considering that node outgoing traffic includes two parts, node outgoing to load and node outgoing to pipe, equation (14) can be generalized to the system using topological equations, which are expressed as:
AoutmT0+dTn＝AinmTM (15)
in the formula, m is a flow vector in the dynamic model, and d is a node injection flow vector. By bringing formula (12) and formula (13) into formula (15) in this order, it is possible to obtain:
MoutTn＝JTn+B (16)
in the formula, MoutThe matrix is an outflow flow matrix of nodes in the dynamic model, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by an initial condition. MouJ and B are respectively:
Mout＝diag(AoutmApn+d),J＝AinmαApn,B＝AinmβΤ0 (17)
definition Vh,srAnd Vh,nsRespectively, the source node and the non-source node of the thermodynamic system, equation (17) can be expanded as follows:
in the formula (I), the compound is shown in the specification,
and
outgoing traffic matrices, T, being a set of source nodes and a set of non-source nodes, respectivelyn,srAnd Tn,nsTemperature vectors for source and non-source nodes, J11、J12、J21And J22Respectively, a block matrix, B, in a coefficient matrix JsrAnd BnsActing on source and non-source nodes for initial conditionsA temperature component.
According to the Gaussian elimination method, a source-load equation in the dynamic model of the heat supply network can be obtained and expressed as follows:
in the above embodiment, the step 30) specifically includes:
step 301), establishing a grid sensitivity calculation formula, including node voltage amplitude, phase angle sensitivity about node injection active power and reactive power, and line transmission power sensitivity to node injection active power and reactive power, as shown in formula (20) and formula (21), respectively:
in the formula, U and theta are the node voltage amplitude and phase angle respectively, P and Q are the node injected active power and reactive power respectively, and PlAnd QlActive and reactive power, S, respectively, for the transmission of network branchesUPSensitivity, S, representing the magnitude of the node voltage with respect to the active power injected into the nodeθPSensitivity, S, representing the phase angle of the node voltage with respect to the active power injected into the nodeUQRepresenting the sensitivity of the node voltage amplitude to the node injected reactive power; sθQRepresenting the sensitivity of the node voltage phase angle to the node injected reactive power.
Step 302) establishing a heat supply network sensitivity calculation formula according to a source-load equation in the heat supply network dynamic model, wherein the heat supply network sensitivity calculation formula comprises the sensitivity of the water supply temperature at a source node relative to the water supply temperature at a non-source node, the sensitivity of the water supply temperature at a node relative to the return water temperature at the node, and the sensitivity of the water supply temperature at the node relative to the thermal power of the node, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the formula (I), the compound is shown in the specification,
indicating the temperature of the water supply at node i at time k,
represents the return water temperature of the node i at the moment k,
represents the thermal power of node i at time k, diInjection flow for node i, CρIs specific heat capacity of working medium, Vh,ldFor a set of load nodes in a heat network, K1,ji,nkIs K1The jth row and the ith column of the block matrix.
Step 303) establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model. The coupling device model may be expressed as:
φCHP＝Ch-cuPCHP (25)
φHP＝ckPHP (26)
in the formula, phiCHPAnd PCHPRespectively representing the thermal power and the electric power output by the cogeneration plant, ChTo the maximum output thermal power of the cogeneration plant, cuIs the heat-to-power ratio of a cogeneration plantHPAnd PHPRespectively representing the thermal power output and the electric power consumed by the thermoelectric conversion device, ckIs heatThe thermoelectric ratio of the electrical conversion device.
According to the formula (22) and the formula (24), the sensitivity of the thermal power of the source node in the heat supply network relative to the thermal power of the load node can be obtained. By combining the coupling equipment model, the sensitivities of the node injection power, the voltage amplitude and the branch transmission power in the power grid to the thermal power of the load node in the heat supply network can be obtained, and the sensitivities are respectively shown as a formula (28) and a formula (29).
In the formula (I), the compound is shown in the specification,
representing the injected active power of node i at time k,
representing the transmission power, V, of the branch at time kcpRepresenting a coupling node set in the electric heating integrated energy system, eta represents the conversion coefficient of the coupling equipment, and eta is-1/c for the electric heating cogeneration equipmentu(ii) a For a thermoelectric conversion device, η is 1/ck。
Step 304) constructing a dynamic safety control strategy according to the coupling relation of the electric heating system, wherein the dynamic safety control strategy comprises adjusting the voltage amplitude of a node injection power control node and the active power transmitted by a branch, adjusting the water supply temperature of a source node to control the water supply temperature of a load node, and adjusting the voltage amplitude of the node water supply temperature control node and the active power transmitted by the branch, which are respectively shown in formulas (30) to (33).
ΔP＝(U-Ulim+dU)/SUP (30)
Wherein, Δ P represents the control quantity of the node injection power,
control quantity, U, representing supply water temperature at node i at time klimAnd Pl,limRepresenting node voltage amplitude threshold and branch transmission active power threshold, dU and dPlIs the margin of the node voltage amplitude.
The control of the pipe end temperature before and after and the comparison of the node injection active power are shown in fig. 3 and fig. 4 respectively.
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.
Claims (4)
1. A dynamic safety control method for an electric heating comprehensive energy system is characterized by comprising the following steps:
step 10) deducing transmission coefficients of boundary conditions and initial conditions based on a heat supply network dynamic model of a central implicit differential format, and establishing a pipeline equivalent model;
step 20) establishing a dynamic source-load equation of the heat supply network based on the pipeline equivalent model and in combination with a network topology equation;
and step 30) establishing a sensitivity calculation formula of the electric power and thermal power network based on the power flow model and the dynamic source-load equation of the heat supply network, establishing a global mapping relation of the electric heating comprehensive energy system by combining the coupling equipment model, and establishing a dynamic safety control method.
2. The dynamic safety control method of the electric-thermal integrated energy system according to claim 1, wherein the step 10) specifically comprises:
step 101) according to the central implicit difference format, establishing a heat supply network dynamic model as follows:
in the formula, i and j represent a space step and a time step of the temperature distribution of the pipeline, respectively, T is the temperature of the pipeline with the ambient temperature as a reference value,
represents the temperature, μ, of the segment at time j at the ith of the pipe1-3The method comprises the following steps of constructing a constant coefficient term for expressing temperature transmission characteristics, specifically:
wherein Δ x and Δ t are the space and time step of the difference; v is the flow velocity of the pipeline, and r is the thermal resistance coefficient of the pipeline;
step 102) expressing the initial conditions, the boundary conditions and the temperature of the end of the pipeline in a vector form, and introducing a permutation and combination calculation formula
Defining reduced coefficient terms simultaneously
Initial condition T0Boundary condition T0And pipe end temperature TMThe vector is shown in equation (5):
in the formula, M and N are the space and time segment number of the pipeline temperature respectively;
permutation and combination calculation formula
Expressed as:
step 103) arbitrary boundary conditions
And status point
There are two types of transmission paths:
(1)
first experiences a spatial transmission loss to
Then experiences M-1 section space transmission loss and k-j section time transmission loss to
(2)
Experiences a loss of transmission in space and time to
Then undergoes M-1 section space transmission loss and k-j-1 section time transmission loss to
αkjThe transmission relationship of the boundary condition at the moment j to the branch terminal temperature at the moment k is shown as follows:
in the formula, the 1 st term and the 2 nd term on the right of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively,
and
denotes alphakjCombined parameter item of middle two kinds of transmission pathRespectively expressed as:
step 104) any initial conditions
And status point
There are two types of transmission paths:
(1)Tj 0the transmission loss is firstly lost to T for a period of timej 1Then experiences the spatial transmission loss of the M-j section and the time transmission loss of the N-j section to
(2)Tj 0Experiences a loss of transmission in space and time to
Then undergoes M-j-1 space transmission loss and k-1 time transmission loss to
βkjAnd (3) representing the transmission relation of the j-th section initial condition to the branch terminal temperature at the k moment, and represented as:
in the formula, the 1 st and 2 nd terms on the right of the equal sign correspond to the transmission coefficients, ψ, of the path (1) and the path (2), respectivelykj1And psikj2Is represented by betakjThe combined parameter items of the two kinds of transmission paths are respectively expressed as:
step 105) combining the transmission matrices of the initial conditions and the boundary conditions, the pipe end temperature at a single moment is expressed as:
expanding the formula (11) from a single moment k to a full period, wherein the temperature vector at the tail end of the pipeline is expressed as a linear combination of an initial condition and a boundary condition, and is expressed as a formula (12):
TM＝αT0+βT0 (12)
where α and β are the transmission matrices for the boundary and initial conditions, respectively.
3. The dynamic safety control method of the electric-thermal integrated energy system according to claim 1, wherein the step 20) is as follows:
step 201) constructing a multi-type network topological equation comprising a pipeline initial end temperature-node temperature incidence matrix ApnPipeline flow-inflow node incidence matrix AinPipeline flow-outflow node incidence matrix Aout(ii) a The block matrix is an N-dimensional square matrix, NbAnd NhThe number of the pipelines and the number of the nodes of the heat supply network are respectively, and each topological equation is specifically expressed as follows:
(1) pipeline initial end temperature-node temperature correlation matrix Apn：ApnContaining Nb×NhA square matrix, if the initial end temperature of the pipeline i is the same as the node j, ApnThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
(2) pipe flow-inflow node incidence matrix Ain：AinContaining Nh×NbA square matrix, if the flow flows into the node i from the pipeline j, AinThe square matrix of the ith row and the jth column is a unit square matrix, otherwise, the unit square matrix is a unit zero square matrix；
(3) Pipe flow-inflow node incidence matrix Aout：AoutContaining Nh×NbA square matrix, if the flow flows into the pipeline j from the node i, AoutThe square matrixes of the ith row and the jth column are unit square matrixes, otherwise, the unit square matrixes are zero square matrixes;
step 202) combining a pipeline equivalent model to deduce a dynamic source-load equation of the heat supply network; based on pipeline initial end temperature vector T0And node temperature TnBy means of the correlation matrix ApnThe relationship, expressed as:
T0＝ApnTn (13)
the tube temperature mixing equation is expressed as:
in the formula, Tn,iDenotes the temperature, T, of node iM,bDenotes the end temperature, m, of the pipe bbDenotes the mass flow of the pipe b, Es,iRepresenting a set of pipes with node i as the head node, Ee,iRepresenting a set of pipes, V, with pipe i as the last nodehThe node set node outflow flow of the heat supply network comprises two parts of node outflow to load and node outflow to pipeline, and the formula (14) is popularized to the system by using a topological equation and expressed as follows:
AoutmT0+dTn＝AinmTM (15)
in the formula, m is a flow vector in the dynamic model, and d is a node injection flow vector;
bringing formula (12) and formula (13) into formula (15) in this order, yields:
MoutTn＝JTn+B (16)
in the formula, MoutThe method comprises the following steps of A, taking an outflow flow matrix of nodes in a dynamic model, J taking a transmission matrix of different node temperatures in the dynamic model, and B taking a temperature component of an initial condition action; mouJ and B are respectively:
Mout＝diag(AoutmApn+d),J＝AinmαApn,B＝AinmβΤ0 (17)
Vh,srand Vh,nsRespectively, the source node and the non-source node of the thermodynamic system, then equation (17) is expanded as:
in the formula (I), the compound is shown in the specification,
and
outgoing traffic matrices, T, being a set of source nodes and a set of non-source nodes, respectivelyn,srAnd Tn,nsTemperature vectors for source and non-source nodes, J11、J12、J21And J22Respectively, a block matrix, B, in a coefficient matrix JsrAnd BnsTemperature components acting on the source node and the non-source node for initial conditions;
according to a Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network is obtained and expressed as:
4. the dynamic safety control method of the electric-thermal integrated energy system according to claim 1, wherein the step 30) is as follows:
step 301), establishing a grid sensitivity calculation formula, including node voltage amplitude, phase angle sensitivity about node injection active power and reactive power, and line transmission power sensitivity to node injection active power and reactive power, as shown in formula (20) and formula (21), respectively:
in the formula, U and theta are the node voltage amplitude and phase angle respectively, P and Q are the node injected active power and reactive power respectively, and PlAnd QlActive and reactive power, S, respectively, for the transmission of network branchesUPSensitivity, S, representing the magnitude of the node voltage with respect to the active power injected into the nodeθPSensitivity, S, representing the phase angle of the node voltage with respect to the active power injected into the nodeUQIndicating the sensitivity of the node voltage amplitude with respect to the node injected reactive power, SθQRepresenting the sensitivity of the node voltage phase angle to the node injected reactive power;
step 302) establishing a heat supply network sensitivity calculation formula according to a source-load equation in the heat supply network dynamic model, wherein the heat supply network sensitivity calculation formula comprises the sensitivity of the water supply temperature at a source node relative to the water supply temperature at a non-source node, the sensitivity of the water supply temperature at a node relative to the return water temperature at the node, and the sensitivity of the water supply temperature at the node relative to the thermal power of the node, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the formula (I), the compound is shown in the specification,
indicating the temperature of the water supply at node i at time k,
represents the return water temperature of the node i at the moment k,
represents the thermal power of node i at time k, diInjection flow for node i, CρIs specific heat capacity of working medium, Vh,ldFor a set of load nodes in a heat network, K1,ji,nkIs K1Elements of the jth row and the ith column in the block matrix of the jth row and the ith column;
step 303) establishing a global mapping relation of the electric heating comprehensive energy system by combining a coupling equipment model, wherein the coupling equipment model is expressed as:
φCHP＝Ch-cuPCHP (25)
φHP＝ckPHP (26)
in the formula, phiCHPAnd PCHPRespectively representing the thermal power and the electric power output by the cogeneration plant, ChTo the maximum output thermal power of the cogeneration plant, cuIs the heat-to-power ratio of a cogeneration plantHPAnd PHPRespectively representing the thermal power output and the electric power consumed by the thermoelectric conversion device, ckIs the thermoelectric ratio of the thermoelectric conversion device;
according to the formula (22) and the formula (24), obtaining the sensitivity of the thermal power of the source node in the heat supply network relative to the thermal power of the load node; and (3) obtaining the sensitivities of the injection power of the node in the power grid, the voltage amplitude and the branch transmission power relative to the thermal power of the load node in the heat supply network by combining a coupling equipment model, wherein the sensitivities are respectively expressed as a formula (28) and a formula (29):
in the formula (I), the compound is shown in the specification,
representing the injected active power, P, of node i at time kl kRepresenting the transmission power, V, of the branch at time kcpRepresenting a coupling node set in the electric heating integrated energy system, eta represents the conversion coefficient of the coupling equipment, and eta is-1/c for the electric heating cogeneration equipmentu(ii) a For a thermoelectric conversion device, η is 1/ck；
Step 304) according to the coupling relation of the electric heating system, a dynamic safety control method is established, which comprises the steps of adjusting the node injection power to control the node voltage amplitude and the branch transmission active power, adjusting the source node water supply temperature to control the load node water supply temperature, adjusting the source node water supply temperature to control the node voltage amplitude and the branch transmission active power, and respectively entering (30) to (33):
ΔP＝(U-Ulim+dU)/SUP (30)
wherein, Δ P represents the control quantity of the node injection power,
control quantity, U, representing supply water temperature at node i at time klimAnd Pl,limRepresenting node voltage amplitude threshold and branch transmission active power threshold, dU and dPlIs the margin of the node voltage amplitude.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110828652.4A CN113486532B (en) | 2021-07-22 | 2021-07-22 | Dynamic safety control method for electric heating comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110828652.4A CN113486532B (en) | 2021-07-22 | 2021-07-22 | Dynamic safety control method for electric heating comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113486532A true CN113486532A (en) | 2021-10-08 |
| CN113486532B CN113486532B (en) | 2023-12-19 |
Family
ID=77942872
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110828652.4A Active CN113486532B (en) | 2021-07-22 | 2021-07-22 | Dynamic safety control method for electric heating comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113486532B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011134432A1 (en) * | 2010-04-30 | 2011-11-03 | 新奥科技发展有限公司 | Smart energy network for achieving optimum utilization of energy and method for providing energy trading and service |
| CN110516951A (en) * | 2019-08-22 | 2019-11-29 | 南京工业大学 | A kind of integrated energy system dispatching method of dynamic interval |
| CN111898816A (en) * | 2020-07-21 | 2020-11-06 | 南京信息工程大学 | Dynamic state estimation method for comprehensive energy system |
| CN112330127A (en) * | 2020-10-29 | 2021-02-05 | 东南大学 | Static safety control method for multi-time-scale electric heating integrated energy system |
| CN113111515A (en) * | 2021-04-14 | 2021-07-13 | 东南大学 | Unified modeling method of comprehensive energy system |
-
2021
- 2021-07-22 CN CN202110828652.4A patent/CN113486532B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011134432A1 (en) * | 2010-04-30 | 2011-11-03 | 新奥科技发展有限公司 | Smart energy network for achieving optimum utilization of energy and method for providing energy trading and service |
| CN110516951A (en) * | 2019-08-22 | 2019-11-29 | 南京工业大学 | A kind of integrated energy system dispatching method of dynamic interval |
| CN111898816A (en) * | 2020-07-21 | 2020-11-06 | 南京信息工程大学 | Dynamic state estimation method for comprehensive energy system |
| CN112330127A (en) * | 2020-10-29 | 2021-02-05 | 东南大学 | Static safety control method for multi-time-scale electric heating integrated energy system |
| CN113111515A (en) * | 2021-04-14 | 2021-07-13 | 东南大学 | Unified modeling method of comprehensive energy system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113486532B (en) | 2023-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020093296A1 (en) | Interval power flow calculation method for power-heat integrated energy system | |
| CN112330127A (en) | Static safety control method for multi-time-scale electric heating integrated energy system | |
| CN104009247B (en) | A kind of Solid Oxide Fuel Cell local temperature method of estimation | |
| CN111428351B (en) | Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method | |
| CN108920866B (en) | Heat supply network dynamic regulation operating parameter estimation method based on moving horizon estimation theory | |
| CN109726483A (en) | A kind of radial heat supply network model of electric heating interconnection integrated energy system and its system | |
| CN110728032B (en) | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network | |
| Marechal et al. | Targeting the optimal integration of steam networks: Mathematical tools and methodology | |
| Huang et al. | A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale | |
| CN113486532A (en) | Dynamic safety control method for electric heating comprehensive energy system | |
| CN112257279B (en) | Method for constructing feasible region of electric heating comprehensive energy system | |
| CN113255105B (en) | Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure | |
| CN114221346A (en) | Load flow calculation method of comprehensive energy system | |
| CN112528482B (en) | Cascading failure simulation method for thermoelectric coupling system under extremely cold disaster | |
| CN113111515A (en) | Unified modeling method of comprehensive energy system | |
| CN113690891B (en) | Analysis-method-based probability power flow determination method for electric-thermal interconnection comprehensive energy system | |
| CN110737993B (en) | Multi-energy complementary system operation boundary analysis method considering load uncertainty | |
| CN113111555B (en) | Quality control thermodynamic system energy flow rapid calculation method based on superposition decoupling method | |
| CN113919754B (en) | Block chain-based distributed state estimation method for comprehensive energy system | |
| CN113221428B (en) | Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system | |
| CN112182905B (en) | Heat supply pipe network simulation method and device for comprehensive energy system | |
| CN114549232A (en) | Hybrid energy flow calculation method for electricity-heat comprehensive energy system | |
| Chen et al. | Steady-state Flow Analysis of District Heating System Using the Holomorphic Embedding | |
| CN112417698B (en) | Dynamic two-port thermodynamic system model based on quality adjustment | |
| Da-Peng et al. | Research on modeling and parameter identification of circulating cooling water 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 |