CN113486532B - Dynamic safety control method for electric heating comprehensive energy system - Google Patents

Dynamic safety control method for electric heating comprehensive energy system Download PDF

Info

Publication number
CN113486532B
CN113486532B CN202110828652.4A CN202110828652A CN113486532B CN 113486532 B CN113486532 B CN 113486532B CN 202110828652 A CN202110828652 A CN 202110828652A CN 113486532 B CN113486532 B CN 113486532B
Authority
CN
China
Prior art keywords
node
temperature
matrix
power
pipeline
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110828652.4A
Other languages
Chinese (zh)
Other versions
CN113486532A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN202110828652.4A priority Critical patent/CN113486532B/en
Publication of CN113486532A publication Critical patent/CN113486532A/en
Application granted granted Critical
Publication of CN113486532B publication Critical patent/CN113486532B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Primary Health Care (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Human Resources & Organizations (AREA)
  • General Health & Medical Sciences (AREA)
  • Operations Research (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Control Of Temperature (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a dynamic safety control method of an electric heating comprehensive energy system, which comprises the following steps: 10 Based on the heat supply network dynamic model of the central implicit differential format, deducing a transmission matrix of boundary conditions and initial conditions, and establishing a pipeline equivalent model; 20 Based on the equivalent model of the pipeline, combining with a network topology equation to establish a source-load equation of the heat supply network dynamic state; 30 Based on the power grid tide model and the heat supply network dynamic source-load equation, establishing a sensitivity calculation formula of the electric power and heat supply network, and combining the coupling equipment model, establishing a global mapping relation of the electric heating comprehensive energy system, and constructing a dynamic safety control method.

Description

Dynamic safety control method for electric heating comprehensive energy system
Technical Field
The invention relates to the field of modeling and operation analysis of energy systems, in particular to a dynamic safety control method of an electric heating comprehensive energy system.
Background
The growing environmental problems and energy crisis promote the transformation of energy technology, and the electric heating comprehensive energy system is used as a main form of energy transformation, is beneficial to improving the complementarity of electric energy and heat energy utilization, and is widely applied in practical engineering. However, due to the multi-time scale characteristic differences of the electric power system and the thermodynamic system, the model of the electrothermal integrated energy group is essentially a complex set of high-dimensional constant and partial differential equations, which are very complex to solve. Furthermore, since electrical energy propagates at the speed of light in electrical power systems, it tends to stabilize in milliseconds; the heat energy is transmitted in the thermodynamic system through the flow of the working medium, and a certain delay is needed to reach stability after the control instruction is issued, so that the control time scale in the electric power system is generally shorter, and the control time scale of the thermodynamic system is generally longer.
The existing research generally establishes a safety control strategy according to the energy flow calculation result at the control moment, but the dynamic process evolution between the control time intervals is often ignored due to the longer control time intervals of the thermodynamic system, so that the establishment of the control strategy is more unilateral. In addition, the existing research generally adopts a steady-state model to carry out modeling analysis, but ignores dynamic and multi-time scale characteristics in the electric heating comprehensive energy system, and the obtained result has larger difference from the actual situation. In addition, because of the coupling of heat and water, 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, and lacks an accurate and quantitative safety control strategy, and adjusts the possible out-of-limit state in the operation process, thereby ensuring the operation of the system.
Disclosure of Invention
In order to solve the defects in the background art, the invention aims to provide a dynamic safety control method of an electric heating comprehensive energy system, which deduces transmission coefficients of initial conditions and boundary conditions in a differential-format heat supply network dynamic model, establishes a pipeline equivalent model of state quantity distribution, constructs a heat supply network dynamic source-charge equation by combining a topological equation, reveals a source charge mapping relation, and provides a dynamic safety control strategy through global sensitivity calculation of the electric heating comprehensive energy system. Provides theoretical guidance for the safe operation of the electric heating comprehensive energy system.
The aim of the invention can be achieved by the following technical scheme:
a dynamic safety control method of an electric heating comprehensive energy system comprises the following steps:
step 10), based on a heat supply network dynamic model in a central implicit differential format, deducing transmission coefficients of boundary conditions and initial conditions, and establishing a pipeline equivalent model;
step 20), establishing a source-load equation of the heat supply network dynamic based on the pipeline equivalent model and combining a network topology equation;
step 30) based on a power grid tide model and a heat supply network dynamic source-load equation, establishing a sensitivity calculation formula of an electric power and heat supply network, and combining a coupling equipment model, establishing a global mapping relation of an electric heating comprehensive energy system, and constructing a dynamic safety control method.
Further, the step 10) specifically includes:
step 101) establishing a dynamic model of the heat supply network according to the central implicit differential format as follows:
wherein i and j respectively represent a space step and a time step of the pipeline temperature distribution, T is the pipeline temperature taking the ambient temperature as a reference value,represents the temperature, μ of the section at the i-th point of the pipe at the instant j 1-3 The constant coefficient term is constructed for expressing the temperature transmission characteristic, specifically:
wherein Δx and Δt are the differential spatial and temporal steps; v is the pipeline flow rate, r is the pipeline thermal coefficient;
step 102) expressing the initial conditions, boundary conditions and pipeline end temperature in vector form and introducing the permutation and combination calculation typeAt the same time define the simplification coefficient item->Initial condition T 0 Boundary condition T 0 Pipe end temperatureT M The vector is shown in formula (5):
wherein M and N are the space and time segment numbers of the pipeline temperature respectively;
permutation and combination calculation formulaExpressed as:
step 103) arbitrary boundary conditionsAnd status point->There are two types of transmission paths:
(1)first undergo a section of space transmission loss to +.>Then the M-1 section space transmission loss and the k-j section time transmission loss are subjected to +.>
(2)First undergo a period of space and time transmission loss to +.>Then experience M-1 section space transmission loss and k-j-1 section time transmission lossConsumption to->
α kj The transmission relation of the boundary condition at the moment j to the tail end temperature of the branch at the moment k is expressed as follows:
wherein items 1 and 2 on the right side of the equal sign correspond to the transmission coefficients of path (1) and path (2), respectively,andrepresenting alpha kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 104) any initial conditionsAnd status point->There are two types of transmission paths:
(1)first, a certain period of time is passed until the transmission loss reaches->Then the space transmission loss of M-j sections and the time transmission loss of N-j sections are subjected to +.>
(2)First undergo a period of space and time transmission loss to +.>Then the space transmission loss of M-j-1 section and the time transmission loss of k-1 section are carried out to +.>
β kj The transmission relation of the j-th section initial condition to the end temperature of the k-moment branch is expressed as follows:
wherein the 1 st and 2 nd items on the right side of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively, ψ kj1 Sum phi kj2 Representing beta kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 105) combining the transmission matrix of the initial conditions and the boundary conditions, the pipe end temperature at a single instant is expressed as:
expanding equation (11) from a single time k to full time period, the pipe end temperature vector is expressed as a linear combination of initial conditions and boundary conditions, expressed as equation (12):
T M ＝αT 0 +βT 0 (12)
where α and β are the transmission matrices of the boundary and initial conditions, respectively.
Further, the step 20) specifically includes the following steps:
step 201) constructing a multi-type network topology equation comprising a pipeline start point temperature-node temperature correlation matrix A pn Pipeline flow-inflow node incidence matrix A in Pipeline flow-outflow node association matrix A out The method comprises the steps of carrying out a first treatment on the surface of the The blocking matrix is an N-dimensional square matrix, N b And N h The number of pipes and the number of nodes of the heat supply network are respectively, and each topological equation is specifically expressed as follows:
(1) Pipeline start end temperature-node temperature correlation matrix A pn ：A pn Comprising N b ×N h A square matrix, if the temperature of the initial end of the pipeline i is the same as the temperature of the node j, A pn The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flow-inflow node incidence matrix A in ：A in Comprising N h ×N b A square matrix, if the flow flows from the pipeline j into the node i, A in The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(3) Pipeline flow-inflow node incidence matrix A out ：A out Comprising N h ×N b A square matrix, if the flow flows from the node i into the pipeline j, A out The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
step 202) combining a pipeline equivalent model, and pushing a dynamic source-load equation of the heat conducting network; based on the temperature vector T at the beginning of the pipeline 0 And a node temperature T n Through the association matrix A pn The association, expressed as:
T 0 ＝A pn T n (13)
the pipe temperature mixing equation is expressed as:
wherein T is n,i Representing the temperature of node i, T M,b Representing the conduit bEnd temperature, m b Representing the mass flow of the pipe b, E s,i Representing a set of pipes with node i as the head node, E e,i Representing a set of pipes ending with pipe i, V h Node outflow traffic for a node set of the heat supply network comprises two parts, namely node outflow to load and node outflow to pipeline, and the formula (14) is generalized to the system by using a topological equation, and is expressed as:
A out mT 0 +dT n ＝A in mT M (15)
wherein m is a flow vector in the dynamic model, and d is a node injection flow vector;
bringing the formula (12) and the formula (13) into the formula (15) in this order gives:
M out T n ＝JT n +B (16)
wherein M is out The method is characterized in that the method comprises the steps that an outflow flow matrix of a node in a dynamic model is obtained, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by initial conditions; m is M ou J and B are respectively:
M out ＝diag(A out mA pn +d),J＝A in mαA pn ,B＝A in mβΤ 0 (17)
V h,sr and V h,ns Respectively a source node and a non-source node set of the thermodynamic system, then the formula (17) expands to:
in the method, in the process of the invention,and->Outflow traffic matrix, T, for the source node set and the non-source node set, respectively n,sr And T n,ns Temperature vector for source node and non-source node, J 11 、J 12 、J 21 And J 22 Respectively a blocking matrix in the coefficient matrix J, B sr And B ns Acting on temperature components of the source node and the non-source node for an initial condition;
according to the Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network is obtained, and the equation is expressed as follows:
further, the step 30) specifically includes the following steps:
step 301) establishing a power grid sensitivity calculation type, including the node voltage amplitude, the sensitivity of the phase angle to the node injected active power and reactive power, and the sensitivity of the line transmission power to the node injected active power and reactive power, as shown in the formula (20) and the formula (21):
wherein U and theta are respectively the node voltage amplitude and phase angle, P and Q are respectively the node injected active power and reactive power, and P l And Q l Active power and reactive power respectively transmitted by power grid branches S UP Representing the sensitivity of the node voltage amplitude to the node injected active power, S θP Representing the sensitivity of the node voltage phase angle with respect to the node injected active power, S UQ Representing the sensitivity of the node voltage amplitude with respect to the node injection of reactive power, S θQ Representing the sensitivity of the node voltage phase angle with respect to the node injected reactive power;
step 302), according to a source-load equation in a heat supply network dynamic model, establishing a heat supply network sensitivity calculation type, wherein the heat supply network sensitivity calculation type comprises the sensitivity of water supply temperature at a source node relative to water supply temperature at a non-source node, the sensitivity of node water supply temperature relative to node backwater temperature, and the sensitivity of node water supply temperature relative to node heat power, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the method, in the process of the invention,represents the water supply temperature of node i at time k +.>Represents the backwater temperature of node i at time k, +.>Represents the thermal power of node i at time k, d i For the injection flow of node i, C ρ Is the specific heat capacity of working medium, V h,ld K is a set of load nodes in a heat supply network 1,ji,nk For K 1 The j-th row and the i-th column of the element of the nth row and the k-th column in the block matrix;
step 303) building 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 ＝C h -c u P CHP (25)
φ HP ＝c k P HP (26)
in phi CHP And P CHP Respectively representing the thermal power and the electric power output by the cogeneration equipment, C h C for maximum output thermal power of the cogeneration plant u Is a thermoelectric powerHeat-power ratio phi of co-production equipment HP And P HP C represents the thermal power output from the thermoelectric conversion device and the consumed electrical power, respectively k A thermoelectric ratio that is a thermoelectric conversion device;
obtaining the sensitivity of the thermal power of the internal node of the heat supply network to the thermal power of the load node according to the formula (22) and the formula (24); and combining the coupling equipment model to obtain the sensitivity 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, wherein the sensitivity is shown in the formula (28) and the formula (29):
in the method, in the process of the invention,represents the injected active power of node i at time k, < >>Representing the transmission power of the branch at time k, V cp Represents a coupling node set in the electric heating integrated energy system, eta represents a conversion coefficient of a coupling device, and eta= -1/c for an electric heating cogeneration device u The method comprises the steps of carrying out a first treatment on the surface of the For the thermoelectric conversion device, η=1/c k ；
Step 304) constructing a dynamic safety control method according to the coupling relation of the electric heating system, wherein the dynamic safety control method comprises the steps of 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, adjusting the voltage amplitude of the water supply temperature control node of the source node and the active power transmitted by the branch, and respectively inputting the voltage amplitude and the active power into formulas (30) to (33):
ΔP＝(U-U lim +dU)/S UP (30)
where Δp represents the control amount of the node injection power,indicating the control quantity of the water supply temperature of the node i at the moment k, U lim And P l,lim Representing node voltage amplitude threshold and branch transmission active power threshold, dU and dP l Is a margin for the node voltage magnitude.
The invention has the beneficial effects that:
the method establishes the analysis relation between the dynamic characteristics of the state quantity of the heat supply network and the initial conditions and boundary conditions, deduces the direct mapping function of the source and the load, quantitatively describes the influence of the dynamic characteristics of the heat supply network on the operation of the power grid, is beneficial to accurately formulating the dynamic safety control strategy of the electric heating comprehensive energy system taking the dynamic characteristics of the heat supply network into account, and ensures the stable operation of the system.
Drawings
The invention is further described below 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 graph showing temperature comparisons of the ends of pipes before and after control in accordance with an embodiment of the present invention;
fig. 4 is a graph showing the comparison of active power injected into the nodes before and after control in an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1: a dynamic safety control method of an electric heating comprehensive energy system comprises the following steps:
step 10), based on a heat supply network dynamic model in a central implicit differential format, deducing transmission coefficients of boundary conditions and initial conditions, and establishing a pipeline equivalent model;
step 20), establishing a source-load equation of the heat supply network dynamic based on the pipeline equivalent model and combining a network topology equation;
step 30) based on a power grid tide model and a heat supply network dynamic source-load equation, establishing a sensitivity calculation formula of an electric power and heat supply network, and establishing a global mapping relation of an electric heating comprehensive energy system by combining a coupling equipment model to construct a dynamic safety control strategy.
The step 10) specifically comprises the following steps:
step 101) establishing a dynamic model of the heat supply network according to the central implicit differential format as follows:
wherein i and j respectively represent a space step and a time step of the pipeline temperature distribution, T is the pipeline temperature taking the ambient temperature as a reference value,represents the temperature, μ of the section at the i-th point of the pipe at the instant j 1-3 A constant coefficient term is constructed for expressing the temperature transmission characteristics, which is specifically:
wherein Δx and Δt are the differential spatial and temporal steps; v is the pipeline flow rate, r is the pipeline thermal coefficient;
step 102) expressing the initial conditions, boundary conditions and pipeline end temperature as vectors and introducing the permutation and combination calculation typeAt the same time define the simplification coefficient item->Initial condition T 0 Boundary condition T 0 Pipe end temperature T M The vector is shown in formula (5). Where M and N are the number of spatial and temporal segments of the pipe temperature, respectively.
Permutation and combination calculation formulaCan be expressed as:
step 103) arbitrary boundary conditionsAnd status point->There are two types of transmission paths: (1)/>First undergo a section of space transmission loss to +.>Then the M-1 section space transmission loss and the k-j section time transmission loss are subjected to +.>(2)/> First undergo a period of space and time transmission loss to +.>Then the M-1 section space transmission loss and the k-j-1 section time transmission loss are subjected to +.>Definition alpha kj The transmission relation of the boundary condition at the moment j to the temperature of the tail end of the branch at the moment k can be expressed as follows:
wherein items 1 and 2 on the right side of the equal sign correspond to the transmission coefficients of path (1) and path (2), respectively,andrepresenting alpha kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 104) any initial conditionsAnd status point->There are two types of transmission paths: (1)/>First, a certain period of time is passed until the transmission loss reaches->Then the space transmission loss of M-j sections and the time transmission loss of N-j sections are subjected to +.>(2)/> First undergo a period of space and time transmission loss to +.>Then the space transmission loss of M-j-1 section and the time transmission loss of k-1 section are carried out to +.>Definition beta kj The transmission relation of the j-th section initial condition to the end temperature of the k-moment branch can be expressed as follows:
wherein the 1 st and 2 nd items on the right side of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively, ψ kj1 Sum phi kj2 Representing beta kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 105) combining the transmission matrix of the initial conditions and the boundary conditions, the pipe end temperature at a single instant can be expressed as:
expanding equation (11) from a single time k to full time, the pipe end temperature vector may be represented as a linear combination of initial and boundary conditions, and may be represented as equation (12), where α and β are the transmission matrices of the boundary and initial conditions, respectively.
T M ＝αT 0 +βT 0 (12)
Said step 20) comprises:
step 201) constructing a multi-type network topology equation comprising a pipeline start point temperature-node temperature correlation matrix A pn Pipeline flow-inflow node incidence matrix A in Pipeline flow-outflow node association matrix A out . Defining a blocking matrix as an N-dimensional square matrix, N b And N h The number of pipes and nodes of the heat supply network respectively, each topological equation can be specifically expressed as:
(1) Pipeline start end temperature-node temperature correlation matrix A pn ：A pn Comprising N b ×N h A square matrix, if the temperature of the initial end of the pipeline i is the same as the temperature of the node j, A pn The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flowQuantity-inflow node correlation matrix A in ：A in Comprising N h ×N b A square matrix, if the flow flows from the pipeline j into the node i, A in The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flow-inflow node incidence matrix A out ：A out Comprising N h ×N b A square matrix, if the flow flows from node i into pipe j, a out The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
step 202) combining the pipeline equivalent model, and pushing a dynamic source-load equation of the heat conducting network. Based on the temperature vector T at the beginning of the pipeline 0 And a node temperature T n Can pass through the incidence matrix A pn The association, expressed as:
T 0 ＝A pn T n (13)
the pipe temperature mixing equation can be expressed as:
wherein T is n,i Representing the temperature of node i, T M,b Represents the end temperature of the pipe b, m b Representing the mass flow of the pipe b, E s,i Representing a set of pipes with node i as the head node, E e,i Representing a set of pipes ending with pipe i, V h Is a set of nodes of a heat supply network. Considering that the node outflow traffic includes both the node outflow to load and the node outflow to pipe, the topology equation can be used to generalize equation (14) to the system, expressed as:
A out mT 0 +dT n ＝A in mT M (15)
where m is the flow vector in the dynamic model and d is the node injection flow vector. By bringing the formula (12) and the formula (13) into the formula (15) in this order, it is possible to obtain:
M out T n ＝JT n +B (16)
wherein M is out The method is characterized in that the method is used for generating an outflow flow matrix of nodes in a dynamic model, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by initial conditions. M is M ou J and B are respectively:
M out ＝diag(A out mA pn +d),J＝A in mαA pn ,B＝A in mβΤ 0 (17)
definition V h,sr And V h,ns The set of source nodes and non-source nodes of the thermodynamic system, respectively, then equation (17) can be expanded to:
in the method, in the process of the invention,and->Outflow traffic matrix, T, for the source node set and the non-source node set, respectively n,sr And T n,ns Temperature vector for source node and non-source node, J 11 、J 12 、J 21 And J 22 Respectively a blocking matrix in the coefficient matrix J, B sr And B ns The temperature components of the source node and the non-source node are acted upon for the initial conditions.
According to the Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network can be obtained, and the equation is expressed as follows:
the step 30) includes:
step 301) establishing a power grid sensitivity calculation type, including the node voltage amplitude, the sensitivity of the phase angle to the node injected active power and reactive power, and the sensitivity of the line transmission power to the node injected active power and reactive power, as shown in the formula (20) and the formula (21):
wherein U and theta are respectively the node voltage amplitude and phase angle, P and Q are respectively the node injected active power and reactive power, and P l And Q l Active power and reactive power respectively transmitted by power grid branches S UP Representing the sensitivity of the node voltage amplitude to the node injected active power, S θP Representing the sensitivity of the node voltage phase angle with respect to the node injected active power, S UQ Representing the sensitivity of the node voltage amplitude with respect to the node injected reactive power; s is S θQ Indicating the sensitivity of the node voltage phase angle with respect to the node injected reactive power.
Step 302), according to a source-load equation in a heat supply network dynamic model, establishing a heat supply network sensitivity calculation type, wherein the heat supply network sensitivity calculation type comprises the sensitivity of water supply temperature at a source node relative to water supply temperature at a non-source node, the sensitivity of node water supply temperature relative to node backwater temperature, and the sensitivity of node water supply temperature relative to node heat power, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the method, in the process of the invention,representation ofWater supply temperature of node i at time k +.>Represents the backwater temperature of node i at time k, +.>Represents the thermal power of node i at time k, d i For the injection flow of node i, C ρ Is the specific heat capacity of working medium, V h,ld K is a set of load nodes in a heat supply network 1,ji,nk For K 1 The j-th row and the i-th column of the partitioning matrix are the elements of the n-th row and the k-th column.
Step 303) combining the coupling equipment model to establish a global mapping relation of the electric heating comprehensive energy system. The coupling device model can be expressed as:
φ CHP ＝C h -c u P CHP (25)
φ HP ＝c k P HP (26)
in phi CHP And P CHP Respectively representing the thermal power and the electric power output by the cogeneration equipment, C h C for maximum output thermal power of the cogeneration plant u For the heat-power ratio, phi, of a cogeneration plant HP And P HP C represents the thermal power output from the thermoelectric conversion device and the consumed electrical power, respectively k Is the thermoelectric ratio of the thermoelectric conversion device.
From equations (22) and (24), the sensitivity of the thermal power of the source node of the thermal network with respect to the thermal power of the load node can be obtained. And combining the coupling equipment model to obtain the sensitivity 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, wherein the sensitivity is shown in a formula (28) and a formula (29) respectively.
In the method, in the process of the invention,represents the injected active power of node i at time k, < >>Representing the transmission power of the branch at time k, V cp Represents a coupling node set in the electric heating integrated energy system, eta represents a conversion coefficient of a coupling device, and eta= -1/c for an electric heating cogeneration device u The method comprises the steps of carrying out a first treatment on the surface of the For the thermoelectric conversion device, η=1/c k 。
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 the steps of 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, adjusting the voltage amplitude of the water supply temperature control node of the source node and the active power transmitted by the branch, and respectively inputting the voltage amplitude and the active power into formulas (30) to (33).
ΔP＝(U-U lim +dU)/S UP (30)
Where Δp represents the control amount of the node injection power,indicating the control quantity of the water supply temperature of the node i at the moment k, U lim And P l,lim Representing node voltage amplitude threshold and branch transmission active power threshold, dU and dP l Is a margin for the node voltage magnitude.
Application examples:
the system shown in fig. 2 is taken as an example.
As shown in fig. 1, the embodiment of the invention provides a dynamic safety control method of an electric heating integrated energy system, which comprises the following steps:
step 10), based on a heat supply network dynamic model in a central implicit differential format, deducing transmission coefficients of boundary conditions and initial conditions, and establishing a pipeline equivalent model;
step 20), establishing a source-load equation of the heat supply network dynamic based on the pipeline equivalent model and combining a network topology equation;
step 30) based on a power grid tide model and a heat supply network dynamic source-load equation, establishing a sensitivity calculation formula of an electric power and heat supply network, and establishing a global mapping relation of an electric heating comprehensive energy system by combining a coupling equipment model to construct a dynamic safety control strategy.
In the above embodiment, the step 10) specifically includes:
step 101) establishing a dynamic model of the heat supply network according to the central implicit differential format as follows:
wherein i and j respectively represent a space step and a time step of the pipeline temperature distribution, T is the pipeline temperature taking the ambient temperature as a reference value,represents the temperature, μ of the section at the i-th point of the pipe at the instant j 1-3 A constant coefficient term is constructed for expressing the temperature transmission characteristics, which is specifically:
/>
wherein Δx and Δt are the differential spatial and temporal steps; v is the pipeline flow rate, r is the pipeline thermal coefficient;
step 102) expressing the initial conditions, boundary conditions and pipeline end temperature as vectors and introducing the permutation and combination calculation typeAt the same time define the simplification coefficient item->Initial condition T 0 Boundary condition T 0 Pipe end temperature T M The vector is shown in formula (5). Where M and N are the number of spatial and temporal segments of the pipe temperature, respectively.
Permutation and combination calculation formulaCan be expressed as:
step 103) arbitrary boundary conditionsAnd status point->There are two types of transmission paths: (1)/>First undergo a section of space transmission loss to +.>Then the M-1 section space transmission loss and the k-j section time transmission loss are subjected to +.>(2) First undergo a period of space and time transmission loss to +.>Then the M-1 section space transmission loss and the k-j-1 section time transmission loss are subjected to +.>Definition alpha kj The transmission relation of the boundary condition at the moment j to the temperature of the tail end of the branch at the moment k can be expressed as follows:
wherein items 1 and 2 on the right side of the equal sign correspond to the transmission coefficients of path (1) and path (2), respectively,andrepresenting alpha kj Combination of two types of transmission pathsParameter items, expressed as:
step 104) any initial conditionsAnd status point->There are two types of transmission paths: (1)/>First, a certain period of time is passed until the transmission loss reaches->Then the space transmission loss of M-j sections and the time transmission loss of N-j sections are subjected to +.>(2) First undergo a period of space and time transmission loss to +.>Then the space transmission loss of M-j-1 section and the time transmission loss of k-1 section are carried out to +.>Definition beta kj The transmission relation of the j-th section initial condition to the end temperature of the k-moment branch can be expressed as follows:
wherein the 1 st and 2 nd items on the right side of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively, ψ kj1 Sum phi kj2 Representing beta kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 105) combining the transmission matrix of the initial conditions and the boundary conditions, the pipe end temperature at a single instant can be expressed as:
expanding equation (11) from a single time k to full time, the pipe end temperature vector may be represented as a linear combination of initial and boundary conditions, and may be represented as equation (12), where α and β are the transmission matrices of the boundary and initial conditions, respectively.
T M ＝αT 0 +βT 0 (12)
In the above embodiment, the step 20) specifically includes:
step 201) constructing a multi-type network topology equation comprising a pipeline start point temperature-node temperature correlation matrix A pn Pipeline flow-inflow node incidence matrix A in Pipeline flow-outflow node association matrix A out . Defining a blocking matrix as an N-dimensional square matrix, N b And N h The number of pipes and nodes of the heat supply network respectively, each topological equation can be specifically expressed as:
(1) Pipeline start end temperature-node temperature correlation matrix A pn ：A pn Comprising N b ×N h A square matrix, if the temperature of the initial end of the pipeline i is the same as the temperature of the node j, A pn The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flow-inflow node incidence matrix A in ：A in Comprising N h ×N b A square matrix, if the flow flows from the pipeline j into the node i, A in The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flow-inflow node incidence matrix A out ：A out Comprising N h ×N b A square matrix, if the flow flows from node i into pipe j, a out The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
step 202) combining the pipeline equivalent model, and pushing a dynamic source-load equation of the heat conducting network. Based on the temperature vector T at the beginning of the pipeline 0 And a node temperature T n Can pass through the incidence matrix A pn The association, expressed as:
T 0 ＝A pn T n (13)
the pipe temperature mixing equation can be expressed as:
wherein T is n,i Representing the temperature of node i, T M,b Represents the end temperature of the pipe b, m b Representing the mass flow of the pipe b, E s,i Representing a set of pipes with node i as the head node, E e,i Representing a set of pipes ending with pipe i, V h Is a set of nodes of a heat supply network. Considering that the node outflow traffic includes both the node outflow to load and the node outflow to pipe, the topology equation can be used to generalize equation (14) to the system, expressed as:
A out mT 0 +dT n ＝A in mT M (15)
where m is the flow vector in the dynamic model and d is the node injection flow vector. By bringing the formula (12) and the formula (13) into the formula (15) in this order, it is possible to obtain:
M out T n ＝JT n +B (16)
wherein M is out Is a dynamic modelAnd the outflow flow matrix of the middle node, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by the initial condition. M is M ou J and B are respectively:
M out ＝diag(A out mA pn +d),J＝A in mαA pn ,B＝A in mβΤ 0 (17)
definition V h,sr And V h,ns The set of source nodes and non-source nodes of the thermodynamic system, respectively, then equation (17) can be expanded to:
in the method, in the process of the invention,and->Outflow traffic matrix, T, for the source node set and the non-source node set, respectively n,sr And T n,ns Temperature vector for source node and non-source node, J 11 、J 12 、J 21 And J 22 Respectively a blocking matrix in the coefficient matrix J, B sr And B ns The temperature components of the source node and the non-source node are acted upon for the initial conditions.
According to the Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network can be obtained, and the equation is expressed as follows:
in the above embodiment, the step 30) specifically includes:
step 301) establishing a power grid sensitivity calculation type, including the node voltage amplitude, the sensitivity of the phase angle to the node injected active power and reactive power, and the sensitivity of the line transmission power to the node injected active power and reactive power, as shown in the formula (20) and the formula (21):
wherein U and theta are respectively the node voltage amplitude and phase angle, P and Q are respectively the node injected active power and reactive power, and P l And Q l Active power and reactive power respectively transmitted by power grid branches S UP Representing the sensitivity of the node voltage amplitude to the node injected active power, S θP Representing the sensitivity of the node voltage phase angle with respect to the node injected active power, S UQ Representing the sensitivity of the node voltage amplitude with respect to the node injected reactive power; s is S θQ Indicating the sensitivity of the node voltage phase angle with respect to the node injected reactive power.
Step 302), according to a source-load equation in a heat supply network dynamic model, establishing a heat supply network sensitivity calculation type, wherein the heat supply network sensitivity calculation type comprises the sensitivity of water supply temperature at a source node relative to water supply temperature at a non-source node, the sensitivity of node water supply temperature relative to node backwater temperature, and the sensitivity of node water supply temperature relative to node heat power, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the method, in the process of the invention,the water supply temperature of node i at time k is indicated,/>represents the backwater temperature of node i at time k, +.>Represents the thermal power of node i at time k, d i For the injection flow of node i, C ρ Is the specific heat capacity of working medium, V h,ld K is a set of load nodes in a heat supply network 1,ji,nk For K 1 The j-th row and the i-th column of the partitioning matrix are the elements of the n-th row and the k-th column.
Step 303) combining the coupling equipment model to establish a global mapping relation of the electric heating comprehensive energy system. The coupling device model can be expressed as:
φ CHP ＝C h -c u P CHP (25)
φ HP ＝c k P HP (26)
in phi CHP And P CHP Respectively representing the thermal power and the electric power output by the cogeneration equipment, C h C for maximum output thermal power of the cogeneration plant u For the heat-power ratio, phi, of a cogeneration plant HP And P HP C represents the thermal power output from the thermoelectric conversion device and the consumed electrical power, respectively k Is the thermoelectric ratio of the thermoelectric conversion device.
From equations (22) and (24), the sensitivity of the thermal power of the source node of the thermal network with respect to the thermal power of the load node can be obtained. And combining the coupling equipment model to obtain the sensitivity 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, wherein the sensitivity is shown in a formula (28) and a formula (29) respectively.
In the method, in the process of the invention,represents the injected active power of node i at time k, < >>Representing the transmission power of the branch at time k, V cp Represents a coupling node set in the electric heating integrated energy system, eta represents a conversion coefficient of a coupling device, and eta= -1/c for an electric heating cogeneration device u The method comprises the steps of carrying out a first treatment on the surface of the For the thermoelectric conversion device, η=1/c k 。
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 the steps of 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, adjusting the voltage amplitude of the water supply temperature control node of the source node and the active power transmitted by the branch, and respectively inputting the voltage amplitude and the active power into formulas (30) to (33).
ΔP＝(U-U lim +dU)/S UP (30)
Where Δp represents the control amount of the node injection power,indicating the control quantity of the water supply temperature of the node i at the moment k, U lim And P l,lim Representing node voltage amplitude threshold and branch transmission active power threshold, dU and dP l Is a margin for the node voltage magnitude.
The temperature of the pipeline end before and after control and the comparison of the node injection active power are shown in fig. 3 and fig. 4 respectively.
In the description of the present specification, the descriptions of the terms "one embodiment," "example," "specific example," and the like, 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 present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. 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 has shown and described the basic principles, principal 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, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims.

