CN109636148B - 基于能量网络方程的多能流系统的工作状态评估方法 - Google Patents

基于能量网络方程的多能流系统的工作状态评估方法 Download PDF

Info

Publication number
CN109636148B
CN109636148B CN201811443473.3A CN201811443473A CN109636148B CN 109636148 B CN109636148 B CN 109636148B CN 201811443473 A CN201811443473 A CN 201811443473A CN 109636148 B CN109636148 B CN 109636148B
Authority
CN
China
Prior art keywords
energy
network
hot water
equation
phi
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811443473.3A
Other languages
English (en)
Other versions
CN109636148A (zh
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.)
Guangdong Boshen Think Tank Energy Technology Development Co ltd
South China University of Technology SCUT
Original Assignee
Guangdong Boshen Think Tank Energy Technology Development Co ltd
South China University of Technology SCUT
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 Guangdong Boshen Think Tank Energy Technology Development Co ltd, South China University of Technology SCUT filed Critical Guangdong Boshen Think Tank Energy Technology Development Co ltd
Priority to CN201811443473.3A priority Critical patent/CN109636148B/zh
Publication of CN109636148A publication Critical patent/CN109636148A/zh
Application granted granted Critical
Publication of CN109636148B publication Critical patent/CN109636148B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Development Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Heat-Pump Type And Storage Water Heaters (AREA)

Abstract

本发明涉及一种基于能量网络方程的多能流系统的工作状态评估方法，包括：S1，利用热能和电能联合给N个独立的空间加热，建立能量网络，S2，利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律对所述能量网络建立能量网络方程组；S3，联立能量网络方程组与能量系统的约束条件，构成完备的方程组；对能量网络方程组采用标幺制计算方法，将能量网络方程组中不同单位的物理量转化为无量纲参数，统一求解，得到各个单元设备的工作状态。对于较复杂的能量网络，利用方案的方法建立统一的模型和统一的能量网络方程并进行求解，能够很方便的得到网络中各个设备的工作状态。

Description

