CN101825867A - 设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 - Google Patents
设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 Download PDFInfo
- Publication number
- CN101825867A CN101825867A CN200910004575A CN200910004575A CN101825867A CN 101825867 A CN101825867 A CN 101825867A CN 200910004575 A CN200910004575 A CN 200910004575A CN 200910004575 A CN200910004575 A CN 200910004575A CN 101825867 A CN101825867 A CN 101825867A
- Authority
- CN
- China
- Prior art keywords
- optimal operation
- operation point
- steam
- input
- equipment
- 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
Links
Images
Landscapes
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
提供一种能够考虑自家用发电设备的蒸气配管网处的蒸气流从而使运转成本最小化的设备最优化运行***。本发明的设备最优化运行***具备用于计算自家用发电设备(10)的最优运转点的最优运转点计算部(3)以及根据最优运转点对设备进行自动控制的自动控制装置(9),最优运转点计算部(3)将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过在内部存储的规定的运算式来进行运算,从而找出最优运转点,并将该最优运转点作为目标值输出至自动控制装置。
Description
技术领域
本发明涉及自家用发电设备的最优化运行***,以及用于该***中的最优运转点计算方法和最优运转点计算程序。
背景技术
一般来讲,工厂内为了使供给电和蒸气的自家用发电设备在最优运转点运转而进行运转控制,在进行该运转控制的设备最优化运行***中,为了使从锅炉燃料或电力***接受的电量的运转成本最小化,确定锅炉输出功率、涡轮主蒸气量、抽气量、发电机输出功率,进行以这些值为目标值的控制。
这样的设备最优化运行***具有:自动控制装置,对设备的蒸气压力、流量和发电机输出进行控制;最优运转点计算装置,利用经由该自动控制装置而输入的设备的状态量、设备的最优化模型及其参数,通过规定的计算式来进行计算,从而计算出使设备的运转成本最小化的最优运转点,并作为上述的自动控制装置的设定值而输出该最优运转点(例如参照专利文献1、专利文献2)。
而且,这样的设备最优化运行***中的最优运转点是通过解析下式的最优化模型来计算出的。
【式1】
目标函数:
∑cj*FFj(燃料成本)…(1)
约束条件:
∑Ei=EL(电力供求平衡)…(2)
∑Sj-∑SHi=SHL(高压蒸气供求平衡)…(3)
∑SMi=SML(中压蒸气供求平衡)…(4)
∑SLi=SLL(低压蒸气供求平衡)…(5)
Sj=fj(FFj),j=1,…,n(锅炉的输入输出特性)…(6)
SHi=gi(Ei,SMi,SLi),i=1,…,m
(涡轮和发电机的输入输出特性)…(7)
(各变量的上下阈值以及等式约束)…(8)
在此,以下是最优化变量。
【式2】
FFj:锅炉j的燃料流量
Sj:锅炉j的蒸气输出功率
SHi:涡轮发电机i的主蒸气流量
SMi:涡轮发电机i的中压蒸气的抽气流量
SLi:涡轮发电机i的低压蒸气的抽气流量
Ei:涡轮发电机i的发电机输出功率
以下是参数。
【式3】
cj:锅炉j的燃料成本
fj:锅炉j的输入输出特性函数
gi:涡轮发电机i的输入输出特性函数
EL,SHL,SML,SLL:电力,高压蒸气,中压蒸气,低压蒸气的各负荷
在这样的以往的设备最优化运行***中,作为在工厂内的生产中所使用的蒸气的约束条件,如(3)~(5)式所示,仅考虑了蒸气的供求平衡。然而,在实际的设备中,由于存在用于向工厂内供给生产用蒸气的配管网,因此,在设备最优化运行***中计算出的各发电机的主蒸气量和抽气量实际上是否能够输出以及实际上是否能够供给蒸气负荷都是未知的。
专利文献1:日本特开2000-78749号公报
专利文献2:日本特开2004-190620号公报
发明内容
本发明是鉴于上述的以往技术中的问题而作出的,目的在于提供一种能够考虑在蒸气配管网处的蒸气流(蒸気潮流)使运转成本最小化的设备最优化运行***以及用于该***中的最优运转点计算方法、程序。
本发明的一个技术特征在于,提供一种设备最优化运行***,其特征在于,具备:最优运转点计算单元,用于计算自家用发电设备的最优运转点;以及自动控制单元,根据最优运转点对设备进行自动控制,上述最优运转点计算单元将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点,并将该最优运转点作为目标值输出至上述自动控制单元。
本发明的其他的特征在于,提供一种最优运转点计算方法,用于计算自家用发电设备的最优运转点,其特征在于,将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点。
本发明的另一个其他的特征在于,提供一种最优运转点计算程序,利用计算机来计算自家用发电设备的最优运转点,其特征在于,将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电健的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点。
发明效果
根据本发明,提供一种能够考虑在蒸气配管网处的蒸气流以实现运转成本最小化的设备运转的设备最优化运行***以及用于该***中的最优运转点计算方法、程序。
附图说明
图1是本发明的一个实施方式的设备最优化运行***的框图。
图2是表示通过上述实施方式的设备最优化运行***来进行运转控制的自家用发电设备的高压蒸气配管网的图。
图3是说明上述实施方式的设备最优化运行***中的设备最优运转控制***的最优运转点计算部所使用的水头损失式的图。
图4是说明上述实施方式的设备最优化运行***中的设备最优运转控制***的最优运转点计算部所使用的节点方程式的图。
图5是表示上述实施方式的设备最优化运行***中的设备最优运转控制***的处理的流程图。
具体实施方式
以下,参照附图对本发明的实施方式进行详细说明。图1示出了本发明的一个实施方式的设备最优化运行***的结构。该设备最优化运行***包括设备最优运转控制***1以及基于最优运转控制***1来进行运转控制的自家用发电设备10。而且设备最优运转控制***1具备数据输入部2、最优运转点计算部3、设备模型存储部4、数据输出部5、自动控制装置9。
数据输入部2用于通过自动控制装置9输入自家用发电设备10的电力负荷和/或蒸气负荷等的状态量,所输入的状态量被发送至最优运转点计算部3。设备模型存储部4用于存储包含作为规定的计算式的蒸气流方程式的设备的最优化模型及其参数,这些模型被发送至最优运转点计算部3。
最优运转点计算部3根据从数据输入部2发送过来的设备的状态量和从设备模型存储部4发送过来的最优化模型,利用后述的数理规划法等的最优化手法,计算使设备的运转成本变为最小的运转点,以此作为最优运转点,并将计算出的最优运转点向自动控制装置9发送,并且通过数据输出部5将该最优运转点显示在显示器(未图示)的画面上,或者通过打印机(未图示)打印出来。
自动控制装置9进行如下控制:读入自家用发电设备10的状态量并发送至数据输入部2,此外,接收在最优运转点计算部3中计算出的最优运转点,根据该最优运转点来跟踪自家用发电设备10的该机器的最优运转点的目标值。
自家用发电设备10包括锅炉11、12,涡轮(T1、T2)13、14,发电机(G1、G2)15、16以及连接它们的蒸气配管和电流母线等,向工厂内供给电力负荷17、高压蒸气负荷18、中压蒸气负荷19、低压蒸气负荷20。
图2示出了工厂内的高压蒸气配管网30。锅炉11、12所供给的高压蒸气除了被涡轮13、14消耗外,工厂内的高压蒸气负荷18在实际的工厂内还通过配管网而作为高压蒸气负荷(1)21、高压蒸气负荷(2)22被消耗。
接着,对在最优运转点计算部3中被计算的最优化模型进行说明。自家用发电设备10的结构如上所述。如图2所述,将来自锅炉11、12的供给点、向涡轮13、14的供给点、分支点和负荷供给点分别作为节点表示为○1~○7,将节点之间作为管路表示为△1~△6。高压蒸气配管网30作为7节点6管路的模型而被示出。图3是说明水头损失式的图,图4是说明节点方程式的图。
以下,使用图2~图4来说明蒸气流方程式。在图3中,假设节点○i以及○j的蒸气压力分别是Pi、Pj,管路△k的蒸气流量是Fk,则管路的水头损失式可以表示如下。
【式4】
Pi-Pj=Rk*Fk 2…(9)
在此,Rk是表示管摩擦的系数。
此外,在图4中,从与节点○k连接的管路流出的蒸气流量Fi的和与向节点○k供给的蒸气供给量Qk相等,因此节点方程式可以表示如下。
【式5】
∑Fi=Qk…(10)
因此,图2的高压蒸气配管网30的蒸气流方程式可以表示如下。
【式6】
即,[F2]=YCTP…(11)
【式7】
即,CF=Q…(12)
在上述的具有蒸气配管网的自家用发电设备中,将锅炉以及涡轮发电机的输入输出特性、电力以及蒸气的供求平衡以及蒸气的潮流方程式作为约束条件,通过解析下式的最优化模型来计算可使锅炉的燃料成本最小化的锅炉以及涡轮发电机的运转点。
【式8】
目标函数:c1*FF1+C2*FF2(燃料成本)…(13)
约束条件:
E1+E2=EL(电力供求平衡)…(14)
[F2]=Y*CT*P(高压蒸气水头损失式)…(15)
Q=C*F(高压蒸气节点方程式)…(16)
SM1+SM2=SML(中压蒸气供求平衡)…(17)
SL1+SL2=SLL(低压蒸气供求平衡)…(18)
S1=f1(FF1)(锅炉11的输入输出特性)…(19)
S2=f2(FF2)(锅炉12的输入输出特性)…(20)
SH1=g1(E1,SM1,SL1)(涡轮13和发电机15的输入输出特性)
…(21)
SH2=g2(E2,SM2,SL2)(涡轮14和发电机16的输入输出特性)
…(22)
(各变量的上下阈值以及等式约束)…(23)
在此,下面是最优化变量。
【式9】
FF1,FF2:锅炉11,12的燃料流量
S1,S2:锅炉11,12的蒸气输出功率
SH1,SH2:涡轮13,14的主蒸气流量
SM1,SM2:涡轮13,14的中压蒸气的抽气流量
SL1,SL2:涡轮13,14的低压蒸气的抽气流量
E1,E2:发电机15,16的发电机输出功率
P=[P1 P2 P3 P4 P5 P6 P7]T:高压蒸气配管的节点压力矢量
F=[F1 F2 F3 F4 F5 F6]T:高压蒸气配管的管路流量矢量
[F2]=[F1 2 F2 2 F3 2 F4 2 F5 2 F6 2]T:将F的各要素的平方作为要素的矢量
Q=[-SH1 S1 0 S2 -SHL1 -SH2 -SHL2]T:节点流入矢量
下面是参数。
【式10】
c1,c2:锅炉11,12的燃料成本
f1,f2:锅炉11,12的输入输出特性函数
g1,g2:涡轮13,14和发电机15,16的输入输出特性函数
EL,SML,SLL:电力负荷17,中压蒸气负荷19,低压蒸气负荷20
SHL=SH1+SH2:高压蒸气负荷18=高压蒸气负荷21+高压蒸气负荷22
C:高压蒸气配管网的节点及管路连接矩阵
Y=diag[1/Rk]:将表示管摩擦的系数的倒数作为对角要素的矩阵
在这些参数中,电力负荷EL、高压蒸气负荷SH1、SH2、中压蒸气负荷SML,低压蒸气负荷SLL是作为状态量而从数据输入部2被供给的,其他的参数是从设备模型存储部4被读入的。
从(13)式到(23)式所表示的最优化问题是非线性规划问题,可以适用非线性规划法作为最优化手法。尤其是,可以将各机器的输入输出特性(19)式~(22)式近似表示为二次函数,以此来适用二次约束二次规划法。此外,可以通过将(19)式~(22)式逐次进行线性近似(逐次П形近似),以此来适用逐次线性规划法。
接着,对基于本实施方式中的设备最优运转控制***1的自家用发电设备10的最优运转控制动作、用于该***的最优运转方法以及为了实施该最优运转方法而存入该***中的最优运转程序进行说明。图5是表示本实施方式的设备最优运转控制***1的处理的流程图。
设备最优运转控制***1首先由数据输入部2通过自动控制装置9输入设备10的状态量,此外,通过从设备模型存储部4读入数据,从而准备在最优运转点计算部3中被处理的数据(S201)。
然后,最优运转点计算部3编制(13)式~(23)式的蒸气流方程式(S202)。然后,使用该蒸气流方程式进行最优化计算,计算最优运转点(S203)。计算出的最优运转点通过数据输出部5而显示在显示器的画面等上(S204),此外,将最优运转点向自动控制装置9输出(S205)。
自动控制装置9将来自最优运转点计算部3的最优运转点作为控制量,对燃料流量进行控制,并控制主蒸气流量、中压蒸气的抽气流量、低压蒸气的抽气流量,此外,通过控制发电机输出功率,持续进行最优运转。
由此,根据本实施方式的设备最优化运行***,能够考虑蒸气配管网处的蒸气流使自家用发电设备10实现运转成本最小化的运转,而且获得的设备的运转点能够可靠地进行运转,还能够可靠地提供蒸气负荷。
Claims (7)
1.一种设备最优化运行***,其特征在于,具备:
最优运转点计算单元,用于计算自家用发电设备的最优运转点;以及
自动控制单元,根据最优运转点对设备进行自动控制,
上述最优运转点计算单元将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点,并将该最优运转点作为目标值输出至上述自动控制单元。
2.如权利要求1所述的设备最优化运行***,其特征在于,
上述蒸气流方程式包含蒸气的压力和流量的关系式。
3.如权利要求1或2所述的设备最优化运行***,其特征在于,
上述最优运转点计算单元能够适用非线性规划法。
4.如权利要求3所述的设备最优化运行***,其特征在于,
上述非线性规划法是二次约束二次规划法。
5.如权利要求3所述的设备最优化运行***,其特征在于,
上述非线性规划法是逐次线性规划法。
6.一种最优运转点计算方法,用于计算自家用发电设备的最优运转点,其特征在于,
将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点。
7.一种最优运转点计算程序,利用计算机来计算自家用发电设备的最优运转点,其特征在于,
将使燃料成本最小化作为目标函数,将电力供求平衡、蒸气流方程式、蒸气供求平衡、锅炉的输入输出特性、涡轮发电机的输入输出特性、变量的上下阈值作为约束条件,通过内部存储的规定的运算式来进行运算,找出最优运转点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910004575A CN101825867A (zh) | 2009-03-06 | 2009-03-06 | 设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910004575A CN101825867A (zh) | 2009-03-06 | 2009-03-06 | 设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101825867A true CN101825867A (zh) | 2010-09-08 |
Family
ID=42689822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910004575A Pending CN101825867A (zh) | 2009-03-06 | 2009-03-06 | 设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101825867A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107735735A (zh) * | 2015-03-13 | 2018-02-23 | 株式会社东芝 | 设备特性模型学习装置、设备特性模型学习方法及存储介质 |
CN107924162A (zh) * | 2015-07-29 | 2018-04-17 | 株式会社东芝 | 最优控制装置、最优控制方法、计算机程序以及最优控制*** |
CN112032800A (zh) * | 2020-09-09 | 2020-12-04 | 山东鲁航智能科技有限公司 | 一种智能化管网平衡***及其调控方法 |
-
2009
- 2009-03-06 CN CN200910004575A patent/CN101825867A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107735735A (zh) * | 2015-03-13 | 2018-02-23 | 株式会社东芝 | 设备特性模型学习装置、设备特性模型学习方法及存储介质 |
CN107924162A (zh) * | 2015-07-29 | 2018-04-17 | 株式会社东芝 | 最优控制装置、最优控制方法、计算机程序以及最优控制*** |
CN112032800A (zh) * | 2020-09-09 | 2020-12-04 | 山东鲁航智能科技有限公司 | 一种智能化管网平衡***及其调控方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Multi-objective optimization for advanced superheater steam temperature control in a 300 MW power plant | |
CN110571789B (zh) | 数据驱动下基于风电不确定性的电热气网三阶段调度方法 | |
US4607325A (en) | Discontinuous optimization procedure modelling the run-idle status of plural process components | |
CN103180666B (zh) | 发电设备的凝结水流量控制装置以及控制方法 | |
CN110503250B (zh) | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 | |
Snarheim et al. | Control design for Gas Turbine Cycle with CO2 Capture Capabilities | |
Ramírez-Elizondo et al. | Scheduling and control framework for distribution-level systems containing multiple energy carrier systems: Theoretical approach and illustrative example | |
US20160169154A1 (en) | Model-based combined cycle power plant load control | |
Chen et al. | Steady-state security regions of electricity-gas integrated energy systems | |
CN112051810A (zh) | 一种供热机组的多能源分配和调度*** | |
CN101825867A (zh) | 设备最优化运行***,最优运转点计算方法以及最优运转点计算程序 | |
CN106446375B (zh) | 一种基于数据驱动的单元机组锅炉汽轮机控制方法及装置 | |
CA1170778A (en) | Method of minimizing the cost of operation of a process | |
Edlund et al. | Simple models for model-based portfolio load balancing controller synthesis | |
Binder et al. | Hybrid energy system modeling in Modelica | |
JP2009146145A (ja) | プラント最適運用システム、最適運転点計算方法及び最適運転点計算プログラム | |
CN107147139B (zh) | 一种计及线路传输损耗的微电网负荷调度方法 | |
CN115187129A (zh) | 考虑热网质量流量调节的微能源网优化调度方法及*** | |
CN114139349A (zh) | 计及管道爬坡的气电联合***日前机组组合优化方法 | |
CN114936440A (zh) | 多时间尺度下多能源耦合***联立潮流仿真方法和*** | |
Skworcow et al. | Model predictive control for energy and leakage management in water distribution systems | |
CN113592149A (zh) | 一种耦合型新能源综合能源配电网优化重构方法 | |
JP2008280972A (ja) | 蒸気系統の最適運用方法 | |
CN106875049A (zh) | 一种计及阀点效应的电力***动态经济调度方法 | |
CN115080907A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20100908 |