Claims (2)

1. The dynamic safety control method of the electric heating comprehensive energy system is characterized by comprising the following steps of:
step 10), based on a heat supply network dynamic model in a central implicit differential format, deducing transmission coefficients of boundary conditions and initial conditions, and establishing a pipeline equivalent model;
step 20), establishing a source-load equation of the heat supply network dynamic based on the pipeline equivalent model and combining a network topology equation;
step 30) based on a power grid tide model and a heat supply network dynamic source-load equation, establishing a sensitivity calculation type of an electric power and heat supply network, and establishing a global mapping relation of an electric heating comprehensive energy system by combining a coupling equipment model to construct a dynamic safety control method;
the specific steps of the step 10) are as follows:
step 101) establishing a dynamic model of the heat supply network according to the central implicit differential format as follows:
wherein i and j respectively represent a space step and a time step of the pipeline temperature distribution, T is the pipeline temperature taking the ambient temperature as a reference value,represents the temperature, μ of the section at the i-th point of the pipe at the instant j 1-3 The constant coefficient term is constructed for expressing the temperature transmission characteristic, specifically:
wherein Δx and Δt are the differential spatial and temporal steps; v is the pipeline flow rate, r is the pipeline thermal coefficient;
step 102) expressing the initial conditions, boundary conditions and pipeline end temperature in vector form and introducing the permutation and combination calculation typeAt the same time define the simplification coefficient item->Initial condition T 0 Boundary condition T 0 Pipe end temperature T M The vector is shown in formula (5):
wherein M and N are the space and time segment numbers of the pipeline temperature respectively;
permutation and combination calculation formulaExpressed as:
step 103) arbitrary boundary conditionsAnd status point->There are two types of transmission paths:
(1)first undergo a section of space transmission loss to +.>Then the M-1 section space transmission loss and the k-j section time transmission loss are subjected to +.>
(2)First undergo a period of space and time transmission loss to +.>Then the M-1 section space transmission loss and the k-j-1 section time transmission loss are subjected to +.>
α kj The transmission relation of the boundary condition at the moment j to the tail end temperature of the branch at the moment k is expressed as follows:
wherein items 1 and 2 on the right side of the equal sign correspond to the transmission coefficients of path (1) and path (2), respectively,and->Representing alpha kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 104) any initial conditionsAnd status point->There are two types of transmission paths:
(1)warp knittingTransmission loss to->Then the space transmission loss of M-j sections and the time transmission loss of N-j sections are subjected to +.>
(2)First undergo a period of space and time transmission loss to +.>Then the space transmission loss of M-j-1 section and the time transmission loss of k-1 section are carried out to +.>
β kj The transmission relation of the j-th section initial condition to the end temperature of the k-moment branch is expressed as follows:
wherein the 1 st and 2 nd items on the right side of the equal sign correspond to the transmission coefficients of the path (1) and the path (2), respectively, ψ kj1 Sum phi kj2 Representing beta kj The combination parameter items of the two types of transmission paths are respectively expressed as:
step 105) combining the transmission matrix of the initial conditions and the boundary conditions, the pipe end temperature at a single instant is expressed as:
expanding equation (11) from a single time k to full time period, the pipe end temperature vector is expressed as a linear combination of initial conditions and boundary conditions, expressed as equation (12):
T M ＝αT 0 +βT 0 (12)
wherein, alpha and beta are respectively the transmission matrixes of the boundary and the initial condition;
the step 20) is specifically as follows:
step 201) constructing a multi-type network topology equation comprising a pipeline start point temperature-node temperature correlation matrix A pn Pipeline flow-inflow node incidence matrix A in Pipeline flow-outflow node association matrix A out The method comprises the steps of carrying out a first treatment on the surface of the The blocking matrix is an N-dimensional square matrix, N b And N h The number of pipes and the number of nodes of the heat supply network are respectively, and each topological equation is specifically expressed as follows:
(1) Pipeline start end temperature-node temperature correlation matrix A pn ：A pn Comprising N b ×N h A square matrix, if the temperature of the initial end of the pipeline i is the same as the temperature of the node j, A pn The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(2) Pipeline flow-inflow node incidence matrix A in ：A in Comprising N h ×N b A square matrix, if the flow flows from the pipeline j into the node i, A in The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
(3) Pipeline flow-inflow node incidence matrix A out ：A out Comprising N h ×N b A square matrix, if the flow flows from the node i into the pipeline j, A out The square matrix of the ith row and the jth column in the matrix is a unit square matrix, otherwise, the square matrix is a unit zero square matrix;
step 202) combining a pipeline equivalent model, and pushing a dynamic source-load equation of the heat conducting network; based on the temperature vector T at the beginning of the pipeline 0 And a node temperature T n Through the association matrix A pn The association, expressed as:
T 0 ＝A pn T n (13)
the pipe temperature mixing equation is expressed as:
wherein T is n,i Representing the temperature of node i, T M,b Represents the end temperature of the pipe b, m b Representing the mass flow of the pipe b, E s,i Representing a set of pipes with node i as the head node, E e,i Representing a set of pipes ending with pipe i, V h Node outflow traffic for a node set of the heat supply network comprises two parts, namely node outflow to load and node outflow to pipeline, and the formula (14) is generalized to the system by using a topological equation, and is expressed as:
A out mT 0 +dT n ＝A in mT M (15)
wherein m is a flow vector in the dynamic model, and d is a node injection flow vector;
bringing the formula (12) and the formula (13) into the formula (15) in this order gives:
M out T n ＝JT n +B (16)
wherein M is out The method is characterized in that the method comprises the steps that an outflow flow matrix of a node in a dynamic model is obtained, J is a transmission matrix of different node temperatures in the dynamic model, and B is a temperature component acted by initial conditions; m is M ou J and B are respectively:
M out ＝diag(A out mA pn +d),J＝A in mαA pn ,B＝A in mβΤ 0 (17)
V h,sr and V h,ns Respectively a source node and a non-source node set of the thermodynamic system, then the formula (17) expands to:
in the method, in the process of the invention,and->Outflow traffic matrix, T, for the source node set and the non-source node set, respectively n,sr And T n,ns Temperature vector for source node and non-source node, J 11 、J 12 、J 21 And J 22 Respectively a blocking matrix in the coefficient matrix J, B sr And B ns Acting on temperature components of the source node and the non-source node for an initial condition;
according to the Gaussian elimination method, a source-load equation in a dynamic model of the heat supply network is obtained, and the equation is expressed as follows:
2. the method for dynamic safety control of an electric heating integrated energy system according to claim 1, wherein the step 30) is specifically as follows:
step 301) establishing a power grid sensitivity calculation type, including the node voltage amplitude, the sensitivity of the phase angle to the node injected active power and reactive power, and the sensitivity of the line transmission power to the node injected active power and reactive power, as shown in the formula (20) and the formula (21):
wherein U and theta are respectively the node voltage amplitude and phase angle, P and Q are respectively the node injected active power and reactive power, and P l And Q l Active power and reactive power respectively transmitted by power grid branches S UP Representing the sensitivity of the node voltage amplitude to the node injected active power, S θP Representing the sensitivity of the node voltage phase angle with respect to the node injected active power, S UQ Representing the sensitivity of the node voltage amplitude with respect to the node injection of reactive power, S θQ Representing the sensitivity of the node voltage phase angle with respect to the node injected reactive power;
step 302), according to a source-load equation in a heat supply network dynamic model, establishing a heat supply network sensitivity calculation type, wherein the heat supply network sensitivity calculation type comprises the sensitivity of water supply temperature at a source node relative to water supply temperature at a non-source node, the sensitivity of node water supply temperature relative to node backwater temperature, and the sensitivity of node water supply temperature relative to node heat power, which are respectively shown as a formula (22), a formula (23) and a formula (24):
in the method, in the process of the invention,represents the water supply temperature of node i at time k +.>Represents the backwater temperature of node i at time k, +.>Represents the thermal power of node i at time k, d i For the injection flow of node i, C ρ Is the specific heat capacity of working medium, V h,ld K is a set of load nodes in a heat supply network 1,ji,nk For K 1 The j-th row and the i-th column of the element of the nth row and the k-th column in the block matrix;
step 303) building 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 ＝C h -c u P CHP (25)
φ HP ＝c k P HP (26)
in phi CHP And P CHP Respectively representing the thermal power and the electric power output by the cogeneration equipment, C h C for maximum output thermal power of the cogeneration plant u For the heat-power ratio, phi, of a cogeneration plant HP And P HP C represents the thermal power output from the thermoelectric conversion device and the consumed electrical power, respectively k A thermoelectric ratio that is a thermoelectric conversion device;
obtaining the sensitivity of the thermal power of the internal node of the heat supply network to the thermal power of the load node according to the formula (22) and the formula (24); and combining the coupling equipment model to obtain the sensitivity 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, wherein the sensitivity is shown in the formula (28) and the formula (29):
in the method, in the process of the invention,represents the injected active power of node i at time k, P l k Representing the transmission power of the branch at time k, V cp Represents a coupling node set in the electric heating integrated energy system, eta represents a conversion coefficient of a coupling device, and eta= -1/c for an electric heating cogeneration device u The method comprises the steps of carrying out a first treatment on the surface of the For the thermoelectric conversion device, η=1/c k ；
Step 304) constructing a dynamic safety control method according to the coupling relation of the electric heating system, wherein the dynamic safety control method comprises the steps of 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, adjusting the voltage amplitude of the water supply temperature control node of the source node and the active power transmitted by the branch, and respectively inputting the voltage amplitude and the active power into formulas (30) to (33):
ΔP＝(U-U lim +dU)/S UP (30)
where Δp represents the control amount of the node injection power,indicating the control quantity of the water supply temperature of the node i at the moment k, U lim And P l,lim Representing node voltage amplitude threshold and branch transmission active power threshold, dU and dP l Is a margin for the node voltage magnitude.
CN202110828652.4A 2021-07-22 2021-07-22 Dynamic safety control method for electric heating comprehensive energy system Active CN113486532B (en)

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 CN113486532A (en) 2021-10-08
CN113486532B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114021490A (en) * 2021-10-25 2022-02-08 东南大学 Dynamic full-analysis method for hot water network