基于能量网络方程的多能流系统的工作状态评估方法
技术领域
本发明涉及综合能源系统分析技术领域，特别是涉及一种基于能量网络方程的多能流系统的工作状态评估方法。
背景技术
风电、光伏等可再生清洁能源的使用成了一种必然趋势,而这些可再生能源不确定性高，发电量不易控制，接入电网容易引起电网波动。因此，仅以电网作为能量传输的载体已渐渐不能满足需求。在这种情况下，以电、热、冷、气多种形式传输能量的能源互联网就有其优越性，目前能源互联网的示范也越来越多。在上述多能网络中，热电联供网络目前发展最为迅速，自20世纪开始，在全球范围内就开始逐步建立热网，目前，热网在丹麦、瑞典、德国、芬兰等欧洲国家都有了一定程度的普及。目前，对于电-热耦合多能流网络的建模和潮流计算已经有了一系列的研究成果。其中，通过评估单元设备的工作状态来对能量系统的分析和优化提供理论指导。现有的评估单元设备的工作状态通常分别采用电网络理论和流体网络理论计算所有设备的工作状况，这种计算方法复杂。而事实上任何形式能量的传递过程都是在基本强度量差的推动下基本广延量的传递，不同类型的能量传输网络具有相同的传递规律。目前行业内还没有研发一种对能源网络内的能源形势建立统一的模型和统一的能量网络方程来求解单元设备的工作状态的方法。
发明内容
针对现有技术存在的上述问题，本发明提供一种基于能量网络方程的多能流系统的工作状态评估方法。
本申请的具体方案如下：
一种基于能量网络方程的多能流系统的工作状态评估方法，包括：
S1，利用热能和电能联合给N个独立的空间加热，建立能量网络，其中，N≥2；所述能量网络包括热水子网络和电力子网络；
S2，利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律对所述能量网络建立能量网络方程组；
S3，联立能量网络方程组与能量系统的约束条件，构成完备的方程组；对能量网络方程组采用标幺制计算方法，将能量网络方程组中不同单位的物理量转化为无量纲参数，统一求解，得到各个单元设备的工作状态。
优选地，N＝2，所述电力子网络包括：电源b1、水泵b2、电网络中线路的等效电阻b5组成的回路C1；电源b1、设置在空间K内的电热器b3、电网络中线路的等效电阻b5组成的回路C2；电源b1、设置在空间L内的电热器b4、电网络中线路的等效电阻b5组成的回路C3；所述热水子网络包括：热源b6、设置在空间K内的散热片b7、热水网络中管网的等效流阻b9和水泵b10组成的回路C4；热源b6、设置在空间L内的的散热片b8、热水网络中管网的等效流阻b9和水泵b10组成的回路C5；
在热水子网络和电力子网络中，选定节点n1、n3、n4、n6、n7利用广义基尔霍夫定律写出5个通量方程，5个通量方程依次为Φ1-Φ2-Φ3-Φ4＝0，Φ5-Φ1＝0，Φ6-Φ7-Φ8＝0，Φ9-Φ10＝0，Φ10-Φ6＝0；其中电力子网络的通量Φ是电流强度，热水子网络的通量Φ是热水体积流量；
将5个通量方程写成能量网络方程，分别为网络关联矩阵A、广延量流量矩阵Φ和矩阵方程A·Φ＝0的形式，其中，
Figure BDA0001885197200000031
Φ＝[Φ1 Φ2 Φ3 Φ4 Φ5 Φ6 Φ7 Φ8 Φ9 Φ10]T；
在热水子网络和电力子网络中，针对回路C1、C2、C3、C4、C5利用广义基尔霍夫定律写出5个强度量方程，5个强度量方程为△X1+△X2+△X5＝0，-△X2+△X3＝0，-△X3+△X4＝0，△X6+△X7+△X9+△X10＝0，-△X7+△X8＝0；其中电力子网络的强度量X是电压，热水子网络的强度量X是压强，
将5个强度量方程写成网络关联矩阵B、强度量差矩阵△X和矩阵方程B·ΔX＝0；其中
Figure BDA0001885197200000032
△X＝[△X1 △X2 △X3 △X4 △X5 △X6 △X7 △X8 △X9 △X10]T；
所述能量网络中第一类约束条件是两个能量子网之间的耦合关系：水泵将电能转化为热水的压能Φ2△X2+Φ10△X10＝0；已知两个空间的热负荷分别是PK和PL，热水网所占热负荷比例分别是αK和αL；则两个散热片的输出功率分别是αKPK＝ρcΦ7△T和αLPL＝ρcΦ8△T，其中ρ和c分别是热水的密度和质量热容，ΔT是散热片进出口的温差；电热器的热负荷分别是(1-αK)PK＝Φ3△X3和(1-αL)PL＝Φ4△X4；
所述能量网络中第二类约束条件已知的电源电动势△X1以及设备的特征方程：由传递公理得设备的特性方程△Xi＝RiΦi(i＝5,6,7,9)，其中的Ri为电阻或流阻，由设备的特性决定；在稳定的工作状态，所有并联散热片的进出口压差相同，由一台散热片的特性可得到推出压差；
利用上述方程组和上述约束条件采用标幺制计算方法求解，得到能量网络中所有设备的工作状况，所述工作状况为广延量通量和设备两端的强度量差。
优选地，其特征在于，能量系统的约束条件如下：
电源电压是△X1＝-220V，线路等效电阻是R5＝1.0Ω。散热片的流阻分别是R6＝1.37×108Pa·s/m3，R7＝1.92×108Pa·s/m3，热水管网的等效流阻是R9＝0.17×108Pa·s/m3。两个空间的热负荷分别是PK＝20kW和PL＝25kW，热水加热所占比例都是α＝0.8，散热片进回水温差都是△T＝15K；常温下水的密度和质量热容量分别取ρ＝980kg/m3和c＝4190J/(kg·K)。由此可知两个散热片的热水体积流量分别是Φ7＝0.26×10-3m3/s和Φ8＝0.325×10-3m3/s，电热器的热负荷Φ3△X3＝4.0kW和Φ4△X4＝5.0kW；
选定能量网络参数的基准值如下，功率Pb＝500W，电压△Xb,e＝200V，电流Φb,e＝Pb/△Xb,e＝2.5A，电阻Rb,e＝△Xb,e/Φb,e＝80Ω，体积流量Φb,f＝0.5×10-3m3/s，压强△Xb,f＝Pb/Φb,f＝106Pa，流阻Rb,f＝2×109Pa·s/m3。
利用所述基准值写出已知参数的无量纲标幺值：
Figure BDA0001885197200000041
Figure BDA0001885197200000042
Figure BDA0001885197200000043
Figure BDA0001885197200000044
利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律建立能量网络方程组，联立能量网络方程与能量系统的约束条件，构成完备的方程组，将标幺化后的参数值代入矩阵方程组并利用约束方程可解得：
Φ*＝[22.025 0.198 9.701 12.126 22.025 1.170 0.520 0.650 1.170 1.170]T
△X*＝[-1.100 0.825 0.825 0.825 0.275 0.080 0.050 0.050 0.010 -0.140]T
将这些无量纲参数乘以基准值，统一求解，得到能量网络中各设备的工作状态：
Φ＝[55.06 0.50 24.25 30.31 55.06 0.585×10-3 0.26×10-3 0.325×10-30.585×10-3 0.585×10-3]T
Figure BDA0001885197200000051
优选地，每个空间内设置有电热器和散热片，所述利用热能和电能联合给N个独立的空间加热包括：
热能通过管网给每个空间内的散热片提供热能，电能通过每个空间内的电热器转化成热能，电能同时向水泵提供电能，由水泵转化成热水的压能。
优选地，基尔霍夫强度量定律为在能量网络的任何一个回路中，各段线路强度量变化量之和为0；基尔霍夫广延量定律为对于能量网中任何一个节点，流入该节点的广延量等于流出该节点的广延量。
与现有技术相比，本发明具有如下有益效果：
根据论证，本方案的能量网络方程组的求解方法和利用电网络理论和流体网络理论分别计算单元设备的工作状态的结果一致，从而确保本发明建立的能量网络方程以及相应的求解方法是可靠的。对于较复杂的能量网络，利用方案的方法建立统一的模型和统一的能量网络方程并进行求解，能够很方便的得到网络中各个设备的工作状态。
附图说明
图1为本发明的基于能量网络方程的多能流系统的工作状态评估方法的示意性流程图；
图2为本发明的利用热能和电能联合给2个独立的空间加热，建立的能量网络图；
图3为本发明的电力子网络图。
图4为本发明的热水子网络图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1-4、一种基于能量网络方程的多能流系统的工作状态评估方法，包括：
S1，利用热能和电能联合给N个独立的空间加热，建立能量网络，其中，N≥2；所述能量网络包括热水子网络和电力子网络；每个空间内设置有电热器和散热片，所述利用热能和电能联合给N个独立的空间加热包括热能通过管网给每个空间内的散热片提供热能，电能通过每个空间内的电热器转化成热能，电能同时向水泵提供电能，由水泵转化成热水的压能。利用热能和电能联合供热是一种常见的用能方式，其中的用能侧可能有多个热负荷。在本实施例，N＝2，选择两个热负荷建立网络模型，但本发明的理论方法可以推广到多个热负荷系统(加热系统)。利用热水和电联合给两个独立的空间K和空间L加热，其中电源给两个电热器提供电能，热源通过管网给两个散热片提供热能，空间所需的热负荷(热能)由电热器和散热片分担。电力子网络将电源的电能传输到两个空间，通过两个电热器转化成热能。电力子网络同时向水泵提供电能，由水泵转化成热水的压能。热水网络将热源的热能通过散热片和管网分别传递到两个空间，同时消耗压能。两个能量子网络通过水泵和加热空间相互耦合。
S2，利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律对所述能量网络建立能量网络方程组；该方程组的参数是各个单元设备的通量和单元设备两端的强度量差。
S3，联立能量网络方程组与能量系统的约束条件，构成完备的方程组；对能量网络方程组采用标幺制计算方法，将能量网络方程组中不同单位的物理量转化为无量纲参数，统一求解，得到各个单元设备的工作状态，即可评估单元设备的工作状态。
热水子网络和电力子网络分别是流体网络和电网络，热水网络和电网络分别是流体网络和电网络。从能量公理和传递公理看，电网络是能量网络的一个特例，能量网络与电网络有相似的连接方式；电网络中的电压即能量网络的强度量(电压和压强)，电流即能量网络中的广延量(电流和体积流量)。在电网络中，基尔霍夫电压定律和基尔霍夫电流定律是电网络方程的基础，概括了电路中电流和电压分别遵循的基本规律，与电路元件的特性无关，只取决于电路的连接方式，二电路元件特性决定了支路两端电压差与电流的关系；因此，可以得到能量网络的基尔霍夫电压定律和基尔霍夫电流定律，定义为：
(1)基尔霍夫强度量定律：在能量网络的任何一个回路中，各段线(管)路强度量变化量之和为0；
(2)基尔霍夫广延量定律：对于能量网中任何一个节点，流入该节点的广延量等于流出该节点的广延量。
对于一个有B条支路的电网络，共有2B个变量，B个电压值和B个电流值；通过基尔霍夫定律可以写出B个方程，再加上B条支路上电路元件的路特性方程，联立可对该电网络实现求解；利用广义基尔霍夫定律可以写出能量网络方程，从而可以对能量网络中各个支路和节点参数进行求解；求解方法如下：
在本实施例，其特征在于，N＝2，所述电力子网络包括：电源b1、水泵b2、电网络中线路的等效电阻b5组成的回路C1；电源b1、设置在空间K内的电热器b3、电网络中线路的等效电阻b5组成的回路C2；电源b1、设置在空间L内的电热器b4、电网络中线路的等效电阻b5组成的回路C3；所述热水子网络包括：热源b6、设置在空间K内的散热片b7、热水网络中管网的等效流阻b9和水泵b10组成的回路C4；热源b6、设置在空间L内的散热片b8、热水网络中管网的等效流阻b9和水泵b10组成的回路C5；
在热水子网络和电力子网络中，选定节点n1、n3、n4、n6、n7利用广义基尔霍夫定律写出5个通量方程，5个通量方程依次为Φ1-Φ2-Φ3-Φ4＝0，Φ5-Φ1＝0，Φ6-Φ7-Φ8＝0，Φ9-Φ10＝0，Φ10-Φ6＝0；其中电力子网络的通量Φ是电流强度，热水子网络的通量Φ是热水体积流量；
将5个通量方程写成能量网络方程，分别为网络关联矩阵A、广延量流量矩阵Φ和矩阵方程A·Φ＝0的形式，其中，
Figure BDA0001885197200000081
Φ＝[Φ1 Φ2 Φ3 Φ4 Φ5 Φ6 Φ7 Φ8 Φ9 Φ10]T；
在热水子网络和电力子网络中，针对回路C1、C2、C3、C4、C5利用广义基尔霍夫定律写出5个强度量方程，5个强度量方程为△X1+△X2+△X5＝0，-△X2+△X3＝0，-△X3+△X4＝0，△X6+△X7+△X9+△X10＝0，-△X7+△X8＝0；其中电力子网络的强度量X是电压，热水子网络的强度量X是压强，
将5个强度量方程写成网络关联矩阵B、强度量差矩阵△X和矩阵方程B·ΔX＝0；其中
Figure BDA0001885197200000091
△X＝[△X1 △X2 △X3 △X4 △X5 △X6 △X7 △X8 △X9 △X10]T；
所述能量网络中第一类约束条件是两个能量子网之间的耦合关系：水泵将电能转化为热水的压能Φ2△X2+Φ10△X10＝0；已知两个空间的热负荷分别是PK和PL，热水网所占热负荷比例分别是αK和αL；则两个散热片的输出功率分别是αKPK＝ρcΦ7△T和αLPL＝ρcΦ8△T，其中ρ和c分别是热水的密度和质量热容，ΔT是散热片进出口的温差；电热器的热负荷分别是(1-αK)PK＝Φ3△X3和(1-αL)PL＝Φ4△X4；
所述能量网络中第二类约束条件已知的电源电动势△X1以及设备的特征方程：由传递公理得设备的特性方程△Xi＝RiΦi(i＝5,6,7,9)，其中的Ri为电阻或流阻，由设备的特性决定；在稳定的工作状态，所有并联散热片的进出口压差相同，由一台散热片的特性可得到推出压差；
利用上述方程组和上述约束条件采用标幺制计算方法求解，得到能量网络中所有设备的工作状况，所述工作状况为广延量通量和设备两端的强度量差。
比如，能量系统的约束条件如下：
电源电压是△X1＝-220V，线路等效电阻是R5＝1.0Ω。散热片的流阻分别是R6＝1.37×108Pa·s/m3，R7＝1.92×108Pa·s/m3，热水管网的等效流阻是R9＝0.17×108Pa·s/m3。两个空间的热负荷分别是PK＝20kW和PL＝25kW，热水加热所占比例都是α＝0.8，散热片进回水温差都是△T＝15K；常温下水的密度和质量热容量分别取ρ＝980kg/m3和c＝4190J/(kg·K)。由此可知两个散热片的热水体积流量分别是Φ7＝0.26×10-3m3/s和Φ8＝0.325×10-3m3/s，电热器的热负荷Φ3△X3＝4.0kW和Φ4△X4＝5.0kW；
选定能量网络参数的基准值如下，功率Pb＝500W，电压△Xb,e＝200V，电流Φb,e＝Pb/△Xb,e＝2.5A，电阻Rb,e＝△Xb,e/Φb,e＝80Ω，体积流量Φb,f＝0.5×10-3m3/s，压强△Xb,f＝Pb/Φb,f＝106Pa，流阻Rb,f＝2×109Pa·s/m3。
利用所述基准值写出已知参数的无量纲标幺值：
Figure BDA0001885197200000101
Figure BDA0001885197200000102
Figure BDA0001885197200000103
Figure BDA0001885197200000104
利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律建立能量网络方程组，联立能量网络方程与能量系统的约束条件，构成完备的方程组，将标幺化后的参数值代入矩阵方程组并利用约束方程可解得：
Φ*＝[22.025 0.198 9.701 12.126 22.025 1.170 0.520 0.650 1.170 1.170]T
△X*＝[-1.100 0.825 0.825 0.825 0.275 0.080 0.050 0.050 0.010 -0.140]T
将这些无量纲参数乘以基准值，统一求解，得到能量网络中各设备的工作状态：
Figure BDA0001885197200000106
Figure BDA0001885197200000105
根据论证，本方案的能量网络方程组的求解方法和利用电网络理论和流体网络理论分别计算单元设备的工作状态的结果一致，从而确保本发明建立的能量网络方程以及相应的求解方法是可靠的。
对于较复杂的能量网络，利用方案的方法建立统一的模型和统一的能量网络方程并进行求解，能够很方便的得到网络中各个设备的工作状态。在本方案，电热器和散热片的型号以及两种加热方式所占比例可以选择，每一种选择对应一种不同的工作状态，每一种工作状态又对应不同的能耗和投资，由此可以实现对该能量网络的评估、比较和优化。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (4)

