精馏塔优化设计.ppt

《精馏塔优化设计.ppt》由会员分享，可在线阅读，更多相关《精馏塔优化设计.ppt（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、精馏塔优化设计,阮 奇（福州大学化学化工学院） 1 目标函数 2 工艺计算 3 目标函数的求解,1 目标函数,式中 J 精馏塔年总费用，元/年； J1 精馏塔主体年投资折旧费及维修费用，元/年； J2 冷凝器的年折旧及操作费用 ，元/年； J3 再沸器的年折旧及操作费用 ，元/年。,（1）,式中 精馏塔年维修费用比率； 年折旧率；一般可取折旧年限 8 10， 塔板费用，元； 塔体费用，元。,1.1 精馏塔主体年投资折旧费及维修费用 J1,（2）,1.1.1 塔板费用CT（碳钢材料）,1.1.2 塔体费用CH,塔质量，kg；（化工系统工程基础） 塔高，m； 塔壁厚，m；常压 ，可取6mm； 碳钢

2、密度， ； 板间距，m； 塔两端高度，m；包括塔顶空间、塔底裙 座。（常用化工单元设备的设计）,1.1.2.1 计算方法一,其中式中,，,（4）,塔体费用，元/ ；可取4800； 塔截面积， ； 。,1.1.2 塔体费用CH,1.1.2.2 计算方法二,式中,（5）,1.2 冷凝器年运转费用J2,冷却水价格，元/kg；可取0.0004； 冷却水用量，kg/h； 年工作时间，h/a； 冷凝器价格回归系数；取487； 冷凝器价格回归指数；取0.72；冷凝器压力校正系数；取1.0； 冷凝器材质校正系数；取6.5； 冷凝器传热面积， 。,式中,（6）,1.3 再沸器年运转费用 J3,1.3.1 间接蒸

3、汽加热,加热蒸汽价格参数，元/kg；可取0.018；加热蒸汽价格参数，元/ ；可取0.0001；加热蒸汽饱和温度， ； 加热蒸汽流量，kg/h； 再沸器价格回归系数；取415； 再沸器价格回归指数；取0.74； 再沸器压力校正系数；取0.94；再沸器材质校正系数；取5.35； 再沸器传热面积，m2 。,式中,（7）,1.3 再沸器年运转费用J3,1.3.2 直接蒸汽加热,（8）,2 工艺计算,2.1 精馏塔主体年投资折旧费及维修费用J1的计算2.1.1 塔径DT的计算,（9）,2.1.1 塔径DT的计算,空塔气速则由液泛气速求得：,2.1.1 塔径DT的计算,注意点：（1）各参数均以第一块塔板

4、（温度，汽液相组成）为基础进行计算；（2）不同的塔板形式， 图也不同；（3）表面张力 的混合规则为： （A 、B为两组分）；（4）由于优化设计结果还未得到， 未得到；所以在计算以上参数时，回流比可以在 区间内选择一适宜值，它对优化设计结果没有影响，因为 的范围较大。,2.1.2 总板效率ET的计算,式中,（10）,（11）,2.1.3 理论塔板数的计算,2.1.3.1 相平衡关系的表示,对理想溶液，其相平衡关系为：,对非理想溶液，其相平衡关系可以从实验数据中利用三次样条插值得到或通过回归实验数据得到相平衡关系表达式，其形式主要有：,（余国琮式）,（阮奇式）,， 分别为液相和汽相的摩尔组成； ，

5、 相平衡关系回归系数。,式中,注意点： （1） 时，更换操作线方程，注意区别直接与间接蒸汽加热时的提馏段操作线方程的不同； （2）若相平衡关系是分段表示的，则必须判断汽相组成所在的区间来选择相平衡关系式； （3）总理论塔板数为 ，不对理论板数进行圆整是为了保证回流比和塔板数一一对应的关系。 由2.1点的工艺分析可知精馏塔年投资折旧费和维修费是回流比的函数，即 。,2.1.3.2 N的计算,全凝器,2.2 冷凝器年运转费用J2的计算,2.2.1 冷却水用量W及冷凝器传热面积AD的计算,式中,（12）,2.2.1 冷却水用量W及冷凝器传热面积AD的计算,由式（12）可得：,将式（13）和（14）代

6、入式（6）得：,（13）,（14）,（15）,2.2.2 冷凝器冷却水最佳出口温度 的确定,，,其中,为 的解，因为 为非线性方程，无法直接求解，故应采用数值方法，如牛顿迭代法（需对 求导）。,2.2.2 冷凝器冷却水最佳出口温度 的确定,2.3 再沸器年运转费用J3的计算,（16）,2.3.1.1 蒸汽流量Z及再沸器传热面积AB的计算,由式（16）得：,将式（17）和（18）代入式（7）得到：,（17）,（18）,（19）,2.3.1.2 再沸器加热蒸汽最佳温度 的确定,2.3.1.2 再沸器加热蒸汽最佳温度 的确定,为 的解，因为 为非线性方程，无法直接求解，估应采用数值方法，如牛顿迭代法

