下载文档原格式
化工系统工程__化工过程系统稳态模拟与分析2 化工过程系统稳态模拟与分析概述通过对化工工艺流程系统进行稳态模拟与分析也就是对过程系统建立模型并对模型进行求解可以解决下述三方面的问题①过程系统的分析与模拟②过程系统设计③过程系统参数优化①过程系统的分析模拟对某个给定的过程系统模型进行模拟求解可得出该系统的全部状态变量从而可以对该过程系统进行工况分析如图21所示②过程系统设计当对某个或某些系统变量提出设计规定要求时通过调整某些决策变量使模拟结果满足设计规定要求如图22所示③过程系统参数优化过程系统模型与最优化模型联解得到一组使工况目标函数最佳的决策变量优化变量从而实施最佳工况如图所示 2 化工过程系统稳态模拟与分析相关的基本概念 1 系统为了某种目标由共同的物料流或信息流联系在一起的单元组合而形成的整体称为系统 2 子系统组成系统的系统下一层次的事物简单系统子系统就是某个单元复杂系统它的子系统又可能包含有子系统基本概念 3 系统的特性由两方面构成 1系统内各个单元的特性复杂系统则是各子系统的特性 2系统流程的结构特性树结构和再循环结构的概念 4 过程拓扑将过程流程图转换为信息流程图再把信息流程图转变为过程矩阵的过程称为过程拓扑过程流程→信息流程用有向线段表示信息流用方框表示设备或节点信息流程→过程矩阵将信息流程数字化使计算机可以识别根据信息流图可以得出过程矩阵 2．1 过程系统模拟的基本方法过程系统模拟计算量大且复杂手工计算难以完成计算机和计算技术的发展为过程系统的整体研究提供了技术手段各种类型的过程系统模拟软件不断出现但就其模拟计算求解方法而言可以归纳为三类序贯模块法 Sequentia1 Modular Method 面向方程法 Equation Oriented Method 联立方程法联立模块法 Stmultaneously Modular Method 2 11过程系统模拟的序贯模块法序贯模块法按照由各种单元模块组成的过程系统的结构序贯的对各单元模块进行计算从而完成该过程系统的模拟计算的方法序贯模块法对过程系统的模拟以单元模块的模拟计算为基础依据单元模块入口的物流信息以及足够的定义单元特性的信息计算出单元出口物流的信息序贯模块法的优点与实际过程的直观联系强模拟系统软件的建立维护和扩充都很方便易于通用化计算出错时易于诊断出错位置序贯模块法的主要缺点计算效率较低尤其是解决设计和优化问题时计算效率更低序贯模块法计算效率低的原因只能根据模块的输入物流信息计算输出物流信息在进行系统模拟的过程中对有再循环物流单元模块的计算需要考虑断裂物流收敛计算使问题复杂 2 12 过程系统模拟的面向方程法面向方程法将描述整个过程系统的数学方程式联立求解从而得出模拟计算结果的方法面向方程法又称联立方程法面向方程法的优点可以根据问题的要求灵活地确定输入输出变量而不受实际物流和流程结构的影响模型中所有的方程可同时计算和同步收敛面向方程法的问题形成通用软件比较困难不能利用现有大量丰富的单元模块缺乏与实际流程的直观联系计算失败之后难于诊断错误所在对初值的要求比较苛刻计算技术难度较大等 2 13 过程系统模拟的联立模块法联立模块法将过程系统的简化模型方程与单元模块严格模型交替求解又被称作双层法 2．2 过程系统模拟的序贯模块法 2．2．1序贯模块法的基本原理单元模块依据相应过程单元的数学模型和求解算法编制而成的子程序如图28 a 中的闪蒸单元可依据闪蒸单元模型和算法编制成闪蒸单元模块单元模块的单向性结定单元模块的输入物流变量及参数可计算出相应的输出物流变量但不能由检出变量计算输入变量也不能由输入输出变量计算模块参数序贯模块法的基本思想从系统入口物流开始经过对该物流变量进入的单元模块的计算得到输出物流变量这个输出物流变量就是下一个相邻单元的输入物流变量依次逐个的计算过程系统中的各个单元最终计算出系统的输出物流计算得出过程系统中所有的物流变量值即状态变量值 2．2．2 再循环物流的断裂当涉及的系统为无再循环流的树形结构时序贯模块法的模拟计算顺序可以按过程单元的排列顺序一一顺利完成用序贯模块法处理具有再循环物流系统的模拟计算时需要用到系统分解断裂 Tearing 和收敛 Convergence 等多项技术 Step1 假定断裂物流S4的变量值然后依次计算单元模块ABC得到物流S4的变量值 Step2利用收敛单元比较S4与S4的相应变量值若不等则改变S4为新的变量值重复Step1过程直到S4与S4两个变量值相等为止问题收敛单元设置在哪个物流处既如何选择断裂物流本问题中不仅可以是物流S4处也可以设置在物流S2或S3处对于复杂系统收敛单元设置的位置不同其效果也将不同究竟设置在何处为好这要通过断裂技术去解决如何得到新的S4变量值如何保证计算收敛如何加快收敛取决于收敛算法还与断裂物流变量的特性有关 2．2．2 再循环物流的断裂 1 断裂的基本概念首先考察方程组的断裂假设有一个由四个方程四个未知变量组成的方程组也可以由另外的方式进行求解例如假设x2的猜值则 f1解出x3 f2解出x4 f3解出x1 最后利用f4来检验最初没定的猜值x2 是否正确如果f4为零则可认为得到了方程组的解若此处的f4 不为零则需修正x2的值再重新进行迭代计算这样可将四维求解问题降阶成了四个一维问题通过迭代计算把高级方程组降阶为低级方程组的办法称为断裂考察过程系统中的不可分隔子系统如图211断裂物流可以选为S10当然也可以选为S11选择不同的断裂物流则其相应的迭代序列也不一样从表面上看上列的两种计算序列似乎没有什么很大的区别但由于系统中各物流及其变量特性的不同在收敛计算上常是有很大差异的如变量个数的多少方程求解的难易程度等如何选择断裂物流确定迭代序列是实施序贯模块法进行过程系统模拟计算过程中必须要解决的问题 2 断裂方法的研究早在20世纪60年代初就有人提出了断裂的思想此后随着流程模拟技术的不断发展有关研究断裂的文章不断出现他们提出判断最佳断裂的准则分为四类 1 断裂的物流数最少 2 断裂物流的变量数最少 3 断裂物流的权重因子之和最少 4 断裂回路的总次数最少另一种归纳 1断裂的流股数目最少 2断裂流股包含的变量数目最少 3对每一流股选定一个权因子该权因子数值反映了断裂该流股时迭代计算的困难程度应当使所有的断裂流股权因子数值总和最小4选择一组断裂流股使直接代入法具有最好的收敛特性四条准则是一般性的原则 3 回路矩阵过程系统中的简单回路可以用回路矩阵 1oop／stream Matrix 表示矩阵中的行代表回路列代表物流若某回路i中包括有物流J则相应的矩阵元素aij＝1否则为空白或零不独立的列 f 1 与 f 值较大的列相比较若某列中的非零元素与 f 值较大列的非零元素同行则该列相对于 f 值大的列不独立如S2的f 值较大与其余小于它的列相比较会发现S2的非零元素为C行和A行而S1列C行非零 S3A行非零其余列中无与S2同行的非零的元素则判别出 S1 S3相对于S2不独立表示为 S1 S3 S2 S5 S6 S4 流股断裂方法一L - R 分解法 L – R分解法遵循的原则断裂流股数目最少且将所有循环路打开例现有一个为最大循环网的不可分割子系统其信息流图如下1 42 53 S4 S3 S2 S1 S6 S5 S7 S8 4流股断裂方法分析在这个信息流程图中有 8个流股S1S2 S8 五个节点12345构成了ABCD四个环路 1 4 2 5 3 S4 S3 S2 S1 S6 S5 S7 S8 A D C B在Lee – Rudd 法中首先分析信息流图再用环路矩阵表示出来 A B C D 环路S1 S2 S3 S4 S5 S6 S7 S8 01 1 0 0 0 0 0 0 00 0 0 0 1 1 1 1 01 0 0 0 0 0 0 0 11 1 1 0 流股 f R 1 42 53 S4 S3 S2 S1 S6 S5 S7 S8A C DB 矩阵做法Si 流股若在 A 环中出现则标 1若不出现则标 0例如 A 环由S2S3 两流股构成其余为零矩阵中还有加和行用f 表示它由每一列中的非零元素加和构成加和列R它将每一行非零元素加和构成 f 称为环路频率代表某流股出现在所有环路中的次数R 称为环路的秩代表某环路中包含的流股总数经运算可得出加和 f 和R值环路矩阵成为下面样子 A B C D S1 S2 S3 S4 S5 S6 S7 S8 0 1 1 00 0 0 0 0 0 0 0 0 0 1 1 1 1 0 1 0 0 0 0 0 0 0 1 1 1 10 R 2 2 3 4 f 1 2 1 2 1 1 2 1 不独立的列 A B C D S1 S2 S3 S4 S5S6 S7 S8 0 1 1 0 0 00 0 0 0 0 0 0 0 11 1 1 0 1 0 0 0 00 0 0 1 1 1 1 0 R 2 2 34 f 1 2 1 2 1 1 2 1 不独立的列基本概念工艺流程图过程流程过程拓扑举例信息流图-13 序贯模块法的基础是单元模块子程序通常单元模块与过程单元是一一对应的过程单元的输入物流变量即为单元模块的输入单元模块的输出即为过程单元的输出物流变量如 A B H G F E C D 系统分解对复杂系统将所有模型方程全部联立求解很困难直接用序贯法又存在相互影响这时可将该系统分成几个相对独立的部分各自联解再序贯求解将大的复杂系统分解为若干个小的子系统的过程称为大系统的分解目的是识别出不可分割子系统 AB H G F ECD 不可分割子系统不相关子系统 A B H G FE C D A B C A B CG F E D 流股断裂 Tearing 一般对于大系统分解得到的子系统已是不可分隔的如ABC构成的当这样的子系统仍很复杂时联立求解仍困难若断开某一个流股则可采用序贯法求解而断开的流股变量则作为迭代变量选择断裂流股是该技术的关键 A B H G F E C D 断裂物流迭代计算步骤如下该方程组可以通过联立求解得到它的解图210 描述了断裂的过程其中流股x2称为断裂流股该流股只有一个变量x2 称为迭代变量流股的收敛性指的就是其中变量x2 的收敛性能问题如果不选择流股x2是否可达到简化的目的。
收稿日期:2024-02-21。
基金项目:(52200146);(2022-PJLH-03)。
作者简介:(1988-),,,,、、X 。
E-mail:***************通信作者:(1984-),,,,。
E-mail:***************.cn绿氢-绿氨工艺流程的多稳态模拟和动态响应李默1,2,郭永坚3,4,董晓莹5(1.(),100176;2.(),100102;3.(),200062;4.,200237;5.,124221)摘要:-、,、、。
-、、、、,AVEVA Process Simulation --,。
,、,,。
关键词:;;;;;中图分类号:TK89;TM615文献标志码:A文章编号:1671-5292(2024)04-0433-070引言,、。
“”“”,,、、[1-3]。
,(、)、、“-”、,[4,5]。
,50,800t [6]。
,,。
-,,,。
,,、,。
,-,、、,、、、、、,,[7,8]。
,-[9,10]、[11,12]、[13,14]、[15-17],,、、,。
-AVEVA Process Simulation (APS)--,、、、、、,,36h 、、,。
1工艺模型的建立与求解1.1工艺流程-1。
可再生能源Renewable Energy Resources第42卷第4期2024年4月Vol.42No.4Apr.2024·433·WF、SF ,、、,,B-1,C()。
,G,、。
EYALK,EYPEM,。
D-1,D-2,K-1,H-1,。
,,,-,,。
,,,H-1K-1。
K-1E-1、D-3,,E-2,D-4,R-1,D-5,K-2,D-6。
1.2求解方式APS ,,-。
-,、、,,,、、、,“”,、、,、。
SRK 。
2基础数据及工艺参数2.1气象数据The National Solar Radia-tion Data Base(NSRDB)[18](120.35°,43.58°)2020529836h WF-;SF-;B-1-;P-1-;H-1-;R-1-;K-1-;K-2-;E-1-;E-2-;E-3-;D-1-;D-2-;D-3-;D-4-;D-5-;D-6-;EYALK-;EYPEM-可再生能源2024,42(4)图1绿氢-绿氨工艺的原则流程Fig.1Principal flowsheet diagram of green hydrogen-green ammonia processP-1D-1WF SFB-1H-1R-1K-1K-2E-1E-2E-3D-2D-3D-4D-5D-6EYALKEYPEMC3动态场景分析3.1绿氢-绿氨生产完全随气象条件变化波动-。
3.7 过程系统模拟的应用实例化工流程模拟软件的种类虽然较多，但针对化工过程模拟的基本思想和方法却相通。
Aspen Plus应用非常广泛，下面我们以乙苯生产过程乙苯精馏塔生产实例为例介绍Aspen Plus稳态模拟的基本方法和基本思想。
某干气制乙苯装置中，干气中的乙烯与苯催化反应生成乙苯，反应产物中同时包括反应副产品丙苯、二乙苯、多乙苯和没有反应完全的苯。
反应物经过苯塔分离脱除苯后进入乙苯塔、丙苯塔获得乙苯、丙苯产品。
已知进入乙苯塔物料含苯0.2%（摩尔分数，下同）、乙苯77%、丙苯10%、二乙苯12.8%，流量8679kg/h，压力1.45MPa（绝压，下同），温度281℃。
乙苯塔为93块浮阀塔盘，进料位置为第40块塔板，操作中，塔顶压力0.52MPa，冷凝器后压力0.5MPa，塔顶温度210℃，塔底温度254℃，回流量19000kg/h，塔顶产品中乙苯含量99.5%，釜液中乙苯含量1.5%，塔顶为全凝器，塔底为热虹吸式再沸器。
生产中对乙苯产品指标有严格限制，要求乙苯产品纯度大于99.6%，同时为较小物耗要求乙苯塔釜液乙苯含量小于1%，试通过模拟分析提出改进方案。
使用Aspen Plus进行模拟的基本步骤：（1）选择模板（2）选择运行类型（3）创建一个流程（4）规定计算的全局信息（5）规定组分（6）选择物性方法（7）输入物流规定（8）输入模型规定（9）运行模拟（10）检查结果（11）灵敏度分析（12）生成报告1. 建立一个新的运行当启动Aspen Plus并建立一个新的模拟时，可以从一个空白模拟着手或者从一个模板着手，见图3-27。
模板设定了特定工业通常使用的缺省项包括 测量单位、 所要报告的物流组成信息和性质、 物流报告格式、 对游离水选项的缺省设置、 性质方法、 其它特定的应用缺省。
Aspen Plus 内置以下列模板：● 空气分离● 化学工艺● 电解质● 气体加工● 一般工艺● 湿法冶金● 石油● 医药● 冶金● 固体● 特种化工对于每个模板，可以选择米制或英制作为缺省测量单位制，其它单位制也可用。
稳态工况法检测流程
稳态工况法（Steady State Test）是一种用于检测发动机排放的方法，通过在特
定的工况下对发动机进行测试，可以准确地测量其排放物的含量，从而评估其环保性能。
本文将介绍稳态工况法的检测流程，以帮助读者更好地了解这一检测方法。
首先，进行预检测准备工作。
在进行稳态工况法检测之前，需要对测试设备进
行检查和校准，确保其正常工作。
同时，还需要对待测发动机进行准备，包括清洁和调整，以确保测试结果的准确性。
接下来，进行稳态工况法测试。
在测试过程中，需要按照预定的工况条件对发
动机进行稳态运行，同时使用适当的测试设备对排放物进行采样和分析。
这一过程需要严格控制各项参数，确保测试结果的准确性和可比性。
测试完成后，进行数据处理和分析。
通过对采集到的数据进行处理和分析，可
以得出发动机排放物的含量，并进行评估。
同时，还可以对测试结果进行统计和比对，以验证测试的可靠性和准确性。
最后，编制测试报告。
根据测试结果，编制详细的测试报告，包括测试方法、
测试过程、测试结果和评估结论等内容。
这一报告将作为评估发动机环保性能的重要依据，因此需要准确、清晰地呈现测试结果和评估结论。
总的来说，稳态工况法检测流程包括预检测准备、测试过程、数据处理和分析、以及编制测试报告等步骤。
通过严格控制每个环节，并确保测试设备和方法的准确性和可靠性，可以得出准确的发动机排放物含量，并对其环保性能进行评估。
希望本文能够帮助读者更好地了解稳态工况法检测流程，提高对发动机环保性能评估的理解和认识。
