Aspen Plus 做氧化法制备 1,3-丙二醇这个题目我从接到任务到跑通全流程前前后后折腾了将近三周。中间踩过的坑、返工重来的次数说实话比过去做过的很多模拟项目都多。但做完之后收获也大这套流程涉及气固催化氧化、水合反应、高压加氢、多塔精制几乎把 Aspen Plus 里常用的反应器模块、塔模块、物性方法和收敛技巧都用上了。这篇就把整个模拟过程的思路梳理一遍。不只是贴图贴步骤更想把每一步选择的“为什么”讲清楚。1,3-丙二醇1,3-PDO是聚对苯二甲酸丙二醇酯PTT的关键单体这个市场懂行的都懂而目前实现工业化的路线主要有环氧乙烷法和丙烯醛法。丙烯醛法的第一步就是丙烯经空气氧化生成丙烯醛后面接水合和加氢“氧化法”这个名字就是这么来的。你会看到真实工业流程里这一步氧化反应恰恰是控制全流程经济性的核心。我自己是先在 Aspen Plus V14 上搭建的模拟框架全部采用稳态模拟物性方法选的 NRTL-RK。这套方案如果你也打算做类似工艺可以直接参考。1. 工艺路线与模拟思路拆解1.1 氧化法制备 1,3-丙二醇的技术背景工业上 1,3-丙二醇的生产主流路线确实就是两条一条是环氧乙烷氢甲酰化法典型代表是 Shell 的工艺另一条是丙烯醛水合加氢法典型代表是 DuPont 和 Degussa 的路线。我们说的“氧化法”实质上就是丙烯醛路线的延伸它从丙烯出发先通过空气氧化得到丙烯醛再经水合生成 3-羟基丙醛HPA最后加氢得到 1,3-丙二醇。我不否认这条路线最大的优势在于原料丙烯来源广泛而丙烯氧化制丙烯醛的技术本身非常成熟工业上大量用于生产丙烯酸和丙烯醛。催化剂体系以钼-铋系混合氧化物为主反应温度控制在 300-350°C压力通常在常压到 1 MPa 之间单程转化率和选择性都相当可观。这条路线也有棘手的地方中间产物丙烯醛和 3-羟基丙醛都不稳定热敏性强容易发生聚合和深度氧化。这就导致下游分离工段和加氢工段的设计需要格外小心而这些正好是流程模拟里最能体现技术含量的环节。我做模拟的目标很明确验证全流程的物料和能量是否平衡考察反应条件和水合工段的稀释比对最终产品纯度的影响同时为后续设备设计和经济评估提供一个比较可靠的工艺基础数据。1.2 全流程模拟框架设计模拟流程整体分四个工段氧化反应工段丙烯和空气按一定配比进入氧化反应器生成丙烯醛、水以及副产物CO、CO2、乙酸、丙烯酸等。急冷吸收工段反应气经过急冷塔用水吸收丙烯醛同时把不凝气N2、CO、CO2分离出去。水合反应工段丙烯醛水溶液进入水合反应器在酸性催化剂作用下生成 3-羟基丙醛。加氢精制工段HPA 水溶液加氢生成 1,3-丙二醇然后通过多塔精制脱水和脱重组分得到产品。在 Aspen Plus 里模块选型上我的做法是工段模块选择选择理由氧化反应器RStoich化学计量反应器有明确的主副反应和工业数据支撑转化率、选择性急冷塔RadFrac无再沸器吸收塔模式水吸收丙烯醛是典型吸收过程用平衡级模型足够水合反应器RPlug活塞流反应器或 RStoich如果只需要物料平衡RStoich 更稳要研究停留时间则需要 RPlug加氢反应器RStoich 或 Ryield加氢反应深度高主反应选择性以实验数据确定精制塔RadFrac严格精馏计算设计规定设置产品纯度这个模块选择的原则就是“能用简单可靠模型的绝不上复杂模块”。氧化反应器的动力学模型不是没有但公开文献里的动力学参数大多基于特定催化剂换一个催化剂体系就不适用了。既然模拟的核心目的是验证工艺集成方案用 RStoich 把反应收率锁定在合理范围内就够用。后面如果要对反应器做详细设计再单独建动力学模型也不迟。1.3 组分体系的搭建要点组分这块是 Aspen 模拟里最基础也最容易出问题的地方。我的组分表包括原料丙烯C3H6、空气按 N2、O2 考虑主产物丙烯醛C3H4O、3-羟基丙醛C3H6O2、1,3-丙二醇C3H8O2副产物CO、CO2、乙酸、丙烯酸、水惰性组分N21,3-丙二醇和丙烯醛在 Aspen 物性数据库里都有内置参数不需要额外定义。需要注意的其实只有 3-羟基丙醛HPA这个物质在有些版本的数据库里没有完整的纯组分参数。我当时查了一下 V14 的数据块发现部分物性参数缺失比如正常沸点、偏心因子和蒸汽压方程系数。这种缺失会在模拟过程中引发莫名其妙的报错尤其是精馏塔算到最后突然警告某个组分饱和蒸汽压超出拟合范围。遇到这种情况常规做法是勾选物性估算Properties Estimation用基团贡献法主要是 Joback 法和 NIST 的 TDE 方法把缺失参数补上。补完之后最好手动核验一下估算的正常沸点是否在合理区间HPA 的正常沸点应该在 160°C 上下常压下如果估算出来偏差太大说明基团拆分有问题需要手动调整。2. 物性方法与反应器细节的硬核讨论2.1 物性方法为什么选 NRTL-RK 而不是其他流程中的关键组分水、丙烯醛、HPA、1,3-丙二醇都带有强极性基团而且水合反应和吸收过程都涉及强非理想体系。所以物性方法的选择直接决定了精馏塔塔板数、回流比和吸收塔用水量这些核心结果。我在前期对比了三种方案NRTL-RK、UNIQUAC-RK、以及 Wilson。Wilson 在液液分相体系上天然受限而这个体系在急冷塔、脱水塔里是可能出现液液两相的所以直接排除。UNIQUAC 和 NRTL 在极性体系上表现接近但 NRTL 对含水体系的相平衡拟合通常更稳定收敛也更顺。最终定了 NRTL-RK气相用 RK 方程修正非理想性适合中低压体系。这里要特别提醒一个易错点Aspen 的 NRTL 默认参数库里水-丙烯醛、水-HPA、丙烯醛-HPA 这些二元交互参数并不都有。如果直接用默认参数吸收塔和精馏塔算出来的结果可能和实际偏差很大甚至出现荒谬的共沸组成。我的处理办法是优先查 Aspen 内置二元参数数据库NRTL 的 Binary Interaction 数据如果内置没有就用 UNIFAC 基团贡献法估算二元参数。算完之后再拿实验数据或文献数据做校验。比如水-丙烯醛体系的汽液平衡数据公开文献里有不少可以对比一下常压下丙烯醛-水的共沸组成估算值和文献值之间的误差控制在 5% 以内才放心。2.2 氧化反应器的建模与参数设定氧化反应器是整个工艺最核心的设备。丙烯和空气在钼-铋催化剂上发生部分氧化主反应C3H6 O2 → C3H4O H2O这个反应强放热ΔH 比较大温升控制是工程关键。我做模拟时主反应的转化率按丙烯单程转化率 90% 设定选择性 85%即转化为丙烯醛的丙烯占反应掉的丙烯的比例。副反应我考虑了三条深度氧化C3H6 4.5O2 → 3CO2 3H2O生成乙酸C3H6 1.5O2 → C2H4O2 H2O简化处理生成丙烯酸C3H6 1.5O2 → C3H4O2 H2O这些副反应的分配比例我根据工业尾气组成大概估算的然后把总氧耗配平。这里有个操作细节RStoich 模块只需要你指定每个反应的计量系数和反应程度或转化率/选择性如果不指定反应程度它会默认按限制反应物完全反应那氧耗会超出实际导致反应器出口组成失真。因此需要把每个副反应的摩尔转化量显式指定出来或者用 Calculator 工具根据选择性比例自动计算。氧化反应器的温度和压力我设定为 330°C、0.2 MPa绝压这是比较典型的丙烯氧化制丙烯醛条件。压力不能太高压力高了选择性会下降深度氧化加剧。温度也不能太高超过 380°C 催化剂会很快失活而且丙烯醛热分解加剧。这些数值虽然来自文献但在模拟里可以作为基准工况后面用灵敏度分析扫描温度区间看对全流程能耗的影响趋势。2.3 急冷吸收塔的“水量权衡”问题从氧化反应器出来的气体温度很高直接进吸收塔肯定不行。我中间加了一级急冷换热器先把反应气冷却到 90°C 左右注意冷却温度不能太低否则丙烯醛和水会在换热器里形成液态聚合物堵塞设备。这是实际装置里常见的操作事故模拟里也要留个心眼把冷却后温度设置在露点以上比较稳妥。吸收塔我用的是 RadFrac 的 Absorber 模式塔板数取 15 块理论板操作压力 0.15 MPa塔顶温度大约 40-50°C。吸收剂是水水量多少需要权衡水量太小丙烯醛吸收不完全尾气排放损失大还污染环境水量太大后面水合反应器的进料浓度太低反应速率和平衡都受影响而且加氢工段和脱水工段的能耗会大幅上升。我做了一个简单估算来定初始水量。假设尾气中允许丙烯醛含量为 500 ppm体积分数用亨利系数估算溶解量。在 40°C 下丙烯醛在水中的亨利常数约为 0.1 MPa·m³/kmol 左右这个值来源不唯一最好用 Aspen 的物性分析直接查由此推出吸收所需的最小水量然后乘以 1.2-1.5 的安全系数作为初值。最后我在模拟里把水吸收剂流量初值设为 40 kmol/h 左右再通过灵敏度分析优化。2.4 水合反应器与加氢反应器的简化处理水合反应C3H4O H2O → C3H6O2HPA这个反应在无催化剂时速率很慢工程上使用酸性阳离子交换树脂或磷酸作催化剂温度一般控制在 50-70°C压力 0.3-0.8 MPa。反应是放热的温度不宜过高否则 HPA 会脱水变回丙烯醛或者生成低聚物。HPA 在水溶液里也不稳定浓度越高越容易聚合所以水合反应器进料中丙烯醛浓度通常控制在 10% 左右。我在模拟里把水合反应器用 RStoich 建模按丙烯醛转化率 80% 设定也就是说出口还有 20% 的丙烯醛未反应。这部分未反应的丙烯醛在后续加氢工段也会被加氢成丙醇之类的东西所以分离上需要考虑。不过我在主流程里做了一步简化水合反应器出口物料直接进加氢反应器未反应的丙烯醛也一块加氢。这时未反应的丙烯醛加氢生成正丙醇作为杂质在精制塔里脱除。这个做法在模拟里是合理简化真实装置可能在加氢前还要加一个脱醛塔但那是设备设计层面的事了。加氢反应器我设定为 80°C、5 MPa采用 RStoichHPA 加氢生成 1,3-丙二醇转化率按 99.5% 计算。这个转化率来自文献中 Ni 基催化剂的数据实际操作中高氢油比、高压力下达到这个水平并不稀奇。副反应加了一条HPA 加氢过度生成正丙醇选择性 1%用来平衡组成、让精馏塔塔釜有真实杂质。3. 精馏序列设计与 Aspen 实操设定3.1 分离序列的总体思路加氢反应器出口的物料组成大概是1,3-丙二醇、水、正丙醇、少量未转化的 HPA、微量丙烯醛、溶解气体。要得到聚合级 1,3-丙二醇纯度 99.5% 以上我设计的分离三步走闪蒸罐降压闪蒸脱除溶解的 H2 和轻组分如果有剩余 H2操作压力降到 0.5 MPa 左右。脱水塔常压或微负压精馏把大部分水从塔顶蒸出。1,3-丙二醇沸点 214°C常压水沸点 100°C两者沸点差足够大分离难度不大。但要注意塔釜温度不能太高1,3-丙二醇长时间处于 200°C 以上会变色和分解所以精制塔要用减压操作。产品精制塔减压精馏塔顶压力 10-20 kPa塔顶采出纯 1,3-丙二醇塔釜排出高沸物HPA 低聚物、重组分。这套序列本质上是“先脱轻、再脱水、最后脱重”的经典顺序在多元醇生产中很常见。好处是每个塔的分离任务单一操作窗口宽对 Aspen 模拟来说也更好收敛。3.2 脱水塔的设计规定与回流比设定脱水塔我设定 20 块理论板常压操作进料在第 10 块板较中间位置回流比初值给 1.5。然后用 Design Spec 设定塔釜水中 1,3-丙二醇含量不超过 0.1%质量分数塔顶水相中 1,3-丙二醇含量不超过 0.5%。这里有一个实际操作上的小问题水-1,3-丙二醇的分离本身容易但由于进料里含有少量正丙醇正丙醇会和水形成共沸物正丙醇-水共沸组成大约是 43% 正丙醇/57% 水共沸温度 87.7°C导致塔顶水相里混有正丙醇很难单靠精馏彻底分开。工厂里这个位置通常会加一个倾析器或者用共沸精馏/渗透汽化来处理。我在模拟中不追求把这部分正丙醇完全除掉而是允许塔顶水相含少量正丙醇后续作为废水处理。如果你要做的工艺要求废水 COD 很低就得在这个位置单独加一个轻组分脱除塔。Aspen 里做设计规定时Vary 选择回流比或者塔顶采出量。一般来说如果目标是产品纯度Vary 回流比如果目标是控制采出量Vary 馏出速率。脱水塔我选的是 Vary 回流比范围设 0.5-3这样收敛比较稳。3.3 产品精制塔的减压操作要点产品精制塔是整个流程里最考验细节的地方。常压下 1,3-丙二醇沸点约 214°C塔釜温度会超过 220°C这个温度下 1,3-丙二醇会缓慢聚合和分解产品色号变差。所以必须减压。我设置塔顶压力 15 kPa绝压用 Vacuum 模式运行此时塔顶温度大约 140-150°C塔釜温度约 180°C这样既保证分离效果又避免热降解。塔板数取 25 块理论板回流比 2因为减压下相对挥发度增大不需要太高的回流比。设计规定塔顶产品纯度 1,3-丙二醇 ≥ 99.5%质量分数塔釜 1,3-丙二醇损失 ≤ 1%。减压塔的模拟要特别注意一个 Aspen 的提示RadFrac 在低压下算水力学Hydraulics时容易警告“Weir loading too low”或者“Flooding”这是因为减压塔气相体积流量大塔径如果按常压塔估会偏小。我做设备尺寸估算时单独校核了塔径模拟阶段不强制收敛水力学计算直接把塔径尺寸初值给大一点等全流程收敛后再做严格水力学核算。4. 循环流股与全流程收敛策略4.1 全流程收敛的核心难点这套流程有一个让 Aspen 新手很头疼的地方急冷吸收塔用的吸收水实际上有一部分来自脱水塔塔顶回收的水相也就是说流程内部存在一个水循环回路。我最初尝试把吸收水设为一股循环流股结果全流程迭代半天不收敛问题就出在这股循环水上。循环流股不收敛的典型症状是每一次迭代循环流股的组成都在漂移尤其是水相中溶解的微量有机物正丙醇、丙烯醛会在回路里累积Aspen 的 Wegstein 法有时候压不住这种振荡。我用的解决手段有三个先断开循环流股用一个 Water 或 Makeup 流股代替等全流程其他部分收敛后再接上循环接上循环后使用 Calculator Design Spec 控制补充水量让整个水循环达到物料平衡迭代收敛参数从默认的 Wegstein 改成 Broyden然后把收敛容差从 0.0001 放宽到 0.001先得到解再收紧容差复算。这个操作思路适用于所有带循环的化工流程模拟是 Aspen 通用的标准做法。4.2 撕裂流股的选择与初值设定如果你非要一开始就带循环算那就要选好撕裂流股Teardown Stream。我的经验是把吸收塔补水的那股水相流股作为撕裂流股给它一个初值流量 50 kmol/h、纯水、温度 40°C。初值越接近真实解收敛越快。这里有一个“看结果判断初值合理性”的技巧第一次收敛后看撕裂流股的实际组成和流量如果与初值差别很大比如实际流量只有 20 kmol/h 而初值给了 50 kmol/h说明水循环里存在净积累或净消耗要么是吸收塔出料带走的水太多要么是水合反应消耗的水没有补偿。这时候要检查水合反应器的反应程度和脱水塔的分离效率而不是盲目加大初值。物料平衡的直觉在这里比盲目调参重要得多。4.3 全流程模拟收敛顺序建议我做整个流程的收敛顺序是严格分步的先单独收敛氧化反应器设定好进料和反应程度跑通查看出口组成。再接急冷吸收塔吸收水先不用循环水直接用一股新鲜水进料跑通吸收塔。再接水合和加氢反应器此时下游没有分离回流反应器出口直接闪蒸出产品和废水流程能跑通。最后再接脱水塔和产品精制塔把闪蒸后的液体进塔分离塔顶水相引回吸收塔补水形成闭环。每一步都在上一步基础上增加单元出问题时定位范围小。你要是直接一上来就把全流程图全连上然后点 Run大概率会得到一堆红色错误和警告然后开始漫长而无序的排查。5. 常见问题与排查技巧实录5.1 物性缺失与估算错误这是最容易踩的坑。我在模拟初期就遇到 3-羟基丙醛物性参数不全结果水合反应器出口流股进入闪蒸罐时Aspen 直接报错“Component HPA is missing properties for equilibrium calculation。”排查过程很简单打开 Properties 下的 Review 页面看哪个组分缺参数然后用 Property Estimation 估算。但估算之后务必重新检查一下模拟结果HPA 估算的饱和蒸汽压如果偏差大会直接导致闪蒸罐和精馏塔的相分率错误进而让塔的收敛反复失败。5.2 反应器热量不平衡RStoich 反应器的反应热配平如果出了问题Aspen 有时候会在 Energy Balance 汇总里给你一个巨大的冷却负荷但实际上完全不合理。我的做法是在理清反应计量方程时认真配平每个反应的原子数然后查看 Heat of Reaction 是否合理。丙烯氧化是强放热反应反应器冷却负荷应该非常显著如果你算出来反应器还需要外界加热那百分百是反应方程式写错了比如漏写产物水。这种情况下要回到 Chemistry 列表逐条检查计量方程和反应程度。5.3 吸收塔水量过大导致收敛振荡吸收塔的水量如果给得过大会造成两个后果一是尾气带水增多塔顶冷凝器负荷飙升二是水循环回路中的水流量在每次迭代时变动剧烈收敛异常慢。我建议用灵敏度分析扫一遍吸收水量对丙烯醛吸收率的曲线找到一个“吸收率已经接近 100% 但水量还没过度增加”的拐点以这个拐点的水量作为操作值。这个拐点通常在丙烯醛吸收率 99%摩尔分数附近水量再增加对吸收率几乎没有贡献只会白白增加能耗。5.4 精馏塔不收敛的通用处理方法RadFrac 不收敛的原因五花八门塔板数不够、进料位置不合理、回流比初值远离实际、设计规定太苛刻等等。我的通用排查顺序是先不设设计规定给定回流比和采出量跑通基础工况看塔内温度分布和组成分布是否合理有没有明显的温度拐点重组分是否在塔釜富集再逐步加设计规定Vary 范围要宽一些否则边界条件卡死导致振荡如果还是不收敛把塔板数多加 30%重新算一遍。这条方法适用于绝大多数常减压塔和水吸收塔。5.5 常见问题速查表症状可能原因解决方向物性缺失报错HPA 等组分无内置参数启用 Property Estimation手动校验估算值合理性反应器能量严重不平衡计量方程未配平或反应程度设置错误检查 Chemistry 列表核对原子数和反应热吸收塔收敛振荡吸收水量过大或过小用灵敏度分析扫水量-吸收率曲线找拐点脱水塔塔顶水相含醇正丙醇-水共沸接受副产物或在模拟中增加共沸分离单元减压塔水力学警告气相体积流量大塔径不够模拟阶段放宽水力学核算设备设计时单独校核循环回路不收敛水循环存在净积累/消耗断开循环先收敛再接回换 Broyden 迭代法产品纯度达不到 99.5%回流比不足或理论板数不足增加回流比和塔板数用 Design Spec 锁定纯度6. 灵敏度分析与结果解读6.1 氧化反应温度对下游的影响我针对氧化反应器的温度做了一次灵敏度分析范围从 300°C 到 360°C。假设转化率随温度升高而升高但选择性随温度升高而下降这是典型的平行反应特征。结果很明显存在一个最优温度区间大约 320-340°C总收率最高。温度超过 350°C 后虽然转化率高了但 CO2 生成量剧增1,3-丙二醇的最终产量反而下降。这个结论和工业装置的实际操作经验一致。6.2 水合反应器中丙烯醛浓度的影响水合反应器进料中丙烯醛浓度低转化率更高因为水过量有利于平衡移动但 HPA 浓度低导致后续加氢反应器进料体积流量增大设备投资和能耗都上升。我扫了丙烯醛浓度 5%、8%、10%、15% 四个水平结果是在 10% 浓度下全流程能耗最低。这为吸收塔的水量设定提供了一个反向约束不是吸收越好越好还要顾及水合反应器的进料浓度要求。6.3 产品塔回流比与能耗的权衡产品精制塔回流比从 1.5 到 3 扫描产品纯度从 98.2% 提升到 99.9%但再沸器热负荷基本线性增长。从经济角度看回到 99.5% 纯度所需的回流比大约 2.0-2.5 之间这个区间再沸器热负荷增加幅度趋缓是性价比比较高的操作区域。7. 从模拟到工程落地前的最后提醒按照我个人的习惯流程模拟跑通之后并不会马上结束以下几个点是你交成果之前最好再核对一遍的第一模拟结果的质量衡算必须严格闭合。尤其是各股物流的质量流量合计误差不能超过 0.1%。如果你发现某处质量不守恒多半是某个模块设定有问题别指望后面会自动修正。第二能量衡算要逐设备核查。Aspen 默认会给出每个换热设备的热负荷你用手算复核一下反应器冷却负荷和精馏塔再沸器热负荷的量级是否合理。如果出现某个再沸器热负荷比整个装置的总能耗还大大概率是你把塔的压力等级选错了。第三物性方法选完以后不要频繁切换。我在做灵敏度分析时为了省时间把 NRTL-RK 临时改成过理想物性Ideal结果所有塔的收敛行为都变了白忙活了半天。流程模拟里物性方法的一致性比什么都重要。第四关键设备吸收塔、加氢反应器的操作参数最好和装置上下游的约束对齐。比如加氢反应器压力高了下游闪蒸罐和精制塔的压力等级要相应调整否则高压液体直接闪蒸到常压会引起严重的汽蚀问题。我再分享一个调试时的真实感受这套流程最难的部分不是氧化反应器也不是精馏塔而是那根看不见摸不着的水循环回路。它牵一发而动全身吸收塔、水合反应器、脱水塔的设备参数都互相耦合。你在做模拟时一定要养成“每次只动一个参数观察全流程响应”的习惯不要同时调三个变量否则出了问题根本不知道是谁干的。把这些细节都处理完你把模拟结果汇总成工艺物料平衡表PFD 上的物料表的时候会发现Aspen 给出的数据完全可以用于后续的设备选型、能耗计算和经济评估。流程模拟本来就是越做越细的过程第一次跑通只是万里长征第一步后面根据设计工况再迭代优化就是顺理成章的事了。 Aspen Plus模拟氧化法制备1,3-丙二醇全流程关键技巧
Aspen Plus 做氧化法制备 1,3-丙二醇这个题目我从接到任务到跑通全流程前前后后折腾了将近三周。中间踩过的坑、返工重来的次数说实话比过去做过的很多模拟项目都多。但做完之后收获也大这套流程涉及气固催化氧化、水合反应、高压加氢、多塔精制几乎把 Aspen Plus 里常用的反应器模块、塔模块、物性方法和收敛技巧都用上了。这篇就把整个模拟过程的思路梳理一遍。不只是贴图贴步骤更想把每一步选择的“为什么”讲清楚。1,3-丙二醇1,3-PDO是聚对苯二甲酸丙二醇酯PTT的关键单体这个市场懂行的都懂而目前实现工业化的路线主要有环氧乙烷法和丙烯醛法。丙烯醛法的第一步就是丙烯经空气氧化生成丙烯醛后面接水合和加氢“氧化法”这个名字就是这么来的。你会看到真实工业流程里这一步氧化反应恰恰是控制全流程经济性的核心。我自己是先在 Aspen Plus V14 上搭建的模拟框架全部采用稳态模拟物性方法选的 NRTL-RK。这套方案如果你也打算做类似工艺可以直接参考。1. 工艺路线与模拟思路拆解1.1 氧化法制备 1,3-丙二醇的技术背景工业上 1,3-丙二醇的生产主流路线确实就是两条一条是环氧乙烷氢甲酰化法典型代表是 Shell 的工艺另一条是丙烯醛水合加氢法典型代表是 DuPont 和 Degussa 的路线。我们说的“氧化法”实质上就是丙烯醛路线的延伸它从丙烯出发先通过空气氧化得到丙烯醛再经水合生成 3-羟基丙醛HPA最后加氢得到 1,3-丙二醇。我不否认这条路线最大的优势在于原料丙烯来源广泛而丙烯氧化制丙烯醛的技术本身非常成熟工业上大量用于生产丙烯酸和丙烯醛。催化剂体系以钼-铋系混合氧化物为主反应温度控制在 300-350°C压力通常在常压到 1 MPa 之间单程转化率和选择性都相当可观。这条路线也有棘手的地方中间产物丙烯醛和 3-羟基丙醛都不稳定热敏性强容易发生聚合和深度氧化。这就导致下游分离工段和加氢工段的设计需要格外小心而这些正好是流程模拟里最能体现技术含量的环节。我做模拟的目标很明确验证全流程的物料和能量是否平衡考察反应条件和水合工段的稀释比对最终产品纯度的影响同时为后续设备设计和经济评估提供一个比较可靠的工艺基础数据。1.2 全流程模拟框架设计模拟流程整体分四个工段氧化反应工段丙烯和空气按一定配比进入氧化反应器生成丙烯醛、水以及副产物CO、CO2、乙酸、丙烯酸等。急冷吸收工段反应气经过急冷塔用水吸收丙烯醛同时把不凝气N2、CO、CO2分离出去。水合反应工段丙烯醛水溶液进入水合反应器在酸性催化剂作用下生成 3-羟基丙醛。加氢精制工段HPA 水溶液加氢生成 1,3-丙二醇然后通过多塔精制脱水和脱重组分得到产品。在 Aspen Plus 里模块选型上我的做法是工段模块选择选择理由氧化反应器RStoich化学计量反应器有明确的主副反应和工业数据支撑转化率、选择性急冷塔RadFrac无再沸器吸收塔模式水吸收丙烯醛是典型吸收过程用平衡级模型足够水合反应器RPlug活塞流反应器或 RStoich如果只需要物料平衡RStoich 更稳要研究停留时间则需要 RPlug加氢反应器RStoich 或 Ryield加氢反应深度高主反应选择性以实验数据确定精制塔RadFrac严格精馏计算设计规定设置产品纯度这个模块选择的原则就是“能用简单可靠模型的绝不上复杂模块”。氧化反应器的动力学模型不是没有但公开文献里的动力学参数大多基于特定催化剂换一个催化剂体系就不适用了。既然模拟的核心目的是验证工艺集成方案用 RStoich 把反应收率锁定在合理范围内就够用。后面如果要对反应器做详细设计再单独建动力学模型也不迟。1.3 组分体系的搭建要点组分这块是 Aspen 模拟里最基础也最容易出问题的地方。我的组分表包括原料丙烯C3H6、空气按 N2、O2 考虑主产物丙烯醛C3H4O、3-羟基丙醛C3H6O2、1,3-丙二醇C3H8O2副产物CO、CO2、乙酸、丙烯酸、水惰性组分N21,3-丙二醇和丙烯醛在 Aspen 物性数据库里都有内置参数不需要额外定义。需要注意的其实只有 3-羟基丙醛HPA这个物质在有些版本的数据库里没有完整的纯组分参数。我当时查了一下 V14 的数据块发现部分物性参数缺失比如正常沸点、偏心因子和蒸汽压方程系数。这种缺失会在模拟过程中引发莫名其妙的报错尤其是精馏塔算到最后突然警告某个组分饱和蒸汽压超出拟合范围。遇到这种情况常规做法是勾选物性估算Properties Estimation用基团贡献法主要是 Joback 法和 NIST 的 TDE 方法把缺失参数补上。补完之后最好手动核验一下估算的正常沸点是否在合理区间HPA 的正常沸点应该在 160°C 上下常压下如果估算出来偏差太大说明基团拆分有问题需要手动调整。2. 物性方法与反应器细节的硬核讨论2.1 物性方法为什么选 NRTL-RK 而不是其他流程中的关键组分水、丙烯醛、HPA、1,3-丙二醇都带有强极性基团而且水合反应和吸收过程都涉及强非理想体系。所以物性方法的选择直接决定了精馏塔塔板数、回流比和吸收塔用水量这些核心结果。我在前期对比了三种方案NRTL-RK、UNIQUAC-RK、以及 Wilson。Wilson 在液液分相体系上天然受限而这个体系在急冷塔、脱水塔里是可能出现液液两相的所以直接排除。UNIQUAC 和 NRTL 在极性体系上表现接近但 NRTL 对含水体系的相平衡拟合通常更稳定收敛也更顺。最终定了 NRTL-RK气相用 RK 方程修正非理想性适合中低压体系。这里要特别提醒一个易错点Aspen 的 NRTL 默认参数库里水-丙烯醛、水-HPA、丙烯醛-HPA 这些二元交互参数并不都有。如果直接用默认参数吸收塔和精馏塔算出来的结果可能和实际偏差很大甚至出现荒谬的共沸组成。我的处理办法是优先查 Aspen 内置二元参数数据库NRTL 的 Binary Interaction 数据如果内置没有就用 UNIFAC 基团贡献法估算二元参数。算完之后再拿实验数据或文献数据做校验。比如水-丙烯醛体系的汽液平衡数据公开文献里有不少可以对比一下常压下丙烯醛-水的共沸组成估算值和文献值之间的误差控制在 5% 以内才放心。2.2 氧化反应器的建模与参数设定氧化反应器是整个工艺最核心的设备。丙烯和空气在钼-铋催化剂上发生部分氧化主反应C3H6 O2 → C3H4O H2O这个反应强放热ΔH 比较大温升控制是工程关键。我做模拟时主反应的转化率按丙烯单程转化率 90% 设定选择性 85%即转化为丙烯醛的丙烯占反应掉的丙烯的比例。副反应我考虑了三条深度氧化C3H6 4.5O2 → 3CO2 3H2O生成乙酸C3H6 1.5O2 → C2H4O2 H2O简化处理生成丙烯酸C3H6 1.5O2 → C3H4O2 H2O这些副反应的分配比例我根据工业尾气组成大概估算的然后把总氧耗配平。这里有个操作细节RStoich 模块只需要你指定每个反应的计量系数和反应程度或转化率/选择性如果不指定反应程度它会默认按限制反应物完全反应那氧耗会超出实际导致反应器出口组成失真。因此需要把每个副反应的摩尔转化量显式指定出来或者用 Calculator 工具根据选择性比例自动计算。氧化反应器的温度和压力我设定为 330°C、0.2 MPa绝压这是比较典型的丙烯氧化制丙烯醛条件。压力不能太高压力高了选择性会下降深度氧化加剧。温度也不能太高超过 380°C 催化剂会很快失活而且丙烯醛热分解加剧。这些数值虽然来自文献但在模拟里可以作为基准工况后面用灵敏度分析扫描温度区间看对全流程能耗的影响趋势。2.3 急冷吸收塔的“水量权衡”问题从氧化反应器出来的气体温度很高直接进吸收塔肯定不行。我中间加了一级急冷换热器先把反应气冷却到 90°C 左右注意冷却温度不能太低否则丙烯醛和水会在换热器里形成液态聚合物堵塞设备。这是实际装置里常见的操作事故模拟里也要留个心眼把冷却后温度设置在露点以上比较稳妥。吸收塔我用的是 RadFrac 的 Absorber 模式塔板数取 15 块理论板操作压力 0.15 MPa塔顶温度大约 40-50°C。吸收剂是水水量多少需要权衡水量太小丙烯醛吸收不完全尾气排放损失大还污染环境水量太大后面水合反应器的进料浓度太低反应速率和平衡都受影响而且加氢工段和脱水工段的能耗会大幅上升。我做了一个简单估算来定初始水量。假设尾气中允许丙烯醛含量为 500 ppm体积分数用亨利系数估算溶解量。在 40°C 下丙烯醛在水中的亨利常数约为 0.1 MPa·m³/kmol 左右这个值来源不唯一最好用 Aspen 的物性分析直接查由此推出吸收所需的最小水量然后乘以 1.2-1.5 的安全系数作为初值。最后我在模拟里把水吸收剂流量初值设为 40 kmol/h 左右再通过灵敏度分析优化。2.4 水合反应器与加氢反应器的简化处理水合反应C3H4O H2O → C3H6O2HPA这个反应在无催化剂时速率很慢工程上使用酸性阳离子交换树脂或磷酸作催化剂温度一般控制在 50-70°C压力 0.3-0.8 MPa。反应是放热的温度不宜过高否则 HPA 会脱水变回丙烯醛或者生成低聚物。HPA 在水溶液里也不稳定浓度越高越容易聚合所以水合反应器进料中丙烯醛浓度通常控制在 10% 左右。我在模拟里把水合反应器用 RStoich 建模按丙烯醛转化率 80% 设定也就是说出口还有 20% 的丙烯醛未反应。这部分未反应的丙烯醛在后续加氢工段也会被加氢成丙醇之类的东西所以分离上需要考虑。不过我在主流程里做了一步简化水合反应器出口物料直接进加氢反应器未反应的丙烯醛也一块加氢。这时未反应的丙烯醛加氢生成正丙醇作为杂质在精制塔里脱除。这个做法在模拟里是合理简化真实装置可能在加氢前还要加一个脱醛塔但那是设备设计层面的事了。加氢反应器我设定为 80°C、5 MPa采用 RStoichHPA 加氢生成 1,3-丙二醇转化率按 99.5% 计算。这个转化率来自文献中 Ni 基催化剂的数据实际操作中高氢油比、高压力下达到这个水平并不稀奇。副反应加了一条HPA 加氢过度生成正丙醇选择性 1%用来平衡组成、让精馏塔塔釜有真实杂质。3. 精馏序列设计与 Aspen 实操设定3.1 分离序列的总体思路加氢反应器出口的物料组成大概是1,3-丙二醇、水、正丙醇、少量未转化的 HPA、微量丙烯醛、溶解气体。要得到聚合级 1,3-丙二醇纯度 99.5% 以上我设计的分离三步走闪蒸罐降压闪蒸脱除溶解的 H2 和轻组分如果有剩余 H2操作压力降到 0.5 MPa 左右。脱水塔常压或微负压精馏把大部分水从塔顶蒸出。1,3-丙二醇沸点 214°C常压水沸点 100°C两者沸点差足够大分离难度不大。但要注意塔釜温度不能太高1,3-丙二醇长时间处于 200°C 以上会变色和分解所以精制塔要用减压操作。产品精制塔减压精馏塔顶压力 10-20 kPa塔顶采出纯 1,3-丙二醇塔釜排出高沸物HPA 低聚物、重组分。这套序列本质上是“先脱轻、再脱水、最后脱重”的经典顺序在多元醇生产中很常见。好处是每个塔的分离任务单一操作窗口宽对 Aspen 模拟来说也更好收敛。3.2 脱水塔的设计规定与回流比设定脱水塔我设定 20 块理论板常压操作进料在第 10 块板较中间位置回流比初值给 1.5。然后用 Design Spec 设定塔釜水中 1,3-丙二醇含量不超过 0.1%质量分数塔顶水相中 1,3-丙二醇含量不超过 0.5%。这里有一个实际操作上的小问题水-1,3-丙二醇的分离本身容易但由于进料里含有少量正丙醇正丙醇会和水形成共沸物正丙醇-水共沸组成大约是 43% 正丙醇/57% 水共沸温度 87.7°C导致塔顶水相里混有正丙醇很难单靠精馏彻底分开。工厂里这个位置通常会加一个倾析器或者用共沸精馏/渗透汽化来处理。我在模拟中不追求把这部分正丙醇完全除掉而是允许塔顶水相含少量正丙醇后续作为废水处理。如果你要做的工艺要求废水 COD 很低就得在这个位置单独加一个轻组分脱除塔。Aspen 里做设计规定时Vary 选择回流比或者塔顶采出量。一般来说如果目标是产品纯度Vary 回流比如果目标是控制采出量Vary 馏出速率。脱水塔我选的是 Vary 回流比范围设 0.5-3这样收敛比较稳。3.3 产品精制塔的减压操作要点产品精制塔是整个流程里最考验细节的地方。常压下 1,3-丙二醇沸点约 214°C塔釜温度会超过 220°C这个温度下 1,3-丙二醇会缓慢聚合和分解产品色号变差。所以必须减压。我设置塔顶压力 15 kPa绝压用 Vacuum 模式运行此时塔顶温度大约 140-150°C塔釜温度约 180°C这样既保证分离效果又避免热降解。塔板数取 25 块理论板回流比 2因为减压下相对挥发度增大不需要太高的回流比。设计规定塔顶产品纯度 1,3-丙二醇 ≥ 99.5%质量分数塔釜 1,3-丙二醇损失 ≤ 1%。减压塔的模拟要特别注意一个 Aspen 的提示RadFrac 在低压下算水力学Hydraulics时容易警告“Weir loading too low”或者“Flooding”这是因为减压塔气相体积流量大塔径如果按常压塔估会偏小。我做设备尺寸估算时单独校核了塔径模拟阶段不强制收敛水力学计算直接把塔径尺寸初值给大一点等全流程收敛后再做严格水力学核算。4. 循环流股与全流程收敛策略4.1 全流程收敛的核心难点这套流程有一个让 Aspen 新手很头疼的地方急冷吸收塔用的吸收水实际上有一部分来自脱水塔塔顶回收的水相也就是说流程内部存在一个水循环回路。我最初尝试把吸收水设为一股循环流股结果全流程迭代半天不收敛问题就出在这股循环水上。循环流股不收敛的典型症状是每一次迭代循环流股的组成都在漂移尤其是水相中溶解的微量有机物正丙醇、丙烯醛会在回路里累积Aspen 的 Wegstein 法有时候压不住这种振荡。我用的解决手段有三个先断开循环流股用一个 Water 或 Makeup 流股代替等全流程其他部分收敛后再接上循环接上循环后使用 Calculator Design Spec 控制补充水量让整个水循环达到物料平衡迭代收敛参数从默认的 Wegstein 改成 Broyden然后把收敛容差从 0.0001 放宽到 0.001先得到解再收紧容差复算。这个操作思路适用于所有带循环的化工流程模拟是 Aspen 通用的标准做法。4.2 撕裂流股的选择与初值设定如果你非要一开始就带循环算那就要选好撕裂流股Teardown Stream。我的经验是把吸收塔补水的那股水相流股作为撕裂流股给它一个初值流量 50 kmol/h、纯水、温度 40°C。初值越接近真实解收敛越快。这里有一个“看结果判断初值合理性”的技巧第一次收敛后看撕裂流股的实际组成和流量如果与初值差别很大比如实际流量只有 20 kmol/h 而初值给了 50 kmol/h说明水循环里存在净积累或净消耗要么是吸收塔出料带走的水太多要么是水合反应消耗的水没有补偿。这时候要检查水合反应器的反应程度和脱水塔的分离效率而不是盲目加大初值。物料平衡的直觉在这里比盲目调参重要得多。4.3 全流程模拟收敛顺序建议我做整个流程的收敛顺序是严格分步的先单独收敛氧化反应器设定好进料和反应程度跑通查看出口组成。再接急冷吸收塔吸收水先不用循环水直接用一股新鲜水进料跑通吸收塔。再接水合和加氢反应器此时下游没有分离回流反应器出口直接闪蒸出产品和废水流程能跑通。最后再接脱水塔和产品精制塔把闪蒸后的液体进塔分离塔顶水相引回吸收塔补水形成闭环。每一步都在上一步基础上增加单元出问题时定位范围小。你要是直接一上来就把全流程图全连上然后点 Run大概率会得到一堆红色错误和警告然后开始漫长而无序的排查。5. 常见问题与排查技巧实录5.1 物性缺失与估算错误这是最容易踩的坑。我在模拟初期就遇到 3-羟基丙醛物性参数不全结果水合反应器出口流股进入闪蒸罐时Aspen 直接报错“Component HPA is missing properties for equilibrium calculation。”排查过程很简单打开 Properties 下的 Review 页面看哪个组分缺参数然后用 Property Estimation 估算。但估算之后务必重新检查一下模拟结果HPA 估算的饱和蒸汽压如果偏差大会直接导致闪蒸罐和精馏塔的相分率错误进而让塔的收敛反复失败。5.2 反应器热量不平衡RStoich 反应器的反应热配平如果出了问题Aspen 有时候会在 Energy Balance 汇总里给你一个巨大的冷却负荷但实际上完全不合理。我的做法是在理清反应计量方程时认真配平每个反应的原子数然后查看 Heat of Reaction 是否合理。丙烯氧化是强放热反应反应器冷却负荷应该非常显著如果你算出来反应器还需要外界加热那百分百是反应方程式写错了比如漏写产物水。这种情况下要回到 Chemistry 列表逐条检查计量方程和反应程度。5.3 吸收塔水量过大导致收敛振荡吸收塔的水量如果给得过大会造成两个后果一是尾气带水增多塔顶冷凝器负荷飙升二是水循环回路中的水流量在每次迭代时变动剧烈收敛异常慢。我建议用灵敏度分析扫一遍吸收水量对丙烯醛吸收率的曲线找到一个“吸收率已经接近 100% 但水量还没过度增加”的拐点以这个拐点的水量作为操作值。这个拐点通常在丙烯醛吸收率 99%摩尔分数附近水量再增加对吸收率几乎没有贡献只会白白增加能耗。5.4 精馏塔不收敛的通用处理方法RadFrac 不收敛的原因五花八门塔板数不够、进料位置不合理、回流比初值远离实际、设计规定太苛刻等等。我的通用排查顺序是先不设设计规定给定回流比和采出量跑通基础工况看塔内温度分布和组成分布是否合理有没有明显的温度拐点重组分是否在塔釜富集再逐步加设计规定Vary 范围要宽一些否则边界条件卡死导致振荡如果还是不收敛把塔板数多加 30%重新算一遍。这条方法适用于绝大多数常减压塔和水吸收塔。5.5 常见问题速查表症状可能原因解决方向物性缺失报错HPA 等组分无内置参数启用 Property Estimation手动校验估算值合理性反应器能量严重不平衡计量方程未配平或反应程度设置错误检查 Chemistry 列表核对原子数和反应热吸收塔收敛振荡吸收水量过大或过小用灵敏度分析扫水量-吸收率曲线找拐点脱水塔塔顶水相含醇正丙醇-水共沸接受副产物或在模拟中增加共沸分离单元减压塔水力学警告气相体积流量大塔径不够模拟阶段放宽水力学核算设备设计时单独校核循环回路不收敛水循环存在净积累/消耗断开循环先收敛再接回换 Broyden 迭代法产品纯度达不到 99.5%回流比不足或理论板数不足增加回流比和塔板数用 Design Spec 锁定纯度6. 灵敏度分析与结果解读6.1 氧化反应温度对下游的影响我针对氧化反应器的温度做了一次灵敏度分析范围从 300°C 到 360°C。假设转化率随温度升高而升高但选择性随温度升高而下降这是典型的平行反应特征。结果很明显存在一个最优温度区间大约 320-340°C总收率最高。温度超过 350°C 后虽然转化率高了但 CO2 生成量剧增1,3-丙二醇的最终产量反而下降。这个结论和工业装置的实际操作经验一致。6.2 水合反应器中丙烯醛浓度的影响水合反应器进料中丙烯醛浓度低转化率更高因为水过量有利于平衡移动但 HPA 浓度低导致后续加氢反应器进料体积流量增大设备投资和能耗都上升。我扫了丙烯醛浓度 5%、8%、10%、15% 四个水平结果是在 10% 浓度下全流程能耗最低。这为吸收塔的水量设定提供了一个反向约束不是吸收越好越好还要顾及水合反应器的进料浓度要求。6.3 产品塔回流比与能耗的权衡产品精制塔回流比从 1.5 到 3 扫描产品纯度从 98.2% 提升到 99.9%但再沸器热负荷基本线性增长。从经济角度看回到 99.5% 纯度所需的回流比大约 2.0-2.5 之间这个区间再沸器热负荷增加幅度趋缓是性价比比较高的操作区域。7. 从模拟到工程落地前的最后提醒按照我个人的习惯流程模拟跑通之后并不会马上结束以下几个点是你交成果之前最好再核对一遍的第一模拟结果的质量衡算必须严格闭合。尤其是各股物流的质量流量合计误差不能超过 0.1%。如果你发现某处质量不守恒多半是某个模块设定有问题别指望后面会自动修正。第二能量衡算要逐设备核查。Aspen 默认会给出每个换热设备的热负荷你用手算复核一下反应器冷却负荷和精馏塔再沸器热负荷的量级是否合理。如果出现某个再沸器热负荷比整个装置的总能耗还大大概率是你把塔的压力等级选错了。第三物性方法选完以后不要频繁切换。我在做灵敏度分析时为了省时间把 NRTL-RK 临时改成过理想物性Ideal结果所有塔的收敛行为都变了白忙活了半天。流程模拟里物性方法的一致性比什么都重要。第四关键设备吸收塔、加氢反应器的操作参数最好和装置上下游的约束对齐。比如加氢反应器压力高了下游闪蒸罐和精制塔的压力等级要相应调整否则高压液体直接闪蒸到常压会引起严重的汽蚀问题。我再分享一个调试时的真实感受这套流程最难的部分不是氧化反应器也不是精馏塔而是那根看不见摸不着的水循环回路。它牵一发而动全身吸收塔、水合反应器、脱水塔的设备参数都互相耦合。你在做模拟时一定要养成“每次只动一个参数观察全流程响应”的习惯不要同时调三个变量否则出了问题根本不知道是谁干的。把这些细节都处理完你把模拟结果汇总成工艺物料平衡表PFD 上的物料表的时候会发现Aspen 给出的数据完全可以用于后续的设备选型、能耗计算和经济评估。流程模拟本来就是越做越细的过程第一次跑通只是万里长征第一步后面根据设计工况再迭代优化就是顺理成章的事了。