CN105656064B - 利用供热时滞性实现热电机组参与系统调峰调度的方法 - Google Patents

利用供热时滞性实现热电机组参与系统调峰调度的方法 Download PDF

Info

Publication number
CN105656064B
CN105656064B CN201511024034.5A CN201511024034A CN105656064B CN 105656064 B CN105656064 B CN 105656064B CN 201511024034 A CN201511024034 A CN 201511024034A CN 105656064 B CN105656064 B CN 105656064B
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
unit
heat supply
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
CN201511024034.5A
Other languages
English (en)
Other versions
CN105656064A (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.)
Northeast Electric Power University
Original Assignee
Northeast Dianli University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN201511024034.5A priority Critical patent/CN105656064B/zh
Publication of CN105656064A publication Critical patent/CN105656064A/zh
Application granted granted Critical
Publication of CN105656064B publication Critical patent/CN105656064B/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/28Arrangements for balancing of the load in a network by storage of energy
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24DDOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
    • F24D19/00Details
    • F24D19/10Arrangement or mounting of control or safety devices
    • F24D19/1006Arrangement or mounting of control or safety devices for water heating systems
    • F24D19/1009Arrangement or mounting of control or safety devices for water heating systems for central heating
    • 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
    • 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/003Load forecast, e.g. methods or systems for forecasting future load demand
    • 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
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B10/00Integration of renewable energy sources in buildings
    • Y02B10/70Hybrid systems, e.g. uninterruptible or back-up power supplies integrating renewable energies
    • 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
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Thermal Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Coloring Foods And Improving Nutritive Qualities (AREA)
  • Wind Motors (AREA)

Abstract

本发明是一种利用供热时滞性热电机组参与调峰调度的方法,其特点是,包括供热时滞性的数学表达形式、总热负荷的优化调度、利用时滞性调节热电机组的出力范围、风‑电‑热综合调度模型和优化模型求解策略等内容。利用供热时滞性,在几乎不影响人民工作生活需要的前提下,一定程度上解决了以往热电机组不参与调峰或参与调峰容量小为电网调峰带来的难题。具有计算简单、结果合理、易于实现且应用价值高等优点。

Description

