CN114004110A - Quantizing event-driven simulation method for electricity-heat comprehensive energy system - Google Patents

Quantizing event-driven simulation method for electricity-heat comprehensive energy system Download PDF

Info

Publication number
CN114004110A
CN114004110A CN202111407795.4A CN202111407795A CN114004110A CN 114004110 A CN114004110 A CN 114004110A CN 202111407795 A CN202111407795 A CN 202111407795A CN 114004110 A CN114004110 A CN 114004110A
Authority
CN
China
Prior art keywords
variable
time
thermal
simulation
algebraic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111407795.4A
Other languages
Chinese (zh)
Other versions
CN114004110B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN202111407795.4A priority Critical patent/CN114004110B/en
Publication of CN114004110A publication Critical patent/CN114004110A/en
Application granted granted Critical
Publication of CN114004110B publication Critical patent/CN114004110B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12—Timing analysis or timing optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a quantization event-driven simulation method for an electro-thermal integrated energy system, which comprises the steps of firstly inputting simulation parameters aiming at the electro-thermal integrated energy system to be simulated, and carrying out simulation initialization; establishing an electric-thermal comprehensive energy system quantization coupling model to be simulated, wherein the electric-thermal comprehensive energy system quantization coupling model comprises a thermal subsystem model, an electric subsystem model and an electric-thermal coupling interface equation; making the current simulation time equal to the minimum value of the next simulation time in all variables, and integrating the corresponding variables to the current simulation time or updating the value of the current simulation time according to an equation and an integration method; and finally judging whether the simulation time reaches the simulation termination time. The invention is suitable for the dynamic simulation of the comprehensive energy system considering the electric-thermal coupling characteristic and the load dispersion characteristic, can simulate the multi-energy-flow coupling characteristic and the dispersion-continuous mixing characteristic, and provides an efficient simulation method for the electric-thermal comprehensive energy system model.

Description

