CN107808218A - 基于热电比调节的城市能源互联网潮流计算方法 - Google Patents

基于热电比调节的城市能源互联网潮流计算方法 Download PDF

Info

Publication number
CN107808218A
CN107808218A CN201711005289.6A CN201711005289A CN107808218A CN 107808218 A CN107808218 A CN 107808218A CN 201711005289 A CN201711005289 A CN 201711005289A CN 107808218 A CN107808218 A CN 107808218A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
msup
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.)
Pending
Application number
CN201711005289.6A
Other languages
English (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.)
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin 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 State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201711005289.6A priority Critical patent/CN107808218A/zh
Publication of CN107808218A publication Critical patent/CN107808218A/zh
Pending legal-status Critical Current

Links

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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

Abstract

本发明涉及一种基于热电比调节的城市能源互联网潮流计算方法，包括步骤有：(1)配气系统数学模型建立；(2)供热系统模型建立；(3)配电系统模型建立；(4)多能源耦合环节建模；(5)城市能源互联网统一潮流模型建立；(6)城市能源互联网优化潮流计算。本发明以城市能源互联网综合能源运营商的运营成本最低为目标，以热电比作为调节手段，建立了考虑城市配电网，配气网，供热系统以及能源集线器四部分约束的优化潮流模型并进行求解，可有效提高城市能源互联网运行的经济性，发明方法应用简单、快速，具有很好的实用价值。

Description