利用供热时滞性实现热电机组参与系统调峰调度的方法
技术领域
本发明涉及电力系统有功调度领域,是一种利用供热时滞性实现热电机组参与系统调峰调度的方法。
背景技术
当前,我国北方地区冬季供暖期热电机组按照“以热定电”约束运行,热电机组调峰范围小,加之风电反调峰特性,弃风现象严重,导致大量清洁能源的浪费。目前电网中,传统机组的结构已经基本固定,机组的组合模式难以改变,所以欲提高系统调峰能力,挖掘热电机组出力显得尤为重要。事实上,热网和建筑物的储热能力很强,考虑供热时滞性后,抽汽式机组可以解耦其“以热定电”约束,扩大出力调整范围,进而提高调峰能力,以接纳更多风电。现有研究对热网和建筑物的供热时滞性进行了描述,采取多种方法以增强系统调峰能力。热电联产与分产供热相结合的方式不足以满足供热质量;用户对供热的需求具有一定时滞性,短时间内对供热量的调整不至于大范围影响用户需求。研究削减供热量及削减时间对用户室温的影响有助于将供热时滞性应用于热电机组调峰中;蓄热装置能储存多余的能源,在需要时释放,目前蓄热装置主要安装在热电厂附近,且其经济效益不高,对于热电系统来说,用户是很好的蓄热载体,且不需要花费额外费用。因此,将供热受端用户作为蓄热载体,利用供热时滞性增大热电机组调峰能力是一种经济可行的方法,有助于减小系统调峰压力。
室温的变化不仅受供热影响,自然影响也占有很大比重。故室内温度波动性不可避免,即使热电机组出力波动范围不大,实际生活中,温度也不是保持不变的,热电机组利用供热时滞性调峰所产生的温度波动是可以控制在此范围内的。如在高峰负荷,增大热电机组出力和供热量,室温升高,此时供热系统作为一个大的储热装置储存一定热量;在低谷负荷,热电机组减少出力和供热量,由供热系统储存的热量配合供热。
发明内容
本发明的目的是,提供一种具有计算简单、结果合理、易于实现且应用价值高的利用供热时滞性热电机组参与调峰调度的方法。
本发明的目的是由以下技术方案来实现的:一种利用供热时滞性热电机组参与调峰调度的方法,其特征是,它包括以下内容:
1)供热时滞性的数学表达形式
用以下模型来表示热网回水温度τh,t、建筑物室温τn,t与热网供水温度τg,t、室外温度τw,t之间的关系,供热时滞性的约束条件表示为
热电厂抽汽供热流量与热网供回水温度之间的关系用下式表示
式中:J、E表示模型的阶次,温度的单位统一为℃;为抽汽机组t时段抽汽供热流量,MW;G为热网水流量,t/h;cp为定压比热容,J/(kg·℃);
2)总热负荷的优化调度
通过当前时刻的采暖建筑物的室内温度、供热流量、室外温度预测下一时刻的室内温度,根据下一时刻室内温度的变化趋势,实时调节热电机组的供热流量,合理分配全网的供热流量,使所有参与调度的抽汽式机组的发电功率调节范围之和最大,目标函数表示为:
约束条件:
ui+vi=1 (7)
0≤vihi≤vihmed,i (9)
目标函数中,cv保持恒定进气量条件下,改变供热抽取两时发电机功率的较少量,cv,1为最大电力出力条件下与之对应的cv值,cv,2为最小电力出力条件下与之对应的cv值;cm为在背压运行条件下热功率和电功率的弹性系数;hmed为在机组发电功率最小条件下的供热功率;为抽气机组的供热功率上限值;分别为抽气机组的发电功率上、下限值;Ki为常数,β是系统内热电机组的数目;约束条件中,供热功率hi要满足功率上下限约束;各发电机组的发热总和要满足供求平衡;ui、vi为0、1变量,分别表示抽气机组的关状态(0)和开状态(1);
3)利用时滞性调节热电机组的出力范围
基于抽汽式机组的运行工况,对于一个给定的供热功率h,其输出电功率可在一定可行的范围内进行调节,改变供热功率h值,会使机组的出力调节范围产生相应的影响;传统的热电机组的调度方式下,抽汽式机组的供热水平基本保持不变,因此抽汽式机组的调峰能力有限,抽汽式机组的电热特性数学描述如下:
其中,分别为抽汽式机组i在凝气工况下最小、最大有功出力;为抽汽式机组i在t时刻的发电出力;式(10)给出了抽汽式机组热功率与发电功率之间的数学关系,由此能够计算出在给定的热功率的条件下抽汽式机组的发电功率的上下限,考虑了热网和采暖建筑物的供热时滞性后,热电机组的供热流量根据采暖室温的变化作相应的调整,通过降低热电机组在负荷低谷时段的供热水平,来增大抽汽式机组的发电功率的上下限,从而使热电机组具备调峰的能力;
4)风-电-热综合调度模型
在满足电、热平衡的前提下,以系统发电成本最小为目标,通过对热电机组、常规机组和风电机组的出力的优化分配,因此,优化目标函数为
其中:F是系统的煤耗总量;是系统内第i台热电机组t时段的煤耗量,且是系统内第i台常规机组t时段的煤耗量,β是系统内热电机组的数目;δ是系统内常规机组的数目;T为运行周期;
约束条件
A.供电平衡
忽略网损条件下的系统供电平衡约束
其中:为系统内第i台抽汽式机组t时段的电出力;是系统内第i台纯凝机组t时段的电出力;是风电场t时段实际出力;是系统在t时段的负荷预测值;
B.供热平衡
忽略传输损耗条件下的供热平衡约束
其中:是第j座热电厂i台机组t时段的热出力;是区域i在t时段的热负荷预测值;
C.各类机组电出力约束
D.抽汽式机组电出力约束:如式(10);
E.纯凝式机组电出力约束:
F.风电场出力约束
其中,分别为纯凝机组电出力最小、最大值;是风电场在t时段的预测出力;
G.抽汽式机组热出力约束
抽汽式机组的热功率约束按照式(4)~(9);
H.机组爬坡率约束
其中,为纯凝式机组爬坡速率下、上限值(MW/min),T60为单段运行时长,通常为60min,机组最大爬坡速率一般处于最大功率4%~5%的水平;
5)优化模型求解策略
首先采集历史时刻的典型负荷数据,确定需要优化的时段和机组类型;在未进行优化调度时,计算室温和供水温度;获得小时级下的室内温度与供热功率之间的关系,确定热电机组的总供热功率;在保证全网热电机组的电出力调节范围最大的前提下,采用0-1整数规划的方法,求解各台热电机组的热出力,由此确定热电机组的发电出力上下限,模型以24小时为一个运行周期,以4小时为一个时段,从第N时段起采集前一时段中t-1、t-2、t-3时刻的室内温度、室外温度、供回水温度作为已知条件,逐时段向前递推室内温度和供热流量,以常规机组、热电机组和风电的发电功率作为决策变量进行优化计算。
本发明的利用供热时滞性实现热电机组参与系统调峰调度的方法,通过对一个含热电机组、常规机组和风电场的电热综合调度,将风—电—热联合调度优化运行视为初始条件的非线性规划问题,模型以24小时为一个运行周期,以4小时为一个时段,从第N时段起采集前一时段中t-1、t-2、t-3时刻的室内温度、室外温度、供回水温度作为已知条件,逐时段向前递推室内温度和供热流量,以常规机组、热电机组和风电的发电功率作为决策变量进行优化计算;利用供热时滞性,在几乎不影响人民工作生活需要的前提下,一定程度上解决了以往热电机组不参与调峰或参与调峰容量小为电网调峰带来的难题。具有计算简单、结果合理、易于实现且应用价值高等优点。
附图说明
图1抽汽式供热机组运行工况示意图;
图2为200MW机组优化运行前后电出力上下限对比图;
图3为优化运行前后室内温度对比图。
具体实施方式
下面利用附图和实施例对本发明作进一步说明。
本发明的一种利用供热时滞性热电机组参与调峰调度的方法,包括以下内容:
1)供热时滞性的数学表达形式
用以下模型来表示热网回水温度τh,t、建筑物室温τn,t与热网供水温度τg,t、室外温度τw,t之间的关系,供热时滞性的约束条件表示为
热电厂抽汽供热流量与热网供回水温度之间的关系用下式表示
式中:J、E表示模型的阶次,温度的单位统一为℃;为抽汽机组t时段抽汽供热流量,MW;G为热网水流量,t/h;cp为定压比热容,J/(kg·℃);
2)总热负荷的优化调度
通过当前时刻的采暖建筑物的室内温度、供热流量、室外温度预测下一时刻的室内温度,根据下一时刻室内温度的变化趋势,实时调节热电机组的供热流量,合理分配全网的供热流量,使所有参与调度的抽汽式机组的发电功率调节范围之和最大,目标函数表示为:
约束条件:
ui+vi=1 (7)
0≤vihi≤vihmed,i (9)
目标函数中,cv保持恒定进气量条件下,改变供热抽取量时发电机功率的较少量,cv,1为最大电力出力条件下与之对应的cv值,cv,2为最小电力出力条件下与之对应的cv值;cm为在背压运行条件下热功率和电功率的弹性系数;hmed为在机组发电功率最小条件下的供热功率;为抽气机组的供热功率上限值;分别为抽气机组的发电功率上、下限值;Ki为常数,β是系统内热电机组的数目;约束条件中,供热功率hi要满足功率上下限约束;各发电机组的发热总和要满足供求平衡;ui、vi为0、1变量,分别表示抽气机组的关状态(0)和开状态(1);
3)利用时滞性调节热电机组的出力范围
基于抽汽式机组的运行工况,对于一个给定的供热功率h,其输出电功率可在一定可行的范围内进行调节,改变供热功率h值,会使机组的出力调节范围产生相应的影响;传统的热电机组的调度方式下,抽汽式机组的供热水平基本保持不变,因此抽汽式机组的调峰能力有限,抽汽式机组的电热特性数学描述如下:
其中,分别为抽汽式机组i在凝气工况下最小、最大有功出力;为抽汽式机组i在t时刻的发电出力;式(10)给出了抽汽式机组热功率与发电功率之间的数学关系,由此能够计算出在给定的热功率的条件下抽汽式机组的发电功率的上下限,考虑了热网和采暖建筑物的供热时滞性后,热电机组的供热流量根据采暖室温的变化作相应的调整,通过降低热电机组在负荷低谷时段的供热水平,来增大抽汽式机组的发电功率的上下限,从而使热电机组具备调峰的能力;
4)风-电-热综合调度模型
在满足电、热平衡的前提下,以系统发电成本最小为目标,通过对热电机组、常规机组和风电机组的出力的优化分配,因此,优化目标函数为
其中:F是系统的煤耗总量;是系统内第i台热电机组t时段的煤耗量,且是系统内第i台常规机组t时段的煤耗量,β是系统内热电机组的数目;δ是系统内常规机组的数目;T为运行周期;
约束条件
A.供电平衡
忽略网损条件下的系统供电平衡约束
其中:为系统内第i台抽汽式机组t时段的电出力;是系统内第i台纯凝机组t时段的电出力;是风电场t时段实际出力;是系统在t时段的负荷预测值;
B.供热平衡
忽略传输损耗条件下的供热平衡约束
其中:是第j座热电厂i台机组t时段的热出力;是区域i在t时段的热负荷预测值;
C.各类机组电出力约束
D.抽汽式机组电出力约束:如式(10);
E.纯凝式机组电出力约束:
F.风电场出力约束:
其中,分别为纯凝机组电出力最小、最大值;是风电场在t时段的预测出力;
G.抽汽式机组热出力约束
抽汽式机组的热功率约束按照式(4)~(9);
H.机组爬坡率约束
其中,为纯凝式机组爬坡速率下、上限值(MW/min),T60为单段运行时长,通常为60min,机组最大爬坡速率一般处于最大功率4%~5%的水平;
5)优化模型求解策略
首先采集历史时刻的典型负荷数据,确定需要优化的时段和机组类型;在未进行优化调度时,计算室温和供水温度;获得小时级下的室内温度与供热功率之间的关系,确定热电机组的总供热功率;在保证全网热电机组的电出力调节范围最大的前提下,采用0-1整数规划的方法,求解各台热电机组的热出力,由此确定热电机组的发电出力上下限,模型以24小时为一个运行周期,以4小时为一个时段,从第N时段起采集前一时段中t-1、t-2、t-3时刻的室内温度、室外温度、供回水温度作为已知条件,逐时段向前递推室内温度和供热流量,以常规机组、热电机组和风电的发电功率作为决策变量进行优化计算。
按照本发明的利用供热时滞性实现热电机组参与系统调峰调度的方法,大幅增大热电机组调峰能力。如图2所示,以200MW机组为例,在优化前,热电机组有功输出的调节范围为37.9MW,且极端情况下仅有25MW;通过利用供热时滞性,热电机组在需要参与调峰时有功输出的调节范围增大到51.9MW,较优化前正常运行状态增大了36.94%,较优化前极端情况增大了107.6%;而优化前后室温的变化却极小,如图3所示,优化前,室温一直在20.8摄氏度至21.8摄氏度之间波动,优化运行后,室温在20.6摄氏度至21.6摄氏度之间波动,室内最高、最低温度均降低0.2摄氏度,这对人民工作生活几乎没有影响。

