CN104657789B - 一种蒸汽动力系统的运行操作方法 - Google Patents
一种蒸汽动力系统的运行操作方法 Download PDFInfo
- Publication number
- CN104657789B CN104657789B CN201510073874.4A CN201510073874A CN104657789B CN 104657789 B CN104657789 B CN 104657789B CN 201510073874 A CN201510073874 A CN 201510073874A CN 104657789 B CN104657789 B CN 104657789B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- fuel
- steam
- mtd
- 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.)
- Expired - Fee Related
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"
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种蒸汽动力系统的运行操作方法，该方法包括：采集蒸汽动力系统中各设备的性能特征参数以及蒸汽动力系统的工艺参数；根据蒸汽动力系统的能量、质量守恒方程建立所述蒸汽动力系统的混合整数线性规划模型，并设定所述混合整数线性规划模型的优化目标函数；以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得目标优化参数；将所述目标优化参数作为蒸汽动力系统的运行参数进行蒸汽动力系统的运行控制。本发明以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得了满足生产成本和排放成本两个目标的优化策略，可帮助企业在保证生产成本的同时，适应节能环保要求趋势，科学调整生产运行策略。
Description
技术领域
本发明涉及蒸汽动力系统技术领域，具体涉及一种蒸汽动力系统的运行操作方法。
背景技术
石化企业是典型的高耗能产业，如何通过改进设备和操作工艺来降低生产成本和环境成本，一直是该行业所面临的严峻问题。在石化企业中，蒸汽动力系统是供能中的重要环节，其通过消耗燃料油、燃料气及新鲜水，为整个生产过程提供不同等级蒸汽、动力乃至电力能源。
蒸汽动力系统流程结构固定的系统参数优化，主要包括现有系统操作条件的优化和新系统设计或旧系统改造时对某一结构方案的设计参数和操作条件的优化。目前，针对这类蒸汽动力系统操作参数优化的方法主要有线性规划法、非线性规划法和带有混合整数的非线性规划法。
申请公布号为CN103544551A的发明专利申请文献公开了一种蒸汽动力系统的操作优化方法，该方法包括以下步骤：确定蒸汽动力系统所包含的各设备的性能特征参数以及所述蒸汽动力系统的工艺参数；根据所述蒸汽动力系统的能量守恒方程、所述蒸汽动力系统的质量守恒方程、所述各设备的能量守恒方程、所述各设备的质量守恒方程以及所述各设备的性能特征参数和所述蒸汽动力系统的工艺参数建立所述蒸汽动力系统的非线性数学模型；对所述非线性数学模型进行模拟求解，得到模拟运算结果；设定所述非线性数学模型中优化变量的取值范围，以及设定所述非线性数学模型的优化目标函数；将所述模拟运算结果作为所述非线性数学模型优化运算的初始可行解，在所述优化变量的取值范围内，采用序列二次规划法和序贯线性规划法混合优化方法对所述非线性数学模块进行优化求解，得到优化求解结果；判断所述优化求解结果是否使得所述优化目标函数的值在所述优化变量的取值范围内达到最小，如果是，则将所述优化求解结果作为所述蒸汽动力模型的运行参数。
传统的控制系统和方法侧重于不同等级蒸汽的生产、消耗平衡关系，使得以最低的生产成本满足生产需求。然而，这种粗放型的生产控制方式并没有考虑到由此带来的环境成本，尤其在当前环境问题日益突出的大背景下，企业在生产中不仅需要考虑“节能”，更应统筹考虑“减排”。章建栋等提出了瓦斯与蒸汽动力系统的集成优化调度策略，但其重载实现瓦斯和蒸汽在全厂的调度平衡关系，而未能考虑环境因素对于实际生产方法的影响。戴文智和罗向龙等在蒸汽调度中考虑了环境因素，但是其实施案例中汽轮机的功能较为单一，同时对于环境因素的考虑也只局限于在目标函数中增加环境成本。
发明内容
本发明提供了一种蒸汽动力系统的运行操作方法，该运行操作方法满足生产成本和排放成本两个目标的优化策略，获得了蒸汽动力系统的最佳生产加工方案。
一种蒸汽动力系统的运行操作方法，包括以下步骤：
(1)采集蒸汽动力系统中各设备的性能特征参数以及蒸汽动力系统的工艺参数；
(2)根据蒸汽动力系统的能量、质量守恒方程建立所述蒸汽动力系统的混合整数线性规划模型，并设定所述混合整数线性规划模型的优化目标函数；
(3)以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得目标优化参数；
(4)将所述目标优化参数作为蒸汽动力系统的运行参数进行蒸汽动力系统的运行控制。
所述的设备为蒸汽动力系统中的产能设备、耗能设备以及运能系统(不同等级蒸汽管网)间的减温减压器；
其中，所述产能设备涉及的性能特征参数为：最大负荷，最大燃料流量，热损失比例和负荷比例回归函数的两个回归参数，蒸汽从过热温度到饱和温度的散热强度、饱和水介于新鲜水入锅炉时温度和蒸汽出锅炉时温度之间的比热以及启停费用；耗能设备涉及的性能特征参数为：最大负荷，最大能量负荷与全负荷效率回归函数的回归参数，蒸汽焓降以及启停费用。
所述蒸汽动力系统的工艺参数包括：燃料性能参数、生产需求参数以及物料成本参数。
其中，燃料性能参数为：热值、含S量、含C量、含N量以及燃料的燃烧产物含量；生产需求参数为：不同温度压力的蒸汽需求量、汽轮机动力需求量以及电力需求量；物料成本参数为：燃料、电和水的单价、废气排放收费标准。
本发明蒸汽动力系统的优化目标是在满足蒸汽、动力和电力需求的前提条件下，调整燃料参数、产能和耗能设备的运行参数、运能系统(蒸汽管网)内部不同等级蒸汽的平衡、设备启停参数以及外购电量和新鲜水购买量，从而在实现生产计划的前提下，同时兼顾生产成本和环境成本。所以，本发明蒸汽动力系统运行操作方法的优化目标为两个：(1)生产成本，(2)环境成本。
所述优化目标函数为：
式中，Ffuel，boiler，t为燃料在周期t内的流量；Cfuel为燃料的单价；Fwat，boiler，t为新鲜水在周期t内的流量；Cwat为新鲜水的单价；Elenet，t为周期t内的外购电量；Cele为外购电量的单价；Cch_bl为产能设备运行状态切换一次对应的成本；Cch_tb为耗能设备运行状态切换一次对应的成本；K为0-1变量，用来表征设备启停可选约束是否采用，K＝1表示采用约束，K＝0表示不采用约束；
以及：
式中，为产能设备中CO2的气体排放量；为产能设备中SO2的气体排放量；为产能设备中NOx的气体排放量；g为燃料燃烧后的排放气体，包括CO2、SO2和NOx；Cg为CO2、SO2和NOx排放费用的单价。
燃料燃烧排放出的污染物含量与燃料的成分组成直接相关，可以根据燃料燃烧产物的排放量，以及排放气体价格标准的变化以及气体污染物的排放量要求，来调整燃料种类，从而实现燃料使用种类和产物排放量的调控。
燃料性能参数方程为：
式中，Ffuel为燃料流量；和为气体污染物的排放量；ratC，fuel、ratN，fuel和ratS，fuel为燃料中C、N、S元素含量比率。
产能设备性能参数方程为：
式中，Tin为新鲜水入产能设备时的温度；为蒸汽出产能设备的温度；Cp为饱和水介于Tin和之间的比热；q为蒸汽从过热温度到饱和温度的散热强度；a和b为产能设备热损失比例和负荷比例回归函数的两个回归参数；Fstm，boiler，t为产能设备在周期t的蒸汽流量；Qfuel为燃料提供的热能。
以及：
式中，fuel为产能设备所用的燃料；Ffuel，boiler，t为燃料fuel在周期t内的流量；Hfuel为燃料fuel热值。
耗能设备性能参数方程为
式中，ΔHturb为蒸汽焓降；为耗能设备最大蒸汽负荷；A和B是耗能设备最大能量负荷与其全负荷效率回归函数的回归参数；Fstm，turb，t为耗能设备在周期t内的耗汽流量；Eturb，t为周期t内耗能设备turb提供的动力。
所述优化目标参数包括：燃料的种类和流量，产能设备、耗能设备和减温减压器的运行负荷，新鲜水流量，外购电量以及设备启停状态。
在炼厂或企业中，设备的平稳运行对于安全生产至关重要，设备频繁启停一方面会减少设备使用寿命，增加维修成本，另一方面还会影响正常的工作质量，降低炼油企业生产效益，所以要对设备运行状态的变化进行约束。
设备启停状态的确定方程为：
式中，Zn，t用来表征相比于周期t-1(t＝{2，3，…T})，设备n在周期t是否发生运行状态的切换；Yn，t为0-1变量，用来表征设备是否运行，Yn，t＝1表示设备处于开启状态，Yn，t＝0表示设备处于关闭状态。
与现有技术相比，本发明具有以下有益效果：
(1)本发明根据蒸汽动力系统的能量、质量守恒方程建立所述蒸汽动力系统的混合整数线性规划模型，并设定优化目标函数，以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得了满足生产成本和排放成本两个目标的优化策略，可帮助企业在保证生产成本的同时，适应节能环保要求趋势，科学调整生产运行策略；
(2)本发明将设备启停状态作为优化目标参数，并将该目标参数设为可选约束，提高了本发明运行操作方法的灵活性和应用性，便于该方法的推广应用。
附图说明
图1为本发明蒸汽动力系统结构示意图；
图2为本发明蒸汽动力系统运行操作方法的运行流程示意图；
图3为本发明实施案例1蒸汽动力系统设备拓扑图；
图4为本发明实施案例1中环境标准ES1的生产调度策略；
图5为本发明实施案例1中环境标准ES2的生产调度策略。
具体实施方式
为了更清晰的表述本发明内容，下面将结合本发明的附图和实施案例进行全面详细的说明。
首先，蒸汽动力系统不仅涉及能量的分配，同时与物料的加工生产有着紧密的联系，因此需要明确蒸汽动力系统在整个企业生产环节中所发挥的作用。
如图1所示，蒸汽动力系统主要由产能设备、不同等级蒸汽管网组成的运能系统以及耗能设备组成；与此同时，不同等级的蒸汽管网间还连接有减温减压器，高等级蒸汽通过减温减压器进行蒸汽降压进入低等级蒸汽管网中，实现不同等级蒸汽管网的生产调度平衡。
产能设备的主要功能是将水加热为高温高压蒸汽，用以满足换热、做功和发电的需求，产能设备做功后蒸汽降级为低等级蒸汽，并汇入低等级蒸汽管网。在具体实施时，产能设备主要考虑两种类型的锅炉，一种是只消耗燃料油的锅炉，另一种是既可以消耗燃料油，又可以消耗燃料气的锅炉，这两种锅炉在最大蒸汽产量上有所不同。
耗能设备的主要功能是通过消耗蒸汽，为生产提供驱动力、电力以及完成换热。在具体实施时，耗能设备主要考虑两种类型的汽轮机：发电汽轮机和动力汽轮机，其中发电汽轮机消耗高压或中压蒸汽以提供生产所需电力，如果自产电力不足以满足生产需要，则可以从外界电网购买进行补充；动力汽轮机消耗中压或低压蒸汽，为泵、压缩机等设备提供动力。其中，不同的耗能设备的产率和特征参数也各不相同。
本发明提供了一种蒸汽动力系统的运行操作方法，该方法的运行流程如图2所示，具体涉及以下步骤：
1、采集蒸汽动力系统中各设备的性能特征参数以及蒸汽动力系统的工艺参数；
该步骤1中提到的各设备为：产能设备、耗能设备和减温减压器；
其中，产能设备涉及的性能特征参数为：最大负荷，最大燃料流量，热损失比例和负荷比例回归函数的两个回归参数，蒸汽从过热温度到饱和温度的散热强度、饱和水介于新鲜水入锅炉时温度和蒸汽出锅炉时温度之间的比热以及启停费用。
耗能设备涉及的性能特征参数为：最大负荷，最大能量负荷与全负荷效率回归函数的回归参数，蒸汽焓降以及启停费用。
蒸汽动力系统的工艺参数主要包括：燃料性能参数、生产需求参数以及物料成本参数。
其中，燃料性能参数为：热值、含S量、含C量、含N量以及燃料的燃烧产物含量。
生产需求参数为：不同温度压力的蒸汽需求量、汽轮机动力需求量以及电力需求量。
物料成本参数为：燃料、电和水的单价、废气排放收费标准。
2、根据蒸汽动力系统的能量、质量守恒方程建立所述蒸汽动力系统的混合整数线性规划模型，并设定所述混合整数线性规划模型的优化目标函数；
为了便于理解，首先对本发明混合整数线性规划模型中的相关符号进行说明：
(1)上下标：drv：动力设备；g：气体污染物类别；n蒸汽动力系统中的设备；i：流入设备n的第i种物流；j：流出设备n的第j种物流；net：外界电网；t：生产周期；
(2)集合BOILER：产能设备集合；DRIVE：动力设备集合；DTURB：驱动动力设备的耗能设备集合；TURB：耗能设备集合；ETURB：驱动发电机的耗能设备集合；FUEL：燃料集合；STM：不同等级蒸汽集合；
(3)参变量D：公用工程需求量；Ele：蒸汽动力系统外购电或发电耗汽设备产电量；F：燃料或蒸汽的流量；G：气体污染物的排放量；H：热焓或热值；Q：外界提供给装置的热能；Y：表示设备是否运行的0-1变量；Z：表示设备是否发生启停转换的0-1变量。
本发明主要遵循以下约束(即守恒定律)，主要分为以下几类：通用约束、装置约束和供需约束。
一、通用约束：主要涉及各周期内设备的物料平衡约束、能量平衡约束以及设备的操作约束等。
(1)物料平衡
设备n(n＝{1，2，…N}在周期t的所有物料流入量等于所有物料的流出量之和，如式(1)所示)，
式中，n为蒸汽动力系统中产能或耗能设备或输送管网；i为流入设备n的第i种物料；j为流出设备n的第j种物料。
(2)能量平衡
设备n在周期t的所有供给能量之和等于流出能量之和，如式(2)所示，
式中，H为周期t流入或流出设备n的物料热焓值；Q为外界提供给设备n的热能。
(3)设备操作约束
设备n在周期t的工作负荷应处于正常范围，如式(3)所示，
式中，和分别指设备n在周期t的操作运行下限和上限。
在炼厂或企业中，设备的平稳运行对于安全生产至关重要，设备频繁启停一方面会减少设备使用寿命，增加维修成本，另一方面还会影响正常的工作质量，降低炼油企业生产效益，所以要对设备运行状态的变化进行约束，如式(4)～(7)所示；
式中，Zn，t用来表征相比于周期t-1(t＝{2，3，…T})，设备n在周期t是否发生运行状态的切换；Yn，t为0-1变量，用来表征设备是否运行，Yn，t＝1则表示设备处于开启状态，Yn，t＝0则表示设备处于关闭状态。
(4)污染物排放约束
Gboiler，gKg≤Gmax，g (26)
式中，g为燃料燃烧后的排放气体，包括CO2、SO2和NOx；G为排放量；K为0-1变量，用来表征污染物可选约束是否采用，K＝1表示采用污染物排放约束，K＝0表示不采用。
二、装置约束：主要是涉及产能设备和耗能设备的能量守恒。
(1)产能设备的约束
基于现有的热动力学原理和能量守恒定律可以建立控制系统中产能设备的线性模型。根据能量守恒定律，产能设备所耗燃料的总燃烧热能不仅可以从蒸汽能量结合锅炉工作效率推导得出，如式(9)所示；也可以由燃料的热值计算得出，如式(10)所示，由上述两式可得到产能设备的线性模型。
式中，Tin为新鲜水入产能设备时的温度；为蒸汽出产能设备的温度；Cp为饱和水介于Tin和之间的比热；q为蒸汽从过热温度到饱和温度的散热强度；a和b为产能设备热损失比例和负荷比例回归函数的两个回归参数；Fstm，boiler，t为产能设备在周期t内的蒸汽流量；Qfuel为燃料提供的热能。
式中，fuel为产能设备所用的燃料；Ffuel，boiler，t为燃料fuel在周期t内的流量；Hfuel为燃料fuel热值。
(2)耗能设备的约束
耗能设备将蕴含在高温高压蒸汽中的热能转化为机械能或电能，以驱动设备做功或者发电来满足生产所需。公式如式(11)所示；
式中，ΔHturb为蒸汽焓降；为耗能设备最大蒸汽负荷；A和B是耗能设备最大能量负荷与其全负荷效率回归函数的回归参数；Fstm，turb，t为耗能设备在周期t内的耗汽流量；Eturb，t为周期t内耗能设备turb提供的动力。Yturb，t为周期t内耗能设备开启或者关闭的状态。
三、供需约束：主要是满足蒸汽动力系统对蒸汽、动力和电力的需求，其约束公式如式(12)-(14)所示，
式中，Fi，stm，t为在周期t内从管路i流入运能系统等级为stm的蒸汽流量；Fj，stm，t为在周期t内从管路j流出运能系统等级为stm的蒸汽流量。
式中，Eturb，t为周期t内耗能设备turb提供的动力；Dturb，t为周期t内耗能设备turb的动力需求量。
式中，Eleturb，t为在周期t内耗能设备turb的发电量；Elenet，t为周期t内炼厂从外界电网所购买的电量；Dele，t为周期t内炼厂的电力需求量。
3、以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得目标优化参数；
本发明蒸汽动力系统的优化目标是在满足蒸汽、动力和电力需求的前提条件下，调整燃料参数、产能和耗能设备的运行参数、运能系统(蒸汽管网)内部不同等级蒸汽的平衡、设备启停参数以及外购电量和新鲜水购买量，从而实现生产成本和环境成本的同时兼顾，完成生产计划。所以，本发明蒸汽动力系统运行操作方法的优化目标为两个：(1)生产成本，(2)环境成本。
目标(1)具体包括：产能设备燃料成本、新鲜水的成本、外购电力成本以及设备的启停成本(启停成本是否作为目标(1)进行考虑，具有可选择性)，具体公式如式(15)所示：
式中，Ffuel，boiler，t为燃料在周期t内的流量；Cfuel为燃料的单价；Fwat，boiler，t为新鲜水在周期t内的流量；Cwat为新鲜水的单价；Elenet，t为周期t内的外购电量；Cele为外购电量的单价；Cch_bl为产能设备运行状态切换一次对应的成本；Cch_tb为耗能设备运行状态切换一次对应的成本；K为0-1变量，用来表征设备启停可选约束是否采用，K＝1表示采用约束，K＝0表示不采用约束。
目标(2)用于实现环境成本的控制，具体包括：CO2气体排放成本、SO2气体排放成本与NOx气体排放成本。
式中，为产能设备中CO2的气体排放量；为产能设备中SO2的气体排放量；为产能设备中NOx的气体排放量；g为燃料燃烧后的排放气体，包括CO2、SO2和NOx；Cg为CO2、SO2和NOx排放费用的单价。
燃料燃烧排放出的污染物含量与燃料的成分组成直接相关，可以根据公式(17)～(19)计算得出燃料燃烧产物的排放量，然后根据排放气体价格标准的变化以及气体污染物的排放量要求调整燃料种类，从而实现燃料燃烧产物种类和排放量的调控。
式中，Ffuel为燃料流量；和为气体污染物的排放量；ratC，fuel、ratN，fuel和ratS，fuel为燃料中C、N、S元素含量比率。
对上述公式(1)～(19)使用增广ε约束法进行求解，从而获取优化目标变量，即：在周期t内，燃料种类和流量Ffuel，boiler，t，产能设备蒸汽流量Fstm，boiler，t，耗能设备耗汽流量Fstm，turb，t，新鲜水流量Fwat，boiler，t，从外界电网所购买的电量Elenet，t，设备的启停选择K，以及减温减压器运行负荷(即周期t内，流入蒸汽管网的蒸汽流量Fi，stm，t和流出该蒸汽管网的蒸汽流量Fj，stm，t)
4、将所述目标优化参数作为蒸汽动力系统的运行参数进行蒸汽动力系统的运行控制。将上述优化目标变量参数数据输入控制系统中，指导蒸汽动力系统各周期的生产，实现本发明生成成本、环境成本同时调控的目的。
实施案例1
针对某化工企业的蒸汽动力系统进行运行操作方法的优化。
如图3所示，该蒸汽动力系统包括三种不同等级的蒸汽管网，上述蒸汽管网中的蒸汽分别为：高压蒸汽(9.20MPa，527℃)、中压蒸汽(3.61MPa，433℃)及低压蒸汽(1.10MPa，333℃)。整个蒸汽动力系统中由4台锅炉(产能设备)供应蒸汽，其中，1#和2#锅炉用来生产高压蒸汽，且只消耗燃料油。3#和4#锅炉用来提供中压蒸汽，不仅能消耗燃料油，还可以使用燃料气。图3中，汽轮机(耗能设备)共有11台，汽轮机1#、2#和6#、7#为发电汽轮机，驱动发电机为系统提供电能。3#-5#、8#-11#为动力汽轮机，用来驱动动力设备。与此同时，换热过程产生的中低压蒸汽也将汇入蒸汽系统中。该实例所研究的多周期模型总共包含6个连续生产周期，每个周期持续时间为1小时，并主要从生产费用和环境费用两方面进行分析。该案例建模使用GAMS 24.2.2，并采用CPLEX求解器。
针对上述蒸汽动力系统，锅炉和汽轮机的参数分别如表1和表2所示，六个周期中蒸汽、动力和电力的需求预测值如表3所示，资源单价如表4所示。
表1锅炉参数
表2蒸汽轮机参数(单位同上)
表3六个周期生产供需预测值(表中内容作了修改)
注：3类公用工程需求单位分别为t·h-1、MW·h-1、MW·h-1。
表4资源单价
燃料燃烧所排污染物含量与燃料的成分组成直接相关，表5为炼油企业两种常用燃料的相关参数。
表5燃料参数
目前，国内有两种较为典型的收费标准，如表6所示，其中：ES1是基于2003年以来的国内现行污染物收费标准；ES2数据源自外部成本计划，一种为研究空气污染物所造成的环境破坏而发展来的量化外部成本的方法。
表6两种环境收费标准
注：收费标准已按照当前汇率进行换算6.2￥/1$
综合上述参数，通过公式(1)～(19)换算，采用增广ε约束法计算方法求解，获得优化目标变量，即系统优化方案。在本发明所述控制系统的生产调控下，以及上述两种排放标准下，蒸汽动力系统中的原料及操作成本和环境成本的组合方案散点图如图4和图5所示，每种费用标准下有多种可选的生产操作策略。
进一步地，若分析图4和图5中满足f1+f2最小的操作策略，其燃料消耗与污染物排放情况如表7所示：
表7两种收费标准下控制系统的资源调配结果(即满足f1+f2最小)
分析表7可知，当排放费用标准提高时，控制系统将调整产能、耗能设备的运行负荷，协调燃料和电力等能源在系统中的分配，并协同CO2、SO2和NOx气体排放量的减少，使企业生产方案更加符合节能减排的要求。
Claims (5)
1.一种蒸汽动力系统的运行操作方法，其特征在于，包括以下步骤：
(1)采集蒸汽动力系统中各设备的性能特征参数以及蒸汽动力系统的工艺参数；
所述蒸汽动力系统的工艺参数包括：燃料性能参数、生产需求参数和物料成本参数；
(2)根据蒸汽动力系统的能量、质量守恒方程建立所述蒸汽动力系统的混合整数线性规划模型，并设定所述混合整数线性规划模型的优化目标函数；
所述优化目标函数为：
<mrow>
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<mo>(</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</munder>
<msub>
<mi>F</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</munder>
<msub>
<mi>F</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Ele</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</munder>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>c</mi>
<mi>h</mi>
<mo>_</mo>
<mi>b</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
</mrow>
</munder>
<msub>
<mi>Z</mi>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>c</mi>
<mi>h</mi>
<mo>_</mo>
<mi>t</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Ffuel,boiler,t为燃料在周期t内的流量；Cfuel为燃料的单价；Fwat,boiler,t为新鲜水在周期t内的流量；Cwat为新鲜水的单价；Elenet,t为周期t内的外购电量；Cele为外购电量的单价；Cch_bl为产能设备运行状态切换一次对应的成本；Cch_tb为耗能设备运行状态切换一次对应的成本；K为0-1变量，用来表征设备启停可选约束是否采用，K＝1表示采用约束，K＝0表示不采用约束；Zboiler,t表示周期t内产能设备是否发生启停转换的0-1变量；Zturb,t表示周期t内耗能设备是否发生启停转换的0-1变量；
以及：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>SO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>SO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>NO</mi>
<mi>X</mi>
</msub>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NO</mi>
<mi>X</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，为产能设备中CO2的气体排放量；为产能设备中SO2的气体排放量；为产能设备中NOx的气体排放量；为CO2排放费用的单价；为SO2排放费用的单价；为NOx排放费用的单价；
(3)以生产成本和废气排放成本为目标，采用增广ε约束法求解优化目标函数，获得目标优化参数；
所述优化目标参数包括：燃料的种类和流量，产能设备、耗能设备和减温减压器的运行负荷，新鲜水流量，外购电量以及设备启停状态；
(4)将所述目标优化参数作为蒸汽动力系统的运行参数进行蒸汽动力系统的运行控制。
2.如权利要求1所述的运行操作方法，其特征在于，燃料性能参数方程为：
<mrow>
<msub>
<mi>G</mi>
<mrow>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mn>3.667</mn>
<mo>&times;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>rat</mi>
<mrow>
<mi>C</mi>
<mo>,</mo>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<msub>
<mi>NO</mi>
<mi>X</mi>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mn>1.63</mn>
<mo>&times;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>rat</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>0.000938</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<msub>
<mi>SO</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mn>2</mn>
<mo>&times;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>rat</mi>
<mrow>
<mi>S</mi>
<mo>,</mo>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Ffuel为燃料流量；和为气体污染物的排放量；ratC,fuel、ratN,fuel和ratS,fuel为燃料中C、N、S元素含量比率。
3.如权利要求1所述的运行操作方法，其特征在于，产能设备性能参数方程为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>aF</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
<mi>max</mi>
</msubsup>
<msub>
<mi>Y</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
</mrow>
<mo>)</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>m</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
<mo>&Element;</mo>
<mi>F</mi>
<mi>U</mi>
<mi>E</mi>
<mi>L</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>&Element;</mo>
<mi>B</mi>
<mi>O</mi>
<mi>I</mi>
<mi>L</mi>
<mi>E</mi>
<mi>R</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Tin为新鲜水入产能设备时的温度；为蒸汽出产能设备的温度；Cp为饱和水介于Tin和之间的比热；q为蒸汽从过热温度到饱和温度的散热强度；a和b为产能设备热损失比例和负荷比例回归函数的两个回归参数；为耗能设备最大蒸汽负荷；Yboiler,t为周期t内产能设备开启或者关闭的状态；Fstm,boiler,t为产能设备在周期t内的蒸汽流量；Qfuel为燃料提供的热能；
以及：
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>H</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
<mo>&Element;</mo>
<mi>F</mi>
<mi>U</mi>
<mi>E</mi>
<mi>L</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mi>e</mi>
<mi>r</mi>
<mo>&Element;</mo>
<mi>B</mi>
<mi>O</mi>
<mi>I</mi>
<mi>L</mi>
<mi>E</mi>
<mi>R</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，fuel为产能设备所用的燃料；Ffuel,boiler,t为燃料fuel在周期t内的流量；Hfuel为燃料fuel热值。
4.如权利要求1所述的运行操作方法，其特征在于，耗能设备性能参数方程为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>6</mn>
<mrow>
<mn>5</mn>
<mi>B</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;H</mi>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mi>A</mi>
<msubsup>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<msubsup>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>max</mi>
</msubsup>
<msub>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>t</mi>
<mi>u</mi>
<mi>r</mi>
<mi>b</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mi>U</mi>
<mi>R</mi>
<mi>B</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，ΔHturb为蒸汽焓降；为耗能设备最大蒸汽负荷；A和B是耗能设备最大能量负荷与其全负荷效率回归函数的回归参数；Fstm,turb,t为耗能设备在周期t内的耗汽流量；Yturb,t为周期t内耗能设备开启或者关闭的状态；Eturb,t为周期t内耗能设备turb提供的动力。
5.如权利要求1所述的运行操作方法，其特征在于，设备启停状态的确定方程为：
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>N</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>N</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>N</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mn>2</mn>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>N</mi>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Zn,t用来表征相比于周期t-1(t＝{2,3,…T})，设备n在周期t是否发生运行状态的切换；Yn,t为0-1变量，用来表征设备是否运行，Yn,t＝1表示设备处于开启状态，Yn,t＝0表示设备处于关闭状态。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510073874.4A CN104657789B (zh) | 2015-02-12 | 2015-02-12 | 一种蒸汽动力系统的运行操作方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510073874.4A CN104657789B (zh) | 2015-02-12 | 2015-02-12 | 一种蒸汽动力系统的运行操作方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104657789A CN104657789A (zh) | 2015-05-27 |
| CN104657789B true CN104657789B (zh) | 2018-01-23 |
Family
ID=53248886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510073874.4A Expired - Fee Related CN104657789B (zh) | 2015-02-12 | 2015-02-12 | 一种蒸汽动力系统的运行操作方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104657789B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111475913A (zh) * | 2020-03-04 | 2020-07-31 | 曹雁青 | 一种蒸汽动力系统的运行操作优化方法及系统 |
| CN112307635B (zh) * | 2020-11-04 | 2023-12-19 | 万华化学集团股份有限公司 | 一种蒸汽系统扩能改造优化方法、电子设备及存储介质 |
| CN114755581B (zh) * | 2022-06-16 | 2022-09-02 | 德电北斗电动汽车有限公司 | 对置活塞磁力线性发电机的性能参数测试方法及相关装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103544551A (zh) * | 2013-11-12 | 2014-01-29 | 北京宜能高科科技有限公司 | 一种蒸汽动力系统的操作优化方法和装置 |
-
2015
- 2015-02-12 CN CN201510073874.4A patent/CN104657789B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103544551A (zh) * | 2013-11-12 | 2014-01-29 | 北京宜能高科科技有限公司 | 一种蒸汽动力系统的操作优化方法和装置 |
Non-Patent Citations (3)
| Title |
|---|
| 大型石化企业蒸汽动力系统运营优化;罗向龙;《石油炼制与化工》;20090512;第40卷(第5期);48-52 * |
| 石化企业蒸汽动力系统运行优化研究综述;戴文智;《计算机与应用化学》;20130528;第30卷(第5期);519-522 * |
| 蒸汽动力系统_经济成本分摊及运行优化;罗向龙;《化学工程》;20130715;第41卷(第7期);29-33 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104657789A (zh) | 2015-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Verrilli et al. | Model predictive control-based optimal operations of district heating system with thermal energy storage and flexible loads | |
| Mirzaei et al. | Evaluating the impact of multi-carrier energy storage systems in optimal operation of integrated electricity, gas and district heating networks | |
| Niu et al. | Implementation of a price-driven demand response in a distributed energy system with multi-energy flexibility measures | |
| Rong et al. | Polygeneration systems in buildings: A survey on optimization approaches | |
| Afzali et al. | Optimal design, operation and analytical criteria for determining optimal operating modes of a CCHP with fired HRSG, boiler, electric chiller and absorption chiller | |
| Cormio et al. | A regional energy planning methodology including renewable energy sources and environmental constraints | |
| Feng et al. | Energy, economy, and environment analysis and optimization on manufacturing plant energy supply system | |
| Li et al. | A reserve capacity model of AA-CAES for power system optimal joint energy and reserve scheduling | |
| CA2745357C (en) | Energy management system | |
| CN103544551B (zh) | 一种蒸汽动力系统的操作优化方法和装置 | |
| Zhang et al. | Optimization of distributed integrated multi-energy system considering industrial process based on energy hub | |
| Wasilewski | Optimisation of multicarrier microgrid layout using selected metaheuristics | |
| CN104657789B (zh) | 一种蒸汽动力系统的运行操作方法 | |
| Li et al. | Optimal design of combined cooling, heating and power multi-energy system based on load tracking performance evaluation of adjustable equipment | |
| CN104318321A (zh) | 一种钢铁企业多能源介质集成调度优化方法 | |
| Klöckner et al. | Is the coherence of coal phase-out and electrolytic hydrogen production the golden path to effective decarbonisation? | |
| CN113850409A (zh) | 一种考虑可再生能源接入的综合能源系统优化运行方法 | |
| Serrano-Arévalo et al. | The impact of circular economy indicators in the optimal planning of energy systems | |
| CN104424372B (zh) | 用于动力驱动设备选型的蒸汽动力系统建模方法和装置 | |
| CN102880912B (zh) | 一种煤气安全与经济评估方法 | |
| Wang et al. | Proactive reconfiguration of heat-exchanger supernetworks | |
| Abolarin et al. | An approach to energy management: A case study of a medium scale printing press in Lagos, Nigeria | |
| de Queiroz Lamas et al. | Thermoeconomic analysis applied to an alternative wastewater treatment | |
| Ayodele et al. | Mathematical methods and software tools for designing and economic analysis of hybrid energy system | |
| CN110619129A (zh) | 基于sos-2约束实现泛能站优化的线性化方法和装置 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180123 Termination date: 20210212 |