常减压蒸馏过程实时动态模拟

符号说明
q F 进料摩尔流量, kmol/ h U 持气量或持液量, kmol Q 对塔板的加热量, kJ/ h H 体系的焓, J/ mol q L 溢流摩尔流量, kmol/ h q 上升蒸汽摩尔流量, kmol/ h q S 液相侧线采出摩尔流量, kmol/ h q G 气相侧线采出摩尔流量, kmol/ h x 液相 摩尔分数 y 气相摩尔分数 z 进料摩尔分数 t 时间, h 下角标 i 组分 j 塔板
qL, j - 1 x i, j - 1 +
qj + 1 yi , j + 1 -
( qL , j + qS , j ) x i , j - ( qj + qG, j ) y i, j
热量平衡方程为:
d(
Uj
H
t j
)
dt
=
qF, jH F, j +
qL , j - 1 hj- 1 +
qj + 1 Hj + 1 -
( qL , j + qS , j ) hj - ( qj + qG, j ) H j - Qj
组分摩尔分数归一方程:
c
c
xi, j=
y i, j = 1 0
i= 1
i= 1
压力为温度、容积、气相滞料量的函数:
Pj = P( T j , V T , UjV )
2. 2 模型求解 对动态模型的计算采用了 跟踪逼近技术 [ 7] ,
流程的实时动态模拟, 其误差基本上控制在 5% 以 内, 完全满足工业装置实时动态模拟的要求。其误 差的产生主要有以下几个原因: 首先石油的组成比 较复杂, 原始实沸点蒸馏数据存在一定误差; 切割实 沸点曲线时选择了 43 个实组分, 虽原油中确实存在 这些组分, 但这 43 个组分仍不能完全代表原油的组 成; 在常减压蒸馏过程模拟时选择不同的状态方程, 也会得出不同的结果, 本文选择了 SRK 方程, 如选 PR 方程或 L ee- Kesler 方程则 结果会 略有 不同, 而 SRK 方程计算速度略快; 从上述图表中可以发现减 压塔的模拟比常压塔的模拟误差大, 这主要是因为 汽液平衡中的相平衡常数与温度和压力关系较为密 切, 而减压塔的温度和压力跨度较大, 故相平衡常数 的计算误差也较大。
第 28 卷 第 1 期 2001 年
北 京化 工大 学学报 JOU RNAL OF BEIJING U NI VERSIT Y O F CHEM ICA L T ECHN OL OGY
Vol. 28, N o. 1 2 00 1
研究简报
常减压蒸馏过程实时动态模拟
李旦杰1) 王健红2) * 方 刚2)
0
105 4 101. 3 4 0 161 5 155 8 3 5
5
144 2 137. 3 4 8 212 3 205 7 3 1
10
154 6 152. 7 1 3 226 7 221 3 2 4
30
ቤተ መጻሕፍቲ ባይዱ
173 7 170 1 2 1 251 5 256 1 1 8
50
187 1 176. 7 5 6 269 1 257 5 4 3
表 2 实沸点蒸馏曲线数据实测值与计算值比较 T able 2 Comparison between the actual and calculated
values of t rue boiling curve
( 实沸点 蒸馏) / %
常一线对应温度 t 1/
常二线对应温度 t2/
实测 计算 误差/ % 实测 计算 误差/ %
1 选定组分的方法
原油是复杂烃类的混合物, 所含组分很多, 而且 往往含有某些非烃组分, 测定其准确组成极为困难。
收稿日期: 2000- 04- 12 第一作者: 男, 1961 年生, 高级工程师, 工学硕士 * 通讯联系人
早期对原油汽液平衡计算一般采用经验方法, 主要 是利用三种蒸馏曲线图( 恩氏蒸馏曲线图、实沸点蒸 馏曲线图和平衡气化蒸馏曲线图) 之间的关系, 得出 原油的泡点、露点以及一系列平衡气化温度与馏出 体积分数之间的关系。有了平衡气化蒸馏数据, 有 关的工艺计算就能进行。但由于这些图表有一定的 局限性, 工艺计算结果不够准确[ 2] 。
于严格模型, 并可使计算效率提高 10 倍至数十倍,
达到超实时的模拟速度。
3 模拟计算结果
常减压蒸馏装置模拟软件的开发以机理模型为 基础, 采用统一的基础物性数据, 通过有效的数值计 算策略求解了由质量衡算、能量衡算、相平衡、流动、 压力平衡等方程构成的微分- 代数方程组, 使数学 模型的计算精度与计算速度得到了很好的统一。精 馏塔及各侧线的模拟结果见表 1、表 2 及图 1。
4 结束语
从模拟结果可见, 利用实组分切割法确定原油 的组分及组成, 并做出以机理模型为基础的精馏塔 模型, 结合常用的状态方程描述原油炼制过程体系 的组成与热力学性质是可靠易行的。如运用跟踪逼 近法, 则可大大提高计算速度, 满足实时动态计算的 需要。该方法的使用, 对开发炼油装置过程模拟系 统, 有较高应用价值。
要求。针对以上问题, 设计了一种新的动态模型求
解策略, 对每一级塔板进行闪蒸计算, 求出每级塔板
的压力和滞料量, 再由水力学计算解出塔内的气液 相流量。实际计算表明该模型能较好地解决上述问
题。因此根据塔板物理模型, 可写出相应的动态模 拟数学方程[ 7] 。
其物料平衡方程为:
d
U i, dt
j
=
qF, j z i, j +
实组分切割法能较好地反映石油馏分组成和性 质。实沸点蒸馏要求使用比工业生产装置分离能力 高的设备进行蒸馏, 可按每馏出 3% ( 体积分数) 或 每隔 10 取一个馏分, 得到馏出温度与馏出量关系 数据。所得到的实沸点 蒸馏数据是一 系列离散的 点, 为了把这些点联成光滑曲线, 采用了插值函数法 来处理这些数据。通过比较多项式插值和样条函数 插值发现, 当节点较多时, 拉格朗日多项式对被插值 函数的逼近较差, 特别是边界点附近波动较大; 而用
93 2
动态模拟数据( 运行至稳态)
p / MPa
qm / ( t h - 1)
0 05
1 755 2
1 23
0 04
28 7
11
- 0 094 8
48 2
qv / ( m3 h- 1) 665 1 332 1
第 1期
李旦杰等: 常减压蒸馏过程实时动 态模拟
83
由表 1、表 2 和图 1 可以看出, 运用 实组分切 割法 处理了石油馏分, 对常减压蒸馏过程实现了全
82
北京化工 大学学报
2001 年
样条函数插值法则能避免这些问题, 使曲线具有较 好的连续性、光滑性[ 3] 。
用三次插值样条函数法对所得的实沸点蒸馏数 据进行处理, 可得到一个二阶光滑的分段三次多项 式。这个多项式表达了一个光滑连续的实沸点蒸馏 曲线, 因此可以利用实组分数据库中的沸点值找到 与之相对应的体积分数。利用构造出的实沸点蒸馏 曲线, 并根据原油的实际情况, 大约 20~ 30 选择 一种实组分作为石油馏分, 进行组分切割, 选择的实 组分主要是正构烷烃和常见芳烃。所做大量数值试 验证明, 不考虑异构等复杂分子作为实组分不至造 成明显误差, 这与文献[ 4] 的结论一致。
参考文献
[ 1] 王健红, 姚 飞, 马润宇. 石油 馏分宏观 性质计算的 实 组分切割法. 见: 中国化工 学会化学工 程专业 委员会. 第七届全国化 学工程 论文 报告 会论 文集. 北京: 化 学 工业出版社, 1994. 1211~ 1214
完整循环后, 再按一定准则调用严格模型进行判敛
并更新简化模型参数。因显著减少了调用严格模型
的次数, 可大大加快运算速度。数值实验表明, 因积
分时间 dt ( 0 5 s) 很小, 简化模型参数的变化远远小
于模型误差、截断误差和舍入误差, 故由此所得积分
曲线与全部使用严格模型所得积分曲线几乎没有差
别。在整个积分过程中, 简化模型动态地逼近、收敛
在石油馏分体系的过程模拟领域, 通常沿用虚 拟组分法来处理密度、焓、逸度、粘度、导热系数等性 质计算的问题。虚组分切割法是把原油的实沸点蒸 馏曲线分割为有限个窄馏分, 每个窄馏分当作纯组 分来处理, 以窄馏分的平均性质计算而得到该组分 的性质。严格地说, 每个窄馏分的平均性质, 在纵坐 标上的位置应按图解积分来确定。但当窄馏分宽度 足够窄, 在其沸程内的蒸馏曲线接近直线时可取该 窄馏分的中沸点。该法在解决蒸馏过程模拟问题时 起到了很大作用。但虚组分是客观上并不存在的组 分, 无论所选馏分怎样窄, 都难与实际组分相对应。 其相对分子质量、蒸气压等物性关联式来自经验, 适 用范围窄、外推性欠佳, 出现重大误差的可能性大。
即采用具有物理意义的逸度和焓简化模型逼近严格
模型, 求解时交替使用严格模型和简化模型, 简化模
型在必要时由严格模型校正其参数, 并且收敛于严 格模型。
实际求解过程为: 在积分计算和迭代计算时使
用解析性好的简化模型计算相平衡和焓、比热, 模型
其它部分仍为严格模型, 故可采用高阶收敛方法, 克 服数值误差, 极大提高计算速度; 在完成一个或数个
t/ 100 4 195 3 104 2 349 1
87 8
实测数据
p / MPa 0 05
qm / ( t h- 1) 1 768 1
1 23 0 04 11 - 0 093 6
29 6 50 1
qv / ( m3 h- 1) 673 9 349 2
t/ 101 5 196 8 105 8 358 3
表 1 为精馏塔主要参数对比值。表 2 根据动态 模拟的结果, 将各侧线油计算结果的实组分组成转 化为实沸点蒸馏数据, 并与实际运行结果比较。