Claims (1)

1.一种利用供热时滞性热电机组参与调峰调度的方法,其特征是,它包括以下内容:
1)供热时滞性的数学表达形式
用以下模型来表示热网回水温度τh,t、建筑物室温τn,t与热网供水温度τg,t、室外温度τw,t之间的关系,供热时滞性的约束条件表示为
<mrow> <msub> <mi>&amp;tau;</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>j</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>J</mi> </munderover> <msub> <mi>&amp;beta;</mi> <mi>j</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>g</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>J</mi> </munderover> <msub> <mi>&amp;gamma;</mi> <mi>j</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>w</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>&amp;tau;</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>e</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>E</mi> </munderover> <msub> <mi>&amp;theta;</mi> <mi>e</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>e</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>E</mi> </munderover> <msub> <mi>&amp;phi;</mi> <mi>e</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>g</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>e</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>E</mi> </munderover> <msub> <mi>&amp;omega;</mi> <mi>e</mi> </msub> <msub> <mi>&amp;tau;</mi> <mrow> <mi>w</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mi>e</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
热电厂抽汽供热流量与热网供回水温度之间的关系用下式表示
<mrow> <msubsup> <mi>P</mi> <mi>h</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msub> <mi>Gc</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;tau;</mi> <mrow> <mi>g</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;tau;</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中:J、E表示模型的阶次,温度的单位统一为℃;为抽汽机组t时段抽汽供热流量,MW;G为热网水流量,t/h;cp为定压比热容,J/(kg·℃);
2)总热负荷的优化调度
通过当前时刻的采暖建筑物的室内温度、供热流量、室外温度预测下一时刻的室内温度,根据下一时刻室内温度的变化趋势,实时调节热电机组的供热流量,合理分配全网的供热流量,使所有参与调度的抽汽式机组的发电功率调节范围之和最大,目标函数表示为:
<mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;beta;</mi> </munderover> <msubsup> <mi>P</mi> <mi>i</mi> <mi>max</mi> </msubsup> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>v</mi> <mn>1</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>&amp;lsqb;</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>c</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mi>i</mi> <mi>min</mi> </msubsup> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>v</mi> <mn>2</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
约束条件:
<mrow> <mn>0</mn> <mo>&amp;le;</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mi>h</mi> <mrow> <mi>T</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;beta;</mi> </munderover> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>H</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
ui+vi=1 (7)
<mrow> <msub> <mi>u</mi> <mi>i</mi> </msub> <msub> <mi>h</mi> <mrow> <mi>m</mi> <mi>e</mi> <mi>d</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <msubsup> <mi>h</mi> <mrow> <mi>T</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
0≤vihi≤vihmed,i (9)
目标函数中,cv保持恒定进气量条件下,改变供热抽取量时发电机功率的较少量,cv,1为最大电力出力条件下与之对应的cv值,cv,2为最小电力出力条件下与之对应的cv值;cm为在背压运行条件下热功率和电功率的弹性系数;hmed为在机组发电功率最小条件下的供热功率;为抽气机组的供热功率上限值;Pi min、Pi max分别为抽气机组的发电功率上、下限值;Ki为常数,β是系统内热电机组的数目;约束条件中,供热功率hi要满足功率上下限约束;各发电机组的发热总和要满足供求平衡;ui、vi为0、1变量,分别表示抽气机组的关状态(0)和开状态(1);
3)利用时滞性调节热电机组的出力范围
基于抽汽式机组的运行工况,对于一个给定的供热功率h,其输出电功率可在一定可行的范围内进行调节,改变供热功率h值,会使机组的出力调节范围产生相应的影响;传统的热电机组的调度方式下,抽汽式机组的供热水平基本保持不变,因此抽汽式机组的调峰能力有限,抽汽式机组的电热特性数学描述如下:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;GreaterEqual;</mo> <mi>min</mi> <mrow> <mo>(</mo> <msub> <mi>c</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>,</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>min</mi> </msubsup> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>v</mi> <mn>2</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mi>v</mi> <mn>1</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中,分别为抽汽式机组i在凝气工况下最小、最大有功出力;为抽汽式机组i在t时刻的发电出力;式(10)给出了抽汽式机组热功率与发电功率之间的数学关系,由此能够计算出在给定的热功率的条件下抽汽式机组的发电功率的上下限,考虑了热网和采暖建筑物的供热时滞性后,热电机组的供热流量根据采暖室温的变化作相应的调整,通过降低热电机组在负荷低谷时段的供热水平,来增大抽汽式机组的发电功率的上下限,从而使热电机组具备调峰的能力;
4)风-电-热综合调度模型
在满足电、热平衡的前提下,以系统发电成本最小为目标,通过对热电机组、常规机组和风电机组的出力的优化分配,因此,优化目标函数为
其中:F是系统的煤耗总量;是系统内第i台热电机组t时段的煤耗量,且 是系统内第i台常规机组t时段的煤耗量,β是系统内热电机组的数目;δ是系统内常规机组的数目;T为运行周期;
约束条件
A.供电平衡
忽略网损条件下的系统供电平衡约束
<mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mi>&amp;beta;</mi> </mrow> </munder> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mi>&amp;delta;</mi> </mrow> </munder> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>f</mi> <mi>s</mi> </mrow> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>P</mi> <mi>D</mi> <mi>t</mi> </msubsup> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
其中:为系统内第i台抽汽式机组t时段的电出力;是系统内第i台纯凝机组t时段的电出力;是风电场t时段实际出力;是系统在t时段的负荷预测值;
B.供热平衡
忽略传输损耗条件下的供热平衡约束
<mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mi>&amp;beta;</mi> </mrow> </munder> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>j</mi> </msubsup> <mo>=</mo> <msubsup> <mi>P</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>h</mi> </mrow> <mi>j</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
其中:是第j座热电厂i台机组t时段的热出力;是区域i在t时段的热负荷预测值;
C.各类机组电出力约束
D.抽汽式机组电出力约束:如式(10);
E.纯凝式机组电出力约束:
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
F.风电场出力约束:
<mrow> <mn>0</mn> <mo>&lt;</mo> <msubsup> <mi>P</mi> <mrow> <mi>f</mi> <mi>s</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>W</mi> <mi>D</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
其中,分别为纯凝机组电出力最小、最大值;是风电场在t时段的预测出力;
G.抽汽式机组热出力约束
抽汽式机组的热功率约束按照式(4)~(9);
H.机组爬坡率约束
<mrow> <mo>-</mo> <msubsup> <mi>r</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>d</mi> <mi>o</mi> <mi>w</mi> <mi>n</mi> </mrow> </msubsup> <msub> <mi>T</mi> <mn>60</mn> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>r</mi> <mrow> <mi>C</mi> <mi>H</mi> <mi>P</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msubsup> <msub> <mi>T</mi> <mn>60</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mo>-</mo> <msubsup> <mi>r</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>d</mi> <mi>o</mi> <mi>w</mi> <mi>n</mi> </mrow> </msubsup> <msub> <mi>T</mi> <mn>60</mn> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>r</mi> <mrow> <mi>C</mi> <mi>O</mi> <mi>N</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msubsup> <msub> <mi>T</mi> <mn>60</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
其中,为第i台纯凝式机组向下/上最大爬坡速率限值,MW/min,为第i台抽汽式机组向下/上最大爬坡速率限值,MW/min;T60为单段运行时长,为60min,机组最大爬坡速率处于最大功率4%~5%的水平;
5)优化模型求解策略
首先采集历史时刻的典型负荷数据,确定需要优化的时段和机组类型;在未进行优化调度时,计算室温和供水温度;获得小时级下的室内温度与供热功率之间的关系,确定热电机组的总供热功率;在保证全网热电机组的电出力调节范围最大的前提下,采用0-1整数规划的方法,求解各台热电机组的热出力,由此确定热电机组的发电出力上下限,模型以24小时为一个运行周期,以4小时为一个时段,从第N时段起采集前一时段中t-1、t-2、t-3时刻的室内温度、室外温度、供回水温度作为已知条件,逐时段向前递推室内温度和供热流量,以常规机组、热电机组和风电的发电功率作为决策变量进行优化计算。
CN201511024034.5A 2015-12-31 2015-12-31 利用供热时滞性实现热电机组参与系统调峰调度的方法 Active CN105656064B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511024034.5A CN105656064B (zh) 2015-12-31 2015-12-31 利用供热时滞性实现热电机组参与系统调峰调度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511024034.5A CN105656064B (zh) 2015-12-31 2015-12-31 利用供热时滞性实现热电机组参与系统调峰调度的方法

