CN115906411A - Electric heating comprehensive energy system optimal energy flow modeling method and system considering full dynamic - Google Patents

Electric heating comprehensive energy system optimal energy flow modeling method and system considering full dynamic Download PDF

Info

Publication number
CN115906411A
CN115906411A CN202211302252.0A CN202211302252A CN115906411A CN 115906411 A CN115906411 A CN 115906411A CN 202211302252 A CN202211302252 A CN 202211302252A CN 115906411 A CN115906411 A CN 115906411A
Authority
CN
China
Prior art keywords
thermal
node
model
power
electric
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
CN202211302252.0A
Other languages
Chinese (zh)
Other versions
CN115906411B (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.)
Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
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 Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical Suzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority to CN202211302252.0A priority Critical patent/CN115906411B/en
Publication of CN115906411A publication Critical patent/CN115906411A/en
Application granted granted Critical
Publication of CN115906411B publication Critical patent/CN115906411B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种考虑全动态的电热综合能源系统最优能流建模方法，包括：步骤1，采集电热综合能源系统数据，包括管道长度、热阻、用户温度舒适度区间、电负荷时序分布；步骤2，结合电热综合能源系统数据建立考虑全动态的电热综合能源系统模型；步骤3，建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程，并简化交流潮流方程、离散的热电联产机组和热力系统方程；步骤4，基于步骤2所述简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件，构造最优能流模型。

Figure 202211302252

The invention discloses a method for modeling the optimal energy flow of an electrothermal integrated energy system considering full dynamics, including: step 1, collecting data of the electrothermal integrated energy system, including pipeline length, thermal resistance, user temperature comfort range, and electric load sequence distribution; step 2, combined with the data of the electric-thermal integrated energy system to establish a fully dynamic electric-thermal integrated energy system model; step 3, to establish the simplified equation of the nonlinear AC power flow, the discrete equation of the dynamic cogeneration unit and the thermal system, and simplify the AC power flow Equations, discrete cogeneration units and thermal system equations; step 4, based on the simplified AC power flow equations described in step 2, discrete cogeneration units and thermal system equations, and aiming at minimizing system operating costs The operating security constraints of the system are used to construct an optimal energy flow model.

Figure 202211302252

Description

考虑全动态的电热综合能源系统最优能流建模方法及系统Optimal energy flow modeling method and system for electric-thermal integrated energy system considering full dynamics

技术领域Technical Field

本发明属于能源系统建模与运行分析领域，具体来说，涉及一种考虑全动态的电热综合能源系统最优能流建模方法。The present invention belongs to the field of energy system modeling and operation analysis, and specifically, relates to an optimal energy flow modeling method for an electric-thermal integrated energy system taking into account a full dynamics approach.

背景技术Background Art

随着气候变化与环境污染的负面影响日益加剧，世界各国正在探索清洁、高效、可持续的能源利用新形式。作为综合能源系统的一种典型形式，电热综合能源系统通过联产机组等设备将供电和供热两个子系统相互耦合。与传统的分立式能源系统不同，热电联产系统能够充分利用发电回收的余热来供应部分工业或者民用热负荷，从而提升系统的综合能效；同时，供热负荷的热惯性还能为供电系统消纳更多可再生能源提供灵活性资源。As the negative impacts of climate change and environmental pollution become increasingly severe, countries around the world are exploring new forms of clean, efficient and sustainable energy utilization. As a typical form of integrated energy system, the electric-thermal integrated energy system couples the power supply and heating subsystems through equipment such as cogeneration units. Unlike traditional discrete energy systems, the cogeneration system can make full use of the waste heat recovered from power generation to supply part of the industrial or civil heat load, thereby improving the overall energy efficiency of the system; at the same time, the thermal inertia of the heating load can also provide flexibility resources for the power supply system to absorb more renewable energy.

电热综合系统的最优能流是在给定边界条件下，通过优化系统中某些状态量并使得系统某一性能指标(如经济性、环保性和能效等)达到最优值时的状态分布。最优能流是一个由非线性偏微分-常微分-代数描述的模型，但现有研究一般忽略系统中的非线性与动态进行建模，而简化了电热综合能源系各状态量的复杂特性，所得结果与真实情况存在较大差异。因此，建了一个合适的、准确的最优能流模型，是电热综合能源系统运行优化的关键。The optimal energy flow of the electric-thermal integrated system is the state distribution when a certain performance index of the system (such as economy, environmental protection and energy efficiency, etc.) reaches the optimal value by optimizing certain state quantities in the system under given boundary conditions. The optimal energy flow is a model described by nonlinear partial differential-ordinary differential-algebra, but existing studies generally ignore the nonlinearity and dynamics in the system for modeling, and simplify the complex characteristics of each state quantity of the electric-thermal integrated energy system. The results obtained are quite different from the actual situation. Therefore, building a suitable and accurate optimal energy flow model is the key to optimizing the operation of the electric-thermal integrated energy system.

现有技术文件1(CN113111555A)公开了一种基于叠加解耦法的质调节热力系统能流快速计算方法，包括：10)以环境温度为参考温度，构建热力系统动态模型，并根据质调节特征进行简化；20)基于简化的支路热传导方程，构建热力系统温度动态映射方程与权重矩阵，根据权重矩阵中的数值，确定热力系统中的温度映射方向；30)根据温度映射方向，将原始热力系统解耦为若干个由单一热源供热的辐射状热力系统；40)分别计算各解耦系统内的能流分布，原始热力系统的能流分布即为多个解耦系统的线性叠加，现有技术文件1不足之处在于该技术主要对热力系统的能流分布进行评估，但忽略了电力系统和热力系统间的耦合关系，且缺少热力系统动态特性对于电热综合能源系统运行优化的影响分析。Prior art document 1 (CN113111555A) discloses a method for quickly calculating the energy flow of a quality-regulated thermal system based on the superposition decoupling method, including: 10) taking the ambient temperature as the reference temperature, constructing a dynamic model of the thermal system, and simplifying it according to the quality regulation characteristics; 20) constructing a thermal system temperature dynamic mapping equation and a weight matrix based on the simplified branch heat conduction equation, and determining the temperature mapping direction in the thermal system according to the numerical value in the weight matrix; 30) decoupling the original thermal system into several radial thermal systems heated by a single heat source according to the temperature mapping direction; 40) calculating the energy flow distribution in each decoupled system respectively, and the energy flow distribution of the original thermal system is the linear superposition of multiple decoupled systems. The shortcoming of prior art document 1 is that this technology mainly evaluates the energy flow distribution of the thermal system, but ignores the coupling relationship between the power system and the thermal system, and lacks analysis of the impact of the dynamic characteristics of the thermal system on the operation optimization of the electric-thermal integrated energy system.

发明内容Summary of the invention

本发明所要解决的技术问题是：提出一种考虑全动态的电热综合能源系统最优能流建模方法，该方法详细的对电热综合能源系统内的非线性和动态特性进行建模，并针对不同的非线性项和动态项进行适当简化，使其便于运算，最后全面地考虑系统内各类型状态量的运行安全约束，建立了最优能流模型，实现运行优化的准确性。The technical problem to be solved by the present invention is: to propose a method for modeling the optimal energy flow of an electric-thermal integrated energy system taking into account the full dynamics. This method models the nonlinear and dynamic characteristics within the electric-thermal integrated energy system in detail, and appropriately simplifies different nonlinear terms and dynamic terms to facilitate calculation. Finally, the operational safety constraints of various types of state quantities in the system are comprehensively considered, and an optimal energy flow model is established to achieve the accuracy of operational optimization.

为解决上述技术问题，本技术方案采用一种考虑全动态的电热综合能源系统最优能流建模方法，该方法包括以下步骤：In order to solve the above technical problems, the present technical solution adopts an optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system, which includes the following steps:

步骤1，采集电热综合能源系统数据，包括管道长度、热阻、用户温度舒适度区间、电负荷时序分布；Step 1: Collect data of the electric and thermal integrated energy system, including pipeline length, thermal resistance, user temperature comfort range, and time series distribution of electric load;

步骤2，结合电热综合能源系统数据建立考虑全动态的电热综合能源系统模型；Step 2, combining the electric-thermal integrated energy system data to establish an electric-thermal integrated energy system model that takes into account the full dynamics;

步骤3，建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程，并简化交流潮流方程、离散的热电联产机组和热力系统方程；Step 3, establishing simplified equations of nonlinear AC power flow, discrete equations of dynamic combined heat and power units and thermal systems, and simplifying AC power flow equations, discrete combined heat and power units and thermal system equations;

步骤4，基于步骤2所述简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件，构造最优能流模型。Step 4, based on the simplified AC power flow equation, discrete cogeneration unit and thermal system equation described in step 2, establish the operation safety constraints of the electric and thermal integrated energy system with the goal of minimizing the system operation cost, and construct the optimal energy flow model.

步骤2具体包括：步骤201，建立交流潮流模型；Step 2 specifically includes: Step 201, establishing an AC power flow model;

步骤202，建立包含压缩机模型、燃烧室模型、汽轮机模型、换热器模型的热电联产机组动态模型；Step 202, establishing a dynamic model of a cogeneration unit including a compressor model, a combustion chamber model, a steam turbine model, and a heat exchanger model;

步骤203，建立热力系统动态模型，包括水力部分模型和热力部分模型；其中水力部分模型为：Step 203, establish a dynamic model of the thermal system, including a hydraulic part model and a thermal part model; wherein the hydraulic part model is:

Am＝dAm＝d

BΔp＝0BΔp＝0

Δp＝Km2 Δp＝Km 2

式中，A和B分别为热力系统的节点-支路关联矩阵和回路-支路关联矩阵，m热水管道的质量流量向量，d为节点注入的质量流量向量，Δp为管道压降向量，K为管道摩阻系数；Where A and B are the node-branch association matrix and loop-branch association matrix of the thermal system, respectively; m is the mass flow vector of the hot water pipeline; d is the mass flow vector of the node injection; Δp is the pipeline pressure drop vector; and K is the pipeline friction coefficient;

热力部分模型为：The thermal model is:

Figure SMS_1
Figure SMS_1

Figure SMS_2
Figure SMS_2

Figure SMS_3
Figure SMS_3

式中，x和t分别为空间和时间变量，v为水流流速，Tp表示管道温度，Ti n表示节点i的温度，mb为管道b的质量流量，Tk p,o为管道k的出口温度，Tk p,i为管道k的入口温度，λ为管道热阻系数，Ta为环境温度，