基于热电比调节的城市能源互联网潮流计算方法
技术领域
本发明属于城市能源互联网优化调控,多能源协同运行优化的技术领域，特别是一种基于热电比调节的城市能源互联网潮流计算方法。
背景技术
城市能源互联网由配电子系统、配气子系统、供冷/热子系统等供能子网络耦合形成，并以冷热电三联供机组为核心耦合环节，完成城市电气冷热多种能源的分配、转换和平衡。城市能源互联网可以分为配气系统，配电系统，供热系统及多能源系统耦合环节。其中，配气系统又包括城镇天然气门站，配气网络，储气装置，调压站，各类城镇燃气用户等等；配电系统包括配电变压器、架空线路或者城市电缆、各类电力用户；供热系统包括供热源、供热管网、用户、热转换设备；耦合环节包括冷热电三联供系统，电锅炉，燃气锅炉，热泵，燃气轮机等。城市能源互联网的发展和应用对于提高能源利用效率具有重要的意义。考虑到城市能源互联网所包含的能源体系多样，涉及的能源环节复杂，针对城市能源互联网的优化潮流分析研究极具价值。
最优潮流(Optimal Power Flow，OPF)是指当城市能源互联网的网络拓扑，源-荷信息，网络参数等信息给定时，通过对控制变量的调整同时满足系统的约束条件，通过优化算法使城市能源互联网的某一个目标函数(比如燃料费用最小，网损最小，用户费用最小，能源供应商收益最大等等)达到最大或者最小时的最优潮流分布。针对能源互联网的最优潮流计算方法，已有诸多前期研究，如文献Martinez-MaresA,Fuerte-Esquivel CR.Aunified gas andpower flow analysis in natural gas and electricity couplednetworks[J].IEEE Transactions on Power Systems,2012,27(4):2156-2166.和文献P,Leung T,Gil M,et al.Gas–electricity coordination in competitivemarkets under renewable energy uncertainty[J].IEEE Transactions on PowerSystems,2014,doi:10.1109/TPWRS.2014.2319588.从宏观角度分析能源互联网中电气与天然气系统之间的互动联系，但忽略了电-气耦合环节的一些细节；如文献Guan Xiaohong，Xu Zhanbo，Jia Qingshan.Energy-Efficient Buildings Facilitated by Microgrid[J].IEEE Transactions on Smart Grid，2010，1(3)：243-252，以运行成本最低为目标进行优化潮流计算，从而降低能源互联网的日前运行成本；如文献光伏与微型燃气轮机混合微网能量管理研究[J].电工技术学报，2012，27(1)：74-84，在分析光伏发电和微型燃气轮机动态模型及运行特性的基础上，研究了光伏与微型燃气轮机组成的能源互联网潮流策略；如文献冷热电联供微网优化调度通用建模方法[J].中国电机工程学报，2013，33(31)：26-33.采用通用母线方式，提出了冷热电联供能源互联系统的优化调度通用建模方法，并基于该方法构建了能源互联网系统的动态经济调度模型框架；如文献Hierar-chicalEnergy Management System for Multi-Source Multi-Product Microgrids[J].Renewable Energy，2015，78：621-630、Hier-archical Management for IntegratedCommunity Energy Systems[J].Applied Energy，2015，160：231-243.则分别针对含有热电联产系统、光伏、空调、电储能以及热储能的区域级能源互联网，设计了一种分层能量管理系统架构，能够实现对区域能源互联网中供能单元的统一潮流优化。
然而，上述研究未充分考虑城市能源互联网不同能源系统之间的复杂交互作用机理，在进行潮流计算时往往采用分立方法，系统整体全局优化能力较弱，导致运行成本依旧居高不下。
发明内容
本发明的目的是针对现有技术的不足，而提出一种基于热电比调节的城市能源互联网潮流计算方法。
本发明解决其技术问题是采取以下技术方案实现的：
一种基于热电比调节的城市能源互联网潮流计算方法，包括步骤如下：
(1)配气系统数学模型
①稳态流动方程
配气系统模型是管道的稳态流动方程，对于k节点到m节点的配气网稳态流量fkm，采用下式描述：
其中：pn为标况下的天然气压力基准；Ta为管道中天然气的平均温度；Dkm为管道直径；L为管道的长度；Tn为温度基准；pk为节点k的天然气压力；pm为节点m的天然气压力；γg，Za，F为和天然气有关的系数，可以查表获取；pk＞pm，sg(pk,pm)＝1，若pk＜pm，sg(pk,pm)＝-1；
根据式子(1)，将配气系统稳态流量f改写为与配气网天然气气压压力差ΔΠ有关的函数，表示为：
f＝φ(ΔΠ) (2)
其中：
②连续方程
天然气流量在某一节点处代数和为零，表示如下：
A1f＝L (4)
其中：A1为天然气管网的支路-节点矩阵，L为负荷；
同时规定符号如下：对于支路而言，如果流量顺着支路流入某一节点，则f为正，如果流量顺着支路流出某一节点，则f为负，对于负荷节点而言，L为正，对气源节点而言，L为负；
根据式(3)，将配气网气压平方改写为与配气网天然气气压压力差ΔΠ有关的函数，表示为：
ΔΠ＝-ATΠ (5)
其中A为天然气管网的支路-节点矩阵，
③回路方程
根据基尔霍夫第二定律，对于一个回路而言，回路中的管道压力降为零，表示为：
BΔΠ＝0 (6)
其中：B为天然气网络中的回路-支路矩阵，
(2)供热系统模型
①水动力模型
供热系统中对于某一节点，流入的热水质量等于流出的热水质量；对于某一回路，热水在网络中运动，由于摩擦的原因，水头压力会下降，因此对于一个热水回路而言，回路的压力降为0，表示为：
Asm＝mq (7)
Bhhf＝0 (8)
式中：As为热力管网的支路-节点矩阵；m为管道流量；mq为热力负荷的负荷需求；Bh为热力管网的支路-回路矩阵；hf为压力降，表示为
hf＝Km|m| (9)
式中：K为管道的阻力系数参量，由管道的材质决定，
②热动力模型
热动力模型描述的是供热系统的热平衡行为，用户的热负荷需求取决于供热的供热温度Ts，供热的输出温度To，以及单位时间内流入该热用户的热水质量mq，热力模型为：
Φ＝Cpmq(Ts-To) (10)
热水在网络的流动过程中，假设管道的始端温度为Tstart，末端温度为Tend，环境温度Ta，管道与温降有关的系数为λ，管道的长度为L，那么温降方程表示为：
热水在网络中发生混合须满足能量守恒定律，因此混合方程表述为：
(∑mout)Tout＝∑minTin (12)
(3)配电系统模型
配电系统的模型为反映节点功率与节点电压之间关系的节点功率方程：
其中：为发电机发电有功功率，为发电机发电无功功率，负荷有功功率，负荷无功功率，
其中：Vi，Vj为节点i,j的电压，Gij为节点导纳矩阵的电导，Bij为节点导纳矩阵的电纳；
(4)多能源耦合环节建模
将城市能源互联网耦合环节概括为一个包括输入-输出二端口的网络，并利用一个耦合矩阵来描述输入能源同输出能源之间的关系，称为能源集线器模型，在能源集线器内部，不同耦合模式的设备根据需要任意组合，输入的能源和输出的能源形式也根据需要进行指定，
能源集线器的输入为P输出为L，整个能源集线器表示为：
L＝f(P) (15)
式中，cij为耦合因子，表示第j种形式能源输出与第i种形式能源输入的比值；
(5)城市能互联网统一潮流模型
城市能源互联网的统一潮流计算模型为：
式中：PSP、QSP、ΦSP和LSP为初始状态下系统给定的参数，分别代表有功功率、无功功率、热功率和天然气负荷；第1、2行代表电力系统的模型，含义为系统的有功和无功差值，第3、4、5、6行代表热力系统的模型，含义为系统的节点热功率、供热力管网络回路压力降、供热温度及回热温度差值，第7行代表配气网模型，含义是天然气系统节点负荷差值；Asl、Agl分别为供热力管网络和天然气网络参数，根据网络拓扑形成的矩阵；Cs、Cr为根据热力管网、回热管网求得的矩阵，bs、br为根据供热温度和输出温度计算的矩阵，x＝[θ，U，m，Ts.load,，Tr.load，Π]T为城市能源互联网的状态量；
将牛顿拉夫逊法运用到城市能源互联网联合潮流计算中，其迭代形式如下：
城市能源互联网潮流方程雅克比矩阵J表示为：
式中：ΔFe＝[ΔP，ΔQ]T表示和配电系统有关的差值；ΔFh＝[(ΔΦ，Δp)，(ΔTs，ΔTr)]T表示和供热系统有关的差值；ΔFg＝Δf表示和配气系统有关的差值；xe＝[θ，U]T表示和配电系统有关的状态变量；xh＝[m，(Ts.load，Tr.load)]T表示和供热系统有关的状态变量；xg＝Π表示和配气系统有关的状态变量；Jee表示配电系统自身潮流同自身状态变量之间的关系；Jhh供热系统自身潮流和自身状态变量之间的关系；Jgg表示配气系统自身潮流与自身状态量之间的关系,
在配气系统中，平衡节点设置在气源节点，因此配气网络状态的变化所产生的供需波动将会由气源节点的波动变化来承担，不会传导到热力系统和电力系统因此Jeg、Jhg均为0；对于配电系统，配电系统内部的网路状态或者负荷情况发生变化时，由于配电系统根节点为变压器，因此其网络内部的波动会由变压器的功率变化来承担，波动所产生的影响不会顺延到供热系统和配气系统，因此Jge、Jhe为0；对于供热系统来说，由于其平衡节点实际上是配电系统和配气系统的负荷，具体依赖于能源集线器内部的设备形式，所以供热系统网络状态的改变以及负荷情况的变化，会导致供热系统平衡节点出力的变化而影响到配电系统和配气系统，因此Jeh、Jgh不为0，Jeh、Jgh的表达式依赖于能源集线器的形态，需要按照具体的情况去分析；
(6)城市能源互联网优化潮流
所述能源集线器由热电联供机组和燃气锅炉组成，通过调整流向燃气锅炉和流向热电联产机组的天然气比例，实现能源集线器热电比的调节，能源集线器的耦合矩阵系数为，
式中，v代表天然气的分配比例；代表燃气锅炉的转化效率；cm代表热电联产机组的热电比；代表热电联产机组中天然气向电能的转化效率；代表热电联产机组中天然气向热能的转化效率；
城市能源互联网的目标函数为综合能源运营商的运营成本最低，运营成本主要包括两块，分别是购电成本和购气成本，设综合能源运营商从大电网购买的电量为Pgrid，从城镇天然气门站购买的天然气总量为Pgas,，定义目标函数为：
min f＝cePgrid+cgPgas (29)
其中，ce表示的是电价；cg表示的是天然气价格，决策变量为能源集线器的分配比v，
①等式约束
对于城市能源互联网，等式约束由城市能源互联网潮流计算模型得知：
②不等式约束
约束条件由配电网，配气网，供热系统以及能源集线器四部分的约束组成，对于配电网而言，节点电压必须在某一范围内；
Umin＜Ui＜Umax (31)
对于配气系统而言，配气系统的节点压力在一个范围内；
Πmin＜Πi＜Πmax (32)
对于供热系统，管道流量在一个范围内；
fheat.min＜fheat.i＜fheat.max (33)
考虑到能源集线器中不同设备的运行约束，能源集线器的分配比存在一个上下限，表示为：
0＜v＜1 (34)
城市能源互联网潮流优化模型由式(27)式-(34)式组成。
而且，所述步骤(4)中能源集线器模型具体由热电联供机组和燃气锅炉组成，能源集线器模型的热电比通过调整流向燃气锅炉和流向热电联产机组的天然气比例实现。
本发明的优点和积极效果是：
1、本发明提出一种基于热电比调节的城市能源互联网统一优化潮流计算方法，以能源互联网综合能源运营商的运营成本最低为目标，以热电比作为调节手段，建立了考虑配电网，配气网，供热系统以及能源集线器四部分约束的优化潮流模型并进行求解，可有效提高城市能源互联网运行的经济性；
2、本发明通过仿真表明，本发明应用简单、快速，具有很好的实用价值。
附图说明
图1是本发明中能源集线器原理示意图；
图2是本发明中一种典型能源集线器；
图3是本发明中城市能源互联网系统结构图；
图4是本发明中配电系统潮流计算结果示意图；
图5是本发明中供热系统各支路稳态流量示意图；
图6是本发明中供热系统各节点稳态温度示意图；
图7是本发明中天然气系统稳态流量示意图；
图8是本发明中天然气系统各节点天然气压力示意图；
图9是本发明中能源集线器输出功率示意图；
图10是本发明中配电系统潮流计算结果示意图；
图11是本发明中供热系统各支路稳态流量示意图；
图12是本发明中供热系统各节点稳态温度示意图；
图13是本发明中天然气系统稳态流量示意图；
图14是本发明中天然气系统各节点天然气压力示意图；
图15是本发明中能源集线器计算结果A示意图；
图16是本发明中能源集线器计算结果B示意图。
具体实施方式
以下对本发明实施例做进一步详述：需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其它实施方式，同样属于本发明保护的范围。
一种基于热电比调节的城市能源互联网潮流计算方法，包括步骤如下：
(1)配气系统数学模型
①稳态流动方程
配气系统模型是管道的稳态流动方程，对于k节点到m节点的配气网稳态流量fkm，可以采用下式描述：
其中：pn为标况下的天然气压力基准；Ta为管道中天然气的平均温度；Dkm为管道直径；L为管道的长度；Tn为温度基准；pk为节点k的天然气压力；pm为节点m的天然气压力；γg，Za，F为和天然气有关的系数，可以查表获取；pk＞pm，sg(pk,pm)＝1，若pk＜pm，sg(pk,pm)＝-1；
根据式子(1)，将配气系统稳态流量f改写为与配气网天然气气压压力差ΔΠ有关的函数，表示为：
f＝φ(ΔΠ) (2)
其中：
②连续方程
天然气流量在某一节点处代数和为零，表示如下：
A1f＝L (4)
其中：A1为天然气管网的支路-节点矩阵，L为负荷；
同时规定符号如下：对于支路而言，如果流量顺着支路流入某一节点，则f为正，如果流量顺着支路流出某一节点，则f为负，对于负荷节点而言，L为正，对气源节点而言，L为负；
根据式(3)，将配气网气压平方π改写为与配气网天然气气压压力差ΔΠ有关的函数，可以表示为：
ΔΠ＝-ATΠ (5)
其中A为天然气管网的支路-节点矩阵
③回路方程
根据基尔霍夫第二定律，对于一个回路而言，回路中的管道压力降为零，表示为：
BΔΠ＝0 (6)
其中：B为天然气网络中的回路-支路矩阵
(2)供热系统模型
①水动力模型
供热系统本质上是流体场，需要满足流体运动规律，主要是节点流量连续和回路压损规律，具体地说，就是对于处于网络中的热水，对于某一节点，流入的热水质量等于流出的热水质量；对于某一回路，热水在网络中运动，由于摩擦的原因，水头压力会下降，因此对于一个热水回路而言，回路的压力降是为0，可以表示为：
Asm＝mq (7)
Bhhf＝0 (8)
式中：As为热力管网的支路-节点矩阵；m为管道流量；mq为热力负荷的负荷需求；Bh为热力管网的支路-回路矩阵；hf为压力降，表示为
hf＝Km|m| (9)
式中：K为管道的阻力系数参量，由管道的材质决定，
②热动力模型
热动力模型描述的是供热系统的热平衡行为，对于热用户，用户的热负荷需求取决于供热的供热温度Ts，供热的输出温度To，以及单位时间内流入该热用户的热水质量mq，热力模型为：
Φ＝Cpmq(Ts-To) (10)
热水在网络的流动过程中由于各类损耗，温度会有一定程度的下降，称之为温降，假设管道的始端温度为Tstart，末端温度为Tend，环境温度Ta，管道与温降有关的系数为λ，管道的长度为L，那么温降方程可以表示为：
热水在网络中发生混合须满足能量守恒定律，因此混合方程表述为：
(∑mout)Tout＝∑minTin (12)
(3)配电系统模型
配电系统的模型为反映节点功率与节点电压之间关系的节点功率方程：
其中：为发电机发电有功功率，为发电机发电无功功率，负荷有功功率，负荷无功功率
其中：Vi，Vj为节点i,j的电压，Gij为节点导纳矩阵的电导，Bij为节点导纳矩阵的电纳；
(4)多能源耦合环节建模
考虑到城市能源互联网内电气冷热之间耦合特性与交互机理相当复杂，一方面耦合设备繁多，能源形式多样，因此通用的建模显得尤为重要，对一个能源互联网而言，无论其中电、热、气之间的耦合关系是多么的复杂，都可以概括为需要各种形式能源的输入，最终有其他形式的能源输出，于是可将耦合环节概括为如图1所示的输入-输出二端口网络，并利用一个耦合矩阵来描述输入能源同输出能源之间的关系，可称为能源集线器模型，在能源集线器内部，不同耦合模式的设备可以根据需要任意组合，输入的能源和输出的能源形式也可以根据需要进行指定，实际上能源集线器成了一个灵活的接口。
能源集线器的输入为P输出为L，整个能源集线器表示为：
L＝f(P) (15)
式中，cij为耦合因子，表示第j种形式能源输出与第i种形式能源输入的比值；
(5)城市能互联网统一潮流模型
城市能源互联网的统一潮流计算模型为：
式中：PSP、QSP、ΦSP和LSP为初始状态下系统给定的参数，分别代表有功功率、无功功率、热功率和天然气负荷；第1、2行代表电力系统的模型，含义为系统的有功和无功差值，第3、4、5、6行代表热力系统的模型，含义为系统的节点热功率、供热力管网络回路压力降、供热温度及回热温度差值，第7行代表配气网模型，含义是天然气系统节点负荷差值；Asl、Agl分别为供热力管网络和天然气网络参数，根据网络拓扑形成的矩阵；Cs、Cr为根据热力管网、回热管网求得的矩阵，bs、br为根据供热温度和输出温度计算的矩阵，x＝[θ，U，m，Ts.load,，Tr.load，Π]T为城市能源互联网的状态量；
本发明将牛顿拉夫逊法运用到城市能源互联网联合潮流计算中，其迭代形式如下：
城市能源互联网潮流方程雅克比矩阵J可表示为：
式中：ΔFe＝[ΔP，ΔQ]T表示和配电系统有关的差值；ΔFh＝[(ΔΦ，Δp)，(ΔTs，ΔTr)]T表示和供热系统有关的差值；ΔFg＝Δf表示和配气系统有关的差值；xe＝[θ，U]T表示和配电系统有关的状态变量；xh＝[m，(Ts.load，Tr.load)]T表示和供热系统有关的状态变量；xg＝Π表示和配气系统有关的状态变量；Jee表示配电系统自身潮流同自身状态变量之间的关系；Jhh供热系统自身潮流和自身状态变量之间的关系；Jgg表示配气系统自身潮流与自身状态量之间的关系,
在配气系统中，平衡节点设置在气源节点，因此天然气网络状态的变化所产生的供需波动将会由气源节点的波动变化来承担，不会传导到热力系统和电力系统因此Jeg、Jhg均为0；对于配电系统，配电系统内部的网路状态或者负荷情况发生变化时，由于配电系统根节点为变压器，因此其网络内部的波动会由变压器(平衡节点)的功率变化来承担，波动所产生的影响不会顺延到供热系统和配气系统，因此Jge、Jhe为0；对于供热系统来说，由于其平衡节点实际上是配电系统和配气系统的负荷，具体依赖于能源集线器内部的设备形式，所以供热系统网络状态的改变以及负荷情况的变化，会导致供热系统平衡节点(供热源)出力的变化而影响到配电系统和配气系统，因此Jeh、Jgh不为0，Jeh、Jgh的表达式依赖于能源集线器的形态，需要按照具体的情况去分析；
(6)城市能源互联网优化潮流模型
以图2所示的典型能源集线器为例，介绍其调节热电比的过程，图中的能源集线器改为由热电联供机组和燃气锅炉组成，调节热电比的方法有两种，第一种方法是通过向热电联产机组通入补燃燃油来改变热电联产机组的热电比进而影响能源集线器的热电比，另外一种方法是直接通过影响天然气的分配系数改变热电联产机组和燃气锅炉的出力而改变能源集线器的热电比，本文采用后一种方法，即通过调整流向燃气锅炉和流向热电联产机组的天然气比例，实现能源集线器热电比的调节，从建模的角度出发，相当于改变了能源集线器的耦合矩阵系数，
式中，v代表天然气的分配比例；代表燃气锅炉的转化效率；cm代表热电联产机组的热电比；代表热电联产机组中天然气向电能的转化效率；代表热电联产机组中天然气向热能的转化效率；
城市能源互联网的目标函数为综合能源运营商的运营成本最低，运营成本主要包括两块，分别是购电成本和购气成本，设综合能源运营商从大电网购买的电量为Pgrid，从城镇天然气门站购买的天然气总量为Pgas,，定义目标函数为：
minf＝cePgrid+cgPgas (29)
其中，ce表示的是电价；cg表示的是天然气价格，决策变量为能源集线器的分配比v，
①等式约束
对于城市能源互联网，等式约束由城市能源互联网潮流计算模型得知：
②不等式约束
约束条件由配电网，配气网，供热系统以及能源集线器四部分的约束组成，对于配电网而言，节点电压必须在某一范围内；
Umin＜Ui＜Umax (31)
对于配气系统而言，配气系统的节点压力在一个范围内；
Πmin＜Πi＜Πmax (32)
对于供热系统，管道流量在一个范围内；
fheat.min＜fheat.i＜fheat.max (33)
考虑到能源集线器中不同设备的运行约束，能源集线器的分配比存在一个上下限，表示为：
0＜v＜1 (34)
城市能源互联网潮流优化模型由式(27)式-(34)式组成。
具体实例的原理分析
如图3所示的城市能源互联网算例中配电系统为一个有13个节点的系统，EB13是配电系统的根节点，因此设节点EB13为配电系统的平衡节点，节点EB1连接了分布式风力发电机，节点EB12与节点GB4、节点EB11与节点GB11同为配电系统与配气系统通过能源集线器连接的节点，因此节点EB11，节点EB12可认为是PV节点，余下的配电系统节点视为PQ节点。两个能源集线器，如图3所示，对于EH1，可以认为处于“以电定热”模式，输入为配电系统的电力和天然气系统的天然气，来自配电系统的电力进入电锅炉中，与来自热电联产机组产生的热量，汇合形成供热力管网络的供热源之一。EH1中发出的电量，优先考虑满足电负荷，根据电负荷确定发出的热量。对于EH2，可以认为处于“以热定电”模式，输入为配天然气系统的天然气，EH2中设备有热电联供机组以及燃气锅炉，作为供热系统的平衡节点。EH2中发出的热量，优先考虑满足热负荷，根据热负荷确定发出的电量。系统整体的结构如图3所示。
潮流计算结果如图4-9图所示：
在进行优化潮流计算时，配电网电压上下限为[0.95,1.06]，天然气节点气压上下限为[0，75]mbar，供热管道流量上下限为[0,10]kg/s，分配比上下限为[0,1]，热电联供机组的热电比为1。气价为2元/方，上游电网购电电价为1元/度，系统结构如图10所示。
通过计算不同分配比(热电比)情况下城市综合能源系统的运营成本如下表1所示：
表1不同分配比情况下系统的运营成本
当分配比v为0.6时，系统稳态分析结果如图11-图14所示。
可见，当分配比为0.6时，整个城市综合能源系统的运营成本最低。如图15-图16所示，对于能源集线器1来说，优化后电功率不变，热功率减少，而能源集线器2的电功率和热功率相应增加，相应的能源集线器1的天然气耗量不变，能源集线器2的天然气耗量上升。能源集线器2所接入的位置，相当于配气网的上游部分，因此整个网络的燃气损耗会相应降低。结果表明，系统的优化潮流你结果严重依赖于天然气价格和电力价格。算例是针对某一城市的，因此天然气的价格是固定的，实际上在我国，天然气的价格是因地不同的。某些地方，天然气的价格能到达1-2块，如西宁，呼和浩特等城市，同时，某些城市的天然气价格也居高不下，比如北京，上海，广州，天津等城市。因此在不同的天然气价格下，运算的结果会有较大的变化，因此热电比的调整，依赖于天然气价格以及当地的热电负荷比例。