Publications (2)

Publication Number Publication Date
CN105656064A CN105656064A (zh) 2016-06-08
CN105656064B true CN105656064B (zh) 2018-02-02

Family

ID=56490041

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511024034.5A Active CN105656064B (zh) 2015-12-31 2015-12-31 利用供热时滞性实现热电机组参与系统调峰调度的方法

Country Status (1)

Country Link
CN (1) CN105656064B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106593771B (zh) * 2016-12-16 2019-06-18 国网辽宁省电力有限公司 一种热电联产系统中的电热解耦系统及其控制方法
CN107067146B (zh) * 2017-01-09 2021-01-22 国网浙江省电力有限公司经济技术研究院 考虑热负荷二维可控性的微网热电协调调度方法
CN107120721B (zh) * 2017-05-25 2019-08-09 石家庄华浩能源科技有限公司 一种集中供热动态气候补偿方法
CN108321839B (zh) * 2018-02-09 2020-01-10 大连理工大学 考虑发电集团利益主体协调的风电消纳调度方法
CN108534113A (zh) * 2018-04-28 2018-09-14 赫普科技发展(北京)有限公司 一种负荷侧电蓄热锅炉调频系统和方法
CN110566295B (zh) * 2019-07-27 2021-11-26 华电电力科学研究院有限公司 基于全厂电、热负荷协同调度的双机组耦合调峰方法及装置
CN110930046B (zh) * 2019-11-29 2023-08-22 国网河北省电力有限公司电力科学研究院 一种供热机组深度调峰控制策略及其系统
CN110991749A (zh) * 2019-12-06 2020-04-10 国网冀北电力有限公司电力科学研究院 供热负荷预测方法及装置
CN113834111B (zh) * 2020-06-24 2023-07-14 国网浙江省电力有限公司电力科学研究院 一种天然气分布式能源系统的控制方法
CN111987716A (zh) * 2020-08-17 2020-11-24 南京工程学院 一种多类别蓄热电采暖用户负荷群组合响应调度方法
CN112577100B (zh) * 2020-11-26 2021-12-10 清华大学 供热补偿方法、装置及可读存储介质
CN113587208A (zh) * 2021-08-03 2021-11-02 中国华能集团清洁能源技术研究院有限公司 一种蒸汽管网储能在线定量计算方法和系统
CN113739296B (zh) * 2021-09-08 2022-09-06 山东佐耀科技有限公司 基于模型预测控制的空气源热泵负荷水温控制方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504608A (zh) * 2014-11-04 2015-04-08 国家电网公司 配置储热的热电机组在分时电价下的最优日计划确定方法
CN104967126B (zh) * 2015-07-14 2016-10-05 华中科技大学 一种面向区域电网的跨流域水电站群多电网联合调峰方法

