CN108155644B - 计及辅助服务的动态供热系统与电力系统耦合调度方法 - Google Patents

计及辅助服务的动态供热系统与电力系统耦合调度方法 Download PDF

Info

Publication number
CN108155644B
CN108155644B CN201711439071.1A CN201711439071A CN108155644B CN 108155644 B CN108155644 B CN 108155644B CN 201711439071 A CN201711439071 A CN 201711439071A CN 108155644 B CN108155644 B CN 108155644B
Authority
CN
China
Prior art keywords
power system
scheduling
heat supply
power
unit
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
CN201711439071.1A
Other languages
English (en)
Other versions
CN108155644A (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.)
Tsinghua University
State Grid Corp of China SGCC
State Grid Jilin Electric Power Corp
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
State Grid Jilin Electric Power Corp
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 Tsinghua University, State Grid Corp of China SGCC, State Grid Jilin Electric Power Corp filed Critical Tsinghua University
Priority to CN201711439071.1A priority Critical patent/CN108155644B/zh
Publication of CN108155644A publication Critical patent/CN108155644A/zh
Application granted granted Critical
Publication of CN108155644B publication Critical patent/CN108155644B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提出一种计及辅助服务的动态供热系统与电力系统耦合调度方法,属于电力系统运行技术领域。本方法首先建立由目标函数和约束条件构成的计及辅助服务的动态供热系统与电力系统,模型的约束条件包括电力系统约束条件和供热系统约束条件;将模型转化为矩阵形式的二次规划模型,通过系统辨识方法,在保证供热系统供热完成情况下,利用供热系统储热性能松弛供热系统供热与用热的强制约关系,使得供热系统能够参与电力系统的调峰,分别得到常规机组、风电机组和热电联产机组的发电出力调度方案和热电联产机组的供热出力调度方案。本发明可充分挖掘供热系统的储热性能,通过电力系统的调峰显著提升风电消纳水平,提高热电联合系统运行灵活性。

Description