7、（需 对 求导）。,2.3.2 直接蒸汽加热,式（8）中的Z 为加热蒸汽的流量，所以：,，kg/h,，kmol/h,（20）,3 目标函数的求解,3.1 求解过程,3.2 函数调用关系,精馏塔优化设计所需的若干拟合式,smith关联图C20,，,式中,精馏塔优化设计所需的若干拟合式,，,规整填料压降通用关联图,式中,，,a1=110.356，a2=215.838，a3=173.321，a4=48.72，a5=81.232，a6=391.67，a7=604.92，a8=288.74，a9=9.4368，a10=5.24，a11=130.45，a12=172.24，a13=20.152，a14=2

8、.5515E-03，a15=5.6332，a16=7.4636E-02，a17=5.4214,精馏塔优化设计所需的若干拟合式,新填料通用关联图,式中,，,a1=374.43，a2=950，a3=553.47，a4=273.817，a5=299.96，a6=3513.6，a7=9.582，a8=1303.9，a9=344.11，a10=0.90973，a11=2060.9，a12=9.9448E-02，a13=2049.5,三次样条插值,已知 、 、 及对应函数值 、 、 （ ），要求 对应的 值。 1 线性插值 当 时，将 间函数视为直线（用直线拟合 函数）。 对两点求直线：,三次样条插值,2

9、 抛物线插值，二次插值 ，用抛物线拟合 、 、 三点，构造二次多项式，而后插值求 。 这两种方法线性插值误差较大，抛物线插值误差相对较小，但对于节点处难以光滑。,三次样条插值,3 三次样条插值 三次样条插值是一种分段插值，既能保证结点处曲线连续，又能保证其光滑的一种方法。 构造的三次样条插值函数必须满足： （1） （ ）； （2） 在插值区间 上具有一阶导数（保证曲线光滑）与二阶导数（一阶导数曲线光滑）； （3）在每个插值区间 上， 均为二次多项式。 三次样条插值适用于三种边界条件： （1）已知函数 在各结点 、 、 上的函数值 、 、 ，且已知端点的一阶导数值 、 ； （2）已知函数 在各结

10、点 、 、 上的函数值 、 、 ，且已知端点的二阶导数值 、 ；,三次样条插值,（3）已知函数 在各结点 、 、 上的函数值 、 、 ，且已知端点的函数值、一阶导数值、二阶导数值相等，即： 对汽液相平衡数据，其端点处导数值未知，故没有特殊要求时令端点处 、 以保证插值函数的光滑，所以应采用第二类边界条件。其插值函数 为 ， 为插值点， 为插值结果即为 。 3.1 插值公式 插值点：,；,；,三次样条插值,插值函数：,三次样条插值,插值点上的一阶导数：,三次样条插值,插值点上的二阶导数：,三次样条插值,其中 、 （ ）为实验点（节点），而 （ ）为：,三次样条插值,；,三次样条插值,3.2 插值

11、步骤 （1）输入 、 、 （2）计算 、 、 （3）计算 、 （4）计算 、 （5）计算 （6）计算 （7）计算,微粒群（PSO）算法,式中：,：表示微粒i在d 维中第k次迭代时的速度。：表示微粒i在d 维中第k次迭代时的位置。：惯性权重，表示微粒沿自己原来方向移动的几率。 ： 是一个常数， 是一个0-1之间的随机数，他们的积表示微粒向其经过的最优位置移动的几率。 ： 是一个常数， 是一个0-1之间的随机数，他们的积表示微粒向整个微粒群最优位置移动的几率。：微粒i所经过的最优位置在d 维中的值。：整个微粒群的最优位置在d 维中的值。,（1）,（2）,微粒群（PSO）算法,标准PSO的算法流程如

12、下: (1)随机选取一个微粒群(群体规模为m)，对每个微粒随机给定一个初始位置和初始速度。 (2)计算每个微粒在初始位置时的适应度，并把它赋给 ，选择适应度最大的一个微粒的适应度为所有微粒的最优适应度，把它赋给 ，并记录该微粒的位置。 (3)根据式（1）和式（2）调整微粒的位置。 (4)重新计算每个微粒在初始位置时的适应度，并把它与微粒的 比较，如果当前值大于 ，则把此值赋给 。把所有微粒的 与 比较，如果有微粒的 大于 ，则把此值赋给 ，并记录该微粒的位置。,微粒群（PSO）算法,(5)如未达到收敛条件(通常以一个最大的迭代代数和一个收敛精度为判断依据，满足其中任一条件即可)则返回第（3）步

13、。,在实际的优化问题中，微粒的位置代表所有变量的有效取值空间，而微粒的适应度代表变量取不同的值所对应的目标函数的值。,微粒群算法中参数取值：一般取 ， 取值在0.4-0.8之间，也可按以下公式进行取值:,其中 表示迭代次数， 表示预定的最大迭代次数。,附录：,真空费用模型： 真空系统的动力费用主要是液环真空泵的运行费用也即电费。对于液环真空泵的运行费用的求法未见报道，故参考空气压缩机的运行费用计算液环真空泵的运行费用，则J61可用下式计算,式中 c1,电费，取c10.4元/（kWh）；,液环真空泵的轴功率，kW；,操作时间，h/a；,附录：,1 谭天恩，麦本熙，丁惠华化工原理第二版下册北京化学工业出版社1999,液环真空泵的轴功率 用下式计算:,