1.一种基于能量网络方程的多能流系统的工作状态评估方法，其特征在于，包括：
S1，利用热能和电能联合给N个独立的空间加热，建立能量网络，其中，N≥2；所述能量网络包括热水子网络和电力子网络；
N＝2，所述电力子网络包括：电源b1、水泵b2和电网络中线路的等效电阻b5组成的回路C1，电源b1、设置在空间K内的电热器b3和电网络中线路的等效电阻b5组成的回路C2以及电源b1、设置在空间L内的电热器b4和电网络中线路的等效电阻b5组成的回路C3；
所述热水子网络包括：热源b6、设置在空间K内的散热片b7、热水网络中管网的等效流阻b9和水泵b10组成的回路C4以及热源b6、设置在空间L内的散热片b8、热水网络中管网的等效流阻b9和水泵b10组成的回路C5；
在热水子网络和电力子网络中，选定节点n1、n3、n4、n6和n7利用广义基尔霍夫定律写出5个通量方程，5个通量方程依次为Φ1-Φ2-Φ3-Φ4＝0，Φ5-Φ1＝0，Φ6-Φ7-Φ8＝0，Φ9-Φ10＝0，Φ10-Φ6＝0；其中电力子网络的通量Φ是电流强度，热水子网络的通量Φ是热水体积流量；
将5个通量方程写成能量网络方程，分别为网络关联矩阵A、广延量流量矩阵Φ和矩阵方程A·Φ＝0的形式，其中，
Figure FDA0003787299200000011
Φ＝[Φ1 Φ2 Φ3 Φ4 Φ5 Φ6 Φ7 Φ8 Φ9 Φ10]T；
在热水子网络和电力子网络中，针对回路C1、C2、C3、C4和C5利用广义基尔霍夫定律写出5个强度量方程，5个强度量方程为ΔX1+ΔX2+ΔX5＝0，-ΔX2+ΔX3＝0，-ΔX3+ΔX4＝0，ΔX6+ΔX7+ΔX9+ΔX10＝0，-ΔX7+ΔX8＝0；其中电力子网络的强度量X是电压，热水子网络的强度量X是压强，
将5个强度量方程写成网络关联矩阵B、强度量差矩阵ΔX和矩阵方程B·ΔX＝0；其中
Figure FDA0003787299200000021
ΔX＝[ΔX1 ΔX2 ΔX3 ΔX4 ΔX5 ΔX6 ΔX7 ΔX8 ΔX9 ΔX10]T；
所述能量网络中第一类约束条件是两个能量子网之间的耦合关系：水泵将电能转化为热水的压能Φ2ΔX2+Φ10ΔX10＝0；已知两个空间的热负荷分别是PK和PL，热水网所占热负荷比例分别是αK和αL；则两个散热片的输出功率分别是αKPK＝ρcΦ7ΔT和αLPL＝ρcΦ8ΔT，其中ρ和c分别是热水的密度和质量热容，ΔT是散热片进出口的温差；电热器的热负荷分别是(1-αK)PK＝Φ3ΔX3和(1-αL)PL＝Φ4ΔX4；
所述能量网络中第二类约束条件已知的电源电动势ΔX1以及设备的特征方程：由传递公理得设备的特性方程ΔXi＝RiΦi(i＝5,6,7,9)，其中的Ri为电阻或流阻，由设备的特性决定；在稳定的工作状态，所有并联散热片的进出口压差相同，由一台散热片的特性得到推出压差；
利用上述方程组和上述约束条件采用标幺制计算方法求解，得到能量网络中所有设备的工作状况，所述工作状况为广延量通量和设备两端的强度量差；
S2，利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律对所述能量网络建立能量网络方程组；
S3，联立能量网络方程组与能量系统的约束条件，构成完备的方程组；对能量网络方程组采用标幺制计算方法，将能量网络方程组中不同单位的物理量转化为无量纲参数，统一求解，得到各个单元设备的工作状态。
2.根据权利要求1所述的基于能量网络方程的多能流系统的工作状态评估方法，其特征在于，能量系统的第一类约束条件和第二类约束条件为：
电源电压是ΔX1＝-220V，线路等效电阻是R5＝1.0Ω，散热片的流阻分别是R6＝1.37×108Pa·s/m3，R7＝1.92×108Pa·s/m3，热水管网的等效流阻是R9＝0.17×108Pa·s/m3；两个空间的热负荷分别是PK＝20kW和PL＝25kW，热水加热所占比例都是α＝0.8，散热片进回水温差都是ΔT＝15K；常温下水的密度和质量热容量分别取ρ＝980kg/m3和c＝4190J/(kg·K)；散热片的热水体积流量分别是Φ7＝0.26×10-3m3/s和Φ8＝0.325×10-3m3/s，电热器的热负荷Φ3ΔX3＝4.0kW和Φ4ΔX4＝5.0kW；
选定能量网络参数的基准值如下，功率Pb＝500W，电压ΔXb,e＝200V，电流Φb,e＝Pb/ΔXb,e＝2.5A，电阻Rb,e＝ΔXb,e/Φb,e＝80Ω，体积流量Φb,f＝0.5×10-3m3/s，压强ΔXb,f＝Pb/Φb,f＝106Pa，流阻Rb,f＝2×109Pa·s/m3；
利用所述基准值写出已知参数的无量纲标幺值：
Figure FDA0003787299200000031
Figure FDA0003787299200000032
Figure FDA0003787299200000033
Figure FDA0003787299200000034
利用广义的基尔霍夫强度量定律和基尔霍夫广延量定律建立能量网络方程组，联立能量网络方程与能量系统的约束条件，构成完备的方程组，将标幺化后的参数值代入矩阵方程组并利用约束方程可解得：
Φ*＝[22.025 0.198 9.701 12.126 22.025 1.170 0.520 0.650 1.170 1.170]T
ΔX*＝[-1.100 0.825 0.825 0.825 0.275 0.080 0.050 0.050 0.010 -0.140]T
将这些无量纲参数乘以基准值，统一求解，得到能量网络中各设备的工作状态：
Φ＝[55.06 0.50 24.25 30.31 55.06 0.585×10-30.26×10-3 0.325×10-3 0.585×10-3 0.585×10-3]T
Figure FDA0003787299200000041
3.根据权利要求1所述的基于能量网络方程的多能流系统的工作状态评估方法，其特征在于，每个空间内设置有电热器和散热片，所述利用热能和电能联合给N个独立的空间加热包括：
热能通过管网给每个空间内的散热片提供热能，电能通过每个空间内的电热器转化成热能，电能同时向水泵提供电能，由水泵转化成热水的压能。
4.根据权利要求1所述的基于能量网络方程的多能流系统的工作状态评估方法，其特征在于，基尔霍夫强度量定律为在能量网络的任何一个回路中，各段线路强度量变化量之和为0；
基尔霍夫广延量定律为对于能量网中任何一个节点，流入该节点的广延量等于流出该节点的广延量。
CN201811443473.3A 2018-11-29 2018-11-29 基于能量网络方程的多能流系统的工作状态评估方法 Active CN109636148B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811443473.3A CN109636148B (zh) 2018-11-29 2018-11-29 基于能量网络方程的多能流系统的工作状态评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811443473.3A CN109636148B (zh) 2018-11-29 2018-11-29 基于能量网络方程的多能流系统的工作状态评估方法