Also Published As

Publication number Publication date
CN105656064A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN105656064B (zh) 利用供热时滞性实现热电机组参与系统调峰调度的方法
CN109063925B (zh) 一种计及负荷聚合商的区域综合能源系统优化运行方法
AU2020100983A4 (en) Multi-energy complementary system two-stage optimization scheduling method and system considering source-storage-load cooperation
CN110535185B (zh) 一种考虑需求响应的含光热电站多源优化运行方法
CN105869075A (zh) 一种冷热电联供型微型能源网经济优化调度方法
CN102520675A (zh) 燃气联合循环与太阳能发电联合制热系统及其调度方法
CN106712033B (zh) 一种热电厂消纳弃风的方法
CN110361969B (zh) 一种冷热电综合能源系统优化运行方法
CN112434915B (zh) 面向弃风消纳的区域综合能源系统灵活性优化配置方法
CN110165665A (zh) 一种基于改进多目标粒子群算法的源-荷-储调度方法
CN113592200B (zh) 一种含水源热泵的区域综合能源系统低碳优化运行方法
CN115857348A (zh) 考虑两联供热泵舒适供能的分布式能量系统容量优化方法
CN106953367A (zh) 用于含燃气‑蒸汽联合循环机组热电系统的风电消纳方法
CN110930046A (zh) 一种供热机组深度调峰控制策略及其系统
CN110286590A (zh) 考虑多类型能源存储的冷热电联供系统日前优化调度方法
CN113555907B (zh) 一种考虑非理想通信情况的分布式多能系统滚动优化控制方法
CN114021361A (zh) 一种含蓄热电锅炉负荷的多时间尺度优化调度方法及系统
Wu et al. Day-ahead optimal dispatch with CHP and wind turbines based on room temperature control
JP5903731B2 (ja) 建物用空調熱源システム
Pasban-Gajan et al. Optimal scheduling of renewable-based energy hubs considering time-of-use pricing scheme
CN109921447A (zh) 一种基于储能装置soc动态约束的微网经济调度方法
CN113420413B (zh) 一种基于负荷可塑性的柔性负荷可调节能力量化方法及系统
CN112529728B (zh) 一种基于非合作博弈的综合能源系统优化调度方法
Li et al. Capacity Configuration of Integrated Energy System Considering Equipment Inertia
CN113762643A (zh) 区域综合能源系统的储能容量优化配置方法

Legal Events

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