CN104268789B - 一种确定钢铁企业能源综合调度可行解的方法 - Google Patents

一种确定钢铁企业能源综合调度可行解的方法 Download PDF

Info

Publication number
CN104268789B
CN104268789B CN201410534373.7A CN201410534373A CN104268789B CN 104268789 B CN104268789 B CN 104268789B CN 201410534373 A CN201410534373 A CN 201410534373A CN 104268789 B CN104268789 B CN 104268789B
Authority
CN
China
Prior art keywords
variable
gas
value
variables
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
CN201410534373.7A
Other languages
English (en)
Other versions
CN104268789A (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.)
Wisdri Engineering and Research Incorporation Ltd
Original Assignee
Wisdri Engineering and Research Incorporation Ltd
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 Wisdri Engineering and Research Incorporation Ltd filed Critical Wisdri Engineering and Research Incorporation Ltd
Priority to CN201410534373.7A priority Critical patent/CN104268789B/zh
Publication of CN104268789A publication Critical patent/CN104268789A/zh
Application granted granted Critical
Publication of CN104268789B publication Critical patent/CN104268789B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Manufacturing & Machinery (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明适用于钢铁企业能源调度技术领域，提供了一种确定钢铁企业能源综合调度可行解的方法，包括：建立能源介质管网信息集Ε、单元设备信息集Ψ和每个单元设备的输入或输出变量对应的变量集X；判断所述变量集中无直接关联的单元设备或者不能根据已知的输入/输出变量计算得到的变量或者不能依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到的变量为自由变量；分别确定煤气能源介质和蒸汽及电力能源介质对应的单元设备的自由变量满足在其对应的上下限范围内的初始值；根据所述变量集中所述自由变量和非自由变量的相关性确定非自由变量的值；判断所述变量集中的变量的值是否满足单元设备的工艺模型约束，是，记录变量值输出可行解，否，调整所述变量集中违反约束的变量的值后再次判断，解决了钢铁企业能源综合调度可行解搜索的技术难题。

Description

一种确定钢铁企业能源综合调度可行解的方法
技术领域
本发明属于钢铁企业能源调度技术领域，尤其涉及一种确定钢铁企业能源综合调度可行解的方法。
背景技术
钢铁工业是国民经济的基础性支柱产业，同时又是资源、能源密集型产业。能源消耗是决定钢铁工业生产成本和利润的重要因素，也是影响环境负荷的主要原因。一方面，钢铁企业生产流程长，工序、设备繁多，各工序间相互衔接，且每种工序、设备都与多种能源介质关联；另一方面，钢铁企业需要用到的能源种类超过20种，这些能源介质不仅各自存在产耗、储存、缓冲和输配等多种形态，而且相互之间有着复杂的转换、替代等关联关系，这都使得整个钢铁企业能源***网络结构紧密耦合、错综复杂。因此，对钢铁企业能源***的研究具有理论和现实两方面的重要意义。
目前，规模在300万t以上钢铁企业的能源调度大都具备了较好的网络和数据等硬件基础条件，不同程度地实现了能源集中监视，以及能源调度与生产调度的信息共享，但仍然存在一些不可忽视的关键问题，名称为“钢铁企业能源中心的现状与发展趋势”(作者：张玉庆，徐化岩；冶金自动化，2011,35(4):15-19.)的文献中明确指出：对单一能源介质的调度模型研究较多，多介质协同调度研究较少，而单介质优化无法解决***优化问题，甚至出现不满足其他介质约束而使优化结果无效的情形。
对于钢铁企业能源综合调度问题而言，针对包含钢铁企业三种主要的能源介质：副产煤气、蒸汽和电力及其他介质的多能源综合调度中，存在着大量的线性和非线性、等式和不等式约束条件，如各种能源产耗、储存、缓冲、输配和转换等单元设备需满足工艺约束，每种能源介质自身还需满足动态平衡约束等。对这种复杂的强约束问题，如果没有一个好的约束处理机制作为支撑，在求解计算过程中，出现违反约束的不可行解的概率是非常大的，可以说，此时获得一个可行解的难度将不亚于搜索到全局最优解的难度。可行解难以获取的窘境不仅会降低能源综合调度问题求解计算的速度和效率，将大量的计算时间和资源消耗在可行解的搜索任务上，还可能导致求解过程不收敛，无法输出可行/最优的调度方案，严重影响能源调度与生产计划的实施。
发明内容
本发明实施例的目的在于提供一种确定钢铁企业能源综合调度可行解的方法，以解决现有技术难以获得钢铁企业能源综合调度可行解的问题。
本发明实施例是这样实现的，一种确定钢铁企业能源综合调度可行解的方法，所述方法包括以下步骤：
步骤1，根据钢铁企业能源***网络拓扑结构建立其对应的能源介质管网信息集E、单元设备信息集Ψ和所述单元设备信息集中每个单元设备的输入或输出变量对应的变量集X；
所述单元设备信息集为***中可调度的关键设备和其他公辅设备，包括：煤气混合站、煤气加压站、煤气柜、煤气放散塔、锅炉、蒸汽轮机和减温减压阀单元设备；
步骤2，判断所述变量集中无直接关联的单元设备或者不能根据已知的输入/输出变量计算得到的变量或者不能依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到的变量为自由变量，构建自由变量集X′；
步骤3，分别确定煤气能源介质和蒸汽及电力能源介质对应的单元设备的自由变量满足在其对应的上下限范围内的初始值；
确定所述煤气能源介质对应的单元设备的自由变量的初始值的方法包括：
确定一个种类的各个单元设备的变量的初始值后，对所述确定的初始值进行限幅处理，使其满足在其对应的上下限范围内，再确定下一个种类的各个单元设备的变量的初始值；
所述确定一个种类的各个单元设备的变量的初始值的方法包括：对该种类各个单元设备依次求初始值；单元设备的初始值为煤气富余量减去已经确定的该种类单元设备的变量初始值后与区间[0，1]范围内的随机数的乘积，其中该种类单元设备中的最后一个单元设备相乘的随机数为1；
蒸汽及电力能源介质对应的蒸汽轮机单元设备的自由变量的初始值为：
锅炉产生的蒸汽量减去主生产工序对蒸汽的预测需求量和已经确定的蒸汽轮机的变量初始值后与区间[0，1]范围内的随机数的乘积，其中最后一个蒸汽轮机相乘的随机数为1；
步骤4，根据所述变量集中所述自由变量和非自由变量的相关性确定非自由变量的值；
步骤5，判断所述变量集中的变量的值是否满足单元设备的工艺模型约束，是，记录变量值输出可行解，否，执行步骤6；
步骤6，调整所述变量集中违反约束的变量的值，执行步骤5。
本发明实施例提供的一种确定钢铁企业能源综合调度可行解的方法的有益效果包括：
1、本发明提供的一种确定钢铁企业能源综合调度可行解的方法，以蒙特卡罗随机模拟和启发式算法结合为特征的可行解设计方法，通过分组对煤气、蒸汽和电力子***中的自由变量进行预分配、限幅、偏差修正，以及检验和调整等处理过程，在无需复杂迭代搜索计算的前提下快速得到自然满足所有约束条件的可行解，解决了钢铁企业能源综合调度可行解搜索的技术难题。
2、所提出的方法可以在无需复杂迭代搜索计算的前提下快速得到自然满足所有约束条件的可行解，极大地节省了计算时间和软硬件资源；
3、应用该方法生成的可行解具有较好的多样性特征，这为钢铁企业能源综合调度问题中后续的最优解求解计算提供良好的初始条件。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的确定钢铁企业能源综合调度可行解的方法的总体流程图；
图2是本发明实施例提供的确定钢铁企业能源综合调度可行解的方法的详细流程图；
图3是典型钢铁企业子***拓扑结构图；
图4是典型钢铁企业蒸汽和电力子***拓扑结构图；
图5是钢铁企业实施例的煤气、蒸汽和电力子***拓扑结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
为了说明本发明所述的技术方案，下面通过具体实施例来进行说明。
如图1所示为本发明提供的确定钢铁企业能源综合调度可行解的方法的总体流程图，所述方法包括以下步骤：
步骤1，根据钢铁企业能源***网络拓扑结构建立其对应的能源介质管网信息集E、单元设备信息集Ψ和该单元设备信息集中每个单元设备的输入或输出变量对应的变量集X。
该单元设备信息集为***中可调度的关键设备和其他公辅设备，包括：煤气混合站、煤气加压站、煤气柜、煤气放散塔、锅炉、蒸汽轮机和减温减压阀单元设备。
步骤2，确定变量集中无直接关联的单元设备或者不能根据已知的输入/输出变量计算得到的变量或者不能依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到的变量为自由变量，构建自由变量集X′。
步骤3，分别确定煤气能源介质和蒸汽及电力能源介质对应的单元设备的自由变量满足在其对应的上下限范围内的初始值。
确定煤气能源介质对应的单元设备的自由变量的初始值的方法包括：
确定一个种类的各个单元设备的变量的初始值后，对该确定的初始值进行限幅处理，使其满足在其对应的上下限范围内，再确定下一个种类的各个单元设备的变量的初始值；
确定一个种类的各个单元设备的变量的初始值的方法包括：对该种类各个单元设备依次求初始值；单元设备的初始值为煤气富余量减去已经确定的该种类单元设备的变量初始值后与区间[0，1]范围内的随机数的乘积，其中该种类单元设备中的最后一个单元设备相乘的随机数为1；
蒸汽及电力能源介质对应的蒸汽轮机单元设备的自由变量的初始值为：
锅炉产生的蒸汽量减去主生产工序对蒸汽的预测需求量和已经确定的蒸汽轮机的变量初始值后与区间[0，1]范围内的随机数的乘积，其中最后一个蒸汽轮机相乘的随机数为1。
步骤4，根据变量集中自由变量和非自由变量的相关性确定非自由变量的值。
以已经完成初始化的自由变量为已知条件，计算得到所有非自由变量的值。
步骤5，判断变量集中的变量的值是否满足单元设备的工艺模型约束，是，记录变量值输出可行解，否，执行步骤6。
步骤6，调整变量集中违反约束的变量的值，执行步骤5。
本发明提供的一种确定钢铁企业能源综合调度可行解的方法，以蒙特卡罗随机模拟和启发式算法结合为特征的可行解设计方法，通过分组对煤气、蒸汽和电力子***中的自由变量进行预分配、限幅、偏差修正，以及检验和调整等处理过程，在无需复杂迭代搜索计算的前提下快速得到自然满足所有约束条件的可行解。
实施例一
如图2所示为本发明实施例提供的确定钢铁企业能源综合调度可行解的方法的详细流程图，由图2可知，本发明提供的确定钢铁企业能源综合调度可行解的方法的实施例包括：
在本发明实施例中，步骤1为根据钢铁企业能源***网络拓扑结构对钢铁企业能源综合调度问题初始化的过程，还包括：
步骤101，建立能源介质管网信息集E＝{GAS,STM,ELE,…}和单元设备信息集Ψ＝{MIX,PRE,HLD,EMI,BOI,TUR,VAL,…}。
GAS、STM和ELE分别为煤气、蒸汽和电力能源介质子集；MIX、PRE、HLD、EMI、BOI、TUR和VAL分别为煤气混合站、煤气加压站、煤气柜、煤气放散塔、锅炉、蒸汽轮机和减温减压阀单元设备子集。
步骤102，确定单元设备的工艺模型和约束条件，并将单元设备的工艺模型和约束条件以及按管网归类的各能源介质的物理平衡约束映射为约束集合C＝{C1,C2}，其中，C1为单元设备的工艺约束子集，C2为能源介质的物理平衡约束子集。
在本发明实施例中，步骤2包括：
步骤201，将步骤101中单元设备信息集Ψ包含的每个单元设备的输入/输出变量构造成变量集X＝(x1,x2,…,xi,…,xn)，其中，xi为第i个单元设备的变量值；结合单元设备的工艺特征对X进行相关性分析，将X中的变量区分为自由变量和非自由变量，将自由变量构造成新的变量集X′。
自由变量和非自由变量区分原则包括：
在能源***网络拓扑结构中无直接关联的单元设备的输入变量均作为自由变量；有直接关联的单元设备的输入变量均作为非自由变量。
单元设备的输出变量中，若能依据工艺模型由已知的输入/输出变量计算得到，或者依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到，则作为非自由变量，其余的作为自由变量。
步骤202：确定变量集X中所有变量的值域集D，满足其中，和分别为变量xi的取值范围下限和上限。
本发明实施例中，确定各个集合后，将钢铁企业能源综合调度可行解的搜索问题search(·)转化为由能源介质管网信息集E、单元设备信息集Ψ、约束集合C、变量集X和值域集D构成的五元组，表示为：search(·)＝<E,Ψ,C,X,D>。
本发明实施例中，步骤3中确定煤气能源介质和蒸汽及电力能源介质的单元设备的自由变量的初始值的过程为按能源介质的物理平衡约束子集C2和变量的值域集D对自由变量集X′中变量的初始化的过程，其中，煤气能源介质子集对应的单元设备包括煤气混合站、煤气加压站、煤气柜、煤气放散塔和锅炉，蒸汽和电力能源介质子集对应的单元设备包括：锅炉、蒸汽轮机和减温减压阀单元设备。
具体的，煤气能源介质子集对应的单元设备的自由变量的初始化过程包括：
步骤301，设煤气混合站和煤气加压站子集中包含的设备数量分别为NA，锅炉子集中包含的设备数量分别为NB；煤气混合站和煤气加压站子集中单元设备吸纳的煤气变量为锅炉燃烧消耗的煤气变量为其中，i为煤气介质序号，j为煤气混合站和煤气加压站的设备序号，t为调度周期序号，m为锅炉设备序号。
步骤302，计算和的值：
对于煤气混合站和煤气加压站和前NB-1个锅炉设备：
对于第NB个锅炉设备：
其中，为副产煤气i在调度周期t内的富余量，random(0,1)为取区间[0,1]内的随机数。
步骤303：根据变量的值域集D对和进行限幅处理，即：
其中，和分别为限幅之后的变量值，ubx和ubc、lbx和lbc分别为变量取值范围的上限和下限。
步骤304：累计已分配的变量值和中j∈MIX∪PRE，将煤气富余量与累计值之间的差值分配给煤气柜子集HLD中单元设备的吞吐量，即煤气柜的吞吐量满足：
对进行限幅处理得到
步骤305：累计已分配的变量值和中j∈MIX∪PRE，中(j∈HLD)，将煤气富余量与累计值之间的差值分配给煤气放散塔子集EMI中单元设备的放散量，即煤气放散塔的放散量满足：
对进行限幅处理得到
蒸汽和电力能源介质子集对应的单元设备的自由变量的初始化过程包括：
步骤311，定义ys,m,t为锅炉m在周期t内蒸汽s的产生量，为汽轮机n在周期t内入口蒸汽s的消耗量，其中，s为蒸汽序号，m为锅炉设备序号，n为汽轮机设备序号；锅炉和汽轮机设备子集中包含的设备数量分别为NB和NT。
步骤312，根据步骤301确定的锅炉消耗的煤气变量的初始化值，由锅炉设备的工艺模型计算得到锅炉产生的蒸汽量ys,m,t。
步骤313，计算的值，包括：
对于前NT-1个汽轮机设备：
对于第NT个汽轮机设备：
其中，为在周期t内主生产工序对蒸汽s的预测需求量之和。
步骤314，根据变量的值域集D对进行限幅处理，即：
其中，为限幅之后的变量值，uby和lby分别为变量取值范围的上限和下限。
进一步的，本发明提供的实施例中，步骤3中确定煤气能源介质子集和蒸汽和电力能源介质子集对应的单元设备的自由变量的初始值确定之后还包括：对初始值进行修正，具体为步骤305之后还包括：
步骤306，统计已分配变量值累计和与煤气富余量之间的偏差，定义为Δ，中j∈MIX∪PRE，中j∈HLD，中j∈EMI，则Δ满足：
判断Δ的大小，若Δ不等于0，则需要根据Δ的正负，同时结合变量的值域集D，对 和进行修正，具体修正方法如下：
若Δ＞0，则：
若Δ＜0，则：
其中，和分别为修正之后的变量值。
步骤314之后还包括：
步骤315，统计已分配的变量值并计算锅炉产生的蒸汽量之和与主生产工序预测值之间的差值，再计算变量值与所述差值之间的偏差，定义为Δ，则Δ满足：
判断Δ的大小，若Δ不等于0，则需要根据Δ的正负，同时结合变量的值域集D，对进行修正，具体修正方法如下：
若Δ＞0，则：
若Δ＜0，则：
其中，为修正之后的变量值。
进一步的，本发明实施例中步骤6中调整变量集中违反约束的变量的值的方法包括：
若变量集X中的变量x1和x2违反了单元设备的工艺约束子集C1中的某项约束g(x1,x2)，需要减小/增大x1和x2的值，则调整算法为：依次减小/增大x1和x2的值为x1·(1±ξ1)和x2·(1±ξ2)，直到满足该项约束或者调整x1和x2的值的次数超过上限，ξ1和ξ2分别为变量x1和x2的调整步长，调整次数上限可以为30次。具体如下：
1:t＝0,tmax＝30
2:while(！g(x1，x2)&&t＜tmax)do
3:x1′＝x1·(1±ξ1),x2′＝x2·(1±ξ2)；
4:x1＝x1′,x2＝x2′；
5:t++；
6:end while
其中，t为调整次数，tmax为调整次数上限。
具体执行时需要遍历单元设备的工艺约束子集C1中所有约束。
本发明实施例中步骤5中判断变量集中的变量的值满足单元设备的工艺模型约束记录变量值输出可行解后，还可以包括：
步骤7，判断是否生成了智能算法种群中要求数目的变量组，若否，则跳转至步骤3搜索下一个可行解，是则结束钢铁企业能源综合调度问题可行解的搜索过程，输出可行解集。
实施例二
本发明提供的实施例二为本发明提供的一种确定钢铁企业能源综合调度可行解的方法的具体应用实施例，如图3所示为典型钢铁企业煤气子***拓扑结构图，图4为典型钢铁企业蒸汽和电力子***拓扑结构图。由图3可知，钢铁企业的主要副产煤气有高炉煤气(BFG)、焦炉煤气(COG)和转炉煤气(LDG)，它们分别由高炉、焦炉和转炉产生。各种煤气介质的主管网上配备有对应的煤气柜和放散塔，主要生产工序用户，如炼铁厂、炼钢厂和宽厚板厂等也分别连接到不同的管网中。除单一煤气介质之外，可能还存在混合煤气，它由两种不同的单一介质在混合站混制而成，混合煤气也供给不同的主生产工序用户或热电厂等。由图4可知，蒸汽和电力子***中包含的单元设备主要有锅炉、蒸汽轮机和减温减压器，其中，燃料锅炉用于燃烧煤气或其他燃料产生蒸汽，蒸汽轮机则由蒸汽驱动产生电力；有些钢铁企业可能还包含燃气-蒸汽联合循环发电机组(CCPP)，在图4中，将CCPP拆解为燃气轮机、余热锅炉和蒸汽轮机，燃气轮机用煤气燃烧后产生的高温高压烟气直接推动透平产生电力，做完功之后的烟气输入余热锅炉中产生蒸汽。
图5为国内某钢铁企业的煤气、蒸汽和电力子***拓扑结构简图，针对该钢铁企业的煤气-蒸汽-电力等能源介质混合优化调度问题的可行解进行搜索。副产煤气有高炉煤气(BFG)、焦炉煤气(COG)和转炉煤气(LDG)，蒸汽也分高、中、低压三种不同品质的蒸汽，它们在各个调度周期内的富余量或需求量等预测结果见表1～2。煤气子***中与每种煤气对应分别设置了1座煤气柜和放散塔，蒸汽和电力子***中有2台动力锅炉和2台蒸汽轮机，其中，动力锅炉为全烧煤气锅炉，各设备参数和约束条件见表3～5。
表1 各调度周期内副产煤气的富余量(单位：kNm3/h)
表2 各调度周期内蒸汽和电力的需求量(单位：t/h或MW.h/h)
表3 煤气柜设备参数和约束条件
表4 锅炉设备参数和约束条件
表5 蒸汽轮机设备参数和约束条件
针对上述实施例，调用本发明提出的钢铁企业能源综合调度可行解的设计方法，搜索实施例在调度周期1内的可行解，构造由能源介质管网信息集E、单元设备信息集Ψ、约束集合C、变量集X和值域集D构成的五元组，其中：
E＝{高炉煤气，焦炉煤气，转炉煤气，高压蒸汽，中压蒸汽，低压蒸汽，内部电网}
Ψ＝{高炉煤气柜，焦炉煤气柜，转炉煤气柜，高炉煤气放散塔，焦炉煤气放散塔，转炉煤气放散塔，1#锅炉，2#锅炉，1#蒸汽轮机，2#蒸汽轮机}
X＝{x1,x2,x3,x4,x5,x6,x7,x8,x9,y1,x10,x11,x12,y2,x13,y3,y4,y5,y6,x14,y7,y8,y9}
其中，各变量对应的工艺含义如表6所示。
表6 变量集中各变量对应的工艺含义
按照自由变量和非自由变量的区分原则确定的自由变量集X′为：
X′＝{x1,x2,x3,x7,x8,x9,x10,x11,x12,x13,y4,y5,x14,y8}
约束集合C、中的单元设备的工艺约束子集C1见表3～5所示，能源介质的物理平衡约束子集C2由能源介质管网信息集E中每种介质在调度周期内的供需平衡共同组成，不具体列出，值域集D也可见表3～5所示，不再赘述。
表7列出了30组在值域集D中随机生成的自由变量集，经检验，所有随机生成的解均为不可行解，表7中的“×”表示为不可行解；表8列出了上述随机生成的自由变量集违反约束的详细情况，注意，表中仅仅列举了部分显见的违反约束情况，对于需经过复杂计算来检验的约束条件，表8中没有列出，但是，这并不影响该自由变量集是否为可行解的判断。
表9中列出了30组调用本发明提出的方法生成的自由变量集，经检验，所有解均为可行解，表8中的“√”表示为可行解。
表7和表9中分别列出的30组自由变量集仅仅是一个示例，实际上，作者进行了多次对比试验，结果表明，随机生成的解全部为不可行解，而应用本发明提出的方法生成的解则均为可行解。
表7 随机生成的自由变量集
表8 随机生成的自由变量集违反约束情况表
序号 违反约束情况(部分)
1 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
2 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
3 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、2#锅炉混合煤气热值范围约束、高压蒸汽平衡
4 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
5 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
6 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡
7 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡、1#锅炉凝气范围约束
8 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡、1#锅炉凝气范围约束
9 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡、1#锅炉凝气范围约束
10 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
11 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡
12 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
13 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、2#锅炉混合煤气热值范围约束、高压蒸汽平衡、1#锅炉凝气范围约束
14 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
15 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
16 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡、1#锅炉凝气范围约束
17 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡、1#锅炉凝气范围约束
18 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
19 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
20 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
21 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡、1#锅炉凝气范围约束
22 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
23 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡
24 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
25 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、2#锅炉混合煤气热值范围约束、高压蒸汽平衡、1#锅炉混合煤气热值范围约束
26 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、高压蒸汽平衡
27 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡
28 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡
29 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、1#锅炉混合煤气热值范围约束、高压蒸汽平衡、1#锅炉混合煤气热值范围约束
30 高炉煤气平衡、焦炉煤气平衡、转炉煤气平衡、2#锅炉混合煤气热值范围约束、高压蒸汽平衡、1#锅炉混合煤气热值范围约束
表9 调用本发明提出的方法生成的自由变量集
综上所述，本发明公开的方法较好的解决了所面临的技术难题，提出了一种以蒙特卡罗随机模拟和启发式算法结合为特征的钢铁企业能源综合调度可行解设计方法，可以在无需复杂迭代搜索计算的前提下快速得到自然满足所有约束条件的可行解，极大地节省了计算时间和软硬件资源；应用该设计方法生成的可行解具有较好的多样性特征，这为钢铁企业能源综合调度问题中后续的最优解求解计算提供良好的初始条件。
本领域普通技术人员还可以理解，实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，所述的程序可以存储于一计算机可读取存储介质中，所述的存储介质，包括ROM/RAM、磁盘、光盘等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种确定钢铁企业能源综合调度可行解的方法，其特征在于，所述方法包括：
步骤1，根据钢铁企业能源***网络拓扑结构建立其对应的能源介质管网信息集E、单元设备信息集Ψ和所述单元设备信息集中每个单元设备的输入或输出变量对应的变量集X；
所述单元设备信息集为***中可调度的关键设备和其他公辅设备，包括：煤气混合站、煤气加压站、煤气柜、煤气放散塔、锅炉、蒸汽轮机和减温减压阀单元设备；
步骤2，判断所述变量集中无直接关联的单元设备或者不能根据已知的输入/输出变量计算得到的变量或者不能依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到的变量为自由变量，构建自由变量集X′；
步骤3，分别确定煤气能源介质和蒸汽及电力能源介质对应的单元设备的自由变量满足在其对应的上下限范围内的初始值；
确定所述煤气能源介质对应的单元设备的自由变量的初始值的方法包括：
确定一个种类的各个单元设备的变量的初始值后，对所述确定的初始值进行限幅处理，使其满足在其对应的上下限范围内，再确定下一个种类的各个单元设备的变量的初始值；
所述确定一个种类的各个单元设备的变量的初始值的方法包括：对该种类各个单元设备依次求初始值；单元设备的初始值为煤气富余量减去已经确定的该种类单元设备的变量初始值后与区间[0，1]范围内的随机数的乘积，其中该种类单元设备中的最后一个单元设备相乘的随机数为1；
蒸汽及电力能源介质对应的蒸汽轮机单元设备的自由变量的初始值为：
锅炉产生的蒸汽量减去主生产工序对蒸汽的预测需求量和已经确定的蒸汽轮机的变量初始值后与区间[0，1]范围内的随机数的乘积，其中最后一个所述蒸汽轮机相乘的随机数为1；
步骤4，根据所述变量集中所述自由变量和非自由变量的相关性确定非自由变量的值；
步骤5，判断所述变量集中的变量的值是否满足单元设备的工艺模型约束，是，记录变量值输出可行解，否，执行步骤6；
步骤6，调整所述变量集中违反约束的变量的值，执行步骤5。
2.如权利要求1所述的方法，其特征在于，所述步骤1为根据所述钢铁企业能源***网络拓扑结构对钢铁企业能源综合调度问题初始化的过程，还包括：
步骤101，建立能源介质管网信息集E＝{GAS,STM,ELE,…}和单元设备信息集Ψ＝{MIX,PRE,HLD,EMI,BOI,TUR,VAL,…}；
GAS、STM和ELE分别为煤气、蒸汽和电力能源介质子集；MIX、PRE、HLD、EMI、BOI、TUR和VAL分别为煤气混合站、煤气加压站、煤气柜、煤气放散塔、锅炉、蒸汽轮机和减温减压阀单元设备子集；
步骤102，确定所述单元设备的工艺模型和约束条件，并将所述单元设备的工艺模型和约束条件以及按管网归类的各能源介质的物理平衡约束映射为约束集合C＝{C1,C2}，其中，C1为单元设备的工艺约束子集，C2为能源介质的物理平衡约束子集。
3.如权利要求2所述的方法，其特征在于，所述步骤2包括：
步骤201，将所述步骤101中所述单元设备信息集Ψ包含的每个单元设备的输入/输出变量构造成变量集X＝(x1,x2,…,xi,…,xn)，其中，xi为第i个单元设备的变量值；结合单元设备的工艺特征对X进行相关性分析，将X中的变量区分为自由变量和非自由变量，将自由变量构造成新的变量集X′；
步骤202：确定所述变量集X中所有所述变量的值域集D，满足其中，和分别为变量xi的取值范围下限和上限；
所述自由变量和非自由变量的区分方法包括：
在能源***网络拓扑结构中无直接关联的单元设备的输入变量均作为自由变量；有直接关联的单元设备的输入变量均作为非自由变量；
单元设备的输出变量中，能够依据工艺模型由已知的输入/输出变量计算得到，或者依据其它的功/能量/热量守恒或物料/质量守恒定律间接计算得到，则为非自由变量，其余的为自由变量。
4.如权利要求1所述的方法，其特征在于，所述步骤3中所述煤气能源介质子集对应的单元设备的自由变量的初始化过程包括：
步骤301，设煤气混合站和煤气加压站子集中包含的设备数量分别为NA，锅炉子集中包含的设备数量分别为NB；煤气混合站和煤气加压站子集中单元设备吸纳的煤气变量为锅炉燃烧消耗的煤气变量为其中，i为煤气介质序号，j为煤气混合站和煤气加压站的设备序号，t为调度周期序号，m为锅炉设备序号；
步骤302，计算和的值：
对于煤气混合站和煤气加压站和前NB-1个锅炉设备：
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mi>o</mi> <mi>m</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>*</mo> <mrow> <mo>(</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> <mo>;</mo> </mrow>
<mrow> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mi>o</mi> <mi>m</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>*</mo> <mrow> <mo>(</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>A</mi> </msub> </munderover> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> <mo>;</mo> </mrow>
对于第NB个锅炉设备：
<mrow> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>A</mi> </msub> </munderover> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>;</mo> </mrow>
其中，为副产煤气i在调度周期t内的富余量，random(0,1)为取区间[0,1]内的随机数；
步骤303：根据变量的值域集D对和进行限幅处理，即：
<mrow> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>lb</mi> <mi>x</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>&lt;</mo> <msub> <mi>lb</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msub> <mi>lb</mi> <mi>x</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>ub</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>ub</mi> <mi>x</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>&gt;</mo> <msub> <mi>ub</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>lb</mi> <mi>c</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>&lt;</mo> <msub> <mi>lb</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msub> <mi>lb</mi> <mi>c</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>ub</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>ub</mi> <mi>c</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>c</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>&gt;</mo> <msub> <mi>ub</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
其中，和分别为限幅之后的变量值，ubx和ubc、lbx和lbc分别为变量取值范围的上限和下限；
步骤304：累计已分配的变量值和 中j∈MIX∪PRE，将煤气富余量与累计值之间的差值分配给煤气柜子集HLD中单元设备的吞吐量，即煤气柜的吞吐量满足：
<mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> </mrow> </munder> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>H</mi> <mi>L</mi> <mi>D</mi> <mo>;</mo> </mrow>
对进行限幅处理得到
步骤305：累计已分配的变量值和 中j∈MIX∪PRE，中j∈HLD，将煤气富余量与累计值之间的差值分配给煤气放散塔子集EMI中单元设备的放散量，即煤气放散塔的放散量满足：
<mrow> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> </mrow> </munder> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>H</mi> <mi>L</mi> <mi>D</mi> </mrow> </munder> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>E</mi> <mi>M</mi> <mi>I</mi> <mo>;</mo> </mrow>
对进行限幅处理得到
5.如权利要求1所述的方法，其特征在于，所述步骤3中所述蒸汽和电力能源介质子集对应的单元设备的自由变量的初始化过程包括：
步骤311，定义ys,m,t为锅炉m在周期t内蒸汽s的产生量，为汽轮机n在周期t内入口蒸汽s的消耗量，其中，s为蒸汽序号，m为锅炉设备序号，n为汽轮机设备序号；锅炉和汽轮机设备子集中包含的设备数量分别为NB和NT；
步骤312，根据步骤301确定的锅炉消耗的煤气变量的初始化值，由锅炉设备的工艺模型计算得到锅炉产生的蒸汽量ys,m,t；
步骤313，计算的值，包括：
对于前NT-1个汽轮机设备：
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mi>o</mi> <mi>m</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>*</mo> <mrow> <mo>(</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msub> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>D</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>t</mi> <mi>m</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mi>U</mi> <mi>R</mi> <mo>;</mo> </mrow>
对于第NT个汽轮机设备：
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msub> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>D</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>t</mi> <mi>m</mi> </mrow> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mi>U</mi> <mi>R</mi> <mo>;</mo> </mrow>
其中，为在周期t内主生产工序对蒸汽s的预测需求量之和；
步骤314，根据变量的值域集D对进行限幅处理，即：
<mrow> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>lb</mi> <mi>y</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>&lt;</mo> <msub> <mi>lb</mi> <mi>y</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msub> <mi>lb</mi> <mi>x</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>ub</mi> <mi>y</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>ub</mi> <mi>y</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>&gt;</mo> <msub> <mi>ub</mi> <mi>y</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
其中，为限幅之后的变量值，uby和lby分别为变量取值范围的上限和下限。
6.如权利要求4所述的方法，其特征在于，所述步骤305之后还包括：
步骤306，统计已分配变量值累计和与煤气富余量之间的偏差，定义为Δ，中j∈MIX∪PRE，中j∈HLD，中j∈EMI，则Δ满足：
<mrow> <mi>&amp;Delta;</mi> <mo>=</mo> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> </mrow> </munder> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>H</mi> <mi>L</mi> <mi>D</mi> </mrow> </munder> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>E</mi> <mi>M</mi> <mi>I</mi> </mrow> </munder> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>;</mo> </mrow>
判断Δ的大小，若Δ不等于0，则需要根据Δ的正负，同时结合变量的值域集D，对 和进行修正，具体修正方法如下：
若Δ＞0，则：
<mrow> <msubsup> <mover> <mi>x</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>ub</mi> <mi>x</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>ub</mi> <mi>x</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>ub</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>ub</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>H</mi> <mi>L</mi> <mi>D</mi> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>ub</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>ub</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>E</mi> <mi>M</mi> <mi>I</mi> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>c</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>ub</mi> <mi>c</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>m</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>ub</mi> <mi>c</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> <mo>;</mo> </mrow>
若Δ＜0，则：
<mrow> <msubsup> <mover> <mi>x</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>x</mi> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msubsup> <mover> <mi>x</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mo>{</mo> <mi>M</mi> <mi>I</mi> <mi>X</mi> <mo>&amp;cup;</mo> <mi>P</mi> <mi>R</mi> <mi>E</mi> <mo>}</mo> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>H</mi> <mi>L</mi> <mi>D</mi> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mrow> <mo>(</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>E</mi> <mi>M</mi> <mi>I</mi> <mo>;</mo> </mrow>
<mrow> <msubsup> <mover> <mi>c</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>c</mi> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>m</mi> </munder> <mrow> <mo>(</mo> <msubsup> <mover> <mi>c</mi> <mo>~</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>g</mi> <mi>a</mi> <mi>s</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> <mo>;</mo> </mrow>
其中，和分别为修正之后的变量值。
7.如权利要求5所述的方法，其特征在于，所述步骤314之后还包括：
步骤315，统计已分配的变量值并计算锅炉产生的蒸汽量之和与主生产工序预测值之间的差值，再计算变量值与所述差值之间的偏差，定义为Δ，则Δ满足：
<mrow> <mi>&amp;Delta;</mi> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mi>O</mi> <mi>I</mi> </mrow> </munder> <msub> <mi>y</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>D</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>s</mi> <mi>t</mi> <mi>m</mi> </mrow> </msubsup> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mi>U</mi> <mi>R</mi> </mrow> </munder> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>;</mo> </mrow>
判断Δ的大小，若Δ不等于0，则需要根据Δ的正负，同时结合变量的值域集D，对进行修正，具体修正方法如下：
若Δ＞0，则：
<mrow> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>ub</mi> <mi>y</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>n</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>ub</mi> <mi>y</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mi>U</mi> <mi>R</mi> <mo>;</mo> </mrow>
若Δ＜0，则：
<mrow> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>y</mi> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>n</mi> </munder> <mrow> <mo>(</mo> <msubsup> <mover> <mi>y</mi> <mo>~</mo> </mover> <mrow> <mi>s</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>lb</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mi>U</mi> <mi>R</mi> <mo>;</mo> </mrow>
其中，为修正之后的变量值。
8.如权利要求1所述的方法，其特征在于，所述步骤6中调整变量集中违反约束的变量的值的方法包括：
若变量集X中的变量x1和x2违反了单元设备的工艺约束子集C1中的某项约束g(x1,x2)，需要减小/增大x1和x2的值，则调整算法为：依次减小/增大x1和x2的值为x1·(1±ξ1)和x2·(1±ξ2)，直到满足该项约束g(x1,x2)或者调整x1和x2的值的次数超过上限，ξ1和ξ2分别为变量x1和x2的调整步长，调整上限次数为30；
执行时需要遍历单元设备的工艺约束子集C1中所有约束。
9.如权利要求1所述的方法，其特征在于，所述步骤5中判断所述变量集中的变量的值满足单元设备的工艺模型约束记录变量值输出可行解后，还包括：
步骤7，判断是否生成了智能算法种群中要求数目的变量组，若否，则跳转至所述步骤3搜索下一个可行解，是则结束钢铁企业能源综合调度问题可行解的搜索过程，输出所述可行解集。
CN201410534373.7A 2014-10-11 2014-10-11 一种确定钢铁企业能源综合调度可行解的方法 Active CN104268789B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410534373.7A CN104268789B (zh) 2014-10-11 2014-10-11 一种确定钢铁企业能源综合调度可行解的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410534373.7A CN104268789B (zh) 2014-10-11 2014-10-11 一种确定钢铁企业能源综合调度可行解的方法

Publications (2)

Publication Number Publication Date
CN104268789A CN104268789A (zh) 2015-01-07
CN104268789B true CN104268789B (zh) 2017-12-22

Family

ID=52160309

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410534373.7A Active CN104268789B (zh) 2014-10-11 2014-10-11 一种确定钢铁企业能源综合调度可行解的方法

Country Status (1)

Country Link
CN (1) CN104268789B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104881713B (zh) * 2015-05-22 2018-10-30 中冶南方工程技术有限公司 实现优化算法与钢铁企业能源综合调度问题解耦的方法
CN104991531B (zh) * 2015-05-22 2017-11-21 中冶南方工程技术有限公司 一种确定钢铁企业副产煤气***优化调度可行解的方法
CN105069533B (zh) * 2015-08-19 2018-08-07 浙江大学 一种基于随机预测模型的钢铁企业多能源优化调度方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10112681B4 (de) * 2000-03-31 2004-02-12 International Business Machines Corp. Computersystem mit Verfahren zum Planen von Vorgängen in einem Stahlwerk und Programmspeichereinrichtung zum Durchführen des Verfahrens

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"钢铁企业富余煤气-蒸汽-电力耦合模型及其应用";张琦等;《化工学报》;20110331;第62卷(第3期);第753-758页 *

Also Published As

Publication number Publication date
CN104268789A (zh) 2015-01-07

Similar Documents

Publication Publication Date Title
CN104239981B (zh) 一种基于多目标优化算法的能源动态平衡与优化调度方法
CN104318321B (zh) 一种钢铁企业多能源介质集成调度优化方法
CN104268712B (zh) 基于改进的混合多种群进化算法的能源平衡与调度方法
Xie et al. The environmental efficiency analysis of China’s power generation sector based on game cross-efficiency approach
CN104573875B (zh) 一种低碳化的电源电网优化规划的方法
CN105005940A (zh) 计及相关性的电-气互联***概率最优潮流计算方法
CN109784554B (zh) 一种基于主从博弈的电力***最优调度方法
He et al. A new coal gas utilization mode in China’s steel industry and its effect on power grid balancing and emission reduction
CN104915728B (zh) 一种确定钢铁企业蒸汽***优化调度可行解的方法
CN104268789B (zh) 一种确定钢铁企业能源综合调度可行解的方法
Wen et al. Regional differences and influencing factors in the CO 2 emissions of China’s power industry based on the panel data models considering power-consuming efficiency factor
Sheng et al. How to efficiently promote distributed energy resources in China: Using a nonparametric econometric method
CN104966156B (zh) 一种钢铁企业能源综合调度问题的双层寻优方法
Yu et al. An improved combined heat and power economic dispatch model for natural gas combined cycle power plants
Liu et al. Improving the economic and environmental benefits of the energy system: A novel hybrid economic emission dispatch considering clean energy power uncertainty
Li et al. Development status evaluation and path analysis of regional clean energy power generation in China
Jia et al. Spatial differences and influencing factors of energy poverty: evidence from provinces in China
Li et al. Decentralized optimization for integrated electricity–heat systems with data center based energy hub considering communication packet loss
CN104991531B (zh) 一种确定钢铁企业副产煤气***优化调度可行解的方法
CN116739404A (zh) 一种综合能源***熵态模型中网络化特征模块的构建方法
Ma et al. Energy saving maximization on combined heat and power units in different scenarios
Zhang et al. A Two‐Stage Optimization Model of Capacity Allocation and Regulation Operation for Virtual Power Plant
Tian et al. Assessing the Efficiency and CO2 Reduction Performance of China’s Regional Wind Power Industry Using an Epsilon-Based Measure Model
Ivanova et al. Cost—Benefit analysis of electric boiler at combined heat and power plants
Liu et al. Influencing factors of CO2 emissions in Chinese power industry: a study from the production and consumption perspectives

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