Publications (2)

Publication Number Publication Date
CN109636148A CN109636148A (zh) 2019-04-16
CN109636148B true CN109636148B (zh) 2022-09-13

Family

ID=66069817

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811443473.3A Active CN109636148B (zh) 2018-11-29 2018-11-29 基于能量网络方程的多能流系统的工作状态评估方法

Country Status (1)

Country Link
CN (1) CN109636148B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111509707A (zh) * 2020-04-24 2020-08-07 贵州电网有限责任公司 一种能源互联网热力网络潮流标幺值模型建模方法
CN111859590B (zh) * 2020-07-21 2022-09-16 大连理工大学 一种构建热量输运热阻模型的方法
CN112580214A (zh) * 2020-12-23 2021-03-30 华北电力大学 广义多能源网络的能量传递统一建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102624812A (zh) * 2011-10-31 2012-08-01 李宗诚 全球价值链市场配置工程的ict 技术支持设计
WO2013026731A1 (en) * 2011-08-22 2013-02-28 Abb Technology Ag System and method to optimize operation of a water network
US9915399B1 (en) * 2017-04-18 2018-03-13 Air Products And Chemicals, Inc. Control system in a gas pipeline network to satisfy demand constraints
CN107808218A (zh) * 2017-10-25 2018-03-16 国网天津市电力公司 基于热电比调节的城市能源互联网潮流计算方法
CN107817681A (zh) * 2017-10-16 2018-03-20 清华大学 一种基于双侧等效模型的热网稳态运行状态估计方法
WO2018130231A1 (zh) * 2017-01-11 2018-07-19 东南大学 一种基于热网和房屋热惯性的综合能源系统优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10013512B2 (en) * 2014-08-29 2018-07-03 Schlumberger Technology Corporation Network flow model

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013026731A1 (en) * 2011-08-22 2013-02-28 Abb Technology Ag System and method to optimize operation of a water network
CN102624812A (zh) * 2011-10-31 2012-08-01 李宗诚 全球价值链市场配置工程的ict 技术支持设计
WO2018130231A1 (zh) * 2017-01-11 2018-07-19 东南大学 一种基于热网和房屋热惯性的综合能源系统优化方法
US9915399B1 (en) * 2017-04-18 2018-03-13 Air Products And Chemicals, Inc. Control system in a gas pipeline network to satisfy demand constraints
CN107817681A (zh) * 2017-10-16 2018-03-20 清华大学 一种基于双侧等效模型的热网稳态运行状态估计方法
CN107808218A (zh) * 2017-10-25 2018-03-16 国网天津市电力公司 基于热电比调节的城市能源互联网潮流计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于时变能量网络理论的综合能源系统动态仿真研究;陈皓勇等;《供用电》;20181105(第11期);全文 *
基于能量网络理论的区域综合能源系统运行优化研究;李明等;《分布式能源》;20180815(第04期);全文 *