Citations (5)

* Cited by examiner, † Cited by third party
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 南京工业大学 Dynamic time interval comprehensive energy system scheduling method
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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 南京工业大学 Dynamic time interval comprehensive energy system scheduling method
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
CN113486532A (en) 2021-10-08

Similar Documents

Publication Publication Date Title
CN109242365B (en) Interval power flow calculation method of electricity-heat interconnection comprehensive energy system
CN109726483B (en) Radial heat supply network model of electric-heat interconnection comprehensive energy system and system thereof
CN112330127B (en) Static safety control method for multi-time-scale electric heating comprehensive energy system
CN111428351B (en) Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method
CN113486532B (en) Dynamic safety control method for electric heating comprehensive energy system
CN110707704A (en) Probability power flow analysis method of electric-thermal interconnection comprehensive energy system based on GMM and multi-point linear semi-invariant method
CN112257279B (en) Method for constructing feasible region of electric heating comprehensive energy system
Huang et al. A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale
CN110866213A (en) Multi-network steady-state energy flow analysis method and device for electricity-gas integrated energy system
CN113283077B (en) Correlation-considered power flow calculation method for electricity-heat-hydrogen comprehensive energy system interval
CN114221346A (en) Load flow calculation method of comprehensive energy system
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN110737993B (en) Multi-energy complementary system operation boundary analysis method considering load uncertainty
CN111310310B (en) Thermodynamic system static state trend rapid decoupling calculation method for quantity adjustment
CN113690891B (en) Analysis-method-based probability power flow determination method for electric-thermal interconnection comprehensive energy system
CN113221428B (en) Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system
CN113515853B (en) Optimal scheduling method of electrothermal interconnection comprehensive energy system based on linear equation
CN113111555B (en) Quality control thermodynamic system energy flow rapid calculation method based on superposition decoupling method
Zhu et al. Data-driven multi-energy flow calculation method for integrated energy system
Chen et al. Steady-state flow analysis of district heating system using the holomorphic embedding
Yao et al. A novel cross iteration method for dynamic energy flow calculation of the hot-water heating network in the integrated energy system
Qin et al. Heating network quasi-dynamic model of multi-energy flow system based on forward method
CN115906587B (en) Method for establishing dynamic coupling model of heat supply pipe network of comprehensive energy system
CN115906411B (en) Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics
CN112417698B (en) Dynamic two-port thermodynamic system model based on quality adjustment

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