Quantizing event-driven simulation method for electricity-heat comprehensive energy system
Technical Field
The invention relates to a quantization event driven simulation method. In particular to a quantization event driven simulation method for an electric-thermal integrated energy system.
Background
The comprehensive energy system integrates refrigeration, heat supply and power generation, improves the utilization rate of primary energy through comprehensive utilization of energy, and has great advantages in the aspect of reducing emission. With the rapid development of micro gas turbines and electric refrigeration/heating technologies, the coupling between electric energy and other energy forms is deepened continuously, and a comprehensive energy system integrating multiple energy forms such as cold, heat, electricity, gas and the like becomes an important development direction of a future smart power grid. The comprehensive energy system relates to deep fusion and close interaction of various energy sources, and a scientific and efficient simulation analysis method is a necessary premise for developing planning design and operation control of the comprehensive energy system and has important theoretical significance and practical value. Compared with the traditional power system simulation, the comprehensive energy system has the advantages that the inertia of energy forms such as cold, heat, gas and the like is strong, the time scale is large, a large number of discrete variables are brought by user behaviors on the load side, the discrete variables are frequently interacted with continuous state variables, and the complexity of the comprehensive energy system modeling simulation problem is greatly increased.
A Quantized State System (QSS) is an integration algorithm for solving a continuous-discrete hybrid system, and the QSS algorithm replaces time discretization of classical numerical integration by quantization of state variables, and the state variables of the system change in units of "quanta", and sequentially calculate the time required for each state variable change, thereby advancing integration. The event-driven framework can integrate a QSS method and a time domain integration algorithm, and has a good application prospect in the simulation of the comprehensive energy system. The method considers the electric-thermal coupling characteristics in the electric-thermal comprehensive energy system, accurately and efficiently describes the source-grid-load coupling influence, and is one of the important problems in the prior art. Therefore, the invention designs a quantization event-driven simulation method for the electric-thermal comprehensive energy system by considering the thermal and electric coupling characteristics and the load dispersion characteristics of the electric-thermal comprehensive energy system, and provides a necessary tool for the multi-energy flow high-efficiency simulation calculation of the comprehensive energy system.
Disclosure of Invention
The invention aims to solve the technical problem of providing a quantization event-driven simulation method for an electro-thermal comprehensive energy system, which can simulate multi-energy flow coupling characteristics and discrete-continuous hybrid characteristics.
The technical scheme adopted by the invention is as follows: a quantization event-driven simulation method for an electric-thermal integrated energy system comprises the following steps:
1) inputting the topological connection relation, element parameters, control parameters and simulation calculation parameters of a system aiming at the electricity-heat comprehensive energy system to be simulated, and carrying out simulation initialization; the control parameters comprise a heat source outlet temperature control interval and a user indoor temperature control interval, and the simulation calculation parameters comprise simulation termination time, a quantized integral threshold, an integral error tolerance, a voltage change threshold, initial values of various simulation variables and simulation variable historical quantities, initial values of first, second and third derivatives of thermal state variables, initial values of quantized thermal state variables and initial values of first and second derivatives; setting the last updating time and the next updating time of all variables as 0, and setting the current simulation time t as 0;
2) establishing a quantization coupling model of the electro-thermal integrated energy system to be simulated, wherein the quantization coupling model comprises a thermal link differential equation set, a thermal control link discrete algebraic equation set, an electric link differential algebraic equation set and a coupling interface algebraic equation set;
3) generating a simulation event schedule Tlist＝[th,s,th,d,te,the,teh]Wherein t ish,sIs the next update time of the thermodynamic state variable, th,dFor the next update time of the thermodynamic discrete algebraic variable, teIs in the form of electricityState variable xeAnd power algebraic variable yeNext update time theFor the next update time of the thermo-electrically coupled interface variable, tehThe next update time of the variable of the electric-thermal coupling interface; taking the simulation event time T which occurs first in the next step as min { T }list}；
4) If the current simulation time t is min { t ═ t }h,sIntegrating the differential equation set of the thermodynamic link to the current simulation time t by adopting a three-order quantization state system integration method, and updating the thermodynamic state variable xhNext update time th,sThermodynamic discrete algebraic variable uhNext update time th,dAnd step 5) is executed; otherwise, directly executing the step 5);
5) if the current simulation time t is equal to th,dCalculating the thermal discrete algebraic variable u of the current simulation time t according to the thermal control link discrete algebraic equation systemh(t) and updating the thermal discrete algebraic variable uhNext update time th,dThermo-electric coupling interface variable uheNext update time theExecuting step 6); otherwise, directly executing the step 6);
6) if the current simulation time t is equal to teIntegrating the differential algebraic equation set of the power link to the current simulation time t by adopting a variable step length integration algorithm, and updating the power state variable xeAnd power algebraic variable yeNext update time teVariable u of electro-thermal coupling interfaceehNext update time tehExecuting step 7); otherwise, directly executing the step 7);
7) if the current simulation time t is min { t ═ t }heCalculating a thermo-electric coupling interface variable u at the current simulation time t according to a coupling interface algebraic equation sethe(t) and updating the electro-thermal coupling interface variable uehNext update time theExecuting step 8); otherwise, directly executing the step 8);
8) if the current simulation time t is min { t ═ t }ehCalculating the variable u of the electro-thermal coupling interface at the current simulation time t according to the algebraic equation set of the coupling interfaceeh(t) and update ueh(t) Next replacementNew time tehExecuting step 9); otherwise, directly executing the step 9);
9) judging whether the current simulation moment T reaches the input simulation termination time T, if T is more than or equal to T, ending the simulation, and outputting a simulation result; otherwise, returning to the step 3).
The quantization event-driven simulation method for the electric-thermal integrated energy system greatly shortens the simulation time, has acceptable simulation precision, can efficiently simulate the electric-thermal integrated energy system with a discrete controller at the source load side, and has advantages in multi-time scale discrete-continuous hybrid system modeling simulation. The invention is suitable for the dynamic simulation of the comprehensive energy system considering the electric-thermal coupling characteristic and the load dispersion characteristic, can simulate the multi-energy-flow coupling characteristic and the dispersion-continuous mixing characteristic, and provides an efficient simulation method for the electric-thermal comprehensive energy system model.
Drawings
FIG. 1 is a flow chart of a quantization event driven simulation method for an electric-thermal integrated energy system according to the present invention;
FIG. 2 is a schematic topological diagram of a thermodynamic system;
FIG. 3 is a schematic view of a heat source and heat load configuration;
FIG. 4 is a graph of heat source number 2 outlet temperature versus time in an example of the invention;
FIG. 5 is a graph of the voltage at node 21 over time in an example of the present invention;
FIG. 6 is a graph of photovoltaic power output versus time for an example of the present invention;
FIG. 7 is a graph of photovoltaic power output versus time for an example of the present invention.
Detailed Description
The following describes a quantization event-driven simulation method for an electric-thermal integrated energy system according to the present invention in detail with reference to the following embodiments and the accompanying drawings.
As shown in fig. 1, the quantization event-driven simulation method for an electric-thermal integrated energy system of the present invention includes the following steps:
1) inputting the topological connection relation, element parameters, control parameters and simulation calculation parameters of a system aiming at the electricity-heat comprehensive energy system to be simulated, and carrying out simulation initialization; the control parameters comprise a heat source outlet temperature control interval and a user indoor temperature control interval, and the simulation calculation parameters comprise simulation termination time, a quantized integral threshold, an integral error tolerance, a voltage change threshold, initial values of various simulation variables and simulation variable historical quantities, initial values of first, second and third derivatives of thermal state variables, initial values of quantized thermal state variables and initial values of first and second derivatives; setting the last updating time and the next updating time of all variables as 0, and setting the current simulation time t as 0;
2) establishing a quantization coupling model of the electro-thermal integrated energy system to be simulated, wherein the quantization coupling model comprises a thermal link differential equation set, a thermal control link discrete algebraic equation set, an electric link differential algebraic equation set and a coupling interface algebraic equation set; wherein, the said:
(1) thermodynamic element differential equation set:
Figure BDA0003373010640000031
in the formula, xhThe heat state variables comprise pipeline temperature variables, heat source outlet temperature variables, user radiator outlet temperature variables and user indoor temperature variables; f. ofhThe function relationship of a thermodynamic link differential equation set; u. ofhThe system is a heat discrete algebraic variable and comprises the number of heat source heating units output by each heat source controller and connected to a power grid and the on-off state of a radiator output by each user temperature controller; u. ofehIs an electro-thermal coupling interface variable, is the thermal power of a single heating unit;
(2) a thermal control link discrete algebraic equation system:
uh＝zh(xh,uh-)
in the formula, zhFor the functional relationship of the discrete algebraic equation system of the thermodynamic control link, uh-Historical quantities of the thermal discrete algebraic variables;
(3) electric power link differential algebraic equation system:
Figure BDA0003373010640000032
in the formula (f)eIs a functional relation of a differential equation set of an electric power link, geIs a functional relation of an electric power link algebraic equation system, xeAs a power state variable, yeIs an algebraic variable of power uheThe impedance of a heat source connected to a power grid is a thermal-electric coupling interface variable;
(4) coupling interface algebraic equations:
Figure BDA0003373010640000033
in the formula, gehIs a function relation of an algebraic equation system of an electric-thermal coupling interface, gheIs a functional relation of a thermoelectric coupling algebraic equation system.
3) Generating a simulation event schedule Tlist＝[th,s,th,d,te,the,teh]Wherein t ish,sIs the next update time of the thermodynamic state variable, th,dFor the next update time of the thermodynamic discrete algebraic variable, teAs a power state variable xeAnd power algebraic variable yeNext update time theFor the next update time of the thermo-electrically coupled interface variable, tehThe next update time of the variable of the electric-thermal coupling interface; taking the simulation event time T which occurs first in the next step as min { T }list}；
4) If the current simulation time t is min { t ═ t }h,sIntegrating the differential equation set of the thermodynamic link to the current simulation time t by adopting a three-order quantization state system integration method, and updating the thermodynamic state variable xhNext update time th,sThermodynamic discrete algebraic variable uhNext update time th,dAnd step 5) is executed; otherwise, directly executing the step 5); wherein the content of the first and second substances,
(1) the method for integrating the thermodynamic element differential equation set to the current simulation time t by adopting a three-order quantization state system integration method specifically comprises the following steps:
(1.1) setting the next updating time t of the thermal state variableh,sThe minimum value is the m-th element and is recorded as
Figure BDA0003373010640000041
Is provided with
Figure BDA0003373010640000042
Representing thermodynamic state variable xhThe m-th element in the equation set of thermodynamic section differential
Figure BDA0003373010640000043
The set of equation sequence numbers for the terms is noted as Ωa(ii) a Order set omegab＝ΩaU{m}；
(1.2) is provided with
Figure BDA0003373010640000044
As a thermodynamic state variable xhThe j-th element in the list, j belongs to omegabUsing the following formula
Figure BDA0003373010640000045
That is to say, the
Figure BDA0003373010640000046
Integration to current simulation time t:
Figure BDA0003373010640000047
in the formula (I), the compound is shown in the specification,
Figure BDA0003373010640000048
is th,s,lThe jth element of (1), th,s,lAs a thermodynamic state variable xhThe time of the last update of the time,
Figure BDA0003373010640000049
Figure BDA00033730106400000410
and
Figure BDA00033730106400000411
are respectively as
Figure BDA00033730106400000412
Time of day
Figure BDA00033730106400000413
Value of and
Figure BDA00033730106400000414
the first, second, third derivative of (a).
(2) The updated thermodynamic state variable xhNext update time th,sThe method specifically comprises the following steps:
(2.1) calculating quantized thermodynamic State variables qhAnd q ishFirst and second derivatives of
Figure BDA00033730106400000415
Value of current simulation time t:
Figure BDA00033730106400000416
in the formula (I), the compound is shown in the specification,
Figure BDA00033730106400000417
are each qh、
Figure BDA00033730106400000418
The (c) th element of (a),
Figure BDA00033730106400000419
and
Figure BDA00033730106400000420
are respectively as
Figure BDA00033730106400000421
Time of day
Figure BDA00033730106400000422
Value of and
Figure BDA00033730106400000423
the values of the first and second derivatives of (a),
Figure BDA00033730106400000424
as a thermodynamic state variable xhLast update time th,s,lThe ith element in (1), NsIs the total number of thermodynamic state variables;
(2.2) setting th,sThe minimum value is the m-th element, and the quantized thermodynamic state variable is updated according to the following formula
Figure BDA00033730106400000425
And
Figure BDA00033730106400000426
derivative:
Figure BDA00033730106400000427
in the formula (I), the compound is shown in the specification,
Figure BDA00033730106400000428
respectively a quantized thermal state variable qhThermodynamic state variable xhThe m-th element of (a) to (b),
Figure BDA00033730106400000429
Figure BDA00033730106400000430
respectively at current simulation time t
Figure BDA00033730106400000431
Value of and
Figure BDA00033730106400000432
the values of the first and second derivatives of (a),
Figure BDA00033730106400000433
respectively at current simulation time t
Figure BDA00033730106400000434
Value of and
Figure BDA00033730106400000435
the value of the first and second derivatives of (a);
(2.3) assuming a thermodynamic discrete algebraic variable uhVariable u of electro-thermal coupling interfaceehKeeping the current value unchanged, and making j equal to omegabCalculating the jth element in the thermodynamic state variable
Figure BDA00033730106400000436
First, second and third derivatives at time t
Figure BDA00033730106400000437
Figure BDA00033730106400000438
In the formula (I), the compound is shown in the specification,
Figure BDA00033730106400000439
is composed of
Figure BDA00033730106400000440
The differential equation functional relationship of the thermodynamic link,
Figure BDA00033730106400000441
is composed of
Figure BDA00033730106400000442
First and second derivative expressions of, qh(t)、
Figure BDA00033730106400000443
Respectively quantized heat powerThe values u of the first and second derivatives of the state variable and the quantized thermal state variable at the current simulation time th(t) is the value of the thermal discrete algebraic variable at the current simulation time t, ueh(t) is the value of the variable of the electro-thermally coupled interface at the current simulation time t, omegabIs a set;
(2.4) applying the formula to j ∈ ΩbObtaining
Figure BDA0003373010640000051
Updating thermodynamic state variable xhNext update time th,s：
Figure BDA0003373010640000052
Wherein:
Figure BDA0003373010640000053
Figure BDA0003373010640000054
in the formula (I), the compound is shown in the specification,
Figure BDA0003373010640000055
is th,sIs the quantization integration threshold of the input,
Figure BDA0003373010640000056
as a thermodynamic state variable xhNext update time th,sThe (c) th element of (a),
Figure BDA0003373010640000057
for quantizing thermal state variables qhThe (j) th element of (a),
Figure BDA0003373010640000058
is composed of
Figure BDA0003373010640000059
Time of day
Figure BDA00033730106400000510
The value of (a) is,
Figure BDA00033730106400000511
respectively at current simulation time t
Figure BDA00033730106400000512
Value of and
Figure BDA00033730106400000513
the values of the first and second derivatives of (a),
Figure BDA00033730106400000514
as a thermodynamic state variable xhThe (j) th element of (a),
Figure BDA00033730106400000515
is composed of
Figure BDA00033730106400000516
Time of day
Figure BDA00033730106400000517
The value of (a) is,
Figure BDA00033730106400000518
respectively at current simulation time t
Figure BDA00033730106400000519
Value of and
Figure BDA00033730106400000520
the first, second, third derivative values of;
(2.5) updating the thermodynamic State variable xhLast update time th,s,lInstant command
Figure BDA00033730106400000521
(3) The updated heat discrete algebraic variable uhNext update time th,dThe method specifically comprises the following steps:
let omegacRepresenting thermodynamic state variable xhThe serial number set of the medium heat source outlet temperature variable and the user indoor temperature variable is in the range of omegab∩ΩcIf not, using the p-th element in the thermodynamic state variable
Figure BDA00033730106400000522
To calculate th,dAnd taking the minimum value obtained by solving as the final th,d：
Upper limit event localization equation:
Figure BDA00033730106400000523
lower bound event localization equation:
Figure BDA00033730106400000524
in the formula, th,dAs a discrete algebraic variable u of heathThe next time of update of the time of day,
Figure BDA00033730106400000525
and
Figure BDA00033730106400000526
respectively being the p-th element in the thermodynamic state variable at the current simulation time t
Figure BDA00033730106400000527
Value of and
Figure BDA00033730106400000528
the values of the first, second and third derivatives of (a),
Figure BDA00033730106400000529
Figure BDA00033730106400000530
to represent
Figure BDA00033730106400000531
The upper limit value and the lower limit value of the control interval of the outlet temperature of the heat source or the indoor temperature of the user; if the upper limit event positioning equation and the lower limit event positioning equation are both solved, t is takenh,dAnd if not, taking the minimum value obtained by solving as Th,d。
5) If the current simulation time t is equal to th,dCalculating the thermal discrete algebraic variable u of the current simulation time t according to the thermal control link discrete algebraic equation systemh(t) and updating the thermal discrete algebraic variable uhNext update time th,dThermo-electric coupling interface variable uheNext update time theExecuting step 6); otherwise, directly executing the step 6); the heat power discrete algebraic variable u of the current simulation time t is calculated according to the heat power control link discrete algebraic equation seth(t) and updating the thermal discrete algebraic variable uhNext update time th,dThermo-electric coupling interface variable uheNext update time theThe method specifically comprises the following steps:
(1) calculating a thermal discrete algebraic variable u at the current simulation time th(t)：
uh(t)＝zh(xh(th,s,l),uh-)
In the formula, xhAs a thermodynamic state variable, th,s,lIs xhLast update time of (x)h(th,s,l) Is xhAt th,s,lValue of time uh-As a discrete algebraic variable u of heathA historical amount of (c);
(2) updating a thermodynamic discrete algebraic variable uhHistory amount u ofh-＝uh(t)；
(3) Updating a thermodynamic discrete algebraic variable uhNext update time th,dT is the input simulation termination time;
(4) updating thermo-electric coupling interface variable uheNext update time the＝t。
6) If the current simulation time t is equal to teIntegrating the differential algebraic equation set of the power link to the current simulation time t by adopting a variable step length integration algorithm, and updating the power state variable xeAnd power algebraic variable yeNext update time teVariable u of electro-thermal coupling interfaceehNext update time tehExecuting step 7); otherwise, directly executing the step 7); the electric power state variable xeAnd power algebraic variable yeNext update time teVariable u of electro-thermal coupling interfaceehNext update time tehThe method specifically comprises the following steps:
(1) calculating to obtain the next integral step length delta t according to the integral error tolerance by adopting a variable step length integral algorithm, and enabling the power state variable xeAnd power algebraic variable yeNext update time te＝t+Δt；
(2) Updating an electro-thermal coupling interface variable uehNext update time teh＝t。
7) If the current simulation time t is min { t ═ t }heCalculating a thermo-electric coupling interface variable u at the current simulation time t according to a coupling interface algebraic equation sethe(t) and updating the electro-thermal coupling interface variable uehNext update time theExecuting step 8); otherwise, directly executing the step 8); calculating the thermo-electric coupling interface variable u of the current simulation moment t according to the coupled interface algebraic equation sethe(t) and updating the electro-thermal coupling interface variable uehNext update time theThe method specifically comprises the following steps:
(1) calculating the thermo-electric interface coupling variable u of the current simulation time t by the following formulahe(t)：
uhe(t)＝ghe(uh(t))
In the formula, gheFor the algebraic equation system functional relation of the thermo-electric coupling interface, uh(t) is the value of the thermal discrete algebraic variable at the current simulation time t;
(2) the following threshold condition judgment is carried out:
Figure BDA0003373010640000061
in the formula (I), the compound is shown in the specification,
Figure BDA0003373010640000062
as a discrete algebraic variable u of heathThe (c) th element of (a),
Figure BDA0003373010640000063
is composed of
Figure BDA0003373010640000064
At the value of the current simulation time t,
Figure BDA0003373010640000065
for last time the threshold condition is satisfied
Figure BDA0003373010640000066
Historical value, ΩdIs uhThe number set of the number variables of the middle heat source heating units is accessed; if the threshold condition is satisfied, the power state variable x is seteAnd power algebraic variable yeNext update time te＝t、
Figure BDA0003373010640000067
Otherwise, keep te、
Figure BDA0003373010640000068
The change is not changed;
(3) put hot-electric coupling interface variable uheNext update time theT is the simulation termination time of the input.
8) If the current simulation time t is min { t ═ t }ehCalculating the variable u of the electro-thermal coupling interface at the current simulation time t according to the algebraic equation set of the coupling interfaceeh(t) and update ueh(t) the next update time tehExecuting step 9); otherwise, directly executing the step 9); the electric-thermal coupling interface at the current simulation moment t is calculated according to the algebraic equation set of the coupling interfaceVariable ueh(t), specifically including:
(1) calculating the variable u of the electro-thermal coupling interface at the current simulation time t by the following formulaeh(t)：
ueh(t)＝geh(ye(t))
In the formula, gehIs a function relation of an algebraic equation system of an electric-thermal coupling interface, ye(t) is the value of the power algebraic variable at the current simulation time t;
(2) the following threshold condition judgment is carried out:
Figure BDA0003373010640000069
in the formula (I), the compound is shown in the specification,
Figure BDA00033730106400000610
for algebraicizing the variable y for the powereThe w-th element of (a) is,
Figure BDA00033730106400000611
is composed of
Figure BDA00033730106400000612
At the value of the current simulation time t,
Figure BDA00033730106400000613
for last time the threshold condition is satisfied
Figure BDA00033730106400000614
Historical value, Δ yeFor the input voltage variation threshold, ΩeA power algebraic variable serial number set corresponding to a power grid node accessed by a heat source; if the threshold condition is satisfied, setting
Figure BDA00033730106400000615
Figure BDA00033730106400000616
As the next update time t of the thermodynamic state variableh,sThe k-th element, k being
Figure BDA0003373010640000071
Corresponding to the outlet temperature variable of the heat source in the electric power state variable xeThe serial number in (1); otherwise, it keeps
Figure BDA0003373010640000072
The change is not changed;
(3) variable u of electro-thermal coupling interfaceehNext update time tehT is the simulation termination time of the input.
9) Judging whether the current simulation moment T reaches the input simulation termination time T, if T is more than or equal to T, ending the simulation, and outputting a simulation result; otherwise, returning to the step 3).
Specific examples are given below:
the example is based on an MATLAB programming language environment, realizes a quantization event-driven simulation method for an electro-thermal comprehensive energy system, and verifies and analyzes the method through a heating system example. The hardware platform of the simulation test is a 4-Core PC machine with Intel Core (TM) i7-8700 CPU @3.20GHz and 8GB RAM; the software environment is a 64-bit Windows10 operating system.
The electricity-heat comprehensive energy thermodynamic system consists of a heat supply network and a power grid. The heat supply system is divided into a water supply network and a water return network, the water supply network and the water return network are completely symmetrical, corresponding pipeline branch parameters are the same, a heat source branch and a heat user branch are connected with the water supply network and the water return network to form a closed hot working medium flowing loop, the water supply network and the water return network are composed of pipeline branches and nodes, the nodes of the water supply network are numbered from 1 to 32, the nodes in the water return network are numbered from 33 to 64, and the sequence is consistent with that of the water supply network; the number of the water supply network pipeline branch is 1-32, the number of the water return network pipeline branch is 33-64, and the sequence is consistent with that of the water supply network; the node and branch numbering sequence of the backwater net is consistent with that of the water supply net, is not listed in the figure, and can be easily deduced. The serial numbers of the heat source branch and the heat user branch are 65-88. The system topology is shown in fig. 1, and the No. 65 branch is the position of the No. 1 heat source; the No. 87 branch is the position where the No. 2 heat source is located; branch 88 is where heat source number 3 is located. No. 1 ~ 21 heat loads are located No. 66 ~ 86 branches respectively, and every branch road in No. 66 ~ 86 branches all is 10 buildings heat supply. The power grid adopts an IEEE-33 node power distribution network calculation example, heat sources 1, 2 and 3 are respectively connected to nodes 18, 21 and 31 of the power distribution network, 8 distributed photovoltaic power sources are also connected to the power grid, and the specific positions are shown in figure 2.
Setting the simulation initial time T as 0, the simulation end time T as 3600 s, and the quantization integration threshold value delta Q as 1 multiplied by 10-6Integral error margin of 1 × 10-3Voltage change threshold of 1 × 10-6Initial value x of thermodynamic state variableh(0) Initial value q of quantized thermal state variableh(0)＝xh(0) First, second and third derivative initial values of thermodynamic state variables
Figure BDA0003373010640000073
Figure BDA0003373010640000074
First-order and second-order derivative initial values of quantized thermal state variable initial values
Figure BDA0003373010640000075
Initial value u of heat discrete algebraic variableh(0) Initial value y of power algebraic variablee(0) Historical values u of variablesh-＝uh(0)，u'h＝uh(0)，y'e＝ye(0) Setting simulation scenes as follows:
the indoor temperature control range of the user is 24-25 ℃, and the temperature control range of the heat source outlet is 80-82 ℃. Three methods are adopted to compare simulation accuracy, and the reference method is as follows: the electric-thermal comprehensive energy system adopts a backward differentiation equation (BDF) integration method to solve simultaneously, and the method of the invention comprises the following steps: a quantized event driven simulation method is characterized in that a thermal link is solved by adopting a third-order quantized state system (QSS 3) integration method, and an electric link is solved by adopting a BDF integration method.
FIG. 3 is a schematic view of a heat source and a heat load configuration, in which a temperature controller of the heat source monitors an outlet temperature thereof and controls the number of connected heating units of the heat source to adjust the outlet temperature of the heat source, and a temperature controller in a user building adjusts an indoor temperature by outputting an on-off state of a radiator; FIG. 4 is a graph of the outlet temperature of heat source No. 2 as a function of time, with heat source No. 2 turning off one heating unit and the temperature dropping when the temperature reaches a set upper limit; fig. 5 is a graph showing the voltage change of the node 21 with time, and after the heat source No. 2 is turned off one heating unit, the consumed electric power is reduced, so that the voltage of the node of the grid No. 21 connected with the heat source No. 2 is increased; FIG. 6 is a graph of output of a photovoltaic power source over time; fig. 7 shows the indoor temperature and the average value thereof, the indoor temperature of all buildings is set to 24-25 ℃, and the switch of the radiator is controlled to keep the indoor temperature within the interval. The simulation efficiency is shown in table 1, and as can be seen from table 1, compared with the BDF algorithm, the quantization event-driven simulation method for the electric-thermal integrated energy system greatly shortens the simulation time, has acceptable simulation accuracy, can efficiently simulate the electric-thermal integrated energy system with a discrete controller on the source charge side, and has advantages in multi-time scale discrete-continuous hybrid system modeling simulation.
TABLE 1 simulation efficiency
Method Maximum relative error (%) Simulation time(s)
The method of the invention 1.2×10-5 56.7
Datum - >1h