计及辅助服务的动态供热系统与电力系统耦合调度方法
技术领域
本发明属于电力系统的运行技术领域,特别涉及一种计及辅助服务的动态供热系统与电力系统耦合调度方法。
背景技术
我国北方城市供热系统具有良好的储热性能,而供热的来源主要来自热电联产机组。
热电联产机组也是我国北方地区火电机组的主要形式,其按照一定的供热与供电关系分别向电网与热网同时供电与供热。为保证热量的有效供应,热电联产机组需要输出大量电力功率,造成北方地区丰富的风电无法被有效利用。
利用供热系统的储热性能,松弛供热与用热间的严格约束,可以促进热电联产机组参与电网调峰,进而促进风电的消纳。
辅助服务是指为维护电力系统的安全稳定运行,保证电能质量,除正常电能生产、输送、使用外,由发电企业、电网经营企业和电力用户提供的服务。包括:一次调频、自动发电控制、调峰、无功调节、备用、黑启动等。
传统电力系统的经济调度中未考虑备用辅助服务成本,调度结果易出现备用设置不合理等情况。另外对于我国北方集中供热期的电力系统,由于热电联产机组要维持基本供热水平,传统调度方法的调度策略弃风严重,该问题急需解决。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种计及辅助服务的动态供热系统和电力系统耦合调度方法。本发明可充分挖掘供热系统的储热性能,通过电力系统的调峰显著提升风电消纳水平。
本发明提出的一种计及辅助服务的动态供热系统与电力系统耦合调度方法,其特征在于,该方法包括以下步骤:
1)建立计及辅助服务的动态供热系统与电力系统耦合调度模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;
计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数为全网运行成本最小化,表达式如下:
Figure BDA0001525209600000021
式中,T为调度时段集合,ICHP、ITU和IWD分别为热电联合系统中的热电联产机组编号集合、常规机组编号集合和风电机组的编号集合,
Figure BDA0001525209600000022
为热电联产机组i在t时段的生产成本函数,
Figure BDA0001525209600000023
为常规机组i在t时段的生产成本函数,
Figure BDA0001525209600000024
为风电机组i在t时段的生产成本函数,
Figure BDA0001525209600000025
为常规机组i在t时段提供的备用辅助服务成本函数,t为调度时段编号;
其中,热电联产机组的生产成本函数为:
Figure BDA0001525209600000026
式中,
Figure BDA0001525209600000027
Figure BDA0001525209600000028
分别为热电联产机组i的生产成本系数,
Figure BDA0001525209600000029
Figure BDA00015252096000000210
分别为热电联产机组i在调度时段t的有功功率和供热功率;
常规机组的生产成本函数为:
Figure BDA00015252096000000211
式中,
Figure BDA00015252096000000212
Figure BDA00015252096000000213
分别为常规机组i的发电成本系数,
Figure BDA00015252096000000214
为常规机组i在调度时段t的有功功率;
风电机组的生产成本函数为:
Figure BDA00015252096000000215
式中,
Figure BDA00015252096000000216
为风电机组i的弃风惩罚因子,
Figure BDA00015252096000000217
为风电机组i在调度时段t的可用有功功率,
Figure BDA00015252096000000218
为风电机组i在调度时段t的实际有功功率;
常规机组提供的备用辅助服务成本函数为:
Figure BDA00015252096000000219
式中,
Figure BDA00015252096000000220
Figure BDA00015252096000000221
分别为常规机组i在调度时段t的向上和向下的备用辅助服务竞价价格,
Figure BDA00015252096000000222
Figure BDA00015252096000000223
分别为常规机组i在调度时段t的向上旋转备用和向下旋转备用;
1-2)确定模型的约束条件,包括电力系统约束条件和供热系统约束条件;具体如下:
1-2-1)电力系统约束条件;具体如下:
1-2-1-1)常规机组约束,表达式如下:
Figure BDA0001525209600000031
式中,
Figure BDA0001525209600000032
Figure BDA0001525209600000033
分别为常规机组i的有功功率上限和有功功率下限;
Figure BDA0001525209600000034
式中,
Figure BDA0001525209600000035
Figure BDA0001525209600000036
分别为常规机组i的向上爬坡能力和向下爬坡能力,△T为每个调度时段的长度;
Figure BDA0001525209600000037
Figure BDA0001525209600000038
1-2-1-2)风电机组约束,表达式如下:
Figure BDA0001525209600000039
1-2-1-3)电力系统运行约束,表达式如下:
Figure BDA00015252096000000310
式中,ILD为电力系统所有负荷的编号集合,Di,t为负荷i在调度时段t的负荷大小;
Figure BDA00015252096000000311
式中,ILN代表电力系统线路编号集合,IEPS代表电力系统中节点编号集合,
Figure BDA00015252096000000312
Figure BDA00015252096000000313
Figure BDA00015252096000000314
分别代表电力系统中与节点i相连接的热电联产机组编号集合、常规机组编号集合、风电机组编号集合与负荷的编号集合,Lj代表电力系统中线路j的最大有功功率,SFi-j为节点i到线路j的转移分布因子;
Figure BDA00015252096000000315
式中,SRUt和SRDt分别代表电力系统在调度时段t的向上旋转备用和向下旋转备用;
1-2-2)供热系统约束条件,具体如下:
1-2-2-1)热电联产机组约束,表达式如下:
Figure BDA00015252096000000316
式中,
Figure BDA00015252096000000317
Figure BDA00015252096000000318
分别为热电联产机组i的机组热电出力特性参数;
Figure BDA00015252096000000319
式中,
Figure BDA0001525209600000041
Figure BDA0001525209600000042
分别为热电联产机组i的向上爬坡能力和向下爬坡能力;
1-2-2-2)供热系统约束,表达式如下:
Figure BDA0001525209600000043
式中,
Figure BDA0001525209600000044
为热电联产机组i所对应的换热站在调度时段t的一次管网入口水温,ti为热电联产机组i传热到其所对应的换热站所需时间,
Figure BDA0001525209600000045
Figure BDA0001525209600000046
分别为热电联产机组i的供热系统热特性参数,
Figure BDA0001525209600000047
为调度时段t-1的外界温度;
Figure BDA0001525209600000048
式中,
Figure BDA0001525209600000049
为热电联产机组i所对应换热站在调度时段t的一次管网入口水温下限值;
2)将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化为标准矩阵形式的二次规划模型;
将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化成矩阵形式:用向量xE表示电力系统变量,电力系统变量包括
Figure BDA00015252096000000410
Figure BDA00015252096000000411
所有电力系统变量顺次排列构成向量xE;用xH表示供热系统变量,供热系统变量包括:
Figure BDA00015252096000000412
Figure BDA00015252096000000413
所有供热系统变量顺次排列构成向量xH;则计及辅助服务的动态供热系统与电力系统耦合调度模型转化为如下的矩阵形式:
Figure BDA00015252096000000414
subject to:
xE∈XE (18)
xH∈XH
DxE=ExH
式中,C表示计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;约束集合XE为电力系统的约束条件集合,即式(6)-(13)与(15)所构成的集合;约束集合XH表示供热系统约束条件集合,即式(16)与(17)所构成的集合;约束DxE=ExH表示电力系统和供热系统耦合约束条件,即变量xE中的
Figure BDA00015252096000000415
和变量xH中的
Figure BDA00015252096000000416
对应相等;
3)对步骤2)得到的二次规划模型求解,将求解得到的
Figure BDA00015252096000000417
分别作为常规机组、风电机组和热电联产机组的发电出力进行调度,将
Figure BDA00015252096000000418
作为热电联产机组的供热出力进行调度。
本发明的技术特点及有益效果在于:
本发明在传统电力系统的经济调度中,考虑了备用辅助服务的成本,避免了传统经济调度中机组备用设置不合理的情况;此外,本发明将动态供热系统引入到电力系统的模型中,利用供热系统的储热性能提升电力系统运行灵活性,调度结果相比与传统方法可显著提升风电消纳水平并降低系统运行成本。
本发明可应用于我国北方冬季集中供热期的电力系统运行调度中。
具体实施方式
本发明提出的计及辅助服务的动态供热系统与电力系统耦合调度方法,下面结合具体实施例进一步详细说明如下。
本发明提出的计及辅助服务的动态供热系统与电力系统耦合调度方法,着重考虑备用辅助服务,首先,针对电力系统建立考虑辅助服务市场的经济调度模型,在经济调度的优化中考虑备用的成本;其次,对供热系统进行系统辨识动态等值建模,基于历史统计信息辨识供热系统的动态特性;然后,结合热电联产机组的特性,建立热电联合优化调度模型;最后,求解热电联合优化调度模型,得到调度策略。该方法包括以下步骤:
1)建立计及辅助服务的动态供热系统与电力系统耦合调度模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;
计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数为全网运行成本最小化,表达式如下:
Figure BDA0001525209600000051
式中,T为调度时段集合(一般可取1~24,即一日24小时,间隔为1小时),ICHP、ITU和IWD分别为热电联合系统中的热电联产机组编号集合、常规机组编号集合和风电机组的编号集合,
Figure BDA0001525209600000052
为热电联合系统中热电联产机组i在t时段的生产成本函数,
Figure BDA0001525209600000053
为热电联合系统中常规机组i在t时段的生产成本函数,
Figure BDA0001525209600000054
为热电联合系统中风电机组i在t时段的生产成本函数,
Figure BDA0001525209600000055
为电力系统中常规机组i在t时段提供的备用辅助服务成本函数,t为调度时段编号;
其中,热电联产机组的生产成本函数为:
Figure BDA0001525209600000056
式中,
Figure BDA0001525209600000057
Figure BDA0001525209600000058
分别为热电联产机组i的生产成本系数,生产成本系数为热电联产机组固有参数(可从热电联产机组设计手册处获取),
Figure BDA0001525209600000061
Figure BDA0001525209600000062
分别为热电联产机组i在调度时段t的有功功率和供热功率;
常规机组的生产成本函数为:
Figure BDA0001525209600000063
式中,
Figure BDA0001525209600000064
Figure BDA0001525209600000065
分别为常规机组i的发电成本系数,该系数为常规机组的固有参数(可从常规机组设计手册处获取),
Figure BDA0001525209600000066
为常规机组i在调度时段t的有功功率;
风电机组的生产成本函数为:
Figure BDA0001525209600000067
式中,
Figure BDA0001525209600000068
为风电机组i的弃风惩罚因子,弃风惩罚因子的取值根据对风电的消纳需求确定,由电力系统调度中心根据调度结果反馈进行调节(常规可取值为
Figure BDA0001525209600000069
Figure BDA00015252096000000610
为风电机组i在调度时段t的可用有功功率,
Figure BDA00015252096000000611
为风电机组i在调度时段t的实际有功功率;
电力系统中常规机组提供的备用辅助服务成本函数为:
Figure BDA00015252096000000612
式中,
Figure BDA00015252096000000613
Figure BDA00015252096000000614
分别为常规机组i在调度时段t的向上和向下的备用辅助服务竞价价格,
Figure BDA00015252096000000615
Figure BDA00015252096000000616
分别为常规机组i在调度时段t的向上旋转备用和向下旋转备用;
1-2)确定模型的约束条件,包括电力系统约束条件和供热系统约束条件;具体如下:
1-2-1)电力系统约束条件;具体如下:
1-2-1-1)常规机组约束,表达式如下:
Figure BDA00015252096000000617
式中,
Figure BDA00015252096000000618
Figure BDA00015252096000000619
分别为常规机组i的有功功率上限和有功功率下限;
Figure BDA00015252096000000620
式中,
Figure BDA00015252096000000621
Figure BDA00015252096000000622
分别为常规机组i的向上爬坡能力和向下爬坡能力,△T为每个调度时段的长度(一般可取值为1小时)。
Figure BDA00015252096000000623
Figure BDA00015252096000000624
1-2-1-2)风电机组约束,表达式如下:
Figure BDA0001525209600000071
1-2-1-3)电力系统运行约束,表达式如下:
Figure BDA0001525209600000072
式中,ILD为电力系统所有负荷的编号集合,Di,t为负荷i在调度时段t的负荷大小。
Figure BDA0001525209600000073
式中,ILN代表电力系统线路编号集合,IEPS代表电力系统中节点编号集合,
Figure BDA0001525209600000074
Figure BDA0001525209600000075
Figure BDA0001525209600000076
分别代表电力系统中与节点i相连接的热电联产机组编号集合、常规机组编号集合、风电机组编号集合与负荷的编号集合,Lj代表电力系统中线路j的最大有功功率,SFi-j为节点i到线路j的转移分布因子,转移分布因子从电力系统运行拓扑图上获取。
Figure BDA0001525209600000077
式中,SRUt和SRDt分别代表电力系统在调度时段t所需的向上旋转备用和向下旋转备用。
1-2-2)供热系统约束条件,具体如下:
1-2-2-1)热电联产机组约束,表达式如下:
Figure BDA0001525209600000078
式中,
Figure BDA0001525209600000079
Figure BDA00015252096000000710
分别为热电联产机组i的机组热电出力特性参数(热电联产机组固有参数,可从热电联产机组设计手册处获取)。
Figure BDA00015252096000000711
式中,
Figure BDA00015252096000000712
Figure BDA00015252096000000713
分别为热电联产机组i的向上爬坡能力和向下爬坡能力;
1-2-2-2)供热系统约束,表达式如下:
Figure BDA00015252096000000714
式中,
Figure BDA00015252096000000715
为热电联产机组i所对应的换热站在调度时段t的一次管网入口水温,ti为热电联产机组i传热到其所对应的换热站所需时间,
Figure BDA00015252096000000716
Figure BDA00015252096000000717
分别为热电联产机组i的供热系统热特性参数(以上参数通过供热系统供热历史运行数据拟合),
Figure BDA00015252096000000718
为调度时段t-1的外界温度。
Figure BDA00015252096000000719
式中,
Figure BDA0001525209600000081
为热电联产机组i所对应的换热站在调度时段t的一次管网入口水温下限值。
2)将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化为标准矩阵形式的二次规划模型;
将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化成矩阵形式:用向量xE表示电力系统变量,电力系统变量包括
Figure BDA0001525209600000082
Figure BDA0001525209600000083
所有电力系统变量顺次排列构成向量xE;用xH表示供热系统变量,供热系统变量包括:
Figure BDA0001525209600000084
Figure BDA0001525209600000085
所有供热系统变量顺次排列构成向量xH;则计及辅助服务的动态供热系统与电力系统耦合调度模型转化为如下的矩阵形式:
Figure BDA0001525209600000086
subject to:
xE∈XE (18)
xH∈XH
DxE=ExH
式中,C表示计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;约束集合XE为电力系统的约束条件集合,即式(6)-(13)与(15)所构成的集合;约束集合XH表示供热系统约束条件集合,即式(16)与(17)所构成的集合;约束DxE=ExH表示电力系统和供热系统耦合约束条件,即变量xE中的
Figure BDA0001525209600000087
和变量xH中的
Figure BDA0001525209600000088
对应相等;D、E的每一行与每一个电力系统和供热系统耦合约束条件一一对应,D的每一列与电力系统中的每一个变量一一对应,E的每一列与供热系统中的每一个变量一一对应,其中D、E的每一个元素为该元素所在行所对应的约束条件中其列所代表的变量的系数;
3)对步骤2)得到的标准矩阵形式的二次规划模型求解,典型的求解算法包括内点法、置信域法、单纯性法等。将求解得到的
Figure BDA0001525209600000089
分别作为常规机组、风电机组和热电联产机组的发电出力进行调度,将
Figure BDA00015252096000000810
作为热电联产机组的供热出力进行调度。

