A Constrained Multi-objective Online Operation Optimization Method of Collaborative Distillation and Heat Exchanger Network
-
摘要: 针对蒸馏装置与换热网络间缺乏协同优化导致的分馏精度差和能耗高的问题，提出了一种基于代理模型的约束多目标在线协同操作优化方法.为了解决蒸馏装置与换热网络操作参数协同优化时存在的计算耗时和约束的问题，构建Kriging代理模型来近似目标函数和约束条件，提出了基于随机欠采样和Adaboost的分类代理模型（RUSBoost）来解决类别不平衡的收敛判定预测问题.提出了基于多阶段自适应约束处理的代理模型的模型管理方法，该方法采用基于参考向量激活状态的最大化改善期望准则和可行概率准则更新机制来平衡优化初始阶段种群的多样性和可行性，采用支配参考点的置信下限准则更新机制加快收敛速度.通过不断与机理模型交互来在线更新代理模型，实现在线操作优化.通过测试函数和仿真实例验证了本文方法的有效性.Abstract: To solve low separation precision and high energy-consuming caused by lacking of collaborative operation optimization between atmospheric distillation unit and heat exchanger network, this paper presents a constrained multi-objective online collaborative operation optimization method based on surrogates model. To solve the problems of time-consuming and constraints in operating parameters collaborative optimization of distillation unit and heat exchanger network, the Kriging surrogate models are built to approximate each objective function and constraint, a classification surrogate model based on random undersampling and Adaboost is presented to solve the class imbalance of convergence prediction problem. A model management method of surrogate models based on multi-stage adaptive constraint handling is proposed. The method uses the maximization expected improvement and probability of feasibility criterion updating mechanism based on the reference vector activation state to balance the diversity and feasibility of the population at the early stage of evolution. The convergence rate is accelerated by using lower confidence bound criterion update mechanism dominating the reference point. By constantly interacting with the mechanism model to online update the surrogate model, the online operation optimization is realized. The efficiency of the proposed algorithm is validated by the results of benchmark functions and simulation example.
-
随着社会经济的迅速发展, 工业生产能力增强, 便利了我们的生活. 但与此同时, 工业生产、汽车尾气、食品加工等所产生的废气也带来了严重的环境污染问题. 其中的挥发性有机化合物(Volatile organic compounds, VOCs)是导致酸雨、雾霾等的主要原因[1], 也给人的生理健康带来了极大的危害. 杨威[2]结合国内外的VOCs研究现状, 对大连市的VOCs排放及治理对策进行了详细阐述. 指出了VOCs污染的治理刻不容缓, 对VOCs种类和浓度的监测成为首要问题. 此外在医疗健康[3]、食品行业[4]和智能家居[5]等行业, VOCs 的痕量检测也起着至关重要的作用.
长期以来, 半导体传感器由于方案成熟、灵敏度高、响应速度快且制作简单, 在VOCs 检测方面得到了广泛的应用[6]. 在半导体传感器领域, 气敏响应是判断气体的种类和浓度的重要指标. 然而, 除了少数仅对某种气体响应的传感器之外, 大多数传感器可以同时对多种气体产生响应. 在不同种类不同浓度气体中有相同的气敏响应, 即具有较差的选择性[7], 这成为制约其发展的瓶颈. 即便是具有良好选择性的只对某种气体响应的传感器, 也无法适应复杂环境多种气体的检测. 为此, 我们尝试使用动态测试方法[8]来提高选择性, 并利用传感器对多种气体都有响应的特性, 增加输出信号的信息量, 从而辨别不同种类不同浓度的气体. 但从目前动态测试的发展水平来看, 由于半导体气体传感器的敏感机理和建模研究一直是一个难点, 而温度调制下的气体传感器的动态响应信号和提取的特征参数的物理意义不明确, 导致气体传感器响应模型的建立和温度调制模式的优化困难[9]. 基于目前的研究现状, 以及对静态性能指标的研究, 阐述了静态性能指标与动态测试信号的对应关系, 提出了在动态测试中选取合适传感器的方法.
目前的动态测试方法主要是周期性循环加热波方法[10-13], 有部分文献在周期性循环加热波方法的基础上研究了多级伪随机序列[14-15]、自适应式温度调制[16-17 ]方法. 在周期性循环加热波方法中, 使用调制波形以矩形波[18-20]为主. 为此本文探究了静态响应时间、最佳工作温度指标、不同温度的气敏响应和传感器的复现性与动态响应信号的关系, 提出了适合于动态测试的半导体传感器的选择方法; 以矩形波为例研究了周期、占空比、工作温度范围对动态响应信号的影响, 在不降低动态响应信号品质的前提下, 缩短在实际应用中的响应时间和功耗; 结合支持向量机(Support vector machine, SVM)算法[21]实现了对不同浓度的丙酮、甲醛、甲酸、乙酸丁酯和乙醇的分类识别.
1. 实验过程
1.1 二氧化锡传感器的制备
用电子天平称取0.267 g SnCl2·2H2O与0.30 g的NaOH放入20 ml的无水乙醇和20 ml的去离子水的混合液中并用磁力搅拌器搅拌20 min, 然后将得到的溶液装入反应釜中, 在170 ℃条件下加热24 h. 待加热结束后自然冷却至室温, 用无水乙醇和去离子水交替清洗3次, 获得沉淀物后在60 ℃条件下干燥. 将干燥后的材料进行研磨, 然后在500 ℃条件下煅烧2 h获得SnO2纳米材料.
传感器主要由底座、Ni-Cr加热丝、陶瓷管构成, 如图1(a)所示. 将少量SnO2纳米材料与无水乙醇3 : 1混合. 然后在超声机中超声处理得到粘稠状膏体. 将带有金电极的陶瓷管蘸取少量膏体后在表面形成一层敏感膜. 再将陶瓷管放入马弗炉中在300 ℃下煅烧2 h以去除乙醇, 在表面形成一层致密的氧化膜. 冷却后将Ni-Cr加热丝插入陶瓷管中焊接, 最终放在老化平台上在400 ℃下老化7天[22].
图 1 传感器结构示意图和气敏测试装置示意图Fig. 1 Schematic diagram of sensor structure and Schematic diagram of gas sensing test device1.2 测试过程
本实验气敏测试装置的示意图如图1(b)所示. 气敏测试装置由干燥空气气瓶、反应气室、可编程直流电源、数字源表和计算机组成, 同时通过Labview开发环境实现对传感器阻值的读取, 以及对可编程电源实施远程控制. 气瓶内的气体为纯净干燥合成空气(79%氮气、21%氧气、0%湿度), 打开气瓶的气阀, 气体通入气室之中, 可以将反应气室内的空气送往屋外; 气室中央底座和六角传感器直插连接; 气室一侧接可编程直流电源, 通过调节电源的电压使传感器达到需要的工作温度, 或者按照规定的波形变化; 另一侧接数字源表, 用于测试工作状态下传感器的阻值; 计算机通过Labview虚拟仪表采集信息, 以数据形式, 记录传感器的阻值情况.
首先进行静态测试, 获得传感器的基础数据. 先定义传感器的气敏响应为
$$ \small S=\frac{{R}_{{\rm{air}}}}{{R}_{{\rm{gas}}}}=\frac{{I}_{{\rm{gas}}}}{{I}_{{\rm{air}}}} $$ (1) 式中, Rair和Rgas分别表示半导体传感器在空气中和待测气体中的电阻值, 相应地, Iair和Igas分别代表空气中和待测气体中的电流值[23]. 定义传感器的响应时间T90为从Iair上升到平衡信号值Igas的90%所需要的时间. 气敏响应测试中, 本文共测试了5种气体, 分别为丙酮(CH3COCH3)、甲醛(HCHO)、甲酸(HCOOH)、乙酸丁酯(CH3COO(CH2)3CH3)和乙醇(C2H5OH). 分别做了最佳工作温度、浓度梯度、长期稳定性测试. 根据之前的文献研究, SnO2传感器的最佳工作温度在200 ~ 280 ℃[24-25]. 将测试回路电压设置为10 V, 在100 ppm的不同气体环境, 150 ~ 350 ℃的温度范围内进行测试. 先在以50 ℃为温度梯度做气敏响应测试, 然后在最佳温度区间内逐渐缩小梯度得出最佳工作温度. 在最佳工作温度区间内, 对5 ~ 400 ppm的不同气体进行了浓度梯度测试, 以获得其选择性. 按照5天的间隔, 进行了为期一个月的长期稳定性测试. 根据静态测试结果评价传感器的性能, 并在周期为80 s, 施加电压3 ~ 7 V的三角波下进行了动态测试分析, 描述静态性能指标与动态响应信号的对应关系.
随后研究了矩形波的周期和占空比改变对动态响应信号的影响, 分别测试了周期30 s占空比50%、周期75 s占空比60%和周期50 s占空比60%的矩形波. 研究了矩形波的温度范围的改变对动态响应信号的影响, 分别使用了90 ~ 350 ℃、120 ~ 350 ℃、150 ~ 350 ℃和180 ~ 350 ℃的矩形波对待测气体进行测试, 研究低温段的改变对动态响应信号的影响; 分别使用了120 ~ 320 ℃、120 ~ 350 ℃和120 ~380 ℃的矩形波对待测气体进行测试, 研究高温段的改变对动态响应信号的影响. 并在最佳波形下收集多组不同种类不同浓度的样本, 通过支持向量机算法定量的分析不同的气体. 说明在之前的测试中, 降低温度范围, 缩短周期并没有影响动态响应信号的品质.
2. 实验结果及讨论
2.1 SnO2传感器的响应机理
SnO2的气敏现象大都采用表面控制型敏感材料的传感机理来解释. 当接触到VOCs时, 会在敏感材料表面发生气体的吸附和脱附现象, 由于在此过程中有电子的交换, 所以纳米材料的电导率也随之发生变化, 其变化大小与传感器温度、气体种类、浓度大小等因素有关. 在SnO2半导体材料表面, 主要存在
${{\rm{O}}}^{2-} $ 、${{\rm{O}}}_2^- $ 和${{\rm{O}}}^- $ 三种吸附态离子. 在室温或温度较低时, 吸附的氧主要以${{\rm{O}}}_2^- $ 为主, 材料表面的氧离子吸附处于动态平衡; 当敏感材料处于温度较高(即处于加热工作状态)的情况下, 吸附态的${{\rm{O}}}_2^- $ 夺得电子转变为${{\rm{O}}}^{2-} $ 和${{\rm{O}}}^- $ . 致使材料表面形成正电荷层, 于是其晶界势垒宽度 L变窄, 敏感材料电导率升高, 电阻下降. 当处于加热工作状态下, 敏感材料表面接触到VOCs时, 气体会与吸附态的氧离子发生反应, 并释放电子使材料表面势垒宽度 L 进一步变窄, 电导率快速上升. 当上述吸附过程达到平衡态时, 敏感材料的表面电阻也随之达到稳定态. 温度调制是打破在静态测试中恒定温度下的氧离子种类的动态平衡, 使其晶界势垒宽度 L成周期性变化, 表面吸附离子形态也会发生周期性变化. 在这个过程中又伴随着测试气体与不同吸附态的氧离子反应, 当氧离子态活性与测试气体匹配时就会产生响应曲线的峰值变化$. $ 2.2 SnO2传感器的静态性能
对传感器进行最佳工作温度、浓度梯度、长期稳定性等的测试是判断传感器性能的优劣以及动态测试前期工作的重要步骤. 因此, 我们首先研究了在150 ~ 350 ℃温度范围内100 ppm不同气体下的传感器传感行为, 以获取最佳工作温度. 图2为丙酮、甲醛、甲酸、乙酸丁酯、乙醇在不同温度下的气敏响应. 得出丙酮、甲醛、甲酸、乙酸丁酯和乙醇的最佳工作温度分别为280 ℃、240 ℃、170 ℃、220 ℃、260 ℃. 上述气体在其最佳工作温度下的静态响应曲线如图3所示. 得出丙酮、甲醛、甲酸、乙酸丁酯和乙醇的响应时间(T90) 分别为116 s、98 s、81 s、111 s、72 s, 因为SnO2没有添加催化剂, 响应时间较长, 但稳定性好, 我们选定250 ℃做了浓度梯度测试, 结果如图4(a)所示. 可以很明显地看出不同的气体在浓度不同的情况下也可以获得相同的气敏响应, 选择性差. 传感器的长期稳定性稳定是传感器的测试以及使用的基础, 我们将100 ppm的不同气体, 以5天为间隔, 在一个月内做了长期稳定性测试结果如图4(b)所示. 在一个月的时间内, 气敏响应差距不大, 显示了SnO2优异的稳定性.
图 2 丙酮、甲醛、甲酸、乙酸丁酯和乙醇在不同温度下的气敏响应Fig. 2 Gas sensing response of acetone, formaldehyde, formic acid, butyl acetate, and ethanol at different temperatures
图 3 丙酮、甲醛、甲酸、乙酸丁酯和乙醇在最佳工作温度下的静态响应曲线Fig. 3 Static response curve of acetone, formaldehyde, formic acid, butyl acetate, and ethanol at the optimum operating temperature
图 4 丙酮、甲醛、甲酸、乙酸丁酯、乙醇在不同浓度下的气敏响应和SnO2的长期稳定性测试Fig. 4 Gas sensitive response of acetone, formaldehyde, formic acid, butyl acetate, and ethanol at different concentrations and long term stability test of SnO22.3 静态响应指标与动态响应信号的关系
从图2和图3中可以看出不同气体的最佳工作温度不同, 这是导致不同气体动态响应信号的最高点出现在不同时刻的重要原因. 静态测试中的气敏响应是在恒定功率作用下经过100 s左右读取的数值. 此时温度已经被恒定的加热到某一固定的值, 敏感材料表面吸附的离子之间的反应也达到一种平衡状态. 但在动态测试中, 实际上我们是通过调节电压来间接调节温度的. 某一时刻的温度并不能利用静态测试中的功率—温度对应关系来衡量, 而是受到前置功率的影响. 例如在三角波中, 当电压增加时, 温度与静态测试中同功率的温度相比偏低, 反之亦然. 对于上述气体, 在施加3 ~ 7 V电压、周期为80 s的三角波时, 动态响应信号如图5所示. 明显地看出对于不同的气体有不同的动态响应信号, 且最高峰出现的时刻不一致. 丙酮、甲醛、甲酸、乙酸丁酯、乙醇的最高峰出现时刻分别为第32 s、28 s、25 s、27 s、29 s. 与之前的最佳工作温度测试相对应, 最佳工作温度越低, 最高峰出现的时刻越早. 在之后的模式识别过程中, 最高点出现时刻不同是定性分类的一个重要依据. 受实验设备的限制, 只能测出静态测试中稳定的加热功率对应的温度值. 但可以确定的是, 在动态测试中, 我们要选取对于待测气体有不同的最佳工作温度的传感器, 然后施加温度缓慢变化的加热波形(例如: 三角波、正弦波等), 就可以读取到最高点出现在不同时刻的动态响应信号.
图 5 在施加3 ~ 7 V电压、周期为80 s的三角波时, 丙酮、甲醛、甲酸、乙酸丁酯和乙醇的动态响应信号Fig. 5 Dynamic response signals of acetone, formaldehyde, formic acid, butyl acetate, and ethanol under the triangular wave with a period of 80 s and a voltage of 3 ~ 7 V图2展示了传感器在不同温度下的气敏响应, 这也与动态响应信号的复杂响应瞬变有着密切的关系. 由于不同气体在不同温度下有不同的气敏响应, 也就是说, 在静态测试时, 不同的温度下吸附气体之后的输出电流不同. 这种现象反应在动态测试中, 由于温度的变化, 导致动态响应信号的波形改变. 如图5所示, 不同气体在动态测试时产生不同的动态响应信号. 但由于静态测试的气敏响应需要一定的反应时间, 动态测试中的响应信号与不同温度的气敏响应值不是一一对应的关系.
为了研究复现性与动态响应信号的关系, 以100 ppm乙醇为例, 在400 ℃下做了复现性测试, 结果如图6所示, 复现性主要的影响因素是传感器的脱附性能. 又知, 对同一个传感器而言, 温度越高, 脱附性越好. 测试高温时而不是最佳工作温度时的复现性, 主要是因为在动态测试时, 周期可重复性的原因除了传感器的稳定之外, 较长的高温区间可以增加吸附脱附速度. 之所以选择矩形波进行动态测试, 也是考虑到了矩形波相对于其他波形而言, 有一个稳定的高温时段, 有利于吸附和脱附. 这也是矩形波经常用来当做加热波形的主要原因.
图 6 100 ppm乙醇在400 ℃下的复现性测试结果Fig. 6 Reproducibility test results of 100 ppm ethanol at 400 ℃2.4 矩形波温度调制参数的选择
在周期和占空比的选择中, 均以100 ppm乙醇为例进行绘图展示结果. 先使用了周期为30 s、占空比为50%、温度范围为150 ~ 350 ℃的矩形波进行温度调制. 动态响应信号如图7(a)所示, 不能发现明显的特征来区分气体. 随后改进加热波形, 每个周期增加30 s的高电平, 15 s的低电平, 总周期为75 s, 占空比为60%, 温度范围仍然为150 ~ 350 ℃. 结果如图7(b)所示, 可以明显地发现特征峰, 但同时也发现由于高温下响应速度快, 动态响应信号在20 s之后便稳定在某一输出电流. 因为在实际应用中, 动态测试的加热周期相当于静态测试的响应时间, 响应时间越短, 可实用性越高, 本着这一原则, 缩短高、低电平中无效的加热时间, 尽可能地缩短周期信号. 将高电平缩短15 s, 低电平缩短10 s, 总周期50 s, 占空比60%, 温度范围为150 ~ 350 ℃的动态响应信号如图7(c)所示. 与图7(b)相比, 特征峰并没有明显的变化. 基于响应时间的考虑, 最终选用周期50 s、占空比60%的矩形波作为加热波形.
图 7 不同周期和占空比的矩形波的动态响应信号Fig. 7 Dynamic response signal of rectangular wave with different period and different duty cycle在工作温度的选择中, 均以100 ppm乙醇为例进行绘图展示结果. 在之前对占空比和周期的调制过程中, 使用的是150 ~ 350 ℃的温度范围. 现研究工作温度范围对动态响应信号的影响, 首先对低温段进行分析. 在周期为50 s、 占空比60%的矩形波的基础上, 分别测试了当低温段为90 ℃、120 ℃、150 ℃、180 ℃, 高温段为350 ℃时的动态响应信号. 结果如图8所示, 90 ℃时低温段特征峰较小, 120 ℃、150 ℃、180 ℃时区别不大. 本着低功耗的原则, 选择120 ℃的低温段进行动态测试. 然后对于高温段的选择, 分别测试了当高温段为320 ℃、350 ℃、380 ℃, 低温段为120 ℃时的动态响应信号. 并且从前文的理论中可知, 高温段的不同, 可能导致周期可重复性的不同, 为此, 我们做了周期重复性的测试. 结果如图9所示, 从波形变化的角度分析, 随着温度的增加, 动态响应信号中的最高峰在增加. 从周期可重复性考虑, 明显地由于温度的降低, 在图9(a)中, 320 ℃时周期可重复性明显不如其他两组. 温度升高, 在图9(c)中, 380 ℃周期可重复性也不会有太大的提高. 基于功耗的考虑, 最终采用的温度范围为120 ~ 350 ℃.
图 8 温度范围分别为90 ~ 350 ℃, 120 ~ 350 ℃, 150 ~ 350 ℃, 180 ~ 350 ℃的矩形波的动态响应信号Fig. 8 Dynamic response signals of rectangular wave with temperature range of 90 ~ 350 ℃, 120 ~ 350 ℃, 150 ~ 350 ℃, and 180 ~ 350 ℃, respectively
图 9 温度范围分别为120 ~ 320 ℃, 120 ~ 350 ℃和120 ~ 380 ℃的矩形波的动态响应周期重复性信号Fig. 9 Dynamic response periodic repetitive signal of rectangular wave with temperature range of 120 ~ 320 ℃, 120 ~ 350 ℃, and 120 ~ 380, respectively2.5 数据获取与模式识别
根据之前的分析, 最终采用了周期为50 s、占空比60%、 温度范围为120 ~ 350 ℃的矩形波进行测试. 5 ppm、10 ppm、50 ppm、100 ppm、200 ppm、300 ppm、400 ppm的丙酮、甲醛、甲酸、乙酸丁酯和乙醇的动态响应梯度信号如图10所示(图中曲线所代表的浓度随箭头方向依次递增). 对于不同的气体之间, 它们响应曲线特征峰出现的位置不同; 对于同种气体在不同浓度下的响应, 它们响应曲线走势相同但有不同的幅值. 特征峰出现的位置以及幅值即为判断气体的种类和浓度的重要特征. 在浓度梯度的测试中, 共获得了831组样本. 由于支持向量机基于不同的数学理论, 对高维数据的识别率较高. 所以并不需要对数据进行特征提取, 降低维度. 将未特征化的数据导入到支持向量机算法中. 为了避免过度拟合, 在支持向量机训练中采用5次交叉验证将数据划分为递归训练集和验证集. 训练识别的准确性可以从平均5个训练结果中获得. 测试结果如图11所示, 识别率达到100%.
图 10 丙酮、甲醛、甲酸、乙酸丁酯和乙醇的动态响应梯度信号Fig. 10 Dynamic response gradient signals of acetone, formaldehyde, formic acid, butyl acetate, and ethanol3. 结束语
采用化学沉淀法制备了SnO2气体传感器. 首先对静态性能指标进行分析, 得出了根据静态测试结果选取适合动态测试传感器的一般特点: 具有较快的静态响应时间、不同的气体有不同的最佳工作温度、不同气体在相同温度下气敏响应相差大、具有良好的复现性等. 并得出了最佳工作温度与动态响应信号的最高峰的出现时刻的隐含关系; 不同温度的气敏响应对动态响应波形的影响; 静态测试中的重复性和动态测试中的周期复现性的对应关系. 通过对矩形波周期、占空比和工作温度范围的调整, 最终选取周期50 s、占空比60%、工作温度范围为120 ~ 350 ℃. 在不降低动态响应信号的品质的情况下, 降低了动态测试在实际应用中的响应时间和功耗. 并使用支持向量机验证了动态响应信号的品质, 当样本量为831 (训练集为731)时, 精度高达100%.
-
图 3 常压蒸馏塔与换热网络协同操作优化框架
Fig. 3 Collaborative operation optimization framework of atmospheric distillation column and heat exchanger network
表 1 IGD指标比较结果
Table 1 The comparison result of IGD
测试函数 RK-CRVEA EI-PoF RVEA CF1 0.0169(0.0016) 0.0231(0.0027) 0.0338(0.0021) CF2 0.0131(0.0032) 0.0305(0.0102) 0.0473(0.0095) CF3 0.1322(0.0128) 0.3909(0.0764) 0.4351(0.0856) CF4 0.0382(0.0290) 0.1445(0.0325) 0.1622(0.0413) CF5 0.1754(0.0457) 0.3517(0.0824) 0.7395(0.0881) 
下载: 导出CSV
表 2 决策变量的取值范围
Table 2 The range limit of decision variables
序号 变量 单位 下限 上限 实际值 1 $F_{prod-1}$, 石脑油流量 bbl/h 430 600 450 2 $F_{prod-2}$, 煤油流量 bbl/h 315.5 500 351.5 3 $F_{prod-3}$, 轻柴油流量 bbl/h 550.1 802.1 789 4 $F_{prod-4}$, 重柴油流量 bbl/h 180.5 387.5 200 5 $F_{stm-1}$, 塔底蒸汽注入量 kg/h 3 000 4 000 3 400 6 $F_{stm-2}$, 第一汽提塔蒸汽注入量 kg/h 1 000 2 000 1 350 7 $F_{stm-3}$, 第二汽提塔蒸汽注入量 kg/h 1 000 2 000 1 150 8 $F_{PA-1}$, 常顶循环回流量 bbl/h 1 883 2 483 2 000 9 $F_{PA-2}$, 常一中循环回流量 bbl/h 950 1 550 1 258 10 $F_{PA-3}$, 常二中循环回流量 bbl/h 950 1 550 1 258 11 $Q_{PA-1}$, 常顶循环回流负荷 MW 10 20 12.5 12 $Q_{PA-2}$, 常一中循环回流负荷 MW 6 16 10.8 13 $Q_{PA-3}$, 常二中循环回流负荷 MW 6 16 10.8 14 $T_{f}$, 常压炉出口温度 ℃ 330 370 350
下载: 导出CSV
表 3 原油、产品及公用工程的价格
Table 3 Prices of crude oil, distillation product and utilities
参数 原油$P_{crude}$ 水蒸气$P_{stm}$ 燃料油$P_{fuel}$ 冷凝水$P_{cw}$ 石脑油$F_{prod-1}$ 煤油$F_{prod-2}$ 轻柴油$F_{prod-3}$ 重柴油$F_{prod-4}$ 价格 79.6 0.0055 0.017 $4.74\times10^{-3}$ 103.5 92.7 99 96.6 单位 ＄/bbl ＄/kg ＄/MJ ＄/MJ ＄/bbl ＄/bbl ＄/bbl ＄/bbl
下载: 导出CSV
表 4 产品TBP 95 %规定范围(℃)
Table 4 The value range of product TBP 95 % (℃)
产品 石脑油$T95_1$ 煤油$T95_2$ 轻柴油$T95_3$ 重柴油$T95_4$ 上限($T95^L$) 95 175 285 345 下限($T95^H$) 120 200 310 370
下载: 导出CSV
表 5 相同机理模型评估次数下优化结果的比较
Table 5 The comparison results under the same mechanism model evaluation times
优化前 机理模型优化 代理模型辅助优化 燃料油消耗量(MW) 49.21 48.20 (-2.0%) 45.12 (-8.3%) 冷凝水消耗量(MW) 41.01 39.23 (-4.4%) 36.52 (-10.9%) 能耗成本(M＄/y) 1.12 0.96 (-14.2%) 0.85 (-24.0%) 产品净收益(M＄/y) 18.51 19.33 (4.3%) 19.92 (7.5%)
下载: 导出CSV
-
[1] Inamdar S V, Gupta K S, Saraf D N. Multi-objective optimization of an industrial crude distillation unit using the elitist non-dominated sorting genetic algorithm. Chemical Engineering Research and Design, 2004, 82(5):611-623 doi: 10.1205/026387604323142667 [2] 丁进良, 杨翠娥, 陈远东, 柴天佑.复杂工业过程智能优化决策系统的现状与展望.自动化学报, 2018, 44(11):1931-1943 http://www.aas.net.cn/CN/abstract/abstract19377.shtmlDing Jin-Liang, Yang Cui-E, Chen Yuan-Dong, Chai Tian-You. Research progress and prospects of intelligent optimization decision making in complex industrial process. Acta Automatica Sinica, 2018, 44(11):1931-1943 http://www.aas.net.cn/CN/abstract/abstract19377.shtml [3] Yang S L, Wang J M, Shi L Y, Tan Y J, Qiao F. Engineering management for high-end equipment intelligent manufacturing. Frontiers of Engineering Management, 2018, 5(4):420-450 doi: 10.15302/J-FEM-2018050 [4] Ochoa-Estopier L M, Jobson M, Smith R. The use of reduced models for design and optimisation of heat-integrated crude oil distillation systems. Energy, 2014, 75(5):5-13 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ad5c0a3e0c9aa507611554b3b641568d [5] Liebmann K, Dhole V R, Jobson M. Integrated design of a conventional crude Oil distillation tower using pinch analysis. Chemical Engineering Research and Design, 1998, 76(3):335-347 doi: 10.1205/026387698524767 [6] Ding J L, Yang C E, Jin Y C, Chai T Y. Generalized multi-tasking for evolutionary optimization of expensive problems. IEEE Transactions on Evolutionary Computation, 2019, 23(1):44-58 doi: 10.1109/TEVC.2017.2785351 [7] 丁进良, 杨翠娥, 陈立鹏, 柴天佑.基于参考点预测的动态多目标优化算法.自动化学报, 2017, 43(2):313-320 http://www.aas.net.cn/CN/abstract/abstract19009.shtmlDing Jin-Liang, Yang Cui-E, Chen Li-Peng, Chai Tian-You. Dynamic multi-objective optimization algorithm based on reference point prediction. Acta Automatica Sinica, 2017, 43(2):313-320 http://www.aas.net.cn/CN/abstract/abstract19009.shtml [8] Ochoa-Estopier L M, Jobson M. Optimization of heat-integrated crude oil distillation systems. Part I:The Distillation Model. Industrial & Engineering Chemistry Research, 2015, 54(18):4988-5000 http://cn.bing.com/academic/profile?id=3fc6aa9e00cdc136b2d7c0587bc51a70&encoded=0&v=paper_preview&mkt=zh-cn [9] Ochoa-Estopier L M, Jobson M, Chen L, Rodr-guez-Forero C A, Smith R. Optimization of heat-integrated crude oil distillation systems. Part Ⅱ:Heat exchanger network retrofit model. Industrial & Engineering Chemistry Research, 2015, 54(18):5001-5017 http://cn.bing.com/academic/profile?id=73b74f4f04c094c1b2ea88d07a388664&encoded=0&v=paper_preview&mkt=zh-cn [10] Ochoa-Estopier L M, Jobson M. Optimization of heat-integrated crude oil distillation systems. Part Ⅲ:Optimization Framework. Industrial & Engineering Chemistry Research, 2015, 54(18):5018-5036 http://cn.bing.com/academic/profile?id=b625acd749d31e770c36c08f14430946&encoded=0&v=paper_preview&mkt=zh-cn [11] Al-Mayyahi M A, Hoadley A F A, Smith N E, Rangaiah G P. Investigating the trade-off between operating revenue and CO2 emissions from crude oil distillation using a blend of two crudes. Fuel, 2011, 90(12):3577-3585 doi: 10.1016/j.fuel.2010.12.043 [12] Ji S, Bagajewicz M. Design of crude distillation plants with vacuum units. I. targeting. Industrial & Engineering Chemistry Research, 2002, 54(24):6094-6099 http://cn.bing.com/academic/profile?id=cd8338396883c39aa0ba3ba3017ae8c6&encoded=0&v=paper_preview&mkt=zh-cn [13] Lopez C D C, Hoyos L J, Mahecha C A, Arellano-Garcia H, Wozny G. Optimization model of crude oil distillation units for optimal crude oil blending and operating conditions. Industrial & Engineering Chemistry Research, 2013, 52(36):12993-13005 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1d2d8046aaf44c1aaae574764daad9e9 [14] Ibrahim D, Jobson M, Li J, Guillén-Gosálbez G. Optimization-based design of crude oil distillation units using surrogate column models and a support vector machine. Chemical Engineering Research and Design, 2018, 134:212-225 doi: 10.1016/j.cherd.2018.03.006 [15] Yao H H, Chu J Z. Operational optimization of a simulated atmospheric distillation column using support vector regression models and information analysis. Chemical Engineering Research and Design, 2012, 90(12):2247-2261 doi: 10.1016/j.cherd.2012.06.001 [16] Cheng R, Jin Y C, Olhofer M, Sendhoff B. A reference vector guided evolutionary algorithm for many-objective optimization. IEEE Transactions on Evolutionary Computation, 2016, 40(5):773-791 http://cn.bing.com/academic/profile?id=a78f88991f475d68f269707178b2f22e&encoded=0&v=paper_preview&mkt=zh-cn [17] Jones D R, Schonlau M, Welch W J. Efficient global optimization of expensive black-box functions. Journal of Global Optimization, 1998, 13(4):455-492 doi: 10.1023/A:1008306431147 [18] Martínez-Frutos J, Herrero-Pérez D. Kriging based infill sampling criterion for constraint handling in multi-objective optimization. Journal of Global Optimization, 2016, 64(1):1-19 doi: 10.1007/s10898-015-0328-x [19] Guo D, Jin Y C, Ding J L, Chai T Y. Heterogeneous ensemble based infill criterion for evolutionary multi-objective optimization of expensive problems. IEEE Transactions on Cybernetics, 2019, 49(3):1012-1025 doi: 10.1109/TCYB.6221036 [20] Seiffert C, Khoshgoftaar T M, Van Hulse J, Napolitano A. RUSBoost:a hybrid approach to alleviating class imbalance. IEEE Transactions on Systems, Man, and Cybernetics-Part A:Systems and Humans, 2010, 40(1):185-197 doi: 10.1109/TSMCA.2009.2029559 [21] Wang J H, Liang G X, Zhang J. Cooperative differential evolution framework for constrained multiobjective optimization. IEEE Transactions on Cybernetics, 2019, 49(6):2060-2072 doi: 10.1109/TCYB.2018.2819208 [22] Zhang Q F, Zhou A M, Jin Y C. RM-MEDA:A regularity model-based multi-objective estimation of distribution algorithm. IEEE Transactions on Evolutionary Computation, 2008, 12(1):41-63 doi: 10.1109/TEVC.2007.894202 [23] Zitzler E, Thiele L, Laumanns M, Fonseca C M, Fonseca V G D. Performance assessment of multiobjective optimizers:an analysis and review. IEEE Transactions on Evolutionary Computation, 2003, 7(2):117-132 doi: 10.1109/TEVC.2003.810758 期刊类型引用(2)
1. 王磊，刘爱国，张涛，马崇伟，苗光尧. 多台电压互感器应急装置协同报警方法研究. 自动化与仪器仪表. 2020(11): 216-219 . 
百度学术2. 周晓君，阳春华，桂卫华. 状态转移算法原理与应用. 自动化学报. 2020(11): 2260-2274 . 本站查看其他类型引用(4)
-
下载:
计量
- 文章访问数: 2010
- HTML全文浏览量: 370
- PDF下载量: 139
- 被引次数: 6