Claims (9)

1. A quantization event-driven simulation method for an electric-thermal integrated energy system is characterized by comprising the following steps:
1) inputting the topological connection relation, element parameters, control parameters and simulation calculation parameters of a system aiming at the electricity-heat comprehensive energy system to be simulated, and carrying out simulation initialization; the control parameters comprise a heat source outlet temperature control interval and a user indoor temperature control interval, and the simulation calculation parameters comprise simulation termination time, a quantized integral threshold, an integral error tolerance, a voltage change threshold, initial values of various simulation variables and simulation variable historical quantities, initial values of first, second and third derivatives of thermal state variables, initial values of quantized thermal state variables and initial values of first and second derivatives; setting the last updating time and the next updating time of all variables as 0, and setting the current simulation time t as 0;
2) establishing a quantization coupling model of the electro-thermal integrated energy system to be simulated, wherein the quantization coupling model comprises a thermal link differential equation set, a thermal control link discrete algebraic equation set, an electric link differential algebraic equation set and a coupling interface algebraic equation set;
3) generating a simulation event schedule Tlist＝[th,s,th,d,te,the,teh]Wherein t ish,sIs the next update time of the thermodynamic state variable, th,dFor the next update time of the thermodynamic discrete algebraic variable, teAs a power state variable xeAnd power algebraic variable yeNext update time theFor the next update time of the thermo-electrically coupled interface variable, tehThe next update time of the variable of the electric-thermal coupling interface; taking the simulation event time T which occurs first in the next step as min { T }list}；
4) If the current simulation time t is min { t ═ t }h,sIntegrating the differential equation set of the thermodynamic link to the current simulation time t by adopting a three-order quantization state system integration method, and updating the thermodynamic state variable xhNext time of update ofTime th,sThermodynamic discrete algebraic variable uhNext update time th,dAnd step 5) is executed; otherwise, directly executing the step 5);
5) if the current simulation time t is equal to th,dCalculating the thermal discrete algebraic variable u of the current simulation time t according to the thermal control link discrete algebraic equation systemh(t) and updating the thermal discrete algebraic variable uhNext update time th,dThermo-electric coupling interface variable uheNext update time theExecuting step 6); otherwise, directly executing the step 6);
6) if the current simulation time t is equal to teIntegrating the differential algebraic equation set of the power link to the current simulation time t by adopting a variable step length integration algorithm, and updating the power state variable xeAnd power algebraic variable yeNext update time teVariable u of electro-thermal coupling interfaceehNext update time tehExecuting step 7); otherwise, directly executing the step 7);
7) if the current simulation time t is min { t ═ t }heCalculating a thermo-electric coupling interface variable u at the current simulation time t according to a coupling interface algebraic equation sethe(t) and updating the electro-thermal coupling interface variable uehNext update time theExecuting step 8); otherwise, directly executing the step 8);
8) if the current simulation time t is min { t ═ t }ehCalculating the variable u of the electro-thermal coupling interface at the current simulation time t according to the algebraic equation set of the coupling interfaceeh(t) and update ueh(t) the next update time tehExecuting step 9); otherwise, directly executing the step 9);
9) judging whether the current simulation moment T reaches the input simulation termination time T, if T is more than or equal to T, ending the simulation, and outputting a simulation result; otherwise, returning to the step 3).
2. The method for the quantized event-driven simulation of the electric-thermal integrated energy system according to claim 1, wherein the step 2) comprises:
(1) thermodynamic element differential equation set:
Figure FDA0003373010630000011
in the formula, xhThe heat state variables comprise pipeline temperature variables, heat source outlet temperature variables, user radiator outlet temperature variables and user indoor temperature variables; f. ofhThe function relationship of a thermodynamic link differential equation set; u. ofhThe system is a heat discrete algebraic variable and comprises the number of heat source heating units output by each heat source controller and connected to a power grid and the on-off state of a radiator output by each user temperature controller; u. ofehIs an electro-thermal coupling interface variable, is the thermal power of a single heating unit;
(2) a thermal control link discrete algebraic equation system:
uh＝zh(xh,uh-)
in the formula, zhFor the functional relationship of the discrete algebraic equation system of the thermodynamic control link, uh-Historical quantities of the thermal discrete algebraic variables;
(3) electric power link differential algebraic equation system:
Figure FDA0003373010630000021
in the formula (f)eIs a functional relation of a differential equation set of an electric power link, geIs a functional relation of an electric power link algebraic equation system, xeAs a power state variable, yeIs an algebraic variable of power uheThe impedance of a heat source connected to a power grid is a thermal-electric coupling interface variable;
(4) coupling interface algebraic equations:
Figure FDA0003373010630000022
in the formula, gehIs a function relation of an algebraic equation system of an electric-thermal coupling interface, gheIs a functional relation of a thermoelectric coupling algebraic equation system.
3. The method for the quantization event-driven simulation of the electric-thermal energy integration system according to claim 1, wherein the step 4) of integrating the thermodynamic element differential equation set to the current simulation time t by using a third-order quantization state system integration method specifically comprises:
(1) setting the next update time t of the thermal state variableh,sThe minimum value is the m-th element and is recorded as
Figure FDA0003373010630000023
Is provided with
Figure FDA0003373010630000024
Representing thermodynamic state variable xhThe m-th element in the equation set of thermodynamic section differential
Figure FDA0003373010630000025
The set of equation sequence numbers for the terms is noted as Ωa(ii) a Order set omegab＝Ωa∪{m}；
(2) Is provided with
Figure FDA0003373010630000026
As a thermodynamic state variable xhThe j-th element in the list, j belongs to omegabUsing the following formula
Figure FDA0003373010630000027
That is to say, the
Figure FDA0003373010630000028
Integration to current simulation time t:
Figure FDA0003373010630000029
in the formula (I), the compound is shown in the specification,
Figure FDA00033730106300000210
is th,s,lThe jth element of (1), th,s,lAs a thermodynamic state variable xhThe time of the last update of the time,
Figure FDA00033730106300000211
Figure FDA00033730106300000212
and
Figure FDA00033730106300000213
are respectively as
Figure FDA00033730106300000214
Time of day
Figure FDA00033730106300000215
Value of and
Figure FDA00033730106300000216
the first, second, third derivative of (a).
4. The electric-thermal integrated energy system oriented quantization event-driven simulation method as claimed in claim 1, wherein the step 4) of updating the thermodynamic state variable xhNext update time th,sThe method specifically comprises the following steps:
(1) calculating quantized thermodynamic state variables qhAnd q ishFirst and second derivatives of
Figure FDA00033730106300000217
Value of current simulation time t:
Figure FDA0003373010630000031
in the formula (I), the compound is shown in the specification,
Figure FDA0003373010630000032
are each qh、
Figure FDA0003373010630000033
The (c) th element of (a),
Figure FDA0003373010630000034
and
Figure FDA0003373010630000035
are respectively as
Figure FDA0003373010630000036
Time of day
Figure FDA0003373010630000037
Value of and
Figure FDA0003373010630000038
the values of the first and second derivatives of (a),
Figure FDA0003373010630000039
as a thermodynamic state variable xhLast update time th,s,lThe ith element in (1), NsIs the total number of thermodynamic state variables;
(2) let th,sThe minimum value is the m-th element, and the quantized thermodynamic state variable is updated according to the following formula
Figure FDA00033730106300000310
And
Figure FDA00033730106300000311
derivative:
Figure FDA00033730106300000312
in the formula (I), the compound is shown in the specification,
Figure FDA00033730106300000313
respectively a quantized thermal state variable qhThermodynamic state variable xhThe m-th element of (a) to (b),
Figure FDA00033730106300000314
Figure FDA00033730106300000315
respectively at current simulation time t
Figure FDA00033730106300000316
Value of and
Figure FDA00033730106300000317
the values of the first and second derivatives of (a),
Figure FDA00033730106300000318
Figure FDA00033730106300000319
respectively at current simulation time t
Figure FDA00033730106300000320
Value of and
Figure FDA00033730106300000321
the value of the first and second derivatives of (a);
(3) assuming a thermodynamic discrete algebraic variable uhVariable u of electro-thermal coupling interfaceehKeeping the current value unchanged, and making j equal to omegabCalculating the jth element in the thermodynamic state variable
Figure FDA00033730106300000322
First, second and third derivatives at time t
Figure FDA00033730106300000323
Figure FDA00033730106300000324
In the formula (I), the compound is shown in the specification,
Figure FDA00033730106300000325
is composed of
Figure FDA00033730106300000326
The differential equation functional relationship of the thermodynamic link,
Figure FDA00033730106300000327
is composed of
Figure FDA00033730106300000328
First and second derivative expressions of, qh(t)、
Figure FDA00033730106300000329
Respectively the values u of the quantized thermal state variable and the first and second derivatives of the quantized thermal state variable at the current simulation time th(t) is the value of the thermal discrete algebraic variable at the current simulation time t, ueh(t) is the value of the variable of the electro-thermally coupled interface at the current simulation time t, omegabIs a set;
(4) is represented by the following formula to j ∈ omegabObtaining
Figure FDA00033730106300000330
Updating thermodynamic state variable xhNext update time th,s：
Figure FDA00033730106300000331
Wherein:
Figure FDA00033730106300000332
Figure FDA00033730106300000333
in the formula (I), the compound is shown in the specification,
Figure FDA00033730106300000334
is th,sIs the quantization integration threshold of the input,
Figure FDA00033730106300000335
as a thermodynamic state variable xhNext update time th,sThe (c) th element of (a),
Figure FDA00033730106300000336
for quantizing thermal state variables qhThe (j) th element of (a),
Figure FDA00033730106300000337
is composed of
Figure FDA00033730106300000338
Time of day
Figure FDA00033730106300000339
The value of (a) is,
Figure FDA00033730106300000340
respectively at current simulation time t
Figure FDA00033730106300000341
Value of and
Figure FDA00033730106300000342
the values of the first and second derivatives of (a),
Figure FDA00033730106300000343
as a thermodynamic state variable xhThe (j) th element of (a),
Figure FDA00033730106300000344
is composed of
Figure FDA00033730106300000345
Time of day
Figure FDA00033730106300000346
The value of (a) is,
Figure FDA00033730106300000347
respectively at current simulation time t
Figure FDA00033730106300000348
Value of and
Figure FDA00033730106300000349
the first, second, third derivative values of;
(5) updating thermodynamic state variable xhLast update time th,s,lInstant command
Figure FDA00033730106300000350
5. The method as claimed in claim 1, wherein the step 4) of updating the discrete algebraic variables u is implemented by using a quantization event-driven simulation method for the electric-thermal integrated energy systemhNext update time th,dThe method specifically comprises the following steps:
let omegacRepresenting thermodynamic state variable xhThe serial number set of the medium heat source outlet temperature variable and the user indoor temperature variable is in the range of omegab∩ΩcIf not, using the p-th element in the thermodynamic state variable
Figure FDA0003373010630000041
To calculate th,dAnd taking the minimum value obtained by solving as the final th,d：
Upper limit event localization equation:
Figure FDA0003373010630000042
lower bound event localization equation:
Figure FDA0003373010630000043
in the formula, th,dAs a discrete algebraic variable u of heathThe next time of update of the time of day,
Figure FDA0003373010630000044
and
Figure FDA0003373010630000045
respectively being the p-th element in the thermodynamic state variable at the current simulation time t
Figure FDA0003373010630000046
Value of and
Figure FDA0003373010630000047
the values of the first, second and third derivatives of (a),
Figure FDA0003373010630000048
Figure FDA0003373010630000049
to represent
Figure FDA00033730106300000410
The upper limit value and the lower limit value of the control interval of the outlet temperature of the heat source or the indoor temperature of the user; if the upper limit event positioning equation and the lower limit event positioning equation are bothIf there is no solution, then get th,dAnd if not, taking the minimum value obtained by solving as Th,d。
6. The method as claimed in claim 1, wherein the step 5) of calculating the discrete algebraic variant u of thermal power at the current simulation time t according to the discrete algebraic equation system of thermal control linkh(t) and updating the thermal discrete algebraic variable uhNext update time th,dThermo-electric coupling interface variable uheNext update time theThe method specifically comprises the following steps:
(1) calculating a thermal discrete algebraic variable u at the current simulation time th(t)：
uh(t)＝zh(xh(th,s,l),uh-)
In the formula, xhAs a thermodynamic state variable, th,s,lIs xhLast update time of (x)h(th,s,l) Is xhAt th,s,lValue of time uh-As a discrete algebraic variable u of heathA historical amount of (c);
(2) updating a thermodynamic discrete algebraic variable uhHistory amount u ofh-＝uh(t)；
(3) Updating a thermodynamic discrete algebraic variable uhNext update time th,dT is the input simulation termination time;
(4) updating thermo-electric coupling interface variable uheNext update time the＝t。
7. The electric-thermal integrated energy system oriented quantization event-driven simulation method as claimed in claim 1, wherein the electric power state variable x in step 6)eAnd power algebraic variable yeNext update time teVariable u of electro-thermal coupling interfaceehNext update time tehThe method specifically comprises the following steps:
(1) by integrating with the applied step-sizeThe algorithm calculates according to the integral error tolerance to obtain the next integral step length delta t and enables the power state variable xeAnd power algebraic variable yeNext update time te＝t+Δt；
(2) Updating an electro-thermal coupling interface variable uehNext update time teh＝t。
8. The method as claimed in claim 1, wherein the step 7) of calculating the thermo-electric coupling interface variable u at the current simulation time t according to the coupled interface algebraic equation systemhe(t) and updating the electro-thermal coupling interface variable uehNext update time theThe method specifically comprises the following steps:
(1) calculating the thermo-electric interface coupling variable u of the current simulation time t by the following formulahe(t)：
uhe(t)＝ghe(uh(t))
In the formula, gheFor the algebraic equation system functional relation of the thermo-electric coupling interface, uh(t) is the value of the thermal discrete algebraic variable at the current simulation time t;
(2) the following threshold condition judgment is carried out:
Figure FDA0003373010630000051
in the formula (I), the compound is shown in the specification,
Figure FDA0003373010630000052
as a discrete algebraic variable u of heathThe (c) th element of (a),
Figure FDA0003373010630000053
is composed of
Figure FDA0003373010630000054
At the value of the current simulation time t,
Figure FDA0003373010630000055
for last time the threshold condition is satisfied
Figure FDA0003373010630000056
Historical value, ΩdIs uhThe number set of the number variables of the middle heat source heating units is accessed; if the threshold condition is satisfied, the power state variable x is seteAnd power algebraic variable yeNext update time te＝t、
Figure FDA0003373010630000057
Otherwise, keep te、
Figure FDA0003373010630000058
The change is not changed;
(3) put hot-electric coupling interface variable uheNext update time theT is the simulation termination time of the input.
9. The method for quantization event-driven simulation of an electric-thermal energy integration system as claimed in claim 1, wherein the step 8) is performed to calculate the variables u of the electric-thermal coupling interface at the current simulation time t according to the algebraic equation system of the coupling interfaceeh(t), specifically including:
(1) calculating the variable u of the electro-thermal coupling interface at the current simulation time t by the following formulaeh(t)：
ueh(t)＝geh(ye(t))
In the formula, gehIs a function relation of an algebraic equation system of an electric-thermal coupling interface, ye(t) is the value of the power algebraic variable at the current simulation time t;
(2) the following threshold condition judgment is carried out:
Figure FDA0003373010630000059
in the formula (I), the compound is shown in the specification,
Figure FDA00033730106300000510
for algebraicizing the variable y for the powereThe w-th element of (a) is,
Figure FDA00033730106300000511
is composed of
Figure FDA00033730106300000512
At the value of the current simulation time t,
Figure FDA00033730106300000513
for last time the threshold condition is satisfied
Figure FDA00033730106300000514
Historical value, Δ yeFor the input voltage variation threshold, ΩeA power algebraic variable serial number set corresponding to a power grid node accessed by a heat source; if the threshold condition is satisfied, setting
Figure FDA00033730106300000515
Figure FDA00033730106300000516
As the next update time t of the thermodynamic state variableh,sThe k-th element, k being
Figure FDA00033730106300000517
Corresponding to the outlet temperature variable of the heat source in the electric power state variable xeThe serial number in (1); otherwise, it keeps
Figure FDA00033730106300000518
The change is not changed;
(3) variable u of electro-thermal coupling interfaceehNext update time tehT is the simulation termination time of the input.
CN202111407795.4A 2021-11-24 2021-11-24 Quantized event driven simulation method for electric-thermal integrated energy system Active CN114004110B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111407795.4A CN114004110B (en) 2021-11-24 2021-11-24 Quantized event driven simulation method for electric-thermal integrated energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111407795.4A CN114004110B (en) 2021-11-24 2021-11-24 Quantized event driven simulation method for electric-thermal integrated energy system