Also Published As

Publication number Publication date
CN109636148A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109636148B (zh) 基于能量网络方程的多能流系统的工作状态评估方法
CN113190942B (zh) 供热/冷系统虚拟储能能力计算方法、装置与电子设备
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN109492325B (zh) 一种基于扩展能源集线器的多能耦合系统潮流分析方法
CN106339794A (zh) 一种电‑热耦合多能流网络节点能价计算方法
Macía et al. Influence parameters on the performance of an experimental solar-assisted ground-coupled absorption heat pump in cooling operation
Alptekin et al. Design and experimental investigation of a thermoelectric self-powered heating system
CN100585881C (zh) 一种晶体硅太阳能电池组件光电光热利用系统
Gou et al. Holistic analysis and optimization of distributed energy system considering different transport characteristics of multi-energy and component efficiency variation
RU2687382C1 (ru) Способ работы тепловой электрической станции и устройство для его реализации
Chen et al. Novel dynamic equivalent circuit model of integrated energy systems
RU2015147057A (ru) Солнечный проточный нагреватель с коллектором-аккумулятором
CN111656098B (zh) 用于能源网的改善的利用率的方法
Liu et al. A tool for integrated analysis of multi-vector district energy networks
Khanjari et al. Perormance evaluation of liquid-cooled photovoltaic system using numerical method
CN114549232A (zh) 电-热综合能源系统混合能流计算方法
Wasik et al. Optimization of PVT geometry in application of ground heat source regeneration
KR101515374B1 (ko) 열병합발전 시스템으로 구축된 개별 열/전기 에너지 및 중앙 열/전기 에너지 공급 연계 시스템, 및 이를 운용하는 방법
Hazi et al. Energy efficiency of the PVT system used in industry
JP5745911B2 (ja) 熱供給システム
CN112070274A (zh) 一种综合能源系统的*效率评估方法及系统
Chen et al. The application of energy network theory in the analysis of district electricity and heating system
Tang et al. Linearized modeling of integrated electricity and district heating systems with VF-VT strategy based on McCormick envelopes
Capocelli et al. Two-Channels Thermal Energy Storage Tank: Experiments and Short-Cut Modelling
Chen et al. Energy flow model analysis and configuration optimization of thermal management system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant