CN104123593B - 基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 - Google Patents
基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 Download PDFInfo
- Publication number
- CN104123593B CN104123593B CN201410339431.0A CN201410339431A CN104123593B CN 104123593 B CN104123593 B CN 104123593B CN 201410339431 A CN201410339431 A CN 201410339431A CN 104123593 B CN104123593 B CN 104123593B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- load
- mtd
- mtr
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,所述步骤包括:步骤1、从厂内各台机组的DCS的实时数据库里读取给定时刻的实时运行数据,包括工质侧、烟气侧全流程的温度、压力、流量等测点;步骤2、结合工质物性参数库、烟气物性参数库,以及全流程能量衡算模型,在线更新每台机组的煤耗特性曲线;步骤3、在当前负荷调度周期内,设置中调总负荷指令及负荷调度模式,实时优化算法根据各台机组的煤耗特性曲线进行机组间负荷分配。本发明具有成本低、计算速度快、适应性强的特点,对燃煤电站的优化运行、节能减排具有重要现实意义。
Description
技术领域
本发明涉及火力发电控制领域的调度方法,具体地,涉及一种基于煤耗特性曲线在线滚动更新的多模式负荷调度方法。
背景技术
自厂网分离政策实施以来,电网以竞价招标的方式向电厂购电,这就对电厂并网供电提出了更高的要求,主要表现在电价、响应速度、稳定性等方面。为了在满足电网各项要求的同时实现经济效益最大化,以提高自身的竞争力,电厂对机组的运行优化与控制有着迫切的需求。厂内不同机组的设备状态、煤种、运行水平等都存在着一定的差异,造成各台机组在同一负荷出力下的煤耗量并不相同,这一总体效率通常用负荷-煤耗特性曲线表征。显然,在不同机组间进行合理的负荷分配,可以使电厂在满足电网中调总负荷指令的同时,尽可能地发挥各台机组的最优性能,从而达到降低总煤耗的目的。
经过对现有技术的检索,中国专利申请号201310194342.7,公开日2013-9-25,记载了一种基于改进多目标粒子群算法的环境经济发电调度方法,以燃料费用最低和污染气体排放量最少为调度目标,该方法采用了多目标粒子群算法实现环境与经济多目标调度,但它采用了固定的机组煤耗特性、污染物排放特性,不能动态反映机组特性的时变性;同时粒子群算法相对复杂,工程上不易实现。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于煤耗特性曲线在线滚动更新的多模式负荷调度方法。该发明利用机组的实时运行数据,得到平稳负荷下的平均煤耗点,进而滚动更新各机组的煤耗特性曲线,并在此基础上实现多模式下的负荷调度,对我国燃煤电站的优化运行、节能减排具有重要现实意义。
为实现上述目的,本发明所述方法具体包括以下步骤:
步骤1、从厂内各台机组的分布式控制系统(DCS)的实时数据库里读取给定时刻的运行工况实时数据,包括工质侧、烟气侧全流程的温度、压力、流量以及机组功率测点。
步骤2、结合工质物性参数库、烟气物性参数库,以及全流程能量衡算模型,在线更新每台机组的煤耗特性曲线;
所述的全流程能量衡算模型包括计算工质在水冷壁、各级过热器/再热器、省煤器中的吸热量,金属壁蓄放热量,以及锅炉热损失;
所述的煤耗特性曲线表征机组负荷与标准煤耗量的关系,采用多项式回归;
优选地,煤耗特性曲线采用最小二乘方法对负荷与煤耗量的关系进行二次多项式回归。
所述的在线更新机组的煤耗特性曲线的方法,具体为:
建立平均煤耗矩阵A1,…,A10,每个矩阵维数均为20行2列。将负荷区间[Pmin,Pmax]分成等长的10段,第i段对应矩阵Ai。矩阵Ai存放对应负荷范围内的20组平均煤耗点(P,C)。Ai矩阵依次叠放组合构成200行2列的矩阵A。其中Pmin,Pmax为机组允许的最小和最大负荷。
建立瞬时煤耗矩阵B,行数可变,列数为2。用于记录在线运行时某一稳定负荷下的瞬时煤耗点(P,C)。瞬时煤耗点的均值即为平均煤耗点。
在每一个时刻,首先判断当前时刻负荷P是否平稳:
若负荷平稳,则根据全流程能量衡算模型,计算当前时刻的入炉煤总能量输出Q,进而得到标煤耗量C(标煤低位发热量29.3MJ/kg),并存入矩阵B。
否则,若上一时刻负荷平稳,则说明负荷刚由平稳负荷状态切换至变负荷状态,对矩阵B中的瞬时煤耗点中的元素P、C分别求均值得到平均煤耗点,并将该新的平均煤耗点根据P的范围滚动存入对应的矩阵Ai,然后清空矩阵B,对A矩阵进行回归得到新的煤耗特性曲线。
若上一时刻负荷不平稳,说明当前时刻仍处于变负荷状态,则不作任何修改。
所述的判断负荷平稳的方法,具体为:
用flag标记负荷状态,flag=0表示平稳负荷,flag=1表示变负荷;同时定义平稳负荷计数器和变负荷计数器;在当前时刻,滚动更新由前N个时刻的负荷数据组成的负荷时间序列P(1),P(2),P(3),…,P(N)。然后分别计算该向量的以下指标:平均斜率、极差、方差。其中,平均斜率是指(N-1)个斜率的算术平均。(N-1)个斜率按如下方法计算:[P(N)-P(N-1)]/Δt,[P(N)-P(N-2)]/(2Δt),[P(N)-P(N-3)]/(3Δt),…,[P(N)-P(1)]/[(N-1)Δt]。
判断负荷状态的条件:
(a)平均斜率绝对值大于阈值TA;(b)极差大于阈值TB;(c)方差大于阈值TC。
负荷状态判断步骤:
根据上一个采样时刻的负荷状态和当前时刻的判断负荷状态条件,判断当前时刻的负荷状态:
若上一时刻为稳定负荷状态(flag=0),则只要(a)、(b)、(c)三个条件中至少一个得到满足,变负荷计数器记数增加1,否则变负荷计数器清零且flag置0;如果变负荷计数器中的数值超过一定的阈值MB,则判断负荷已处于变负荷状态,将flag置1并将稳定负荷计数器清零;
若上一时刻为变负荷状态(flag=1),则当(a)、(b)、(c)三个条件均不满足时,稳定负荷计数器记数增加1,否则稳定负荷计数器清零且flag置1;如果稳定负荷计数器中的数值超过一定的阈值MW,则判断负荷已处于稳定负荷状态,将flag置0并将变负荷计数器清零。
阈值TA,TB,TC分别根据历史运行数据中的稳定负荷与变负荷状态下的平均斜率绝对值、极差、方差的统计值确定,MB、MW根据变负荷切换时间的历史统计特征加以确定。
步骤3、在当前的负荷调度周期内,设置中调总负荷指令及负荷调度模式,实时优化算法根据各台机组的煤耗特性曲线进行机组间负荷分配。
所述的多模式厂级负荷调度方法,包括简单经济负荷调度、允许启停经济负荷调度、最速响应负荷调度、经济与快速多目标负荷调度4种模式;
所述4种负荷调度模式根据运行需要选择其中之一;
所述实时优化算法可以采用非线性优化方法,如单纯形算法,或启发式算法,如模拟退火算法;
所述4种负荷调度模式下的最优化问题为:
(1)简单经济负荷调度
上式中,Fi为第i台机组的煤耗量;F为各台机组总煤耗量;Pi为分配到第i台机组的负荷;n为厂内参与调度的机组总数;fi为第i台机组的煤耗特性曲线;P为中调总负荷指令;Pimin、Pimax为第i台机组允许的最小和最大负荷。
(2)允许启停经济负荷调度
其中,
上式中,Ui为第i台机组的启(以1表示)停(以0表示)状态;Ui0为上一个调度周期第i台机组的启停状态;Si为机组启动或停机耗费折算标准煤耗;Ai为机组停机耗煤;Bi为机组启动耗煤;ΔT为调度周期;Tiqt为第i台机组的启停时间。
(3)最速响应负荷调度
minT=min{max(Ti(Pi))}
上式中,T为机组总负荷达到中调负荷指令的过渡用时;Ti为第i台机组达到负荷Pi的过渡用时;Vi_up、Vi_down为第i台机组升、降负荷速率;Vimax_up、Vimax_down为第i台机组升、降负荷速率的最大值。
(4)经济与快速多目标负荷调度
上式中,G为经济与快速综合指标;Fmin为简单经济负荷调度的煤耗量;Tmin为最速响应调度的煤耗量;α为经济与快速多目标权重系数;W为归一化因子,通过离线实验确定。
与现有技术相比,本发明具有以下有益效果:
本发明实现了厂级负荷的多模式调度,充分利用了机组运行的DCS数据,滚动更新各台机组的煤耗特性,可以及时反映机组的特性变化,同时多模式的负荷调度方法为实际运行提供很大的选择性和灵活性。该发明具有成本低、计算速度快、适应性强的特点,对燃煤电站的优化运行、节能减排具有重要现实意义。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1本发明一实施例多模式负荷调度流程示意图;
图2本发明一实施例在线更新煤耗特性曲线示意图;
图3本发明一实施例经济与快速多目标负荷调度模式不同α的煤耗与负荷过渡用时。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本实施例涉及某电厂4台328.5MW亚临界机组,提供一种基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,见图1,具体包括以下步骤:
步骤1、从厂内各台机组的DCS的实时数据库里读取给定时刻的运行工况实时数据,包括工质侧、烟气侧全流程的温度、压力、流量以及机组功率测点。
步骤2、结合工质物性参数库、烟气物性参数库,以及全流程能量衡算模型,在线更新每台机组的煤耗特性曲线;
所述工质物性参数库,是指根据水和水蒸气热力性质工业公式(IAPWS-IF97)开发的具有可并行调用的、区域自动判别、批处理运算等特点的用于在线计算的工质物性参数库,可参考文献:王旭辉,于彤,惠兆宇,袁景淇,用于火电全范围仿真的工质物性参数数据库,控制工程,2011;18:131-133。
所述烟气物性参数库,是指通过烟气压力和温度实时数据在线计算空气的比热和密度的物性数据库。可参考文献:蔡惟,于彤,惠兆宇,袁景淇,张锐峰,陈宇,火电锅炉排烟热损失的在线估计,控制工程,2011;18:149-151。
所述的全流程能量衡算模型包括计算工质在水冷壁、各级过热器/再热器、省煤器中的吸热量,金属壁蓄放热量,以及锅炉热损失;此部分可以采用发明专利《一种燃煤电站入炉煤低位发热量的实时辨识方法(专利申请号201310697798.5)》中的技术实现。
所述的煤耗特性曲线表征机组负荷与标准煤耗量之间的关系,采用最小二乘方法对负荷与煤耗量的关系进行二次多项式回归。
所述的在线更新机组的煤耗特性曲线的方法,流程图见图2,具体为:
建立平均煤耗矩阵A1,…,A10,每个矩阵维数均为20行2列。将负荷区间[160,328.5]分成等长的10段,第i段对应矩阵Ai。矩阵Ai存放对应负荷范围内的20组平均煤耗点(P,C)。Ai矩阵依次叠放组合构成200行2列的矩阵A。
建立瞬时煤耗矩阵B,行数可变,列数为2。用于记录在线运行时某一稳定负荷下的瞬时煤耗点(P,C)。瞬时煤耗点的均值即为平均煤耗点。
在每一个时刻,首先判断当前时刻负荷P是否平稳:
若负荷平稳,则根据全流程能量衡算模型,计算当前时刻的入炉煤总能量输出Q,进而得到标煤耗量C(标煤低位发热量29.3MJ/kg),并存入矩阵B。
否则,若上一时刻负荷平稳,则说明负荷刚由平稳负荷状态切换至变负荷状态,对矩阵B中的瞬时煤耗点中的元素P、C分别求均值得到平均煤耗点,并将该新的平均煤耗点根据P的范围滚动存入对应的矩阵Ai,然后清空矩阵B,对A矩阵进行回归得到新的煤耗特性曲线。
若上一时刻负荷不平稳,说明当前时刻仍处于变负荷状态,则不作任何修改。
所述的判断负荷平稳的方法,具体为:
用flag标记负荷状态,flag=0表示平稳负荷,flag=1表示变负荷;同时定义稳定负荷计数器和变负荷计数器;在当前时刻,滚动更新由前N个时刻的负荷数据组成的负荷时间序列P(1),P(2),P(3),…,P(N)。然后分别计算该向量的以下指标:平均斜率、极差、方差。其中,平均斜率是指(N-1)个斜率的算术平均。(N-1)个斜率按如下方法计算:[P(N)-P(N-1)]/Δt,[P(N)-P(N-2)]/(2Δt),[P(N)-P(N-3)]/(3Δt),…,[P(N)-P(1)]/[(N-1)Δt]。
判断负荷状态的条件:
(a)平均斜率绝对值大于阈值TA;(b)极差大于阈值TB;(c)方差大于阈值TC。
所述负荷状态判断步骤:
根据上一个采样时刻的负荷状态和当前时刻的判断负荷状态条件,判断当前时刻的负荷状态:
若上一时刻为稳定负荷状态(flag=0),则只要(a)、(b)、(c)三个条件中至少一个得到满足,变负荷计数器记数增加1,否则变负荷计数器清零且flag置0;如果变负荷计数器中的数值超过一定的阈值MB,则判断负荷已处于变负荷状态,将flag置1并将稳定负荷计数器清零;
若上一时刻为变负荷状态(flag=1),则当(a)、(b)、(c)三个条件均不满足时,稳定负荷计数器记数增加1,否则稳定负荷计数器清零且flag置1;如果稳定负荷计数器中的数值超过一定的阈值MW,则判断负荷已处于稳定负荷状态,将flag置0并将变负荷计数器清零。
阈值TA,TB,TC分别根据历史运行数据中的稳定负荷与变负荷状态下的平均斜率绝对值、极差、方差的统计值确定,MB、MW根据历史运行数据中变负荷切换时间统计值进行确定。
本实施例中,N=24,Δt=5,TA=0.75,TB=3,TC=0.6,MB=12,MW=120。当然,在其他实施例中也可以根据实际需要采用其他数值。
步骤3、设置中调总负荷指令及负荷调度模式,实时优化算法根据各台机组的煤耗特性曲线进行机组间负荷分配。
所述实时优化算法采用单纯形算法;
所述的多模式厂级负荷调度方法,包括简单经济负荷调度、允许启停经济负荷调度、最速响应负荷调度、经济与快速多目标负荷调度4种模式;
所述4种负荷调度模式由运行人员根据运行需要选择其中之一;
设在某调度周期时刻,四台机组的煤耗特性曲线为:
F1=5.8×10-4P2+4.89×10-3P+67.81
F2=4.3×10-4P2+1.95×10-1P+40.62
F3=4.1×10-4P2+2.29×10-1P+47.31
F4=1.0×10-3P2-6.56×10-2P+84.02
且Pimin=160MW,Pimax=328.5MW,Vimax_up=5MW/min、Vimax_down=3MW/min.
所述4种负荷调度模式下的最优化问题为:
(1)简单经济负荷调度
上式中,Fi为第i台机组的煤耗量;F为各台机组总煤耗量;Pi为分配到第i台机组的负荷;n为厂内参与调度的机组总数;fi为第i台机组的煤耗特性曲线;P为中调总负荷指令;Pimin、Pimax为第i台机组允许的最小和最大负荷。
表1简单经济负荷调度结果
*用Pi组成的4维向量表示分配结果,下同
(2)允许启停经济负荷调度
其中,
上式中,Ui为第i台机组的启(以1表示)停(以0表示)状态;Ui0为上一个调度周期第i台机组的启停状态;Si为机组启动或停机耗费折算标准煤耗;Ai为机组停机耗煤;Bi为机组启动耗煤;ΔT为调度周期,取2h;Tiqt为第i台机组的启停时间,取0.5h。
表2允许启停的经济负荷调度结果*
*机组启动折算耗煤71.4t,机组停机折算耗煤42.8t
(3)最速响应负荷调度
minT=min{max(Ti(Pi))}
上式中,T为机组总负荷达到中调负荷指令的过渡用时;Ti为第i台机组达到负荷Pi的过渡用时;Vi_up、Vi_down为第i台机组升、降负荷速率;Vimax_up、Vimax_down为第i台机组升、降负荷速率的最大值。
表3最速响应的负荷调度结果*
*上一个调度周期的总负荷指令为P0=800MW,各机组负荷分配为[285,160,160,195],本调度周期的总负荷指令P=1000MW
(4)经济与快速多目标负荷调度
上式中,G为经济与快速综合指标;Fmin为简单经济负荷调度的煤耗量;Tmin为最速响应调度的煤耗量;α为经济与快速多目标权重系数;W为归一化因子,通过离线实验确定。
表4经济与快速多目标负荷调度结果*
*上一个调度周期的总负荷指令为P0=800MW,各机组负荷分配为[285,160,160,195],本调度周期的总负荷指令P=1000MW,W=500
本实施例中,选择α=0.5,则调度结果为[328.5,223.6,216.2,231.6]。
本发明利用了机组运行的DCS数据,滚动更新各台机组的煤耗特性,及时反映机组的特性变化,同时多模式的负荷调度方法为实际运行提供很大的选择性和灵活性。该发明具有成本低、计算速度快、适应性强的特点,对燃煤电站的优化运行、节能减排具有重要现实意义。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (4)
1.一种基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,其特征在于,该方法具体包括以下步骤:
步骤1、从厂内各台机组的分布式控制系统DCS的实时数据库里读取给定时刻的实时运行数据,包括工质侧、烟气侧全流程的温度、压力、流量以及机组功率测点;
步骤2、结合工质物性参数库、烟气物性参数库,以及全流程能量衡算模型,在线更新每台机组的煤耗特性曲线;
所述工质物性参数库,是指根据水和水蒸气热力性质工业公式IAPWS-IF97开发的具有可并行调用的、区域自动判别、批处理运算特点的用于在线计算的工质物性参数库;
所述烟气物性参数库,是指通过烟气压力和温度实时数据在线计算空气的比热和密度的物性数据库;
步骤2所述的全流程能量衡算模型包括计算工质在水冷壁、各级过热器/再热器、省煤器中的吸热量,金属壁蓄放热量,以及锅炉热损失;所述的煤耗特性曲线表征机组负荷与标准煤耗量的关系,采用多项式回归;
步骤2所述的在线更新机组的煤耗特性曲线,具体为:
建立平均煤耗矩阵A1,…,A10,每个矩阵维数均为20行2列,将负荷区间[Pmin,Pmax]分成等长的10段,第i段对应矩阵Ai,矩阵Ai存放对应负荷范围内的20组平均煤耗点(P,C),Ai矩阵依次叠放组合构成200行2列的矩阵A,其中Pmin,Pmax为机组允许的最小和最大负荷;
建立瞬时煤耗矩阵B,行数可变,列数为2,用于记录在线运行时某一稳定负荷下的瞬时煤耗点(P,C),瞬时煤耗点的均值即为平均煤耗点;
在每一时刻,首先判断当前时刻负荷P是否平稳,负荷状态分为平稳负荷和变负荷:
若负荷平稳,则根据全流程能量衡算模型,计算当前时刻的入炉煤总能量输出Q,进而得到标煤耗量C,并存入矩阵B;
否则,若上一时刻负荷平稳,说明负荷刚由平稳负荷状态切换至变负荷状态,则对矩阵B中的瞬时煤耗点中的元素P、C分别求均值得到平均煤耗点,并将该新的平均煤耗点根据P的范围滚动存入对应的矩阵Ai,然后清空矩阵B,对A矩阵进行回归得到新的煤耗特性曲线;
若上一时刻负荷不平稳,说明当前时刻仍处于变负荷状态,则不作任何修改;
步骤3、在当前负荷调度周期内,设置中调总负荷指令及负荷调度模式,实时优化算法根据各台机组的煤耗特性曲线进行机组间负荷分配。
2.根据权利要求1所述的基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,其特征在于,所述的负荷平稳,判断方法具体为:
用flag标记负荷状态,flag=0表示平稳负荷,flag=1表示变负荷;同时定义平稳负荷计数器和变负荷计数器;在当前时刻,滚动更新由前N个时刻的负荷数据组成的负荷时间序列P(1),P(2),P(3),…,P(N),然后分别计算向量的以下指标:平均斜率、极差、方差,其中,平均斜率是指(N-1)个斜率的算术平均,(N-1)个斜率按如下方法计算:[P(N)-P(N-1)]/Δt,[P(N)-P(N-2)]/(2Δt),[P(N)-P(N-3)]/(3Δt),…,[P(N)-P(1)]/[(N-1)Δt];
判断负荷状态的条件:
(a)平均斜率绝对值大于阈值TA;(b)极差大于阈值TB;(c)方差大于阈值TC;
所述负荷状态判断步骤:
根据上一个采样时刻的负荷状态和当前时刻的判断负荷状态条件,判断当前时刻的负荷状态:
若上一时刻为平稳负荷状态(flag=0),则只要(a)、(b)、(c)三个条件中至少一个得到满足,变负荷计数器记数增加1,否则变负荷计数器清零且flag置0;如果变负荷计数器中的数值超过一定的阈值MB,则判断负荷已处于变负荷状态,将flag置1并将平稳负荷计数器清零;
若上一时刻为变负荷状态(flag=1),则当(a)、(b)、(c)三个条件均不满足时,平稳负荷计数器记数增加1,否则平稳负荷计数器清零且flag置1;如果平稳负荷计数器中的数值超过一定的阈值MW,则判断负荷已处于平稳负荷状态,将flag置0并将变负荷计数器清零;
阈值TA,TB,TC分别根据历史运行数据中的平稳负荷与变负荷状态下的平均斜率绝对值、极差、方差的统计值确定,MB、MW根据变负荷切换时间的历史统计特征加以确定。
3.根据权利要求1-2任一项所述的基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,其特征在于,所述的负荷调度模式,包括简单经济负荷调度、允许启停经济负荷调度、最速响应负荷调度、经济与快速多目标负荷调度4种模式,所述4种负荷调度模式根据运行需要选择其中之一;所述实时优化算法采用非线性优化方法。
4.根据权利要求3所述的基于煤耗特性曲线在线滚动更新的多模式负荷调度方法,其特征在于,4种负荷调度模式下的最优化问题为:
(1)简单经济负荷调度
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>F</mi>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
上式中,Fi为第i台机组的煤耗量;F为各台机组总煤耗量;Pi为分配到第i台机组的负荷指令;n为厂内参与调度的机组总数;fi为第i台机组的煤耗特性曲线;P为中调总负荷指令;Pimin、Pimax为第i台机组允许的最小和最大负荷;
(2)允许启停经济负荷调度
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>F</mi>
<mo>=</mo>
<mi>min</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>q</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>min</mi>
<mrow>
<mo>{</mo>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>q</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,
上式中,Ui为第i台机组的启(以1表示)停(以0表示)状态;Ui0为上一个调度周期第i台机组的启停状态;Si为机组启动或停机耗费折算标准煤耗;Ai为机组停机折算耗煤;Bi为机组启动折算耗煤;ΔT为调度周期;Tiqt为第i台机组的启停时间;
(3)最速响应负荷调度
minT=min{max(Ti(Pi))}
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mo>_</mo>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msub>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>_</mo>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mo>_</mo>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msub>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>_</mo>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
上式中,T为机组总负荷达到中调负荷指令的过渡用时;Ti为第i台机组达到负荷指令Pi的过渡用时;Vi_up、Vi_down为第i台机组升、降负荷速率;Vimax_up、Vimax_down为第i台机组升、降负荷速率的最大值;
(4)经济与快速多目标负荷调度
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>G</mi>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>{</mo>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<mi>W</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>F</mi>
<mi>min</mi>
</msub>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>}</mo>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mo>_</mo>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msub>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>_</mo>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mo>_</mo>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msub>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>_</mo>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
上式中,G为经济与快速综合指标;Fmin为简单经济负荷调度的煤耗量;Tmin为最速响应调度的煤耗量;α为经济与快速多目标权重系数;W为归一化因子,通过离线实验确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410339431.0A CN104123593B (zh) | 2014-07-16 | 2014-07-16 | 基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410339431.0A CN104123593B (zh) | 2014-07-16 | 2014-07-16 | 基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104123593A CN104123593A (zh) | 2014-10-29 |
CN104123593B true CN104123593B (zh) | 2017-10-31 |
Family
ID=51768995
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410339431.0A Active CN104123593B (zh) | 2014-07-16 | 2014-07-16 | 基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104123593B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104571018A (zh) * | 2014-12-04 | 2015-04-29 | 上海交通大学 | 一种燃煤电站制粉系统在线负荷分配优化方法 |
CN104901341B (zh) * | 2015-06-26 | 2017-08-04 | 华北电力大学 | 基于机组出力及变负荷速率的火力发电煤耗量计算方法 |
CN107612045A (zh) * | 2017-10-09 | 2018-01-19 | 国网安徽省电力公司 | 一种电厂机组发电量智能分配系统及方法 |
CN107834609B (zh) * | 2017-11-06 | 2021-02-05 | 安徽立卓智能电网科技有限公司 | 一种提高机组稳定性的厂级agc系统分配策略优化方法 |
CN107895209A (zh) * | 2017-11-17 | 2018-04-10 | 上海交通大学 | 混杂燃料机组电厂负荷燃料优化分配方法及系统 |
CN112491049B (zh) * | 2020-11-26 | 2022-10-11 | 贵州电网有限责任公司 | 一种考虑在线煤耗曲线的多能源接入电网优化调度方法 |
CN112510704B (zh) * | 2020-11-26 | 2022-10-11 | 贵州电网有限责任公司 | 一种在线煤耗曲线实时生成方法及系统 |
CN112510703B (zh) * | 2020-11-26 | 2022-10-04 | 贵州电网有限责任公司 | 一种考虑煤耗曲线修正的多能源接入电网优化调度方法 |
CN112835950B (zh) * | 2020-12-09 | 2023-03-28 | 华能陕西发电有限公司 | 一种dcs数据挖掘的湿法脱硫系统达标排放运行曲线的获取系统及方法 |
CN113532894A (zh) * | 2021-06-09 | 2021-10-22 | 苏州西热节能环保技术有限公司 | 一种燃煤电厂热平衡监测方法 |
CN114398779A (zh) * | 2022-01-07 | 2022-04-26 | 华北电力科学研究院有限责任公司 | 火电机组煤耗特性曲线确定方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101206754A (zh) * | 2006-12-21 | 2008-06-25 | 北京华电天仁电力控制技术有限公司 | 基于多种约束规则的火电厂厂级负荷优化分配方法 |
CN101593979A (zh) * | 2009-07-06 | 2009-12-02 | 贵阳高新金辰软件有限公司 | 一种火力发电厂厂级负荷优化分配方法及装置 |
JP2012195990A (ja) * | 2011-03-14 | 2012-10-11 | Omron Corp | 負荷制御装置およびその制御方法、並びに制御プログラム |
CN103326353A (zh) * | 2013-05-21 | 2013-09-25 | 武汉大学 | 基于改进多目标粒子群算法的环境经济发电调度求解方法 |
-
2014
- 2014-07-16 CN CN201410339431.0A patent/CN104123593B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101206754A (zh) * | 2006-12-21 | 2008-06-25 | 北京华电天仁电力控制技术有限公司 | 基于多种约束规则的火电厂厂级负荷优化分配方法 |
CN101593979A (zh) * | 2009-07-06 | 2009-12-02 | 贵阳高新金辰软件有限公司 | 一种火力发电厂厂级负荷优化分配方法及装置 |
JP2012195990A (ja) * | 2011-03-14 | 2012-10-11 | Omron Corp | 負荷制御装置およびその制御方法、並びに制御プログラム |
CN103326353A (zh) * | 2013-05-21 | 2013-09-25 | 武汉大学 | 基于改进多目标粒子群算法的环境经济发电调度求解方法 |
Non-Patent Citations (2)
Title |
---|
对火电厂机组负荷优化分配算法的研究;曹文亮等;《汽轮机技术》;20040228;第46卷(第1期);第43-46页 * |
综合经济目标下电厂机组负荷优化调度方法研究;张强;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20101115(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104123593A (zh) | 2014-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104123593B (zh) | 基于煤耗特性曲线在线滚动更新的多模式负荷调度方法 | |
CN103151803B (zh) | 一种含风电系统机组及备用配置的优化方法 | |
Lee et al. | Determination of optimal contract capacities and optimal sizes of battery energy storage systems for time-of-use rates industrial customers | |
CN105279707B (zh) | 一种综合考虑负荷与风电时序特性的随机生产模拟方法 | |
CN108090666B (zh) | 一种含aa-caes的电网电能与备用容量的协同调度方法 | |
CN111030094B (zh) | 一种多能互补园区供能可靠性评估方法 | |
CN107292766A (zh) | 面向风电消纳的电力系统调峰手段经济性评估方法与系统 | |
CN112270433B (zh) | 考虑可再生能源不确定性和用户满意度的微电网优化方法 | |
CN104268800A (zh) | 基于场景库的风电并网调峰平衡判定方法 | |
CN106684928B (zh) | 一种基于调峰成本的电网调峰裕度计算方法 | |
CN113256045A (zh) | 考虑风光不确定性的园区综合能源系统日前经济调度方法 | |
CN107274069A (zh) | 基于可选多因素条件下火电机组的负荷特性确定方法 | |
CN115640982A (zh) | 基于抽水蓄能优先调节的多能互补系统日前优化调度方法 | |
CN108335232A (zh) | 基于备用整定的双层鲁棒机组组合方法 | |
Zhang et al. | Research on economic optimal dispatching of microgrid cluster based on improved butterfly optimization algorithm | |
CN110689156B (zh) | 一种泛能站优化方法及装置 | |
CN113300361B (zh) | 基于改进多目标法的电热联合系统风电接纳能力评估方法 | |
CN109871999A (zh) | 一种综合能源系统的随机生产模拟方法及系统 | |
CN114188942A (zh) | 一种包含大规模新能源基地的电网调度方法 | |
CN114399162A (zh) | 一种基于能源调度时间自适应变化的滚动优化调度方法 | |
Stark et al. | Design parameters of steam accumulators for the utilization in industrial solid biomass-fuelled CHP plants | |
CN113919676A (zh) | 计及需求响应与电制氢系统的虚拟电厂运行效果评估方法 | |
CN105868920A (zh) | 一种确定火电厂机组小指标竞赛指标考核权重的方法 | |
Zhang et al. | Algorithm on optimal wind power capacity using peak load regulation restraints | |
CN111555270A (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 |