Publications (2)

Publication Number Publication Date
CN114004110A true CN114004110A (en) 2022-02-01
CN114004110B CN114004110B (en) 2024-04-12

Family

ID=79930210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111407795.4A Active CN114004110B (en) 2021-11-24 2021-11-24 Quantized event driven simulation method for electric-thermal integrated energy system

Country Status (1)

Country Link
CN (1) CN114004110B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113949059A (en) * 2021-09-28 2022-01-18 昆明理工大学 Generalized discrete equivalent model building method of comprehensive load

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010089253A1 (en) * 2009-02-03 2010-08-12 Dong Energy Power A/S Distributed electrical power production system and method of control thereof
CN109193626A (en) * 2018-09-15 2019-01-11 天津大学 Unified solving method for model selection and constant volume optimization planning of distributed energy station
CN109541961A (en) * 2018-10-23 2019-03-29 清华大学 For the discrete state event-driven simulation method of power electronics hybrid system emulation
CN110555264A (en) * 2019-08-30 2019-12-10 南方电网科学研究院有限责任公司 Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
CN113095538A (en) * 2020-06-08 2021-07-09 华北电力大学 Flexible operation-oriented wide-load operation dynamic characteristic modeling method for cogeneration unit

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010089253A1 (en) * 2009-02-03 2010-08-12 Dong Energy Power A/S Distributed electrical power production system and method of control thereof
CN109193626A (en) * 2018-09-15 2019-01-11 天津大学 Unified solving method for model selection and constant volume optimization planning of distributed energy station
CN109541961A (en) * 2018-10-23 2019-03-29 清华大学 For the discrete state event-driven simulation method of power electronics hybrid system emulation
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
CN110555264A (en) * 2019-08-30 2019-12-10 南方电网科学研究院有限责任公司 Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method
CN113095538A (en) * 2020-06-08 2021-07-09 华北电力大学 Flexible operation-oriented wide-load operation dynamic characteristic modeling method for cogeneration unit

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113949059A (en) * 2021-09-28 2022-01-18 昆明理工大学 Generalized discrete equivalent model building method of comprehensive load