Claims (2)

1.一种基于热电比调节的城市能源互联网潮流计算方法，其特征在于包括步骤如下：
(1)配气系统数学模型
①稳态流动方程
配气系统模型是管道的稳态流动方程，对于k节点到m节点的配气网稳态流量fkm，采用下式描述：
<mrow> <msub> <mi>f</mi> <mrow> <mi>k</mi> <mi>m</mi> </mrow> </msub> <mo>=</mo> <mi>s</mi> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>p</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>7.6</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </msup> <mo>&amp;CenterDot;</mo> <msub> <mi>T</mi> <mi>n</mi> </msub> </mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> </mfrac> <mo>)</mo> </mrow> <msubsup> <mi>D</mi> <mrow> <mi>k</mi> <mi>m</mi> </mrow> <mrow> <mn>5</mn> <mo>/</mo> <mn>2</mn> </mrow> </msubsup> <mo>&amp;times;</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msubsup> <mi>p</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>p</mi> <mi>m</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <msup> <mi>L</mi> <mrow> <mi>k</mi> <mi>m</mi> </mrow> </msup> <msub> <mi>&amp;gamma;</mi> <mi>g</mi> </msub> <msub> <mi>T</mi> <mi>a</mi> </msub> <msub> <mi>Z</mi> <mi>a</mi> </msub> <mi>F</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中：pn为标况下的天然气压力基准；Ta为管道中天然气的平均温度；Dkm为管道直径；L为管道的长度；Tn为温度基准；pk为节点k的天然气压力；pm为节点m的天然气压力；γg，Za，F为和天然气有关的系数，可以查表获取；pk＞pm，sg(pk,pm)＝1，若pk＜pm，sg(pk,pm)＝-1；
根据式子(1)，将配气系统稳态流量f改写为与配气网天然气气压压力差ΔΠ有关的函数，表示为：
f＝φ(ΔΠ) (2)
其中：
<mrow> <msub> <mi>&amp;Delta;&amp;Pi;</mi> <mrow> <mi>k</mi> <mi>m</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>p</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>p</mi> <mi>m</mi> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
②连续方程
天然气流量在某一节点处代数和为零，表示如下：
A1f＝L (4)
其中：A1为天然气管网的支路-节点矩阵，L为负荷；
同时规定符号如下：对于支路而言，如果流量顺着支路流入某一节点，则f为正，如果流量顺着支路流出某一节点，则f为负，对于负荷节点而言，L为正，对气源节点而言，L为负；
根据式(3)，将配气网气压平方改写为与配气网天然气气压压力差ΔΠ有关的函数，表示为：
ΔΠ＝-ATΠ (5)
其中A为天然气管网的支路-节点矩阵，
③回路方程
根据基尔霍夫第二定律，对于一个回路而言，回路中的管道压力降为零，表示为：
BΔΠ＝0 (6)
其中：B为天然气网络中的回路-支路矩阵，
(2)供热系统模型
①水动力模型
供热系统中对于某一节点，流入的热水质量等于流出的热水质量；对于某一回路，热水在网络中运动，由于摩擦的原因，水头压力会下降，因此对于一个热水回路而言，回路的压力降为0，表示为：
Asm＝mq (7)
Bhhf＝0 (8)
式中：As为热力管网的支路-节点矩阵；m为管道流量；mq为热力负荷的负荷需求；Bh为热力管网的支路-回路矩阵；hf为压力降，表示为
hf＝Km|m| (9)
式中：K为管道的阻力系数参量，由管道的材质决定，
②热动力模型
热动力模型描述的是供热系统的热平衡行为，用户的热负荷需求取决于供热的供热温度Ts，供热的输出温度To，以及单位时间内流入该热用户的热水质量mq，热力模型为：
Φ＝Cpmq(Ts-To) (10)
热水在网络的流动过程中，假设管道的始端温度为Tstart，末端温度为Tend，环境温度Ta，管道与温降有关的系数为λ，管道的长度为L，那么温降方程表示为：
<mrow> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>s</mi> <mi>t</mi> <mi>a</mi> <mi>r</mi> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>a</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>&amp;lambda;</mi> <mi>L</mi> </mrow> <mrow> <msub> <mi>C</mi> <mi>p</mi> </msub> <mi>m</mi> </mrow> </mfrac> </mrow> </msup> <mo>+</mo> <msub> <mi>T</mi> <mi>a</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
热水在网络中发生混合须满足能量守恒定律，因此混合方程表述为：
(∑mout)Tout＝∑minTin (12)
(3)配电系统模型
配电系统的模型为反映节点功率与节点电压之间关系的节点功率方程：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mi>P</mi> <mi>g</mi> <mi>i</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mi>l</mi> <mi>i</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mi>g</mi> <mi>i</mi> </msubsup> <mo>-</mo> <msubsup> <mi>Q</mi> <mi>l</mi> <mi>i</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>e</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
其中：为发电机发电有功功率，为发电机发电无功功率，Pl i负荷有功功率，负荷无功功率，
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>P</mi> <mi>i</mi> </msup> <mo>=</mo> <msub> <mi>V</mi> <mi>i</mi> </msub> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <msub> <mi>V</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mi>cos</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mi>cos</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>,</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>Q</mi> <mi>i</mi> </msup> <mo>=</mo> <msub> <mi>V</mi> <mi>i</mi> </msub> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <msub> <mi>V</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mi>sin</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mi>cos</mi> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>,</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
其中：Vi，Vj为节点i,j的电压，Gij为节点导纳矩阵的电导，Bij为节点导纳矩阵的电纳；
(4)多能源耦合环节建模
将城市能源互联网耦合环节概括为一个包括输入-输出二端口的网络，并利用一个耦合矩阵来描述输入能源同输出能源之间的关系，称为能源集线器模型，在能源集线器内部，不同耦合模式的设备根据需要任意组合，输入的能源和输出的能源形式也根据需要进行指定，
能源集线器的输入为P输出为L，整个能源集线器表示为：
L＝f(P) (15)
式中，cij为耦合因子，表示第j种形式能源输出与第i种形式能源输入的比值；
(5)城市能互联网统一潮流模型
城市能源互联网的统一潮流计算模型为：
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>P</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>Q</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;Phi;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>p</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;T</mi> <mi>s</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;T</mi> <mi>r</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>f</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mi>P</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>-</mo> <mi>Re</mi> <mo>{</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>Y</mi> <mi>U</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>Q</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>-</mo> <mi>Im</mi> <mo>{</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>Y</mi> <mi>U</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>p</mi> </msub> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mi>m</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>o</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>B</mi> <mi>h</mi> </msub> <mi>K</mi> <mi>m</mi> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>s</mi> </msub> <msub> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>s</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>CrT</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>r</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>A</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <mi>&amp;phi;</mi> <mrow> <mo>(</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>g</mi> <mi>T</mi> </msubsup> <mi>&amp;Pi;</mi> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>L</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
式中：PSP、QSP、ΦSP和LSP为初始状态下系统给定的参数，分别代表有功功率、无功功率、热功率和天然气负荷；第1、2行代表电力系统的模型，含义为系统的有功和无功差值，第3、4、5、6行代表热力系统的模型，含义为系统的节点热功率、供热力管网络回路压力降、供热温度及回热温度差值，第7行代表配气网模型，含义是天然气系统节点负荷差值；Asl、Agl分别为供热力管网络和天然气网络参数，根据网络拓扑形成的矩阵；Cs、Cr为根据热力管网、回热管网求得的矩阵，bs、br为根据供热温度和输出温度计算的矩阵，x＝[θ，U，m，Ts.load,，Tr.load，Π]T为城市能源互联网的状态量；
将牛顿拉夫逊法运用到城市能源互联网联合潮流计算中，其迭代形式如下：
<mrow> <mo>{</mo> <mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>&amp;Delta;x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>J</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>&amp;Delta;F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msup> <mi>x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <mi>&amp;Delta;x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> </mrow> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
城市能源互联网潮流方程雅克比矩阵J表示为：
<mrow> <mi>J</mi> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>J</mi> <mrow> <mi>e</mi> <mi>e</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>e</mi> <mi>h</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>e</mi> <mi>g</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mi>e</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mi>h</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mi>g</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mrow> <mi>g</mi> <mi>e</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>g</mi> <mi>h</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>g</mi> <mi>g</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>e</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>e</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>e</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>h</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>e</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>g</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>e</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>h</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>g</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>e</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>h</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mi>g</mi> </msub> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
式中：ΔFe＝[ΔP，ΔQ]T表示和配电系统有关的差值；ΔFh＝[(ΔΦ，Δp)，(ΔTs，ΔTr)]T表示和供热系统有关的差值；ΔFg＝Δf表示和配气系统有关的差值；xe＝[θ，U]T表示和配电系统有关的状态变量；xh＝[m，(Ts.load，Tr.load)]T表示和供热系统有关的状态变量；xg＝Π表示和配气系统有关的状态变量；Jee表示配电系统自身潮流同自身状态变量之间的关系；Jhh供热系统自身潮流和自身状态变量之间的关系；Jgg表示配气系统自身潮流与自身状态量之间的关系,
<mrow> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mi>h</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>&amp;Phi;</mi> <mo>,</mo> <mi>p</mi> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>&amp;Phi;</mi> <mo>,</mo> <mi>p</mi> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>T</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>,</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>,</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>T</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>11</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>12</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>21</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>22</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>11</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <mi>&amp;Phi;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;Delta;</mi> <mi>p</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>C</mi> <mi>p</mi> </msub> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>{</mo> <mo>(</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>)</mo> <mo>}</mo> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mi>B</mi> <mi>K</mi> <mo>|</mo> <mi>m</mi> <mo>|</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>12</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>p</mi> </msub> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>{</mo> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mi>m</mi> <mo>}</mo> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>21</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mo>-</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>b</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>b</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>m</mi> <mi>T</mi> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>J</mi> <mrow> <mi>h</mi> <mn>22</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> <mi>T</mi> </msubsup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>T</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> <mi>T</mi> </msubsup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> <mi>T</mi> </msubsup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;Delta;T</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>T</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> <mi>T</mi> </msubsup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>C</mi> <mi>s</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>C</mi> <mi>r</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>J</mi> <mrow> <mi>g</mi> <mi>g</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <msubsup> <mi>DA</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>d</mi> <mrow> <mi>i</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>f</mi> <mi>i</mi> </msub> <mrow> <mn>2</mn> <msub> <mi>&amp;Delta;&amp;Pi;</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>
在配气系统中，平衡节点设置在气源节点，因此配气网络状态的变化所产生的供需波动将会由气源节点的波动变化来承担，不会传导到热力系统和电力系统因此Jeg、Jhg均为0；对于配电系统，配电系统内部的网路状态或者负荷情况发生变化时，由于配电系统根节点为变压器，因此其网络内部的波动会由变压器的功率变化来承担，波动所产生的影响不会顺延到供热系统和配气系统，因此Jge、Jhe为0；对于供热系统来说，由于其平衡节点实际上是配电系统和配气系统的负荷，具体依赖于能源集线器内部的设备形式，所以供热系统网络状态的改变以及负荷情况的变化，会导致供热系统平衡节点出力的变化而影响到配电系统和配气系统，因此Jeh、Jgh不为0，Jeh、Jgh的表达式依赖于能源集线器的形态，需要按照具体的情况去分析；
(6)城市能源互联网优化潮流
所述能源集线器由热电联供机组和燃气锅炉组成，通过调整流向燃气锅炉和流向热电联产机组的天然气比例，实现能源集线器热电比的调节，能源集线器的耦合矩阵系数为，
<mrow> <msub> <mi>L</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>+</mo> <mi>c</mi> <mi>m</mi> </mrow> </mfrac> <msubsup> <mi>vd</mi> <mrow> <mi>g</mi> <mi>e</mi> </mrow> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> </mrow> </msubsup> <mi>P</mi> <mi>g</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>L</mi> <mi>h</mi> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>v</mi> <mo>)</mo> </mrow> <msubsup> <mi>d</mi> <mrow> <mi>g</mi> <mi>h</mi> </mrow> <mrow> <mi>G</mi> <mi>B</mi> </mrow> </msubsup> <mo>+</mo> <mfrac> <mrow> <mi>c</mi> <mi>m</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>c</mi> <mi>m</mi> </mrow> </mfrac> <msubsup> <mi>vd</mi> <mrow> <mi>g</mi> <mi>h</mi> </mrow> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> </mrow> </msubsup> <mo>&amp;rsqb;</mo> <mi>P</mi> <mi>g</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>
式中，v代表天然气的分配比例；代表燃气锅炉的转化效率；cm代表热电联产机组的热电比；代表热电联产机组中天然气向电能的转化效率；代表热电联产机组中天然气向热能的转化效率；
城市能源互联网的目标函数为综合能源运营商的运营成本最低，运营成本主要包括两块，分别是购电成本和购气成本，设综合能源运营商从大电网购买的电量为Pgrid，从城镇天然气门站购买的天然气总量为Pgas,，定义目标函数为：
minf＝cePgrid+cgPgas (29)
其中，ce表示的是电价；cg表示的是天然气价格，决策变量为能源集线器的分配比v，
①等式约束
对于城市能源互联网，等式约束由城市能源互联网潮流计算模型得知：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mi>P</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>-</mo> <mi>Re</mi> <mo>{</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>Y</mi> <mi>U</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>Q</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>-</mo> <mi>Im</mi> <mo>{</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>Y</mi> <mi>U</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>p</mi> </msub> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mi>m</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>o</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>B</mi> <mi>h</mi> </msub> <mi>K</mi> <mi>m</mi> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>s</mi> </msub> <msub> <mi>T</mi> <mrow> <mi>s</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>s</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>CrT</mi> <mrow> <mi>r</mi> <mo>.</mo> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>r</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>A</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <mi>&amp;phi;</mi> <mrow> <mo>(</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>g</mi> <mi>T</mi> </msubsup> <mi>&amp;Pi;</mi> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>L</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msup> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>
②不等式约束
约束条件由配电网，配气网，供热系统以及能源集线器四部分的约束组成，对于配电网而言，节点电压必须在某一范围内；
Umin＜Ui＜Umax (31)
对于配气系统而言，配气系统的节点压力在一个范围内；
Πmin＜Πi＜Πmax (32)
对于供热系统，管道流量在一个范围内；
fheat.min＜fheat.i＜fheat.max (33)
考虑到能源集线器中不同设备的运行约束，能源集线器的分配比存在一个上下限，表示为：
0＜v＜1 (34)
城市能源互联网潮流优化模型由式(27)式-(34)式组成。
2.根据权利要求1所述的基于热电比调节的城市能源互联网潮流计算方法，其特征在于：所述步骤(4)中能源集线器模型具体由热电联供机组和燃气锅炉组成，能源集线器模型的热电比通过调整流向燃气锅炉和流向热电联产机组的天然气比例实现。
CN201711005289.6A 2017-10-25 2017-10-25 基于热电比调节的城市能源互联网潮流计算方法 Pending CN107808218A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711005289.6A CN107808218A (zh) 2017-10-25 2017-10-25 基于热电比调节的城市能源互联网潮流计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711005289.6A CN107808218A (zh) 2017-10-25 2017-10-25 基于热电比调节的城市能源互联网潮流计算方法

Publications (1)

Publication Number Publication Date
CN107808218A true CN107808218A (zh) 2018-03-16

Family

ID=61591736

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711005289.6A Pending CN107808218A (zh) 2017-10-25 2017-10-25 基于热电比调节的城市能源互联网潮流计算方法

Country Status (1)

Country Link
CN (1) CN107808218A (zh)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596453A (zh) * 2018-04-10 2018-09-28 山东大学 考虑网络动态特性的综合能源系统日前优化调度方法及系统
CN108921727A (zh) * 2018-06-30 2018-11-30 天津大学 考虑热负荷动态特性的区域综合能源系统可靠性评估方法
CN109145441A (zh) * 2018-08-21 2019-01-04 清华-伯克利深圳学院筹备办公室 一种热电供能的计算方法、装置、电子设备及存储介质
CN109344480A (zh) * 2018-09-20 2019-02-15 国网四川省电力公司经济技术研究院 考虑蒸汽网络的电热联合系统潮流计算方法
CN109347107A (zh) * 2018-09-29 2019-02-15 河海大学 一种并网运行电热互联综合能源系统潮流计算方法
CN109492325A (zh) * 2018-11-29 2019-03-19 东北大学 一种基于扩展能源集线器的多能耦合系统潮流分析方法
CN109636148A (zh) * 2018-11-29 2019-04-16 华南理工大学 基于能量网络方程的多能流系统的工作状态评估方法
CN109726483A (zh) * 2018-12-29 2019-05-07 国网江苏省电力有限公司南京供电分公司 一种电热互联综合能源系统辐射状热网模型及其系统
CN109783986A (zh) * 2019-02-26 2019-05-21 上海交通大学 基于功率守恒的孤岛型综合能源系统联合潮流获取方法
CN109888790A (zh) * 2019-03-28 2019-06-14 国网福建省电力有限公司经济技术研究院 一种不同运行模式下的区域综合能源系统多能潮流计算方法
CN109901637A (zh) * 2019-04-04 2019-06-18 阳江职业技术学院 一种热电植三联产的温室控制方法
CN110263966A (zh) * 2019-05-06 2019-09-20 天津大学 考虑动态传热过程的电-热综合能源系统优化调度方法
CN110348602A (zh) * 2019-06-06 2019-10-18 国网浙江省电力有限公司经济技术研究院 计及天然气管网和热力管网特性的综合能源系统优化方法
CN110502791A (zh) * 2019-07-22 2019-11-26 清华大学 基于能源集线器的综合能源系统稳态建模方法
CN110689206A (zh) * 2019-10-09 2020-01-14 国电南瑞南京控制系统有限公司 一种能源互联网多能主体能源购买及转换运营方法
CN111030120A (zh) * 2019-12-05 2020-04-17 国网辽宁省电力有限公司大连供电公司 城市综合能源网统一平台潮流计算方法
CN111061986A (zh) * 2019-11-25 2020-04-24 广东电网有限责任公司 一种多运行模式的热电综合能源系统潮流计算方法
CN111414721A (zh) * 2020-02-22 2020-07-14 清华大学 一种用于综合能源系统调度的供热网络水路建模方法
CN111711206A (zh) * 2020-06-30 2020-09-25 天津大学 一种计及热网动态特性的城市热电综合能源系统调度方法
CN111815042A (zh) * 2020-06-30 2020-10-23 天津大学 计及精细化热网模型的电热综合能源系统优化调度方法
CN113487342A (zh) * 2021-06-03 2021-10-08 华中科技大学 一种热力系统作为利益整体参与电力市场的决策方法
CN113837577A (zh) * 2021-09-14 2021-12-24 国网河南省电力公司经济技术研究院 一种农村电热联合系统耦合元件规划方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106532709A (zh) * 2016-10-25 2017-03-22 天津大学 一种含配电网重构的区域综合能源系统最优潮流计算方法
CN106786603A (zh) * 2017-02-22 2017-05-31 天津大学 一种区域综合能源系统多目标最优混合潮流算法
CN106845671A (zh) * 2016-12-12 2017-06-13 长沙理工大学 一种多能流系统多目标最优潮流模型及其求解方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106532709A (zh) * 2016-10-25 2017-03-22 天津大学 一种含配电网重构的区域综合能源系统最优潮流计算方法
CN106845671A (zh) * 2016-12-12 2017-06-13 长沙理工大学 一种多能流系统多目标最优潮流模型及其求解方法
CN106786603A (zh) * 2017-02-22 2017-05-31 天津大学 一种区域综合能源系统多目标最优混合潮流算法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘聪 等: "电/热/气综合能源系统混合潮流计算方法", 《山东工业技术》 *
张以斌: "天然气-电力混合系统分析方法研究", 《中国优秀博士学位论文全文数据库 工程科技II辑》 *
张国衡: "电力系统动态潮流计算及网络拓扑分析", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
徐宪东 等: "区域综合能源系统电/气/热混合潮流算法研究", 《中国电机工程学》 *

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596453B (zh) * 2018-04-10 2020-12-11 山东大学 考虑网络动态特性的综合能源系统日前优化调度方法及系统
CN108596453A (zh) * 2018-04-10 2018-09-28 山东大学 考虑网络动态特性的综合能源系统日前优化调度方法及系统
CN108921727A (zh) * 2018-06-30 2018-11-30 天津大学 考虑热负荷动态特性的区域综合能源系统可靠性评估方法
CN108921727B (zh) * 2018-06-30 2021-04-27 天津大学 考虑热负荷动态特性的区域综合能源系统可靠性评估方法
CN109145441B (zh) * 2018-08-21 2023-06-02 清华大学深圳国际研究生院 一种热电供能的计算方法、装置、电子设备及存储介质
CN109145441A (zh) * 2018-08-21 2019-01-04 清华-伯克利深圳学院筹备办公室 一种热电供能的计算方法、装置、电子设备及存储介质
CN109344480A (zh) * 2018-09-20 2019-02-15 国网四川省电力公司经济技术研究院 考虑蒸汽网络的电热联合系统潮流计算方法
CN109347107A (zh) * 2018-09-29 2019-02-15 河海大学 一种并网运行电热互联综合能源系统潮流计算方法
CN109636148B (zh) * 2018-11-29 2022-09-13 华南理工大学 基于能量网络方程的多能流系统的工作状态评估方法
CN109636148A (zh) * 2018-11-29 2019-04-16 华南理工大学 基于能量网络方程的多能流系统的工作状态评估方法
CN109492325A (zh) * 2018-11-29 2019-03-19 东北大学 一种基于扩展能源集线器的多能耦合系统潮流分析方法
CN109726483A (zh) * 2018-12-29 2019-05-07 国网江苏省电力有限公司南京供电分公司 一种电热互联综合能源系统辐射状热网模型及其系统
CN109783986A (zh) * 2019-02-26 2019-05-21 上海交通大学 基于功率守恒的孤岛型综合能源系统联合潮流获取方法
CN109888790A (zh) * 2019-03-28 2019-06-14 国网福建省电力有限公司经济技术研究院 一种不同运行模式下的区域综合能源系统多能潮流计算方法
CN109888790B (zh) * 2019-03-28 2022-07-05 国网福建省电力有限公司经济技术研究院 一种不同运行模式下的区域综合能源系统多能潮流计算方法
CN109901637A (zh) * 2019-04-04 2019-06-18 阳江职业技术学院 一种热电植三联产的温室控制方法
CN109901637B (zh) * 2019-04-04 2021-03-16 阳江职业技术学院 一种热电植三联产的温室控制方法
CN110263966A (zh) * 2019-05-06 2019-09-20 天津大学 考虑动态传热过程的电-热综合能源系统优化调度方法
CN110348602B (zh) * 2019-06-06 2021-09-21 国网浙江省电力有限公司经济技术研究院 计及天然气管网和热力管网特性的综合能源系统优化方法
CN110348602A (zh) * 2019-06-06 2019-10-18 国网浙江省电力有限公司经济技术研究院 计及天然气管网和热力管网特性的综合能源系统优化方法
CN110502791B (zh) * 2019-07-22 2021-05-11 清华大学 基于能源集线器的综合能源系统稳态建模方法
CN110502791A (zh) * 2019-07-22 2019-11-26 清华大学 基于能源集线器的综合能源系统稳态建模方法
CN110689206A (zh) * 2019-10-09 2020-01-14 国电南瑞南京控制系统有限公司 一种能源互联网多能主体能源购买及转换运营方法
CN110689206B (zh) * 2019-10-09 2022-11-01 国电南瑞南京控制系统有限公司 一种能源互联网多能主体能源购买及转换运营方法
CN111061986A (zh) * 2019-11-25 2020-04-24 广东电网有限责任公司 一种多运行模式的热电综合能源系统潮流计算方法
CN111061986B (zh) * 2019-11-25 2023-12-26 广东电网有限责任公司 一种多运行模式的热电综合能源系统潮流计算方法
CN111030120A (zh) * 2019-12-05 2020-04-17 国网辽宁省电力有限公司大连供电公司 城市综合能源网统一平台潮流计算方法
CN111030120B (zh) * 2019-12-05 2023-07-28 国网辽宁省电力有限公司大连供电公司 城市综合能源网统一平台潮流计算方法
CN111414721A (zh) * 2020-02-22 2020-07-14 清华大学 一种用于综合能源系统调度的供热网络水路建模方法
CN111711206B (zh) * 2020-06-30 2023-01-24 天津大学 一种计及热网动态特性的城市热电综合能源系统调度方法
CN111711206A (zh) * 2020-06-30 2020-09-25 天津大学 一种计及热网动态特性的城市热电综合能源系统调度方法
CN111815042A (zh) * 2020-06-30 2020-10-23 天津大学 计及精细化热网模型的电热综合能源系统优化调度方法
CN113487342A (zh) * 2021-06-03 2021-10-08 华中科技大学 一种热力系统作为利益整体参与电力市场的决策方法
CN113837577A (zh) * 2021-09-14 2021-12-24 国网河南省电力公司经济技术研究院 一种农村电热联合系统耦合元件规划方法
CN113837577B (zh) * 2021-09-14 2024-03-19 国网河南省电力公司经济技术研究院 一种农村电热联合系统耦合元件规划方法

Similar Documents

Publication Publication Date Title
CN107808218A (zh) 基于热电比调节的城市能源互联网潮流计算方法
Wang et al. Optimal scheduling strategy of district integrated heat and power system with wind power and multiple energy stations considering thermal inertia of buildings under different heating regulation modes
Wang et al. Operation optimization of regional integrated energy system based on the modeling of electricity-thermal-natural gas network
Jiang et al. Coordinated operation of gas-electricity integrated distribution system with multi-CCHP and distributed renewable energy sources
Zheng et al. Integrated heat and power dispatch truly utilizing thermal inertia of district heating network for wind power integration
CN106532709B (zh) 一种含配电网重构的区域综合能源系统最优潮流计算方法
CN110163443B (zh) 考虑电-气综合需求响应的天然气调压站微能网优化调度方法
Yang et al. Economic-emission dispatch problem in integrated electricity and heat system considering multi-energy demand response and carbon capture Technologies
CN107769215A (zh) 基于能源集线器的园区混合能源系统优化调度方法
CN105048468B (zh) 基于分布式计算的输配电网一体化电压稳定评估方法
Wang et al. Decentralized optimization of coordinated electrical and thermal generations in hierarchical integrated energy systems considering competitive individuals
CN109523076A (zh) 一种储能参与下的综合能源系统日前稳态优化分析方法
Li et al. Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors
CN107947182A (zh) 一种动态潮流分析方法及动态潮流分析系统
CN112035984A (zh) 一种电-气-储区域综合能源系统协同规划方法
Wang et al. Stackelberg game-based optimal scheduling of integrated energy systems considering differences in heat demand across multi-functional areas
Wang et al. Integrated modeling of regional and park-level multi-heterogeneous energy systems
Zhang et al. Bi-layer economic scheduling for integrated energy system based on source-load coordinated carbon reduction
CN113141005B (zh) 一种面向新能源消纳的综合能源系统多时间尺度调度方法
Man et al. State estimation for integrated energy system containing electricity, heat and gas
Liu et al. A new power flow model for combined heat and electricity analysis in an integrated energy system
Zhou et al. Digital Twin Modeling and Simulation of Distributed and Centralized Integrated Energy System on CloudPSS-IESLab
CN113221428B (zh) 电-热综合能源系统动态能流计算的快速分解方法
Luo et al. Optimal operation strategy of combined heat and power system considering demand response and household thermal inertia
CN113487188A (zh) 考虑电气联合价格引导机制的综合能源系统优化调度方法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180316

RJ01 Rejection of invention patent application after publication