Figure SMS_4
和
Figure SMS_5
分别表示以节点i为首节点和末节点管道集。Where x and t are spatial and temporal variables, respectively; v is the water velocity; Tp is the pipe temperature; Tin is the temperature of node i ; mb is the mass flow rate of pipe b; Tkp ,o is the outlet temperature of pipe k; Tkp ,i is the inlet temperature of pipe k; λ is the thermal resistance coefficient of the pipe; Ta is the ambient temperature;
Figure SMS_4
and
Figure SMS_5
They represent the pipeline sets with node i as the first node and the last node respectively.

压缩机模型为：The compressor model is:

p2,t＝CPR1×p1,t p2,t ＝ CPR1 ×p1 ,t

Figure SMS_6
Figure SMS_6

Pc,t＝Cama,t(T2,t-T1,t)P c,t = C a m a,t (T 2,t -T 1,t )

式中，t为时间标志，p1,t和p2,t分别表示t时刻压缩机的入口和出口压力，T1,t和T2,t分别表示t时刻压缩机的入口和出口温度，CPR1表示压缩机的压比，β1表示空气绝热系数，η1表示压缩机效率，Pc,t表示t时刻压缩机耗功，Ca表示空气比热容，ma,t表示流入空气的质量流量。Wherein, t is the time mark, p1,t and p2,t represent the inlet and outlet pressures of the compressor at time t respectively, T1,t and T2 ,t represent the inlet and outlet temperatures of the compressor at time t respectively, CPR1 represents the pressure ratio of the compressor, β1 represents the air adiabatic coefficient, η1 represents the compressor efficiency, Pc,t represents the power consumption of the compressor at time t, Ca represents the specific heat capacity of air, and ma ,t represents the mass flow rate of the inflowing air.

燃烧室模型为：The combustion chamber model is:

Figure SMS_7
Figure SMS_7

式中，T3,t分别表示t时刻的燃烧室温度，β2表示燃烧室的储热系数，Hg表示燃烧室热值，LHV表示燃料的低位热值，mf,t表示流入燃料的质量流量，Cs为混合烟气的比热容。Wherein, T3 ,t represents the combustion chamber temperature at time t, β2 represents the heat storage coefficient of the combustion chamber, Hg represents the combustion chamber calorific value, LHV represents the lower heating value of the fuel, mf ,t represents the mass flow rate of the inflowing fuel, and Cs is the specific heat capacity of the mixed flue gas.

汽轮机模型为：The turbine model is:

p3,t＝CPR2p2,t p 3,t = CPR 2 p 2,t

Figure SMS_8
Figure SMS_8

Pb,t＝Cs(ma,t+mf,t)(T4,t-T3,t)P b,t =C s (m a,t +m f,t )(T 4,t -T 3,t )

式中，T4,t分别表示t时刻的汽轮机出口温度，β3表示燃烧室的绝热效率，Pb,t为t时刻汽轮机生产的总功率，CPR2表示汽轮机的压比。Where T4 ,t represents the turbine outlet temperature at time t, β3 represents the adiabatic efficiency of the combustion chamber, Pb,t is the total power produced by the turbine at time t, and CPR2 represents the pressure ratio of the turbine.

联产机组生产的用于发电和供热的功率分别为：The power produced by the combined unit for electricity and heat supply is:

Pg,t＝η2η3(Pb,t-Pc,t)P g,t = η 2 η 3 (P b,t -P c,t )

φg,t＝η2(1-η3)(Pb,t-Pc,t)φ g,t =η 2 (1-η 3 )(P b,t -P c,t )

式中，Pg,t和φg,t分别为联产机组用于发电和供热的功率，η2表示汽轮机机械效率，η3表示热电比。Where P g,t and φ g,t are the power of the cogeneration unit for power generation and heat supply, respectively, η 2 represents the mechanical efficiency of the steam turbine, and η 3 represents the heat-to-electricity ratio.

换热器模型为：The heat exchanger model is:

Figure SMS_9
Figure SMS_9

式中，β4表示换热器的储热系数，Cw表示水的比热容，mw,t为换热器中水流的质量流量，T5,t和T6,t分别表示换热器的出口和入口温度。Where β4 represents the heat storage coefficient of the heat exchanger, Cw represents the specific heat capacity of water, mw ,t is the mass flow rate of water in the heat exchanger, and T5,t and T6,t represent the outlet and inlet temperatures of the heat exchanger, respectively.

优选地，步骤3包括：Preferably, step 3 comprises:

步骤301，在输网级系统中，令节点电压的相角差为0，令节点电压幅值为1对非线性交流潮流模型进行简化得到电力系统潮流模型：Step 301, in the transmission grid system, the phase angle difference of the node voltage is set to 0, and the node voltage amplitude is set to 1 to simplify the nonlinear AC power flow model to obtain the power system power flow model:

Figure SMS_10
Figure SMS_10

Figure SMS_11
Figure SMS_11

Figure SMS_12
Figure SMS_12

Figure SMS_13
Figure SMS_13

式中，i和j分别表示节点编号，

Figure SMS_14
表示节点i的电压幅值，PG,i和PL,i表示节点i的发电机有功功率和负荷有功功率，QG,i和QL,i表示节点i的发电机无功功率和负荷无功功率，Gij和Bij表示节点i和节点j之间的电导与电纳，θij表示节点i和节点j之间的相角差，Pij和Qij表示节点i和节点j之间的传输的有功功率和无功功率；In the formula, i and j represent the node numbers,
Figure SMS_14
represents the voltage amplitude of node i, PG,i and PL,i represent the generator active power and load active power of node i, QG,i and QL,i represent the generator reactive power and load reactive power of node i, Gij and Bij represent the conductance and susceptance between node i and node j, θij represents the phase angle difference between node i and node j, Pij and Qij represent the active power and reactive power transmitted between node i and node j;

步骤202，采用离散的时间和空间步长解法将热力系统模型中的热力部分模型离散化，得到离散化热力部分模型：Step 202, using a discrete time and space step solution method to discretize the thermal part model in the thermal system model to obtain a discretized thermal part model:

Figure SMS_15
Figure SMS_15

Figure SMS_16
Figure SMS_16

式中，J1，J2，J3和J4分别为传递系数；

Figure SMS_17
表示j时刻i+1处的管道温度，
Figure SMS_18
表示j+1时刻i+1处的管道温度，
Figure SMS_19
表示j时刻i处的管道温度，
Figure SMS_20
表示j+1时刻i处的管道温度；Δx和Δt分别为离散的空间和时间步长，L和Γ分别为管道长度和时间区间，Nx和Nt分别为空间和时间步数，
Figure SMS_21
Where, J 1 , J 2 , J 3 and J 4 are the transfer coefficients respectively;
Figure SMS_17
represents the pipe temperature at time i+1,
Figure SMS_18
represents the pipe temperature at i+1 at time j+1,
Figure SMS_19
represents the pipe temperature at time i,
Figure SMS_20
represents the pipe temperature at time i j+1; Δx and Δt are the discrete space and time steps, L and Γ are the pipe length and time interval, Nx and Nt are the space and time steps, respectively.
Figure SMS_21

步骤203，采用向后欧拉格式将燃烧室模型和换热器模型离散化；Step 203, discretizing the combustion chamber model and the heat exchanger model using a backward Euler format;

其中则离散燃烧室模型为：The discrete combustion chamber model is:

Figure SMS_22
Figure SMS_22

离散换热器模型为：The discrete heat exchanger model is:

Figure SMS_23
Figure SMS_23

步骤4包括：Step 4 includes:

步骤401，建立电力系统运行安全约束，包括：Step 401, establishing power system operation safety constraints, including:

节点电压幅值约束：

Figure SMS_24
相角约束：
Figure SMS_25
Node voltage amplitude constraint:
Figure SMS_24
Phase Angle Constraint:
Figure SMS_25

支路传输功率约束：

Figure SMS_26
Branch transmission power constraints:
Figure SMS_26

发电机有功功率约束：

Figure SMS_27
Generator active power constraint:
Figure SMS_27

发电机无功功率约束：

Figure SMS_28
Generator reactive power constraints:
Figure SMS_28

式中，

Figure SMS_29
表示节点i和节点j之间传输的视在功率上限，
Figure SMS_30
和
Figure SMS_31
为发电机i生产有功功率的下限和上限，
Figure SMS_32
和
Figure SMS_33
为发电机i生产无功功率的下限和上限；In the formula,
Figure SMS_29
represents the apparent upper limit of the power transmitted between node i and node j,
Figure SMS_30
and
Figure SMS_31
are the lower and upper limits of active power produced by generator i,
Figure SMS_32
and
Figure SMS_33
The lower and upper limits of reactive power produced by generator i;

步骤402，建立热力系统运行安全约束，包括：Step 402, establishing thermal system operation safety constraints, including:

节点压力约束：

Figure SMS_34
Nodal pressure constraints:
Figure SMS_34

管道和节点的质量流量的容量约束：

Figure SMS_35
Capacity constraints for mass flow rates at pipes and nodes:
Figure SMS_35

管道质量流量变化率约束：γminmi,t-1≤mi,t≤γmaxmi,t-1,γmindi,t-1≤di,t≤γmaxdi,t-1；Pipeline mass flow rate change rate constraints: γ min m i,t-1 ≤m i,t ≤γ max m i,t-1 ,γ min d i,t-1 ≤d i,t ≤γ max d i,t-1 ;

节点供水温度约束：

Figure SMS_36
Node water supply temperature constraints:
Figure SMS_36

节点回水温度约束：

Figure SMS_37
Node return water temperature constraint:
Figure SMS_37

式中，

Figure SMS_40
和
Figure SMS_42
分别为节点i的水压的下限和上限，
Figure SMS_45
和
Figure SMS_39
分别为管道i的质量流量下限和上限，
Figure SMS_41
和
Figure SMS_43
分别为节点i的质量流量的下限和上限，γmin和γmax分别为质量流量变化率的下限和上限，
Figure SMS_44
和
Figure SMS_38
分别为t时刻节点i供水温度的下限和上限，
Figure SMS_46
和
Figure SMS_47
分别为t时刻节点i回水温度的下限和上限；In the formula,
Figure SMS_40
and
Figure SMS_42
are the lower and upper limits of the water pressure at node i,
Figure SMS_45
and
Figure SMS_39
are the lower and upper limits of the mass flow rate of pipeline i, respectively.
Figure SMS_41
and
Figure SMS_43
are the lower and upper limits of the mass flow rate of node i, respectively; γ min and γ max are the lower and upper limits of the mass flow rate change rate, respectively;
Figure SMS_44
and
Figure SMS_38
are the lower and upper limits of the water supply temperature at node i at time t,
Figure SMS_46
and
Figure SMS_47
are the lower and upper limits of the return water temperature at node i at time t respectively;

步骤403，建立热电联产机组运行安全约束，包括：Step 403, establishing safety constraints for the operation of the cogeneration unit, including:

压缩机入口温度和压力约束：T1 min≤T1,t≤T1 max，

Figure SMS_48
Compressor inlet temperature and pressure constraints: T 1 min ≤ T 1, t ≤ T 1 max ,
Figure SMS_48

燃烧室输入燃料的质量流量约束：

Figure SMS_49
The mass flow rate constraint of the fuel input to the combustion chamber is:
Figure SMS_49

汽轮机出口压力约束：

Figure SMS_50
Steam turbine outlet pressure constraint:
Figure SMS_50

燃料和空气混合比约束：αminmf,t≤ma,t≤αmaxmf,t；Fuel and air mixture ratio constraints: α min m f,t ≤m a,t ≤α max m f,t ;

燃烧室温度约束：T3 min≤T3,t≤T3 max Combustion chamber temperature constraint: T 3 min ≤T 3,t ≤T 3 max

汽轮机出口温度约束：

Figure SMS_51
Turbine outlet temperature constraint:
Figure SMS_51

热电比约束：

Figure SMS_52
Thermoelectric ratio constraints:
Figure SMS_52

式中，T1 min和T1 max分别为压缩机入口温度的下限和上限，

Figure SMS_54
和
Figure SMS_59
分别为压缩机入口压力的下限和上限，
Figure SMS_60
和
Figure SMS_55
分别为燃料质量流量的下限和上限，
Figure SMS_56
和
Figure SMS_57
分别为汽轮机出口压力的下限和上限，αmin和αmax分布为空气燃料比的下限和上限，T3 min和T3 max分别为燃烧室温度的下限和上限，和T4 max分别为汽轮机出口温度的下限和上限，
Figure SMS_58
和
Figure SMS_53
分别为热电比的下限和上限；Where T 1 min and T 1 max are the lower and upper limits of the compressor inlet temperature, respectively.
Figure SMS_54
and
Figure SMS_59
are the lower and upper limits of the compressor inlet pressure, respectively.
Figure SMS_60
and
Figure SMS_55
are the lower and upper limits of the fuel mass flow rate, respectively.
Figure SMS_56
and
Figure SMS_57
are the lower and upper limits of the turbine outlet pressure, α min and α max are the lower and upper limits of the air-fuel ratio, T 3 min and T 3 max are the lower and upper limits of the combustion chamber temperature, and T 4 max are the lower and upper limits of the turbine outlet temperature,
Figure SMS_58
and
Figure SMS_53
are the lower and upper limits of the thermoelectric ratio, respectively;

步骤404，以最小化系统运行成本为目标，建立最优能流模型的目标函数：Step 404, with the goal of minimizing the system operation cost, establish the objective function of the optimal energy flow model:

Figure SMS_61
Figure SMS_61

式中，F表示电热综合能源系统的总运行成本，minF表示电热综合能源系统的最小总运行成本，c1fi,t表示电热综合能源系统内第i个热电联产机组在t时间段的总运行成本；c2gj,t表示电热综合能源系统内第j个普通发电机在t时间段的总运行成本；c1和c2分别为热电联产机组的燃料单价和普通发电机的燃煤单价，

Figure SMS_62
为普通发电机集合，
Figure SMS_63
为热电联产机组集合，fi,t和gj,t分别为t时间段联产机组i和普通发电机j的成本函数；Nt为设定计算周期内的优化采样时刻集分别表示为：Where F represents the total operating cost of the electric-thermal integrated energy system, minF represents the minimum total operating cost of the electric-thermal integrated energy system, c 1 fi,t represents the total operating cost of the ith cogeneration unit in the electric-thermal integrated energy system in time period t; c 2 g j,t represents the total operating cost of the jth ordinary generator in the electric-thermal integrated energy system in time period t; c 1 and c 2 are the fuel unit price of the cogeneration unit and the coal unit price of the ordinary generator, respectively.
Figure SMS_62
For a common generator set,
Figure SMS_63
is the set of cogeneration units, fi ,t and gj ,t are the cost functions of cogeneration unit i and ordinary generator j in time period t respectively; Nt is the set of optimized sampling time in the set calculation period and is expressed as:

Figure SMS_64
Figure SMS_64

Figure SMS_65
Figure SMS_65

式中，μ11，μ12和μ13分别为联产机组发电功率和燃料流量之间的折算系数，μ21，μ22和μ23分别为普通发电机发电功率和燃煤量之间的折算系数；PG,j表示节点j的发电机有功功率；PG,i表示节点i的发电机有功功率。In the formula, μ11 , μ12 and μ13 are the conversion coefficients between the power generation of the cogeneration unit and the fuel flow, μ21 , μ22 and μ23 are the conversion coefficients between the power generation of the ordinary generator and the coal consumption, PG,j represents the active power of the generator at node j, and PG,i represents the active power of the generator at node i.

考虑全动态的电热综合能源系统最优能流建模系统，基于考虑全动态的电热综合能源系统最优能流建模方法，系统包括：The optimal energy flow modeling system of the electric-thermal integrated energy system considering the full dynamics is based on the optimal energy flow modeling method of the electric-thermal integrated energy system considering the full dynamics. The system includes:

数据采集模块、能源管理分析模块、逻辑计算模块与供热模型模块；Data acquisition module, energy management analysis module, logic calculation module and heating model module;

数据采集模块用于采集电热综合能源系统数据；The data acquisition module is used to collect data of the electric and thermal integrated energy system;

能源管理分析模块用于建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程；The energy management analysis module is used to establish simplified equations for nonlinear AC power flow, discrete equations for dynamic combined heat and power units and thermal systems;

逻辑计算模块用于基于简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件；The logic calculation module is used to establish the operation safety constraints of the electric and thermal integrated energy system based on the simplified AC power flow equation, discrete cogeneration unit and thermal system equation, with the goal of minimizing the system operation cost;

最优能流建模模块用于根据运行安全约束条件和电热综合能源系统模型建立考虑全动态的最优能流模型。The optimal energy flow modeling module is used to establish an optimal energy flow model that takes into account the full dynamics based on the operational safety constraints and the electric and thermal integrated energy system model.

本发明的有益效果在于：该方法全面的对电热综合能源系统的全动态和非线性特征减小建模，有利于精确刻画系统状态，从而准确制定电热综合能源系统运行控制策略，提高系统经济性和安全性。The beneficial effect of the present invention is that the method comprehensively reduces the modeling of the full dynamic and nonlinear characteristics of the electric-thermal integrated energy system, which is conducive to accurately characterizing the system state, thereby accurately formulating the operation control strategy of the electric-thermal integrated energy system and improving the system economy and safety.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明实施例考虑全动态的电热综合能源系统最优能流建模方法流程示意图；FIG1 is a schematic flow chart of an optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to an embodiment of the present invention;

图2为本发明实施例中热电联产机组结构图；FIG2 is a structural diagram of a cogeneration unit in an embodiment of the present invention;

图3为本发明实施例中采用的热力系统结构图；FIG3 is a structural diagram of a thermal system used in an embodiment of the present invention;

图4为本发明实施例中热力系统中节点1供热功率的时序分布。FIG. 4 is a time series distribution of heating power of node 1 in the thermal system in an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明的技术方案进行清楚、完整地描述。本申请所描述的实施例仅仅是本发明一部分的实施例，而不是全部实施例。基于本发明精神，本领域普通技术人员在没有做出创造性劳动前提下所获得的有所其它实施例，都属于本发明的保护范围。In order to make the purpose, technical scheme and advantages of the present invention clearer, the technical scheme of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. The embodiments described in this application are only embodiments of a part of the present invention, not all embodiments. Based on the spirit of the present invention, other embodiments obtained by ordinary technicians in this field without making creative work are all within the scope of protection of the present invention.

实施例1：一种考虑全动态的电热综合能源系统最优能流建模方法，如图1所示。该方法包括以下步骤：Embodiment 1: A method for modeling optimal energy flow of a fully dynamic electric-thermal integrated energy system, as shown in FIG1 . The method comprises the following steps:

步骤1，建立考虑全动态的电热综合能源系统模型，包括交流潮流模型、动态热电联产机组模型和动态热力系统模型；Step 1: Establish a fully dynamic electric-thermal integrated energy system model, including an AC power flow model, a dynamic cogeneration unit model, and a dynamic thermal system model;

步骤101，建立交流潮流模型，包含节点功率守恒和支路功率守恒方程：Step 101, establish an AC power flow model, including node power conservation and branch power conservation equations:

Figure SMS_66
Figure SMS_66

Figure SMS_67
Figure SMS_67

Figure SMS_68
Figure SMS_68

Figure SMS_69
Figure SMS_69

式中，i和j分别表示节点编号，

Figure SMS_70
表示节点i的电压幅值，PG,i和PL,i表示节点i的发电机有功功率和负荷有功功率，QG,i和QL,i表示节点i的发电机无功功率和负荷无功功率，Gij和Bij表示节点i和节点j之间的电导与电纳，θij表示节点i和节点j之间的相角差，Pij和Qij表示节点i和节点j之间的传输的有功功率和无功功率。In the formula, i and j represent the node numbers,
Figure SMS_70
represents the voltage amplitude of node i, PG,i and PL,i represent the generator active power and load active power of node i, QG,i and QL,i represent the generator reactive power and load reactive power of node i, Gij and Bij represent the conductance and susceptance between node i and node j, θij represents the phase angle difference between node i and node j, and Pij and Qij represent the active power and reactive power transmitted between node i and node j.

步骤102，建立热电联产机组动态模型，包含压缩机模型、燃烧室模型、汽轮机模型、换热器模型，热电联产机组的拓扑及各状态量分布如图2所示。压缩机模型描述了出入口温度、压力关系及消耗电功率计算式，分别表示为：Step 102, establish a dynamic model of the cogeneration unit, including a compressor model, a combustion chamber model, a turbine model, and a heat exchanger model. The topology and state quantity distribution of the cogeneration unit are shown in Figure 2. The compressor model describes the inlet and outlet temperature, pressure relationship and the power consumption calculation formula, which are respectively expressed as:

p2,t＝CPR1×p1,t (5)p 2,t = CPR 1 × p 1,t (5)

Figure SMS_71
Figure SMS_71

Pc,t＝Cama,t(T2,t-T1,t) (7)P c,t = C a m a,t (T 2,t -T 1,t ) (7)

式中，t为时间标志，p1,t和p2,t分别表示t时刻压缩机的入口和出口压力，T1,t和T2,t分别表示t时刻压缩机的入口和出口温度，CPR1表示压缩机的压比，β1表示空气绝热系数，η1表示压缩机效率，Pc,t表示t时刻压缩机耗功，Ca表示空气比热容，ma,t表示流入空气的质量流量。Wherein, t is the time mark, p1,t and p2,t represent the inlet and outlet pressures of the compressor at time t respectively, T1,t and T2 ,t represent the inlet and outlet temperatures of the compressor at time t respectively, CPR1 represents the pressure ratio of the compressor, β1 represents the air adiabatic coefficient, η1 represents the compressor efficiency, Pc,t represents the power consumption of the compressor at time t, Ca represents the specific heat capacity of air, and ma ,t represents the mass flow rate of the inflowing air.

燃烧室模型描述了燃烧室内的能量守恒关系，可表示为：The combustion chamber model describes the energy conservation relationship in the combustion chamber and can be expressed as:

Figure SMS_72
Figure SMS_72

式中，T3,t分别表示t时刻的燃烧室温度，β2表示燃烧室的储热系数，Hg表示燃烧室热值，LHV表示燃料的低位热值，mf,t表示流入燃料的质量流量，Cs为混合烟气的比热容。Wherein, T3 ,t represents the combustion chamber temperature at time t, β2 represents the heat storage coefficient of the combustion chamber, Hg represents the combustion chamber calorific value, LHV represents the lower heating value of the fuel, mf ,t represents the mass flow rate of the inflowing fuel, and Cs is the specific heat capacity of the mixed flue gas.

汽轮机模型包括出入口温度、压力以及产生功率的计算式，可分别表示为：The steam turbine model includes the calculation formulas for inlet and outlet temperature, pressure, and generated power, which can be expressed as:

p3,t＝CPR2p2,t (9)p 3,t =CPR 2 p 2,t (9)

Figure SMS_73
Figure SMS_73

Pb,t＝Cs(ma,t+mf,t)(T4,t-T3,t) (11)P b,t =C s (m a,t +m f,t )(T 4,t -T 3,t ) (11)

式中，T4,t分别表示t时刻的汽轮机出口温度，β3表示燃烧室的绝热效率，Pb,t为t时刻汽轮机生产的总功率，CPR2表示汽轮机的压比。Where T4 ,t represents the turbine outlet temperature at time t, β3 represents the adiabatic efficiency of the combustion chamber, Pb,t is the total power produced by the turbine at time t, and CPR2 represents the pressure ratio of the turbine.

联产机组生产的用于发电和供热的功率可分别表示为：The power produced by the cogeneration unit for electricity generation and heat generation can be expressed as:

Pg,t＝η2η3(Pb,t-Pc,t) (12)P g,t = η 2 η 3 (P b,t -P c,t ) (12)

φg,t＝η2(1-η3)(Pb,t-Pc,t) (13)φ g,t =η 2 (1-η 3 )(P b,t -P c,t ) (13)

式中，Pg,t和φg,t分别为联产机组用于发电和供热的功率，η2表示汽轮机机械效率，η3表示热电比。Where P g,t and φ g,t are the power of the cogeneration unit for power generation and heat supply, respectively, η 2 represents the mechanical efficiency of the steam turbine, and η 3 represents the heat-to-electricity ratio.

换热器模型描述了换热器内部的能量守恒，可表示为：The heat exchanger model describes the energy conservation inside the heat exchanger and can be expressed as:

Figure SMS_74
Figure SMS_74

式中，β4表示换热器的储热系数，Cw表示水的比热容，mw,t为换热器中水流的质量流量，T5,t和T6,t分别表示换热器的出口和入口温度。Where β4 represents the heat storage coefficient of the heat exchanger, Cw represents the specific heat capacity of water, mw ,t is the mass flow rate of water in the heat exchanger, and T5,t and T6,t represent the outlet and inlet temperatures of the heat exchanger, respectively.

步骤103，建立热力系统动态模型，包含水力部分和热力部分。水力部分包括：Step 103, establish a dynamic model of the thermal system, including a hydraulic part and a thermal part. The hydraulic part includes:

Am＝d (15)Am＝d (15)

BΔp＝0 (16)BΔp＝0 (16)

Δp＝Km2 (17)Δp＝Km 2 (17)

式中，A和B分别为热力系统的节点-支路关联矩阵和回路-支路关联矩阵，m热水管道的质量流量向量，d为节点注入的质量流量向量，Δp为管道压降向量，K为管道摩阻系数。其中，式(15)描述了节点处的质量守恒，式(16)描述了回路的压降平衡，式(17)描述了管道压降和管道质量流量的关系。Where A and B are the node-branch association matrix and loop-branch association matrix of the thermal system, respectively, m is the mass flow vector of the hot water pipe, d is the mass flow vector of the node injection, Δp is the pipeline pressure drop vector, and K is the pipeline friction coefficient. Among them, equation (15) describes the mass conservation at the node, equation (16) describes the pressure drop balance of the loop, and equation (17) describes the relationship between the pipeline pressure drop and the pipeline mass flow.

热力部分包括：The thermal part includes:

Figure SMS_75
Figure SMS_75

Figure SMS_76
Figure SMS_76

Figure SMS_77
Figure SMS_77

式中，x和t分别为空间和时间变量，v为水流流速，Tp表示管道温度，Ti n表示节点i的温度，mb为管道b的质量流量，Tk p,o为管道k的出口温度，Tk p,i为管道k的入口温度，λ为管道热阻系数，Ta为环境温度，

Figure SMS_78
和
Figure SMS_79
分别表示以节点i为首节点和末节点管道集。Where x and t are spatial and temporal variables, respectively; v is the water velocity; Tp is the pipe temperature; Tin is the temperature of node i ; mb is the mass flow rate of pipe b; Tkp ,o is the outlet temperature of pipe k; Tkp ,i is the inlet temperature of pipe k; λ is the thermal resistance coefficient of the pipe; Ta is the ambient temperature;
Figure SMS_78
and
Figure SMS_79
They represent the pipeline sets with node i as the first node and the last node respectively.

步骤2，基于考虑全动态的电热综合能源系统模型，建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程；Step 2: Based on the fully dynamic electric-thermal integrated energy system model, the simplified equations of nonlinear AC power flow and the discrete equations of dynamic cogeneration units and thermal systems are established;

步骤2包括：Step 2 includes:

步骤201，对非线性交流潮流模型进行简化。在输网级系统中，节点电压的相角差趋近于0，可得到：Step 201, simplify the nonlinear AC power flow model. In the transmission grid system, the phase angle difference of the node voltage approaches 0, and it can be obtained:

Figure SMS_80
Figure SMS_80

此外，输网级系统中的节点电压幅值趋近于1，可得到：In addition, the node voltage amplitude in the transmission grid system approaches 1, and we can get:

Figure SMS_81
Figure SMS_81

将式(21)和式(22)代入式(1)至式(4)，电力系统潮流模型可重写为：Substituting equations (21) and (22) into equations (1) to (4), the power system flow model can be rewritten as:

Figure SMS_82
Figure SMS_82

Figure SMS_83
Figure SMS_83

Figure SMS_84
Figure SMS_84

Figure SMS_85
Figure SMS_85

步骤202，采用中心隐式格式将热力系统模型中的偏微分方程离散化。首先将研究的时间区间离散化，可得到离散的时间和空间步长，分别为：Step 202, discretize the partial differential equations in the thermal system model using the central implicit format. First, discretize the time interval of study to obtain discrete time and space steps, which are:

Figure SMS_86
Figure SMS_86

式中，Δx和Δt分别为离散的空间和时间步长，L和Γ分别为管道长度和时间区间，Nx和Nt分别为空间和时间步数。式(18)中的各偏微分项可表示为：In the formula, Δx and Δt are the discrete space and time steps, L and Γ are the pipeline length and time interval, Nx and Nt are the space and time steps, respectively. The partial differential terms in formula (18) can be expressed as:

Figure SMS_87
Figure SMS_87

Figure SMS_88
Figure SMS_88

Figure SMS_89
Figure SMS_89

式中，J1，J2，J3和J4分别为传递系数；

Figure SMS_90
表示j时刻i+1处的管道温度，
Figure SMS_91
表示j+1时刻i+1处的管道温度，
Figure SMS_92
表示j时刻i处的管道温度，
Figure SMS_93
表示j+1时刻i处的管道温度；。将式(28)至式(30)代入式(18)，可得到：Where, J 1 , J 2 , J 3 and J 4 are the transfer coefficients respectively;
Figure SMS_90
represents the pipe temperature at time i+1,
Figure SMS_91
represents the pipe temperature at i+1 at time j+1,
Figure SMS_92
represents the pipe temperature at time i,
Figure SMS_93
represents the pipe temperature at time i at j+1; Substituting equations (28) to (30) into equation (18), we can obtain:

Figure SMS_94
Figure SMS_94

Figure SMS_95
Figure SMS_95

式中，J1，J2，J3和J4分别为传递系数。Where J 1 , J 2 , J 3 and J 4 are the transfer coefficients respectively.

步骤203，采用向后欧拉格式将热电联产机组模型中的常微分方程离散化。式(8)和式(14)中的常微分可离散为：Step 203, discretize the ordinary differential equation in the cogeneration unit model using the backward Euler format. The ordinary differential in equation (8) and equation (14) can be discretized as:

Figure SMS_96
Figure SMS_96

Figure SMS_97
Figure SMS_97

将式(33)和式(34)分别代入式(8)和式(14)，可得到：Substituting equation (33) and equation (34) into equation (8) and equation (14) respectively, we can obtain:

Figure SMS_98
Figure SMS_98

Figure SMS_99
Figure SMS_99

步骤3，结合简化的交流潮流方程、离散的热电联产机组和热力系统方程，建立电热综合能源系统的运行安全约束，构造最优能流模型。Step 3, combining the simplified AC power flow equation, discrete cogeneration unit and thermal system equation, establish the operation safety constraints of the electric and thermal integrated energy system and construct the optimal energy flow model.

步骤3包括：Step 3 includes:

步骤301建立电力系统运行安全约束，包括节点电压幅值约束、相角约束、支路传输功率约束，发电机有功功率和无功功率约束，分别如式(37)至(41)所示：Step 301 establishes power system operation safety constraints, including node voltage amplitude constraints, phase angle constraints, branch transmission power constraints, generator active power and reactive power constraints, as shown in equations (37) to (41), respectively:

Figure SMS_100
Figure SMS_100

Figure SMS_101
Figure SMS_101

Figure SMS_102
Figure SMS_102

Figure SMS_103
Figure SMS_103

Figure SMS_104
Figure SMS_104

式中，

Figure SMS_105
表示节点i和节点j之间传输的视在功率上限，
Figure SMS_106
和
Figure SMS_107
为发电机i生产有功功率的下限和上限，
Figure SMS_108
和
Figure SMS_109
为发电机i生产无功功率的下限和上限。In the formula,
Figure SMS_105
represents the apparent upper limit of the power transmitted between node i and node j,
Figure SMS_106
and
Figure SMS_107
are the lower and upper limits of active power produced by generator i,
Figure SMS_108
and
Figure SMS_109
The lower and upper limits for reactive power produced by generator i.

步骤302，建立热力系统运行安全约束，包括节点压力约束，管道和节点的质量流量的容量约束，管道质量流量变化率约束，节点供水温度约束和节点回水温度约束，其中节点供水温度约束和节点回水温度约束对应用户舒适温度区间，分别如式(42)至式(46)所示：Step 302, establish the safety constraints of thermal system operation, including node pressure constraints, capacity constraints of mass flow of pipelines and nodes, pipeline mass flow rate change constraints, node water supply temperature constraints and node return water temperature constraints, where the node water supply temperature constraints and node return water temperature constraints correspond to the user comfort temperature range, as shown in equations (42) to (46), respectively:

Figure SMS_110
Figure SMS_110

Figure SMS_111
Figure SMS_111

γminmi,t-1≤mi,t≤γmaxmi,t-1,γmindi,t-1≤di,t≤γmaxdi,t-1 (44)γ min m i,t-1 ≤m i,t ≤γ max m i,t-1 ,γ min d i,t-1 ≤d i,t ≤γ max d i,t-1 (44)

Figure SMS_112
Figure SMS_112

Figure SMS_113
Figure SMS_113

式中，

Figure SMS_116
和
Figure SMS_118
分别为节点i的水压的下限和上限，
Figure SMS_121
和
Figure SMS_115
分别为管道i的质量流量下限和上限，
Figure SMS_117
和
Figure SMS_122
分别为节点i的质量流量的下限和上限，γmin和γmax分别为质量流量变化率的下限和上限，
Figure SMS_123
和
Figure SMS_114
分别为t时刻节点i供水温度的下限和上限，
Figure SMS_119
和
Figure SMS_120
分别为t时刻节点i回水温度的下限和上限。In the formula,
Figure SMS_116
and
Figure SMS_118
are the lower and upper limits of the water pressure at node i,
Figure SMS_121
and
Figure SMS_115
are the lower and upper limits of the mass flow rate of pipeline i, respectively.
Figure SMS_117
and
Figure SMS_122
are the lower and upper limits of the mass flow rate of node i, respectively; γ min and γ max are the lower and upper limits of the mass flow rate change rate, respectively;
Figure SMS_123
and
Figure SMS_114
are the lower and upper limits of the water supply temperature at node i at time t,
Figure SMS_119
and
Figure SMS_120
are the lower and upper limits of the return water temperature at node i at time t respectively.

步骤303，建立热电联产机组运行安全约束，包括压缩机入口温度和压力约束，燃烧室输入燃料的质量流量约束，汽轮机出口压力约束，燃料和空气混合比约束，燃烧室温度约束，汽轮机出口温度约束，热电比约束，分别表示如式(47)Step 303, establish the safety constraints for the operation of the cogeneration unit, including the compressor inlet temperature and pressure constraints, the mass flow constraint of the fuel input to the combustion chamber, the turbine outlet pressure constraint, the fuel and air mixing ratio constraint, the combustion chamber temperature constraint, the turbine outlet temperature constraint, and the heat-to-power ratio constraint, which are respectively expressed as formula (47):

至式(54)所示。As shown in formula (54).

T1 min≤T1,t≤T1 max (47)T 1 min ≤T 1,t ≤T 1 max (47)

Figure SMS_124
Figure SMS_124

Figure SMS_125
Figure SMS_125

Figure SMS_126
Figure SMS_126

αminmf,t≤ma,t≤αmaxmf,t (51)α min m f,t ≤m a,t ≤α max m f,t (51)

T3 min≤T3,t≤T3 max (52)T 3 min ≤T 3,t ≤T 3 max (52)

Figure SMS_127
Figure SMS_127

Figure SMS_128
Figure SMS_128

式中，T1 min和T1 max分别为压缩机入口温度的下限和上限，

Figure SMS_131
和
Figure SMS_132
分别为压缩机入口压力的下限和上限，
Figure SMS_133
和
Figure SMS_130
分别为燃料质量流量的下限和上限，，
Figure SMS_134
和
Figure SMS_135
分别为汽轮机出口压力的下限和上限，αmin和αmax分布为空气燃料比的下限和上限，T3 min和T3 max分别为燃烧室温度的下限和上限，和T4 max分别为汽轮机出口温度的下限和上限，
Figure SMS_136
和
Figure SMS_129
分别为热电比的下限和上限。Where T 1 min and T 1 max are the lower and upper limits of the compressor inlet temperature, respectively.
Figure SMS_131
and
Figure SMS_132
are the lower and upper limits of the compressor inlet pressure, respectively.
Figure SMS_133
and
Figure SMS_130
are the lower and upper limits of the fuel mass flow rate, respectively,
Figure SMS_134
and
Figure SMS_135
are the lower and upper limits of the turbine outlet pressure, α min and α max are the lower and upper limits of the air-fuel ratio, T 3 min and T 3 max are the lower and upper limits of the combustion chamber temperature, and T 4 max are the lower and upper limits of the turbine outlet temperature,
Figure SMS_136
and
Figure SMS_129
are the lower and upper limits of the thermoelectric ratio respectively.

步骤304，以最小化系统运行成本为目标，建立最优能流模型的目标函数，可表示为：Step 304, with the goal of minimizing the system operation cost, establish the objective function of the optimal energy flow model, which can be expressed as:

Figure SMS_137
Figure SMS_137

式中，c1和c2分别为热电联产机组的燃料单价和普通发电机的燃煤单价，

Figure SMS_138
为普通发电机集合，
Figure SMS_139
为热电联产机组集合，fi,t和gj,t分别为t时间段联产机组i和普通发电机j的成本函数；Nt为设定计算周期内的优化采样时刻集分别表示为：In the formula, c1 and c2 are the fuel unit price of the cogeneration unit and the coal unit price of the ordinary generator, respectively.
Figure SMS_138
For a common generator set,
Figure SMS_139
is the set of cogeneration units, fi ,t and gj ,t are the cost functions of cogeneration unit i and ordinary generator j in time period t respectively; Nt is the set of optimized sampling time in the set calculation period and is expressed as:

Figure SMS_140
Figure SMS_140

Figure SMS_141
Figure SMS_141

式中，In the formula,

μ11，μ12和μ13分别为联产机组发电功率和燃料流量之间的折算系数，μ 11 , μ 12 and μ 13 are the conversion coefficients between the power generation and fuel flow of the cogeneration unit,

μ21，μ22和μ23分别为普通发电机发电功率和燃煤量之间的折算系数。基于此，考虑全动态的电热综合能源系统最优能流模型可表示为：μ 21 , μ 22 and μ 23 are the conversion coefficients between the power generation of ordinary generators and the amount of coal burned. Based on this, the optimal energy flow model of the fully dynamic electric-thermal integrated energy system can be expressed as:

Figure SMS_142
Figure SMS_142

s.t.式(5)-(7)，(9)-(13)，(31)，(37)-(41)，(47)-(54)s.t. formula (5)-(7), (9)-(13), (31), (37)-(41), (47)-(54)

式(15)-(17)，(19)-(20)，(35)-(36)，(42)-(46)(58)Formula (15)-(17), (19)-(20), (35)-(36), (42)-(46)(58)

式(23)-(26)，(37)-(41)，(56)-(57)Formula (23)-(26), (37)-(41), (56)-(57)

以图3所示系统为例，计算周期为24小时，优化时间间隔为20分钟，空间步长为250米，优化所得的热电联产机组的输出热功率如图4所示。Taking the system shown in FIG3 as an example, the calculation cycle is 24 hours, the optimization time interval is 20 minutes, the spatial step length is 250 meters, and the output thermal power of the optimized cogeneration unit is shown in FIG4.

考虑全动态的电热综合能源系统最优能流建模系统，基于考虑全动态的电热综合能源系统最优能流建模方法，系统包括：The optimal energy flow modeling system of the electric-thermal integrated energy system considering the full dynamics is based on the optimal energy flow modeling method of the electric-thermal integrated energy system considering the full dynamics. The system includes:

数据采集模块、能源管理分析模块、逻辑计算模块与供热模型模块；Data acquisition module, energy management analysis module, logic calculation module and heating model module;

数据采集模块用于采集电热综合能源系统数据；The data acquisition module is used to collect data of the electric and thermal integrated energy system;

能源管理分析模块用于建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程；The energy management analysis module is used to establish simplified equations for nonlinear AC power flow, discrete equations for dynamic combined heat and power units and thermal systems;

逻辑计算模块用于基于简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件；The logic calculation module is used to establish the operation safety constraints of the electric-thermal integrated energy system based on the simplified AC power flow equation, discrete cogeneration unit and thermal system equation, with the goal of minimizing the system operation cost;

最优能流建模模块用于根据运行安全约束条件和电热综合能源系统模型建立考虑全动态的最优能流模型。The optimal energy flow modeling module is used to establish an optimal energy flow model that takes into account the full dynamics based on the operational safety constraints and the electric and thermal integrated energy system model.

本发明的有益效果在于，与现有技术相比，该方法全面的对电热综合能源系统的全动态和非线性特征减小建模，有利于精确刻画系统状态，从而准确制定电热综合能源系统运行控制策略，提高系统经济性和安全性。The beneficial effect of the present invention is that, compared with the prior art, this method comprehensively reduces the modeling of the full dynamic and nonlinear characteristics of the electric-thermal integrated energy system, which is conducive to accurately characterizing the system state, thereby accurately formulating the operation control strategy of the electric-thermal integrated energy system and improving the system economy and safety.

本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。The present disclosure may be a system, a method and/or a computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其它自由传播的电磁波、通过波导或其它传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。A computer-readable storage medium may be a tangible device that can hold and store instructions used by an instruction execution device. A computer-readable storage medium may be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples of computer-readable storage media (a non-exhaustive list) include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a static random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanical encoding device, such as a punch card or a raised structure in a groove on which instructions are stored, and any suitable combination of the foregoing. As used herein, a computer-readable storage medium is not to be interpreted as a transient signal per se, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagating through a waveguide or other transmission medium (e.g., a light pulse through a fiber optic cable), or an electrical signal transmitted through a wire.

这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing/processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and/or a wireless network. The network can include copper transmission cables, optical fiber transmissions, wireless transmissions, routers, firewalls, switches, gateway computers, and/or edge servers. The network adapter card or network interface in each computing/processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing/processing device.

用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。The computer program instructions for performing the operation of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages such as Smalltalk, C++, etc., and conventional procedural programming languages such as "C" language or similar programming languages. Computer-readable program instructions may be executed completely on a user's computer, partially on a user's computer, as an independent software package, partially on a user's computer, partially on a remote computer, or completely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer via any type of network including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., using an Internet service provider to connect via the Internet). In some embodiments, an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), may be personalized by utilizing the state information of the computer-readable program instructions, and the electronic circuit may execute the computer-readable program instructions, thereby realizing various aspects of the present disclosure.

这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。Various aspects of the present disclosure are described herein with reference to the flowcharts and/or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and/or block diagram and the combination of each box in the flowchart and/or block diagram can be implemented by computer-readable program instructions.

这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其它设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device that implements the functions/actions specified in one or more boxes in the flowchart and/or block diagram is generated. These computer-readable program instructions can also be stored in a computer-readable storage medium, and these instructions cause the computer, programmable data processing device, and/or other equipment to work in a specific manner, so that the computer-readable medium storing the instructions includes a manufactured product, which includes instructions for implementing various aspects of the functions/actions specified in one or more boxes in the flowchart and/or block diagram.

也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operating steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions/actions specified in one or more boxes in the flowchart and/or block diagram.

附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system, method and computer program product according to multiple embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a part of a module, program segment or instruction, and the part of the module, program segment or instruction contains one or more executable instructions for realizing the specified logical function. In some alternative implementations, the function marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and/or flow chart, and the combination of the square boxes in the block diagram and/or flow chart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of special hardware and computer instructions.

最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents, and any modifications or equivalent replacements that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims (12)

1.考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，该方法包括以下步骤：1. An optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system is characterized in that the method comprises the following steps: 步骤1，采集电热综合能源系统数据，包括管道长度、热阻、用户温度舒适度区间、电负荷时序分布；Step 1: Collect data of the electric and thermal integrated energy system, including pipeline length, thermal resistance, user temperature comfort range, and time series distribution of electric load; 步骤2，结合电热综合能源系统数据建立考虑全动态的电热综合能源系统模型；Step 2, combining the electric-thermal integrated energy system data to establish an electric-thermal integrated energy system model that takes into account the full dynamics; 步骤3，建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程，并简化交流潮流方程、离散的热电联产机组和热力系统方程；Step 3, establishing simplified equations of nonlinear AC power flow, discrete equations of dynamic combined heat and power units and thermal systems, and simplifying AC power flow equations, discrete combined heat and power units and thermal system equations; 步骤4，基于步骤2所述简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件，构造最优能流模型。Step 4, based on the simplified AC power flow equation, discrete cogeneration unit and thermal system equation described in step 2, establish the operation safety constraints of the electric and thermal integrated energy system with the goal of minimizing the system operation cost, and construct the optimal energy flow model. 2.根据权利要求1所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，2. The optimal energy flow modeling method for the fully dynamic electric-thermal integrated energy system according to claim 1 is characterized in that: 步骤2具体包括：步骤201，建立交流潮流模型：Step 2 specifically includes: Step 201, establishing an AC power flow model: 步骤202，建立包含压缩机模型、燃烧室模型、汽轮机模型、换热器模型的热电联产机组动态模型；Step 202, establishing a dynamic model of a cogeneration unit including a compressor model, a combustion chamber model, a steam turbine model, and a heat exchanger model; 步骤203，建立热力系统动态模型，包括水力部分模型和热力部分模型；其中水力部分模型为：Step 203, establish a dynamic model of the thermal system, including a hydraulic part model and a thermal part model; wherein the hydraulic part model is: Am＝dAm＝d BΔp＝0BΔp＝0 Δp＝Km2 Δp＝Km 2 式中，A和B分别为热力系统的节点-支路关联矩阵和回路-支路关联矩阵，m热水管道的质量流量向量，d为节点注入的质量流量向量，Δp为管道压降向量，K为管道摩阻系数；Where A and B are the node-branch association matrix and loop-branch association matrix of the thermal system, respectively; m is the mass flow vector of the hot water pipeline; d is the mass flow vector of the node injection; Δp is the pipeline pressure drop vector; and K is the pipeline friction coefficient; 热力部分模型为：The thermal model is:
Figure FDA0003905277850000021
Figure FDA0003905277850000021
Figure FDA0003905277850000022
Figure FDA0003905277850000022
Figure FDA0003905277850000023
Figure FDA0003905277850000023
 式中，x和t分别为空间和时间变量，v为水流流速，Tp表示管道温度，Ti n表示节点i的温度，mb为管道b的质量流量，
Figure FDA0003905277850000024
为管道k的出口温度，
Figure FDA0003905277850000025
为管道k的入口温度，λ为管道热阻系数，Ta为环境温度，
Figure FDA0003905277850000026
和
Figure FDA0003905277850000027
分别表示以节点i为首节点和末节点管道集。Where x and t are spatial and temporal variables, respectively; v is the water velocity; Tp represents the pipe temperature ; Tin represents the temperature of node i; mb is the mass flow rate of pipe b;
Figure FDA0003905277850000024
is the outlet temperature of pipe k,
Figure FDA0003905277850000025
is the inlet temperature of pipe k, λ is the thermal resistance coefficient of the pipe, Ta is the ambient temperature,
Figure FDA0003905277850000026
and
Figure FDA0003905277850000027
They represent the pipeline sets with node i as the first node and the last node respectively. 3.根据权利要求2所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，3. The optimal energy flow modeling method for the fully dynamic electric-thermal integrated energy system according to claim 2 is characterized in that: 压缩机模型为：The compressor model is: p2,t＝CPR1×p1,t p2,t ＝ CPR1 ×p1 ,t
Figure FDA0003905277850000028
Figure FDA0003905277850000028
 Pc,t＝Cama,t(T2,t-T1,t)P c,t = C a m a,t (T 2,t -T 1,t ) 式中，t为时间标志，p1,t和p2,t分别表示t时刻压缩机的入口和出口压力，T1,t和T2,t分别表示t时刻压缩机的入口和出口温度，CPR1表示压缩机的压比，β1表示空气绝热系数，η1表示压缩机效率，Pc,t表示t时刻压缩机耗功，Ca表示空气比热容，ma,t表示流入空气的质量流量。Wherein, t is the time mark, p1,t and p2,t represent the inlet and outlet pressures of the compressor at time t respectively, T1,t and T2 ,t represent the inlet and outlet temperatures of the compressor at time t respectively, CPR1 represents the pressure ratio of the compressor, β1 represents the air adiabatic coefficient, η1 represents the compressor efficiency, Pc,t represents the power consumption of the compressor at time t, Ca represents the specific heat capacity of air, and ma ,t represents the mass flow rate of the inflowing air. 4.根据权利要求2所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，4. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 2 is characterized in that: 燃烧室模型为：The combustion chamber model is:
Figure FDA0003905277850000031
Figure FDA0003905277850000031
 式中，T3,t分别表示t时刻的燃烧室温度，β2表示燃烧室的储热系数，Hg表示燃烧室热值，LHV表示燃料的低位热值，mf,t表示流入燃料的质量流量，Cs为混合烟气的比热容。Wherein, T3 ,t represents the combustion chamber temperature at time t, β2 represents the heat storage coefficient of the combustion chamber, Hg represents the combustion chamber calorific value, LHV represents the lower heating value of the fuel, mf ,t represents the mass flow rate of the inflowing fuel, and Cs is the specific heat capacity of the mixed flue gas. 5.根据权利要求2所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，5. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 2 is characterized in that: 汽轮机模型为：The turbine model is: p3,t＝CPR2p2,t p 3,t = CPR 2 p 2,t
Figure FDA0003905277850000032
Figure FDA0003905277850000032
 Pb,t＝Cs(ma,t+mf,t)(T4,t-T3,t)P b,t =C s (m a,t +m f,t )(T 4,t -T 3,t ) 式中，T4,t分别表示t时刻的汽轮机出口温度，β3表示燃烧室的绝热效率，Pb,t为t时刻汽轮机生产的总功率，CPR2表示汽轮机的压比。Where T4 ,t represents the turbine outlet temperature at time t, β3 represents the adiabatic efficiency of the combustion chamber, Pb,t is the total power produced by the turbine at time t, and CPR2 represents the pressure ratio of the turbine. 6.根据权利要求2所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，6. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 2 is characterized in that: 联产机组生产的用于发电和供热的功率分别为：The power produced by the combined unit for electricity and heat supply is: Pg,t＝η2η3(Pb,t-Pc,t)P g,t = η 2 η 3 (P b,t -P c,t ) φg,t＝η2(1-η3)(Pb,t-Pc,t)φ g,t =η 2 (1-η 3 )(P b,t -P c,t ) 式中，Pg,t和φg,t分别为联产机组用于发电和供热的功率，η2表示汽轮机机械效率，η3表示热电比。Where P g,t and φ g,t are the power of the cogeneration unit for power generation and heat supply, respectively, η 2 represents the mechanical efficiency of the steam turbine, and η 3 represents the heat-to-electricity ratio. 7.根据权利要求2所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，换热器模型为：7. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 2 is characterized in that the heat exchanger model is:
Figure FDA0003905277850000033
Figure FDA0003905277850000033
 式中，β4表示换热器的储热系数，Cw表示水的比热容，mw,t为换热器中水流的质量流量，T5,t和T6,t分别表示换热器的出口和入口温度。Where β4 represents the heat storage coefficient of the heat exchanger, Cw represents the specific heat capacity of water, mw ,t is the mass flow rate of water in the heat exchanger, and T5,t and T6,t represent the outlet and inlet temperatures of the heat exchanger, respectively. 8.根据权利要求1所述的一种考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，所述步骤3包括：8. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 1, wherein step 3 comprises: 步骤301，在输网级系统中，令节点电压的相角差为0，令节点电压幅值为1对非线性交流潮流模型进行简化得到电力系统潮流模型：Step 301, in the transmission grid system, the phase angle difference of the node voltage is set to 0, and the node voltage amplitude is set to 1 to simplify the nonlinear AC power flow model to obtain the power system power flow model:
Figure FDA0003905277850000041
Figure FDA0003905277850000041
Figure FDA0003905277850000042
Figure FDA0003905277850000042
Figure FDA0003905277850000043
Figure FDA0003905277850000043
Figure FDA0003905277850000044
Figure FDA0003905277850000044
 式中，i和j分别表示节点编号，
Figure FDA0003905277850000045
表示节点i的电压幅值，PG,i和PL,i表示节点i的发电机有功功率和负荷有功功率，QG,i和QL,i表示节点i的发电机无功功率和负荷无功功率，Gij和Bij表示节点i和节点j之间的电导与电纳，θij表示节点i和节点j之间的相角差，Pij和Qij表示节点i和节点j之间的传输的有功功率和无功功率；In the formula, i and j represent the node numbers,
Figure FDA0003905277850000045
represents the voltage amplitude of node i, PG,i and PL,i represent the generator active power and load active power of node i, QG,i and QL,i represent the generator reactive power and load reactive power of node i, Gij and Bij represent the conductance and susceptance between node i and node j, θij represents the phase angle difference between node i and node j, Pij and Qij represent the active power and reactive power transmitted between node i and node j; 步骤202，采用离散的时间和空间步长解法将热力系统模型中的热力部分模型离散化，得到离散化热力部分模型：Step 202, using a discrete time and space step solution method to discretize the thermal part model in the thermal system model to obtain a discretized thermal part model:
Figure FDA0003905277850000046
Figure FDA0003905277850000046
Figure FDA0003905277850000051
Figure FDA0003905277850000051
 式中，J1，J2，J3和J4分别为传递系数；
Figure FDA0003905277850000052
表示j时刻i+1处的管道温度，
Figure FDA0003905277850000053
表示j+1时刻i+1处的管道温度，
Figure FDA0003905277850000054
表示j时刻i处的管道温度，
Figure FDA0003905277850000055
表示j+1时刻i处的管道温度；Δx和Δt分别为离散的空间和时间步长，L和Γ分别为管道长度和时间区间，Nx和Nt分别为空间和时间步数，
Figure FDA0003905277850000056
Where, J 1 , J 2 , J 3 and J 4 are the transfer coefficients respectively;
Figure FDA0003905277850000052
represents the pipe temperature at time i+1,
Figure FDA0003905277850000053
represents the pipe temperature at i+1 at time j+1,
Figure FDA0003905277850000054
represents the pipe temperature at time i,
Figure FDA0003905277850000055
represents the pipe temperature at time i j+1; Δx and Δt are the discrete space and time steps, L and Γ are the pipe length and time interval, Nx and Nt are the space and time steps, respectively.
Figure FDA0003905277850000056
 步骤203，采用向后欧拉格式将燃烧室模型和换热器模型离散化；Step 203, discretizing the combustion chamber model and the heat exchanger model using a backward Euler format; 其中离散燃烧室模型为：The discrete combustion chamber model is:
Figure FDA0003905277850000057
Figure FDA0003905277850000057
 离散换热器模型为：The discrete heat exchanger model is:
Figure FDA0003905277850000058
Figure FDA0003905277850000058
 9.根据权利要求8所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于，9. The optimal energy flow modeling method for a fully dynamic electric-thermal integrated energy system according to claim 8 is characterized in that: 步骤4包括：Step 4 includes: 步骤401，建立电力系统运行安全约束，包括：Step 401, establishing power system operation safety constraints, including: 节点电压幅值约束：
Figure FDA0003905277850000059
相角约束：
Figure FDA00039052778500000510
Node voltage amplitude constraint:
Figure FDA0003905277850000059
Phase Angle Constraint:
Figure FDA00039052778500000510
 支路传输功率约束：
Figure FDA00039052778500000511
Branch transmission power constraints:
Figure FDA00039052778500000511
 发电机有功功率约束：
Figure FDA00039052778500000512
Generator active power constraint:
Figure FDA00039052778500000512
 发电机无功功率约束：
Figure FDA0003905277850000061
Generator reactive power constraints:
Figure FDA0003905277850000061
 式中，
Figure FDA0003905277850000062
表示节点i和节点j之间传输的视在功率上限，
Figure FDA0003905277850000063
和
Figure FDA0003905277850000064
为发电机i生产有功功率的下限和上限，
Figure FDA0003905277850000065
和
Figure FDA0003905277850000066
为发电机i生产无功功率的下限和上限；In the formula,
Figure FDA0003905277850000062
represents the apparent upper limit of the power transmitted between node i and node j,
Figure FDA0003905277850000063
and
Figure FDA0003905277850000064
are the lower and upper limits of active power produced by generator i,
Figure FDA0003905277850000065
and
Figure FDA0003905277850000066
The lower and upper limits of reactive power produced by generator i; 步骤402，建立热力系统运行安全约束，包括：Step 402, establishing thermal system operation safety constraints, including: 节点压力约束：
Figure FDA0003905277850000067
Nodal pressure constraints:
Figure FDA0003905277850000067
 管道和节点的质量流量的容量约束：
Figure FDA0003905277850000068
Capacity constraints for mass flow rates at pipes and nodes:
Figure FDA0003905277850000068
 管道质量流量变化率约束：γminmi,t-1≤mi,t≤γmaxmi,t-1,γmindi,t-1≤di,t≤γmaxdi,t-1；Pipeline mass flow rate change rate constraints: γ min m i,t-1 ≤m i,t ≤γ max m i,t-1 ,γ min d i,t-1 ≤d i,t ≤γ max d i,t-1 ; 节点供水温度约束：
Figure FDA0003905277850000069
Node water supply temperature constraints:
Figure FDA0003905277850000069
 节点回水温度约束：
Figure FDA00039052778500000610
Node return water temperature constraint:
Figure FDA00039052778500000610
 式中，
Figure FDA00039052778500000611
和
Figure FDA00039052778500000612
分别为节点i的水压的下限和上限，
Figure FDA00039052778500000613
和
Figure FDA00039052778500000614
分别为管道i的质量流量下限和上限，
Figure FDA00039052778500000615
和
Figure FDA00039052778500000616
分别为节点i的质量流量的下限和上限，γmin和γmax分别为质量流量变化率的下限和上限，
Figure FDA00039052778500000617
和
Figure FDA00039052778500000618
分别为t时刻节点i供水温度的下限和上限，
Figure FDA00039052778500000619
和
Figure FDA00039052778500000620
分别为t时刻节点i回水温度的下限和上限；In the formula,
Figure FDA00039052778500000611
and
Figure FDA00039052778500000612
are the lower and upper limits of the water pressure at node i,
Figure FDA00039052778500000613
and
Figure FDA00039052778500000614
are the lower and upper limits of the mass flow rate of pipeline i, respectively.
Figure FDA00039052778500000615
and
Figure FDA00039052778500000616
are the lower and upper limits of the mass flow rate of node i, respectively; γ min and γ max are the lower and upper limits of the mass flow rate change rate, respectively;
Figure FDA00039052778500000617
and
Figure FDA00039052778500000618
are the lower and upper limits of the water supply temperature at node i at time t,
Figure FDA00039052778500000619
and
Figure FDA00039052778500000620
are the lower and upper limits of the return water temperature at node i at time t respectively; 步骤403，建立热电联产机组运行安全约束，包括：Step 403, establishing safety constraints for the operation of the cogeneration unit, including: 压缩机入口温度和压力约束：T1 min≤T1,t≤T1 max，
Figure FDA00039052778500000621
Compressor inlet temperature and pressure constraints: T 1 min ≤ T 1, t ≤ T 1 max ,
Figure FDA00039052778500000621
 燃烧室输入燃料的质量流量约束：
Figure FDA00039052778500000622
The mass flow rate constraint of the fuel input to the combustion chamber is:
Figure FDA00039052778500000622
 汽轮机出口压力约束：
Figure FDA00039052778500000623
Steam turbine outlet pressure constraint:
Figure FDA00039052778500000623
 燃料和空气混合比约束：αminmf,t≤ma,t≤αmaxmf,t；Fuel and air mixture ratio constraints: α min m f,t ≤m a,t ≤α max m f,t ; 燃烧室温度约束：
Figure FDA00039052778500000624
Combustion chamber temperature constraints:
Figure FDA00039052778500000624
 汽轮机出口温度约束：
Figure FDA0003905277850000071
Turbine outlet temperature constraint:
Figure FDA0003905277850000071
 热电比约束：
Figure FDA0003905277850000072
Thermoelectric ratio constraints:
Figure FDA0003905277850000072
 式中，T1 min和T1 max分别为压缩机入口温度的下限和上限，
Figure FDA0003905277850000073
和
Figure FDA0003905277850000074
分别为压缩机入口压力的下限和上限，
Figure FDA0003905277850000075
和
Figure FDA0003905277850000076
分别为燃料质量流量的下限和上限，
Figure FDA0003905277850000077
和
Figure FDA0003905277850000078
分别为汽轮机出口压力的下限和上限，αmin和αmax分布为空气燃料比的下限和上限，T3 min和T3 max分别为燃烧室温度的下限和上限，和T4 max分别为汽轮机出口温度的下限和上限，
Figure FDA0003905277850000079
和
Figure FDA00039052778500000710
分别为热电比的下限和上限；Where T 1 min and T 1 max are the lower and upper limits of the compressor inlet temperature, respectively.
Figure FDA0003905277850000073
and
Figure FDA0003905277850000074
are the lower and upper limits of the compressor inlet pressure, respectively.
Figure FDA0003905277850000075
and
Figure FDA0003905277850000076
are the lower and upper limits of the fuel mass flow rate, respectively.
Figure FDA0003905277850000077
and
Figure FDA0003905277850000078
are the lower and upper limits of the turbine outlet pressure, α min and α max are the lower and upper limits of the air-fuel ratio, T 3 min and T 3 max are the lower and upper limits of the combustion chamber temperature, and T 4 max are the lower and upper limits of the turbine outlet temperature,
Figure FDA0003905277850000079
and
Figure FDA00039052778500000710
are the lower and upper limits of the thermoelectric ratio, respectively; 步骤404，以最小化系统运行成本为目标，建立最优能流模型的目标函数：Step 404, with the goal of minimizing the system operation cost, establish the objective function of the optimal energy flow model:
Figure FDA00039052778500000711
Figure FDA00039052778500000711
 式中，F表示电热综合能源系统的总运行成本，minF表示电热综合能源系统的最小总运行成本，c1fi,t表示电热综合能源系统内第i个热电联产机组在t时间段的总运行成本；c2gj,t表示电热综合能源系统内第j个普通发电机在t时间段的总运行成本；c1和c2分别为热电联产机组的燃料单价和普通发电机的燃煤单价，
