灵敏度分析报告aspen
Aspen教程-4高级功能
ASPEN PLUS 高级功能报告人：宋维仁导师：漆志文教授华东理工大学联合化学反应工程研究所2010年6月30日星期三1主要内容灵敏度分析设计规定工况研究Windows应用程序间的交互操作功能2本页已使用福昕阅读器进行编辑。
福昕软件（Ｃ）２００５－２０１０，版权所有，仅供试用。
灵敏度分析灵敏度分析模块是用来分析和研究一个或多个流程操作变量的改变对其它流程变量的影响，它是做工况研究的一个最有用的工具，你可以用灵敏度分析模块来生成随进料物流、模块输入参数或其它输入参数的变化的模拟结果表和图，还可以用它来验证一个设计规定的解是否在操作变量的变化范围内，你甚至还可以用它来做简单的过程优化。
下面我们就以回收废水中的DMF为例来说明相关操作。
3灵敏度分析例子单击此按钮，进入Data Browser窗口。
下面这个例子将演示如何使用灵敏度分析功能来考察精馏塔column的进料位置对塔釜物流W中DMF的质量分数的影响。
4新建灵敏度分析ID点击Model Analysis Tools中的Sensitivity，点击右边的New按钮，出现如右图的对话框（即新建一个灵敏度分析任务的名称，可以采用默认名称，也可重新命名）。
5命名因变量在S-1-Input-Define页中，单击New，新建一个因变量，并且输入这个因变量的名称。
在本例中因变量就是塔釜出料W中DMF的质量分数。
6定义因变量3.Type选择质量分数Mass Frac2.选择物流变量4.物流选择W，组分选择DMF。
1.此页用于定义因变量的类型和内容。
5.单击，进入下一操作。
7定义自变量在S-1-Input-Vary页中，定义自变量的类型和内容。
本例中自变量为精馏塔Column的进料位置。
输入自变量的考察范围和步长。
8分析结果表输入在此页中定义将要得到的灵敏度分析结果表格中要显示的内容。
在第一列中输入数字1，后面输入要列表的因变量名称或因变量的表达式。
至此，灵敏度分析模块已输入完毕，点击运行按钮。
Aspen精馏模拟灵敏度分析与优化实例use
（3）设计实例
常压操作连续筛板精馏塔设计，设计参数如下[1]： 进料组份：水 63.2%、甲醇 38.6%（质量分率）； 处理量：水甲醇混合液 55t/h； 进料热状态：饱和液相进料； 进料压力：125 kPa； 操作压力：110 kPa； 单板压降：≤0.7 kPa； 塔顶馏出液：甲醇量大于 99.5 %（质量分率） 塔底釜液： 水量大于 99.5 %；（质量分率）. 回流比：自选； 全塔效率：ET=52% 热源：低压饱和水蒸汽；
Aspen plus 精馏模拟 灵敏度分析与优化 实例教程 (甲醇-水系统)
仇汝臣 青岛科技大学化工学院
13376489639 8978122@ QQ群：182452945
北京中技培咨询服务中心
第 1 页共 37 页
Aspen plus 精馏模拟灵敏度分析与优化实例教程
1. Aspen Plus 简介
快速访问按钮类似. ·Next 按钮（N->）：设计过程的任意时刻点击它，系统都会自动跳转到当前应当进行的工
作位置，这为我们输入数据提供了极大的方便.
第 2 页共 37 页
2 Aspen Plus 模拟精馏简介
（1）塔模型分类
做塔新流程模拟分析必须先进行简捷塔计算--- 塔的初步设计. 计算结果为理论板数、 进料位置、最小回流比、塔顶/釜热负荷. 然后进行塔精确模拟分析，简捷塔计算结果做为精 确计算的输入依据. 本文以甲醇-水混合物系分离为例，首先介绍初步设计方法，然后介绍复 杂塔模拟计算。为初学者提供帮助。
Dstwu 要求定义组分的份的回收率. 计算得到两种组分的回收率为： 轻关键组分的回收率为 0.9983 重关键组分的回收率为 0.0029 ·定义再沸器和冷凝器的压力（Pressure） ·定义冷凝器类型（Condenser specifications） 选择全冷器（Total Condenser）
ASPEN
塔顶压力 进料压力 塔底压力 压降
10kpa 30 kpa 38.8kpa 0.8kpa
Page 13
回顾最初遇到的问题
输入：操作回流比为４ ** ERROR SPECIFIED REFLUX RATIO (4.0000) IS LESS THAN THE MINIMUM REFLUX RATIO (4.7118). 2*MINIMUM USED TO CALCULATE NUMBER OF STAGES.
1.5045e-09 3.1655e-05
ETHYL-01 STYRE-01
0.5 0.4
0.73326308 0.06473383
0.27139001 0.72857833
未能达到要求，需优化
Page 7
对压降进行灵敏度分析
0 得表格： ok 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 0.75516479 0.7551648 0.75516478 0.7551648 0.5 0.7551646 0.0428014 2 0.0428012 3 0.0428013 5 0.0428012 3 0.0428013 5
4.71181478 7 19.377901 30.969345 19.6928302 18.6928302 0.35789982 MMkcal/hr 0.35661627 MMkcal/hr 125.6367 C 147.574033 C 0.54499987
Actual reflux ratio: Minimum number of stages: Number of actual stages: Feed stage: Number of actual stages above feed: Reboiler heating required: Condenser cooling required: Distillate temperature: Bottom temperature: Distillate to feed fraction:
甲醇装置预精馏塔Aspen模拟任务书
甲醇装置预精馏塔Aspen模拟任务书一、模拟计算依据：1、原料处理量：学号后三位XXX × 100 kg/h；2、粗甲醇液进料组成如表1所示(质量分数)；进料条件为：液相进料温度60℃，进料压力140kPa，塔顶（分凝器气相出料）冷凝器压力130kPa，再沸器压力150kPa；3、分离要求：塔顶甲酸甲酯摩尔回收率为99.99%，塔顶甲醇摩尔回收率为0.7%。
4、物性方法：BWRS表1 进料组成表二、任务1、按计算依据，用简捷法（DSTWU模块）模拟计算预精馏塔以分离粗甲醇中的轻组分（建议实际回流比取最小回流比的1.5倍）。
2、在简捷模拟计算中，通过回流比随理论板数变化曲线，确定适宜回流比、理论板数。
及相应的进料位置、塔顶产品与进料的摩尔流量比（D/F)、最小回流比、最小理论板数、实际理论板数、进料位置以及塔顶温度。
3、根据简捷计算的结果，利用严格法（RadFrac模块）对预精馏塔进行严格计算，进料条件、冷凝器形式、冷凝器压力、再沸器压力、再沸器采用釜式再沸器、产品纯度要求以及物性方法与简捷法相同，用严格法核算任务2中的结果（简捷计算结果）是否达到回收率要求。
4、通过严格法（RadFrac模块）设计规定功能，调整回流比、馏出与进料量比以达到分离要求；5、通过Aspen灵敏度分析功能，在严格法中求取回流比随理论板数据的变化曲线，重新确定适宜回流比、理论板数。
6、绘制塔内温度分布曲线、塔内液相质量组成分布曲线、塔内的气相组成分布曲线。
7、书写模拟报告。
以下为选做部分（评优学生必做）6-1、通过Aspen灵敏度分析功能，在严格法中求取进料板位置与再沸器热负荷的关系曲线，重新确定进料板位置。
6-2、设实际塔板的塔板默弗里效率为60%，在严格法中重新设定塔板数、进料板位置；然后在严格法中初步设定塔板类型为浮阀，查看塔板的水力学性质；6-3、对塔进行校核计算，确定塔的结构尺寸、水力学性能、负荷性能。
Aspen精馏模拟灵敏度分析与优化实例use
第 3 页共 37 页
3. 精馏塔的简捷计算
·设计任务 确定理论塔板数 确定合适的回流比
·DSTWU 精馏模型简介
本例选择 DSTWU 简捷精馏计算模型. DSTWU 可对一个带有分凝器或全凝器一股进料和两种产品的蒸馏塔进行简捷精馏 计算. DSTWU 假设恒定的摩尔溢流量和恒定的相对挥发度
进入 Aspen Plus 后，出现图 1 所示的 Aspen Plus 软件操作界面.
图1
操作界面构成 ·标题条：在该栏目中显示运行标识. 在你给出运行名字之前，Simulation1 是缺省的标识. ·拉式菜单：Aspen Plus 的功能菜单. 这些下拉式菜单与 Windows 的标准菜单类似. ·工艺流程窗口：在该窗口中可以建立及连接所要模拟的工艺流程. ·模式选择按钮：按下此按钮你可以关闭插入对象的插入模式，并返回到选择模式. ·模型库：在这里列出建立模型可用的任何单元操作的模型.. ·状态域：显示当前有关运行的状态信息. ·快速访问按钮：快速执行 Aspen Plus 相应的命令。这些快捷按钮与其它 Windows 程序的
图 3.2-1
图 3.2-2
3.3 定义组份
本节任务： ·输入物料化学成份
单击 N-> 快捷键直到进入进料参数输入页，如图 3.3-1 所示.
第 7 页共 37 页
图 3.3-1
在此窗口中，我们可以定义流程中所涉及的化学组分. 定义方法有两种： 1) 可以在 component ID 或 component name 中直接输入组分的英文名称. 其中 Component ID 是该组分的代号，用户可以进行定义和修改. 2) 可以使用 Aspen plus 提供的 find 工具，查找 Aspen plus 提供的组份. 单击 find 按钮， 进入组份查找页，在对话框中输入组分的英文名称或分子式，也可以输入其部分字符串. 这 里我们输入甲醇分子式 CH4O（注意不能输入 CH3OH），点击 find now 按钮，查找结果出现 在下面列表中，如图 3.3-2.
(完整版)ASPEN错误及解决方法
ZERO FEED - BLOCK BYPASSED.Block:B17 Model: SEP* WARNINGZERO FEED TO THE BLOCK. BLOCK BYPASSEDBlock:B21 Model: HEATER* WARNINGZERO FEED TO THE BLOCK. BLOCK BYPASSEDBlock:B2 Model: MIXER答：此错误为该模块的进料为零，原因为前几个模块中出现零物流输出，建议检查流程中没有物料的流股。
若没有的话，设定撕裂流股，run- reconcile steam -tear steamWARNING IN THE "STREAM" PARAGRAPH WHICH BEGINS ONLINE 45 STREAM NAME: 1 COMPONENT MASS FLOWS OF SUBSTREAM: "MIXED" ARENORMALIZEDTO THE TOTAL MASS FLOW VALUE.答：在定义“Streams”（流股）的时候，原本应该是想按照质量分数或者摩尔分数来设置的，但是没有将“Composition”里默认的“Mole Flow”改为“Mole-Frac”或者“Mass-Frac”。
改过来就好了！ERRORDESIGN SPEC IS NOT SATISFIED BECAUSE ONE OR MORE MANIPULATEDVARIABLE IS AT ITS BOUND.答：有可能按你设置的distillate RATE无法满足设计规定的目标值两种情况：1.设计规定设置的参数有问题导致内部迭代计算振荡发散无法收敛2.塔规定的参数无法达到设计规定的目标值（由两物质的物质、塔操作参数决定）建议：在塔规定的参数下，采用灵敏度分析，这样就可以看出参数A对参数B的变化趋势COMPONENT "YIXI" HAS ZERO FLOWRATE IN THE OUTLET STREAM. FRACTIONALCONVERSION HAS BEEN MODIFIED BY A FACTOR OF "0.7407407" 答：你设定的转化率有问题，出口没流量，好像软件自己调整了，这应该是个警告，不是错误，可以根据需要自己调Block: B1 Model: RADFRAC***SEVERE ERRORALIQUID FEED/PUMPAROUND TO THE TOP STAGE IS REQUIRED WHEN"Q1=0"IS SPECIFIED IN "COL-SPECS"; REQUIRED FEED/PUMPAROUNDHASZERO FLOW.RESTOF BLOCK BYPASSED DUE TO SEVERE ERROR答：当在"冷凝器设定"中，塔顶冷凝器热负荷指定为0时，要求必须有顶层塔板液相进料或中段回流。
2(2)灵敏度分析
c j→ CB 0 2 1 0 2 基 x3 x1 x2 cj-zj x3 x1 15 5 b 35/2 11/2 -1/2
2 x1 0 1 0 0 0 1
1 x2 0 0 1 0 5 1
0 x3 1 0 0 0 1 0
0 x4 5/4 1/4 [-1/4] -1/4 0 0
0 x5 -15/2 -1/2 3/2 -1/2 0 1
-7
0 [2] 1 0 0 1 0
-1/2 0
最优生产计划应为每天生产7/2件家电Ⅰ， 51/4件家电Ⅲ。
分析参数aij的变化
灵 敏 度 分 析 举 例
例 在美佳公司的例子中，若家电Ⅱ每件需设备A,B和 调试工时变为8h、4h、1h，该产品的利润变为3元/件， 试重新确定该公司最优生产计划。
设生产工时变化后的新家电Ⅱ的生产量为x2′，其中：
（2）若家电Ⅰ的利润不变，则家电Ⅱ的利润在什 么范围内变化时，该公司的最优生产计划将不发 生变化？ 设家电Ⅱ的利润为（1+λ）元，如下
项目 CB 基 b 2 x1 1+λ x2 0 x3 0 x4 0 x5
0
2 1+λ
x3
x1 x2 cj-zj
15/2
7/2 3/2
0
1 0 0
0
0 1 0
1
0 0 0
15 / 2 1 / 2 3/2 3 7 4 0 2 2
1 P 6 0 0
5/4 1/ 4 1 / 4
cj→ CB 基 b
2 x1
1 x2
0 x3
0 x4
0 x5
3 x6
灵 敏 度 分 析 举 例
灵敏度分析 aspen
第八章灵敏度分析目得：介绍灵敏度分析得用法，研究过程变量之间得关系。
（1）灵敏度分析●可使用户研究输入变量得变化对过程输出得影响●在灵敏度模块文件夹得Results表上能够查瞧结果●可以把结果绘制成曲线，使不同变量之间得关系更加形象化●在灵敏度模块中对流程输入量所做得改变不会影响模拟，灵敏度研究独立于基础工况模拟而运行●位于/Data/Model Analysis Tools/Sensitivity下（2）灵敏度分析得用法●研究输入变量得变化对过程（模型）得影响●用图表表示输入变量得影响●核实设计规定得解就是否可行●初步优化●用准稳态方法研究时间变化变量（3）灵敏度分析应用步骤a）定义被测量（采集）变量-它们就是在模拟中计算得参量，在第4步将要用到（Sensitivity Input Define页）b）定义被操作（改变得）变量-它们就是要改变得流程变量（Sensitivity Input Vary页）c）定义被操作（改变得）变量范围-被操作变量得变化可以按在一个间隔内等距点或变量值列表来规定（Sensitivity Input Vary页）d）规定要计算得或要制成表得参量-制表参量可以就是任何合法得Fortran表达式，表达式含有步骤1中定义得变量(Sensitivity Input Tabulate页)（4）绘图a）选择包括X轴变量得列，然后选择从Plot菜单下选择X-Axis变量b）选择包括Y轴变量得列，然后选择从Plot菜单下选择Y-Axis变量c）（可选得）选择含有参数变量得列，然后从Plot菜单下选择参数变量d）从Plot菜单下选择Display Plot»要选择一列，用鼠标左键点击列标题（5）注意●只有被输入到流程中得参量才可以被改变或操作●可以改变多个输入●对于每一个被操作（改变得）变量得组合都运行一次模拟（6）灵敏度分析举例以第二章中苯与丙烯为原料合成异丙基苯为例，如下图：冷却器出口温度怎样影响产品物流纯度得？●被调节(被改变)变量就是什么？冷凝器出口温度●被测量(采集)变量就是什么？产品物流中异丙基苯纯度（摩尔分率）打开文件cumene、bkp，另存为cumene-s、bkp，如下图所示：在数据浏览窗口中，点击Model Analysis Tools/Sensitivity，如下图，点击N ew…创建一个新得灵敏度分析：点击New、、、按扭输入创建得灵敏度分析得ID，可以自己指定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第八章灵敏度分析
目的：介绍灵敏度分析的用法，研究过程变量之间的关系。
（1）灵敏度分析
●可使用户研究输入变量的变化对过程输出的影响
●在灵敏度模块文件夹的Results表上能够查看结果
●可以把结果绘制成曲线，使不同变量之间的关系更加形象化
●在灵敏度模块中对流程输入量所做的改变不会影响模拟，灵敏度研究独立于基
础工况模拟而运行
●位于/Data/Model Analysis Tools/Sensitivity下
（2）灵敏度分析的用法
●研究输入变量的变化对过程（模型）的影响
●用图表表示输入变量的影响
●核实设计规定的解是否可行
●初步优化
●用准稳态方法研究时间变化变量
（3）灵敏度分析应用步骤
a）定义被测量（采集）变量
-它们是在模拟中计算的参量，在第4步将要用到（Sensitivity Input Define页）b）定义被操作（改变的）变量
-它们是要改变的流程变量（Sensitivity Input Vary页）
c）定义被操作（改变的）变量围
-被操作变量的变化可以按在一个间隔等距点或变量值列表来规定（Sensitivity
Input Vary页）
d）规定要计算的或要制成表的参量
-制表参量可以是任何合法的Fortran表达式，表达式含有步骤1中定义的变量
(Sensitivity Input Tabulate页)
（4）绘图
a）选择包括X轴变量的列，然后选择从Plot菜单下选择X-Axis变量
b）选择包括Y轴变量的列，然后选择从Plot菜单下选择Y-Axis变量
c）（可选的）选择含有参数变量的列，然后从Plot菜单下选择参数变量
d）从Plot菜单下选择Display Plot
»要选择一列，用鼠标左键点击列标题
（5）注意
●只有被输入到流程中的参量才可以被改变或操作
●可以改变多个输入
●对于每一个被操作（改变的）变量的组合都运行一次模拟（6）灵敏度分析举例
以第二章中苯和丙烯为原料合成异丙基苯为例，如下图：
冷却器出口温度怎样影响产品物流纯度的？
●被调节(被改变)变量是什么？
冷凝器出口温度
●被测量(采集)变量是什么？
产品物流中异丙基苯纯度（摩尔分率）
打开文件cumene.bkp，另存为cumene-s.bkp，如下图所示：
在数据浏览窗口中，点击Model Analysis Tools/Sensitivity，如下图，点击N ew…创建一个新的灵敏度分析：
点击New...按扭
输入创建的灵敏度分析的ID，可以自己指定。
点击N ew…添加新的Fortran variable（Fortran变量）。
点击New...按扭
输入变量的名称，本例中输入PURITY，意指产品物流中异丙基苯纯度。
变量定义选项。
种类（Category）、类型（Type）、物流（Stream）、组分（Component）分别选物流（Streams）、摩尔分率（Mole-Frac）、产品（PRODUCT）、异丙基苯（CUMENE）。
点击Close按扭
设定Vary选项。
块变量（Block-Var）类型，冷凝器模块，温度变量，围100-300，步长5。
设定Tabulate选项，如下：
打开控制面板，运行如下：
查看结果如下：
选中冷凝器温度一列（第二列），点击下拉菜单Plot/X-Axis Variable，如下图：
选中产品纯度一列（第三列），点击下拉菜单Plot/Y-Axis Variable，如下图：
点击下拉菜单Plot/Display Plot，如下图：
得到不同冷凝器温度下产品物流中异丙基苯纯度曲线。