Claims (1)

1.一种计及辅助服务的动态供热系统与电力系统耦合调度方法,其特征在于,该方法包括以下步骤:
1)建立计及辅助服务的动态供热系统与电力系统耦合调度模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;
计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数为全网运行成本最小化,表达式如下:
Figure FDA0001525209590000011
式中,T为调度时段集合,ICHP、ITU和IWD分别为热电联合系统中的热电联产机组编号集合、常规机组编号集合和风电机组的编号集合,
Figure FDA0001525209590000012
为热电联产机组i在t时段的生产成本函数,
Figure FDA0001525209590000013
为常规机组i在t时段的生产成本函数,
Figure FDA0001525209590000014
为风电机组i在t时段的生产成本函数,
Figure FDA0001525209590000015
为常规机组i在t时段提供的备用辅助服务成本函数,t为调度时段编号;
其中,热电联产机组的生产成本函数为:
Figure FDA0001525209590000016
式中,
Figure FDA0001525209590000017
Figure FDA0001525209590000018
分别为热电联产机组i的生产成本系数,
Figure FDA0001525209590000019
Figure FDA00015252095900000110
分别为热电联产机组i在调度时段t的有功功率和供热功率;
常规机组的生产成本函数为:
Figure FDA00015252095900000111
式中,
Figure FDA00015252095900000112
Figure FDA00015252095900000113
分别为常规机组i的发电成本系数,
Figure FDA00015252095900000114
为常规机组i在调度时段t的有功功率;
风电机组的生产成本函数为:
Figure FDA00015252095900000115
式中,
Figure FDA0001525209590000021
为风电机组i的弃风惩罚因子,
Figure FDA0001525209590000022
为风电机组i在调度时段t的可用有功功率,
Figure FDA0001525209590000023
为风电机组i在调度时段t的实际有功功率;
常规机组提供的备用辅助服务成本函数为:
Figure FDA0001525209590000024
式中,
Figure FDA0001525209590000025
Figure FDA0001525209590000026
分别为常规机组i在调度时段t的向上和向下的备用辅助服务竞价价格,
Figure FDA0001525209590000027
Figure FDA0001525209590000028
分别为常规机组i在调度时段t的向上旋转备用和向下旋转备用;
1-2)确定模型的约束条件,包括电力系统约束条件和供热系统约束条件;具体如下:
1-2-1)电力系统约束条件;具体如下:
1-2-1-1)常规机组约束,表达式如下:
Figure FDA0001525209590000029
式中,
Figure FDA00015252095900000210
Pi TU 分别为常规机组i的有功功率上限和有功功率下限;
Figure FDA00015252095900000211
式中,
Figure FDA00015252095900000212
Figure FDA00015252095900000213
分别为常规机组i的向上爬坡能力和向下爬坡能力,△T为每个调度时段的长度;
Figure FDA00015252095900000214
Figure FDA00015252095900000215
1-2-1-2)风电机组约束,表达式如下:
Figure FDA00015252095900000216
1-2-1-3)电力系统运行约束,表达式如下:
Figure FDA00015252095900000217
式中,ILD为电力系统所有负荷的编号集合,Di,t为负荷i在调度时段t的负荷大小;
Figure FDA00015252095900000218
式中,ILN代表电力系统线路编号集合,IEPS代表电力系统中节点编号集合,
Figure FDA00015252095900000219
Figure FDA00015252095900000220
Figure FDA00015252095900000221
分别代表电力系统中与节点i相连接的热电联产机组编号集合、常规机组编号集合、风电机组编号集合与负荷的编号集合,Lj代表电力系统中线路j的最大有功功率,SFi-j为节点i到线路j的转移分布因子;
Figure FDA0001525209590000031
式中,SRUt和SRDt分别代表电力系统在调度时段t的向上旋转备用和向下旋转备用;
1-2-2)供热系统约束条件,具体如下:
1-2-2-1)热电联产机组约束,表达式如下:
Figure FDA0001525209590000032
式中,
Figure FDA0001525209590000033
Figure FDA0001525209590000034
分别为热电联产机组i的机组热电出力特性参数;
Figure FDA0001525209590000035
式中,
Figure FDA0001525209590000036
Figure FDA0001525209590000037
分别为热电联产机组i的向上爬坡能力和向下爬坡能力;
1-2-2-2)供热系统约束,表达式如下:
Figure FDA0001525209590000038
式中,
Figure FDA0001525209590000039
为热电联产机组i所对应的换热站在调度时段t的一次管网入口水温,ti为热电联产机组i传热到其所对应的换热站所需时间,
Figure FDA00015252095900000310
Figure FDA00015252095900000311
分别为热电联产机组i的供热系统热特性参数,
Figure FDA00015252095900000312
为调度时段t-1的外界温度;
Figure FDA00015252095900000313
式中,
Figure FDA00015252095900000314
为热电联产机组i所对应换热站在调度时段t的一次管网入口水温下限值;
2)将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化为标准矩阵形式的二次规划模型;
将步骤1)建立的计及辅助服务的动态供热系统与电力系统耦合调度模型转化成矩阵形式:用向量xE表示电力系统变量,电力系统变量包括
Figure FDA00015252095900000315
Figure FDA00015252095900000316
所有电力系统变量顺次排列构成向量xE;用xH表示供热系统变量,供热系统变量包括:
Figure FDA00015252095900000317
Figure FDA00015252095900000318
所有供热系统变量顺次排列构成向量xH;则计及辅助服务的动态供热系统与电力系统耦合调度模型转化为如下的矩阵形式:
Figure FDA00015252095900000319
式中,C表示计及辅助服务的动态供热系统与电力系统耦合调度模型的目标函数;约束集合XE为电力系统的约束条件集合,即式(6)-(13)与(15)所构成的集合;约束集合XH表示供热系统约束条件集合,即式(16)与(17)所构成的集合;约束DxE=ExH表示电力系统和供热系统耦合约束条件,即变量xE中的
Figure FDA0001525209590000041
和变量xH中的
Figure FDA0001525209590000042
对应相等;
3)对步骤2)得到的二次规划模型求解,将求解得到的
Figure FDA0001525209590000043
分别作为常规机组、风电机组和热电联产机组的发电出力进行调度,将
Figure FDA0001525209590000044
作为热电联产机组的供热出力进行调度。
CN201711439071.1A 2017-12-26 2017-12-26 计及辅助服务的动态供热系统与电力系统耦合调度方法 Active CN108155644B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711439071.1A CN108155644B (zh) 2017-12-26 2017-12-26 计及辅助服务的动态供热系统与电力系统耦合调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711439071.1A CN108155644B (zh) 2017-12-26 2017-12-26 计及辅助服务的动态供热系统与电力系统耦合调度方法