Figure FDA00039052778500000712
为普通发电机集合，
Figure FDA00039052778500000713
为热电联产机组集合，fi,t和gj,t分别为t时间段联产机组i和普通发电机j的成本函数；Nt为设定计算周期内的优化采样时刻集分别表示为：Where F represents the total operating cost of the electric-thermal integrated energy system, minF represents the minimum total operating cost of the electric-thermal integrated energy system, c 1 fi,t represents the total operating cost of the i-th cogeneration unit in the electric-thermal integrated energy system in time period t; c 2 g j,t represents the total operating cost of the j-th ordinary generator in the electric-thermal integrated energy system in time period t; c 1 and c 2 are the fuel unit price of the cogeneration unit and the coal unit price of the ordinary generator, respectively.
Figure FDA00039052778500000712
For a common generator set,
Figure FDA00039052778500000713
is the set of cogeneration units, fi ,t and gj ,t are the cost functions of cogeneration unit i and ordinary generator j in time period t respectively; Nt is the set of optimized sampling time in the set calculation period and is expressed as:
Figure FDA00039052778500000714
Figure FDA00039052778500000714
Figure FDA00039052778500000715
Figure FDA00039052778500000715
 式中，μ11，μ12和μ13分别为联产机组发电功率和燃料流量之间的折算系数，μ21，μ22和μ23分别为普通发电机发电功率和燃煤量之间的折算系数；PG,j表示节点j的发电机有功功率；PG,i表示节点i的发电机有功功率。In the formula, μ11 , μ12 and μ13 are the conversion coefficients between the power generation of the cogeneration unit and the fuel flow, μ21 , μ22 and μ23 are the conversion coefficients between the power generation of the ordinary generator and the coal consumption, PG,j represents the active power of the generator at node j, and PG,i represents the active power of the generator at node i. 10.考虑全动态的电热综合能源系统最优能流建模系统，基于权利要求1-9中任一项所述的考虑全动态的电热综合能源系统最优能流建模方法，其特征在于所述系统包括：10. An optimal energy flow modeling system for an electric-thermal integrated energy system considering full dynamics, based on an optimal energy flow modeling method for an electric-thermal integrated energy system considering full dynamics as described in any one of claims 1 to 9, characterized in that the system comprises: 数据采集模块、能源管理分析模块、逻辑计算模块与供热模型模块；Data acquisition module, energy management analysis module, logic calculation module and heating model module; 数据采集模块用于采集电热综合能源系统数据；The data acquisition module is used to collect data of the electric and thermal integrated energy system; 能源管理分析模块用于建立非线性交流潮流的简化方程、动态热电联产机组和热力系统的离散方程；The energy management analysis module is used to establish simplified equations for nonlinear AC power flow, discrete equations for dynamic combined heat and power units and thermal systems; 逻辑计算模块用于基于简化的交流潮流方程、离散的热电联产机组和热力系统方程，以最小化系统运行成本为目标建立电热综合能源系统的运行安全约束条件；The logic calculation module is used to establish the operation safety constraints of the electric and thermal integrated energy system based on the simplified AC power flow equation, discrete cogeneration unit and thermal system equation, with the goal of minimizing the system operation cost; 最优能流建模模块用于根据运行安全约束条件和电热综合能源系统模型建立考虑全动态的最优能流模型。The optimal energy flow modeling module is used to establish an optimal energy flow model that takes into account the full dynamics based on the operational safety constraints and the electric and thermal integrated energy system model. 11.一种终端，包括处理器及存储介质；其特征在于：11. A terminal comprising a processor and a storage medium; characterized in that: 所述存储介质用于存储指令；The storage medium is used to store instructions; 所述处理器用于根据所述指令进行操作以执行根据权利要求1-9任一项所述方法的步骤。The processor is configured to operate according to the instructions to execute the steps of the method according to any one of claims 1-9. 12.计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-9任一项所述方法的步骤。12. A computer-readable storage medium having a computer program stored thereon, wherein when the program is executed by a processor, the steps of the method according to any one of claims 1 to 9 are implemented.
CN202211302252.0A 2022-10-24 2022-10-24 Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics Active CN115906411B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211302252.0A CN115906411B (en) 2022-10-24 2022-10-24 Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211302252.0A CN115906411B (en) 2022-10-24 2022-10-24 Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics

Publications (2)

Publication Number Publication Date
CN115906411A true CN115906411A (en) 2023-04-04
CN115906411B CN115906411B (en) 2024-06-04

Family

ID=86480336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211302252.0A Active CN115906411B (en) 2022-10-24 2022-10-24 Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics

Country Status (1)

Country Link
CN (1) CN115906411B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150032278A1 (en) * 2013-07-25 2015-01-29 International Business Machines Corporation Managing devices within micro-grids
CN110752605A (en) * 2019-10-29 2020-02-04 国家电网有限公司 An optimal power flow calculation method for an electrothermal coupled integrated energy system
CN110866627A (en) * 2019-08-16 2020-03-06 东南大学 Multi-region electricity-gas coupling comprehensive energy system optimal scheduling method considering step gas price
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
CN111815042A (en) * 2020-06-30 2020-10-23 天津大学 Optimal scheduling method of electric-heat integrated energy system considering refined heat network model
CN112989612A (en) * 2021-03-18 2021-06-18 贵州电网有限责任公司 Linear modeling method for electric heating comprehensive energy system based on Mackek envelope
CN113111555A (en) * 2021-04-14 2021-07-13 东南大学 Mass-adjustment thermodynamic system energy flow rapid calculation method based on superposition decoupling method
CN113725915A (en) * 2021-08-26 2021-11-30 国网浙江省电力有限公司经济技术研究院 Rural electric heating comprehensive energy system operation optimization method considering renewable energy uncertainty and thermal inertia
CN113837589A (en) * 2021-09-18 2021-12-24 清华大学 N-1 scheduling method and device for comprehensive energy system
CN114077934A (en) * 2022-01-19 2022-02-22 国网湖北省电力有限公司经济技术研究院 An integrated energy microgrid interconnection system and its scheduling method
CN114971147A (en) * 2022-03-23 2022-08-30 南京工程学院 Electric-heat-gas comprehensive energy system distribution robust optimization scheduling method based on Hausdorff distance

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150032278A1 (en) * 2013-07-25 2015-01-29 International Business Machines Corporation Managing devices within micro-grids
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
CN110866627A (en) * 2019-08-16 2020-03-06 东南大学 Multi-region electricity-gas coupling comprehensive energy system optimal scheduling method considering step gas price
CN110752605A (en) * 2019-10-29 2020-02-04 国家电网有限公司 An optimal power flow calculation method for an electrothermal coupled integrated energy system
CN111815042A (en) * 2020-06-30 2020-10-23 天津大学 Optimal scheduling method of electric-heat integrated energy system considering refined heat network model
CN112989612A (en) * 2021-03-18 2021-06-18 贵州电网有限责任公司 Linear modeling method for electric heating comprehensive energy system based on Mackek envelope
CN113111555A (en) * 2021-04-14 2021-07-13 东南大学 Mass-adjustment thermodynamic system energy flow rapid calculation method based on superposition decoupling method
CN113725915A (en) * 2021-08-26 2021-11-30 国网浙江省电力有限公司经济技术研究院 Rural electric heating comprehensive energy system operation optimization method considering renewable energy uncertainty and thermal inertia
CN113837589A (en) * 2021-09-18 2021-12-24 清华大学 N-1 scheduling method and device for comprehensive energy system
CN114077934A (en) * 2022-01-19 2022-02-22 国网湖北省电力有限公司经济技术研究院 An integrated energy microgrid interconnection system and its scheduling method
CN114971147A (en) * 2022-03-23 2022-08-30 南京工程学院 Electric-heat-gas comprehensive energy system distribution robust optimization scheduling method based on Hausdorff distance

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BO ZHAO等: "Optimal Scheduling Method for Electrical-Thermal Integrated Energy System Considering Heat Storage Characteristics of Heating Network", IEEE *
LIU, L.等: "Regionmodel and application of regional integrated energy system security analysis", APPL. ENERGY 260, 114268 *
李杭;王晗;侯恺;: "考虑多种能源耦合运行方式的综合能源潮流计算方法", 分布式能源, no. 04 *
黎静华;黄玉金;张鹏;: "综合能源系统多能流潮流计算模型与方法综述", 电力建设, no. 03 *

Also Published As

Publication number Publication date
CN115906411B (en) 2024-06-04

Similar Documents

Publication Publication Date Title
CN108494015B (en) An integrated energy system design method for coordinated interaction of source-load-storage
Xu et al. Quantification of flexibility of a district heating system for the power grid
Tian et al. System modeling and optimal dispatching of multi-energy microgrid with energy storage
CN108960503B (en) Multi-scenario optimization analysis method of integrated energy system based on interior point method
CN110175311B (en) An Optimal Power Flow Calculation Method Based on Multi-energy Coupling Model
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN108241782B (en) A hybrid transient analysis method and hybrid transient analysis system
CN111428351B (en) Power flow calculation method for electric-thermal integrated energy system based on forward push-back method
CN112989612B (en) Electric heating comprehensive energy system linear modeling method based on Mickey envelope
CN110443467A (en) A kind of regional complex energy resource system solar energy digestion capability appraisal procedure
CN110322068A (en) A kind of integrated energy system energy transition equipment configuration method and configuration system
Mu et al. Optimal operation model of integrated energy system for industrial plants considering cascade utilisation of heat energy
CN113919754A (en) A distributed state estimation method for integrated energy systems based on blockchain
CN112926835A (en) Comprehensive energy system optimization scheduling method considering dynamic characteristics of heat supply network
CN110020506B (en) Differential format selection method based on operation optimization of electric heating type comprehensive energy system
CN111310310A (en) A fast decoupling calculation method for static power flow of thermal system for quantity regulation
Zhou et al. Digital twin modeling and simulation of distributed and centralized integrated energy system on CloudPSS-IESLab
CN115906411B (en) Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics
Lei et al. A novel methodology for electric-thermal mixed power flow simulation and transmission loss analysis in multi-energy micro-grids
CN119129884A (en) A method, device and medium for optimizing planning of electric-thermal integrated energy system considering seasonal heat storage
CN110516868B (en) Comprehensive energy system optimization operation model considering network constraints
WO2024108642A1 (en) Operation optimization method for electrical integrated energy system, and related device
CN118199023A (en) Construction method of urban comprehensive energy system station network collaborative planning model
CN109726906B (en) Heat and power cogeneration system day-ahead scheduling method based on heat supply network partial differential equation constraint
CN110737993A (en) Operational Boundary Analysis Method of Multi-energy Complementary System Considering Load Uncertainty

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