Also Published As

Publication number Publication date
CN114004110B (en) 2024-04-12

Similar Documents

Publication Publication Date Title
CN113902040B (en) Method, system, equipment and storage medium for coordinating and optimizing electricity-heat comprehensive energy system
CN109270842B (en) Bayesian network-based regional heat supply model prediction control system and method
Bashash et al. Modeling and control of aggregate air conditioning loads for robust renewable power management
CN103591637B (en) A kind of central heating secondary network runing adjustment method
Homod Analysis and optimization of HVAC control systems based on energy and performance considerations for smart buildings
Zong et al. Application of model predictive control for active load management in a distributed power system with high wind penetration
Li et al. Modeling for indoor temperature prediction based on time-delay and Elman neural network in air conditioning system
Scherer et al. Efficient building energy management using distributed model predictive control
CN104573851B (en) A kind of building hourly load prediction technique based on meteorological temperature forecast
CN111046594B (en) Hot water heating network dynamic simulation method based on cross iteration principle
CN110571865B (en) Active power distribution network source load interaction decision method based on consistency algorithm
CN113203187B (en) Building heating ventilation air conditioning load optimization control method based on partial linear model
CN113190999A (en) Electric heating coordination method and device for regulating flow of heat supply pipe network to improve wind power consumption
Zhang et al. Experimental study on a novel fuzzy control method for static pressure reset based on the maximum damper position feedback
CN114004110A (en) Quantizing event-driven simulation method for electricity-heat comprehensive energy system
Gouadria et al. Comparison between self-tuning fuzzy PID and classic PID controllers for greenhouse system
Xu et al. Optimal building energy management using intelligent optimization
Berouine et al. A predictive control approach for thermal energy management in buildings
CN114282779A (en) Intelligent power grid dispatching method based on economic dispatching target and power parameter
CN111783309A (en) Dynamic simulation method of steam heating network based on internal conservation
CN113065249A (en) Method and device for predicting supply and return water temperature of heating system
CN114114958B (en) Discrete simulation method for state of heating system by considering control characteristics of source load temperature
Li et al. Adaptive predictive control method for improving control stability of air-conditioning terminal in public buildings
de Oliveira et al. Dynamic online optimization of a house heating system in a fluctuating energy price scenario
Li et al. Thermal comfort control based on MEC algorithm for HVAC systems

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