Publications (2)

Publication Number Publication Date
CN108155644A CN108155644A (zh) 2018-06-12
CN108155644B true CN108155644B (zh) 2020-04-24

Family

ID=62462200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711439071.1A Active CN108155644B (zh) 2017-12-26 2017-12-26 计及辅助服务的动态供热系统与电力系统耦合调度方法

Country Status (1)

Country Link
CN (1) CN108155644B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109256787B (zh) * 2018-11-07 2021-06-22 清华-伯克利深圳学院筹备办公室 一种供电系统调频方法、装置、控制设备及存储介质
CN111146819B (zh) * 2019-12-26 2021-03-02 上海电力大学 一种考虑热网特性的电热联合系统辅助服务调度方法
CN117634682B (zh) * 2023-11-28 2024-06-18 国网吉林省电力有限公司经济技术研究院 一种寒冷地区电-热联供型虚拟电厂优化调控方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105761109A (zh) * 2016-02-19 2016-07-13 刘隽琦 虚拟电厂能源管理和电力交易的智能管理系统及其优化运行方法
CN106339772B (zh) * 2016-08-11 2019-06-18 清华大学 基于供热管网储热效益的热-电联合优化调度方法
CN106815661B (zh) * 2017-02-22 2020-10-20 清华大学 一种热电联合系统的分解协调调度方法
CN106992555B (zh) * 2017-04-25 2019-06-04 国网吉林省电力有限公司 一种辅助服务市场环境下的热电联合系统分布式调度方法

Also Published As

Publication number Publication date
CN108155644A (zh) 2018-06-12

Similar Documents

Publication Publication Date Title
CN109063925B (zh) 一种计及负荷聚合商的区域综合能源系统优化运行方法
WO2022100091A1 (zh) 一种广义源储系统调度的集中控制方法
CN103296682B (zh) 一种多时空尺度渐进趋优的负荷调度模式设计方法
CN111738502A (zh) 促进富余风电消纳的多能互补系统需求响应运行优化方法
CN108155644B (zh) 计及辅助服务的动态供热系统与电力系统耦合调度方法
CN105470982A (zh) 一种含介质储能的智能微电网发电功率控制系统及控制方法
CN110829504B (zh) 弃风参与电网调频的电转气-储气-燃气轮机容量优化配置方法及系统
CN106712033A (zh) 一种热电厂消纳弃风的方法
CN114154328A (zh) 提升灵活性的电热综合能源系统两阶段随机优化调度方法
CN113240166A (zh) 考虑高比例新能源消纳的微能源网日前经济调度方法
CN113128844A (zh) 一种基于供电设备容量限制下分布式电源规划方法
CN109376406B (zh) 供能系统超结构模型、建模方法、计算机设备和存储介质
CN109615125B (zh) 一种考虑特高压调峰的多区域随机生产模拟方法及应用
CN109583753B (zh) 一种基于区域内部交易的智能电网监控系统及其控制方法
CN104239960A (zh) 考虑抽水蓄能机组的发电计划优化方法
CN112531785B (zh) 多时间尺度下光热-储热的主动配电网主辅联合调度方法
CN111125611B (zh) 面向多场景的冷-热-电微能网群两阶段优化调度方法
CN116191577A (zh) 一种涉及多电源类型的地区电网多时间尺度无功优化方法
CN106960254B (zh) 一种考虑风电消纳的电转气设备容量优化配置方法
CN115566731A (zh) 一种基于系统多能互补的灵活性电源规划方法
CN112736950B (zh) 一种面向微电网群的公共储能电站配置方法及系统
Fang et al. Optimal operation strategy considering wind power accommodation in heating district
CN113762643A (zh) 区域综合能源系统的储能容量优化配置方法
CN113627720A (zh) 基于电能替代的综合园区风光储系统配置方法
CN110415138B (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
GR01 Patent grant
GR01 Patent grant