分子筛吸附脱水装置再生工艺计算与设计要点
2026/9/19 15:16:06 网站建设 项目流程

简介:这是一份针对天然气分子筛吸附脱水工艺的课程设计报告,聚焦再生工艺的设计计算,适用于油气储运、石油工程等专业学生及相关工艺设计人员。报告从分子筛性质与吸附原理入手,明确双塔交替操作周期,系统给出吸附水量、再生气与冷凝气气量以及加热炉热负荷的计算过程,为深度脱水装置的再生方案提供完整参考。压缩包内为1个doc格式文档,约467KB,内容包含任务书、天然气物性基础、工艺流程选择、再生加热与冷却操作说明,以及再生工艺计算原稿,结构清晰,便于直接对照学习和修改使用。目前已有71人在线学习,对正在完成《油气集输工程》课程设计或从事天然气脱水初步设计的读者,具有较高的实用参考价值。

1. 分子筛吸附脱水装置能不能长期运行,先看再生工艺计算

一套 12 万方/天的分子筛吸附脱水装置,设计压力 5 MPa,进料温度 26 ℃,要求干气水含量低于 1 ppm。多数人盯住的是分子筛装填量和吸附塔高度,真正决定这套装置能连续跑多久的,反而是再生工艺计算:再生气流量给多大、加热炉热负荷按多少选、冷吹到什么温度切换。算小了,床层再生不彻底,下一个吸附周期露点提前击穿;算大了,加热炉和再生气系统选型整体上浮,运行费用每年多出几十万。这份设计把再生过程拆成了一条可手算的完整链条,从双塔 8 小时周期排产、再生热负荷拆分,到再生气量、冷却气量和空塔速度校核,每一步都有明确参数和检验逻辑,值得按工程习惯复算一遍。

2. 4A分子筛选型与双塔循环时间表:孔径匹配决定再生负荷

2.1 从分子直径表看为什么选4A不选13X

分子筛吸附脱水选型,最先看的不是吸附容量,而是孔径和原料气里各组分分子直径的匹配关系。天然气脱水场景中常见分子的公称直径如下:

分子公称直径(×10⁻¹⁰ m)4A分子筛(4.2~4.7 Å)能否进入
H₂O3.2可进入
CO₂2.8可进入
H₂S3.6可进入
CH₄4.0临界,少量共吸附
C₂H₆4.4基本不进入
C₃H₈4.9不进入

水的公称直径 3.2 Å,CO₂ 是 2.8 Å,H₂S 是 3.6 Å,三者都小于 4A 分子筛的孔径上限,因此 4A 能把水、CO₂、H₂S 一并吸附。这里有一个容易被忽略的细节:3A 分子筛孔径只有 3.2~3.3 Å,H₂S 直径 3.6 Å 进不去,CO₂ 的 2.8 Å 能进。如果选择 3A,脱水没问题,但装置就失去了对酸性气体的缓冲能力。设计里明确提到,上游脱碳装置短时间波动导致天然气不达标时,分子筛还能起到一定净化作用,这正是选 4A 而非 3A 的工程理由。

13X 分子筛孔径约 10 Å,几乎能把 C₃⁺重烃也吸进去,再生时重烃在高温下脱附困难,床层容易结焦失活,所以在单纯深度脱水的场合,13X 不是首选。4A 的 4.0 Å 孔径临界点刚好卡住甲烷,原料气中甲烷共吸附量有限,干气损失可控,这是它成为天然气脱水主流吸附剂的核心原因。

2.2 8小时吸附周期的粒度划分与双塔轮换逻辑

吸附和再生是同一个塔在时间上的两种状态。吸附时床层温度会升高,设计里进料 26 ℃、吸附终了床层 35 ℃,温升 5 ℃就是吸附放热的直接体现;再生利用升温脱附原理,让热干气把床层加热到 200 ℃以上,把水从分子筛孔道里赶出来。这种“低温高压吸附、高温低压脱附”的循环,决定了操作周期只能按小时级来排。

下面是双塔交替运行的时间轴,一个周期内 A 塔吸附,B 塔同时完成加热再生、冷吹和切换:

# 双塔8小时循环排产:A塔吸附,B塔再生+冷却 schedule = [ ("0.0~4.5h", "A塔吸附", "B塔加热再生"), ("4.5~7.7h", "A塔吸附", "B塔冷吹降温"), ("7.7~8.0h", "A塔吸附", "B塔升压切换"), ] for slot, a_tower, b_tower in schedule: print(f"{slot} | {a_tower:8s} | {b_tower}")

这段排产逻辑明确了 8 小时周期的三个时间窗口:加热再生 4.5 小时,冷吹 3.2 小时,升压切换 0.3 小时,合计刚好 8 小时。设计里 24 小时周期的分配比例不同,加热占 65%~68%,冷却占 30%,切换占 2%~5%;8 小时周期则压缩了加热和冷却的相对占比,对吸附塔切换阀的动作时间和压力平衡要求更高。实际操作中,如果原料气含水量波动大,为保证干气露点,可以把操作周期从 8 小时缩短到 6 小时甚至 4 小时,但再生加热和冷却时间必须同步压缩,否则再生不彻底。

3. 再生热负荷拆解:钢材、分子筛、瓷球与脱附热的工程配比

3.1 热负荷四件套公式与关键数据

再生热负荷不是只算把水蒸干需要的热量,而是四部分相加:加热分子筛本体、加热吸附器钢材、提供吸附水解吸热、加热床层底部铺垫的瓷球。设计中的参数取值如下:

  • 吸附终了床层温度 t₁ = 35 ℃,再生终点温度 t₂ = 230 ℃;
  • 分子筛装填量 2265 kg,比热 0.96 kJ/(kg·℃);
  • 吸附器筒体 1529.5 kg,吸附件等钢材 2555.5 kg,钢材合计 4085 kg,比热 0.5 kJ/(kg·℃);
  • 瓷球 401.9 kg,比热 0.88 kJ/(kg·℃);
  • 水的脱附热按 4186.8 kJ/kg;
  • 总热负荷考虑 10% 热量损失。

吸附器各部分质量参数汇总如下:

部件/介质质量(kg)比热(kJ/(kg·℃))35→230 ℃热负荷(MJ)占比
分子筛2265.00.9642423.5%
筒体+吸附件钢材4085.00.5039822.1%
吸附水脱附218.084186.8(kJ/kg)91350.6%
瓷球401.90.88693.8%
合计1804100%

这里有一个必须注意的数据问题。原设计表格里吸附水质量写成 2180.8 kg,按 12 万方/天处理量、含水 0.5636%(摩尔分数)反推:8 小时处理原料气约 40000 m³,其中水汽约 225 m³,折合液态水约 181 kg。218.08 kg 与这个量级吻合,2180.8 kg 明显是手误多敲了一位。参与工程计算时,表上数据和气质分析矛盾,优先以物料衡算为准。

3.2 用Python把热负荷算到可复核

把上述参数写成计算脚本,每一步都能对照检查:

# 再生热负荷计算:分子筛 + 钢材 + 吸附水 + 瓷球 m_zeolite = 2265.0 # 分子筛装填量,kg m_steel = 1529.5 + 2555.5 # 筒体+吸附件钢材,kg m_ceramic = 401.9 # 瓷球,kg m_water = 218.08 # 单周期吸附水量,kg(按物料衡算修正) water_ratio = 0.005636 # 原料气含水摩尔分数 q_gas = 120000 / 24 * 8 # 8小时原料气量,m3 m_water_check = q_gas * water_ratio / 22.4 * 18 print(f"物料衡算水量: {m_water_check:.1f} kg") t1, t2 = 35.0, 230.0 dt = t2 - t1 cp_zeolite, cp_steel, cp_ceramic = 0.96, 0.50, 0.88 q_des = 4186.8 # 水脱附热,kJ/kg Q1 = m_zeolite * cp_zeolite * dt Q2 = m_steel * cp_steel * dt Q3 = m_water * q_des Q4 = m_ceramic * cp_ceramic * dt Q = (Q1 + Q2 + Q3 + Q4) * 1.10 # 10%热损失修正 Qh = Q / 4.5 # 加热时间4.5h print(f"分子筛加热: {Q1/1e3:.1f} MJ") print(f"钢材加热: {Q2/1e3:.1f} MJ") print(f"水解吸热: {Q3/1e3:.1f} MJ") print(f"瓷球加热: {Q4/1e3:.1f} MJ") print(f"总热负荷: {Q/1e6:.3f} GJ,每小时 {Qh/1e3:.1f} MJ/h")

代码的逻辑是先把四个分项分别按“质量 × 比热 × 温升”计算,吸附水解吸热单独按“质量 × 脱附热”计算,再加 10% 热损失,最后除以 4.5 小时得到每小时热负荷。物料衡算那一步把水蒸气的标况体积折算成质量,用来验证吸附水量取值是否合理。注意 22.4 是 0 ℃标准状态下的摩尔体积,如果按 20 ℃标态应取 24.06,差异不大,不影响数量级判断。

3.3 再生气量为什么优先用焓差法

再生气量计算是再生工艺里最容易出分歧的地方。设计给出两个依据:再生气在 230 ℃时的平均比热取 3.14 kJ/(kg·℃),入口 260 ℃、出口 200 ℃,温降 60 ℃;同时又给出了再生气在 260 ℃和 115 ℃下的热焓值,分别是 -3776.58 kJ/kg 和 -4167.30 kJ/kg。

# 再生气量对比:平均比热法 vs 热焓差法 Qh = 441.0 * 1000 # 上面算出的每小时热负荷,kJ/h(按修正水量) # 方法一:平均比热法,温差60℃ G_cp = Qh / (3.14 * (260 - 200)) print(f"平均比热法: {G_cp:.0f} kg/h") # 方法二:热焓差法,取260℃与115℃焓值 dh = 4167.30 - 3776.58 # 390.72 kJ/kg G_enthalpy = Qh / dh print(f"焓差法: {G_enthalpy:.0f} kg/h")

两种方法算出来的结果会差约一倍:平均比热法倾向算大,焓差法则更贴近原设计给出的 1083.3 kg/h。原因在于 3.14 kJ/(kg·℃) 是在 230 ℃附近取的平均定压比热,而再生气从 260 ℃降到 200 ℃的实际焓变并不是“平均比热 × 温降”这么线性。工程上我的习惯是:只要工艺包或物性包能给出操作温度区间的焓值,就用焓差法;平均比热法只在初步估算时使用,而且要留 1.5~2 倍的余量再选设备。

4. 冷却气量与空塔速度校核:冷吹热吹为什么取同一流量

4.1 冷却负荷计算与冷吹气量确定

加热再生结束后,床层温度约 230 ℃,需要冷吹降温到 30 ℃才能切换回吸附状态。冷却负荷和加热负荷不同,不需要再计算水解吸热,因为水已经在加热阶段脱除,冷吹只是把床层固体材料的热量带走。

# 冷却负荷:床层230℃降到30℃,不含脱附热 Q_cool = (m_zeolite * 0.96 + m_steel * 0.50 + m_ceramic * 0.88) * (230 - 30) Q_cool *= 1.10 # 10%冷却裕量 Q_cool_h = Q_cool / 3.2 # 冷却时间3.2h # 冷却气平均比热2.9 kJ/(kg·℃),进出口温差按100℃估算 G_cool = Q_cool_h / (2.9 * 100) print(f"冷却气量: {G_cool:.0f} kg/h") # 换算20℃标态体积流量,再生气分子量按17 G_cool_nm3 = G_cool / 17 * 24.06 print(f"标态体积流量: {G_cool_nm3:.0f} Nm3/h")

计算逻辑与加热负荷一致,只是把温差从 195 ℃换成 200 ℃,并去掉吸附水项。冷吹气量的结果和再生气量几乎是同一个数量级,原设计给出的冷吹气量约为 1427.4 m³/h(20 ℃标态),换算成质量流量约 1083 kg/h,与再生气量基本相等。设计最后明确了一个工程做法:冷吹气量和热吹气量虽然理论上不等,但实际操作中为了装置平稳,统一采用 1427.4 m³/h。这个取法很实用,切换时不需要大幅调整流量,可以避免双塔切换瞬间系统压力波动。

4.2 空塔速度C值陷阱:0.167与0.25~0.32的方向差异

空塔速度校核公式是经典的分子筛吸附器流速公式:

允许质量流速 = C × √(ρg × ρs × dp)

其中 ρg 是操作条件下气体密度,ρs 是分子筛堆密度(660 kg/m³),dp 是分子筛直径(3.2 mm 球形),C 是方向系数。这里最容易被新手上手就搞错的是 C 值:气体自上而下流动时 C 在 0.25~0.32,自下而上流动时 C 只取 0.167。原因是向上的气流对床层有举升作用,流速过大会导致分子筛流化或磨损。

设计里再生加热气和冷却气都是自下而上流动,所以 C 统一取 0.167。代入参数校核:

工况气体密度(kg/m³)C值允许质量流速(kg/(m²·s))允许流量(kg/h)实际流量(kg/h)
加热气 260 ℃17.410.167约 1.01约 9300约 1083
冷却气 30 ℃33.880.167约 1.41约 13000约 1083
# 空塔速度校核:按允许质量流速反算最小床层截面积 rho_g = 17.41 # 加热气操作密度,kg/m3 rho_s = 660.0 # 分子筛堆密度,kg/m3 dp = 0.0032 # 球形分子筛直径,m C = 0.167 # 自下而上流动 G_allow = C * (rho_g * rho_s * dp) ** 0.5 A_need = (1083.3 / 3600) / G_allow print(f"允许质量流速: {G_allow:.3f} kg/(m2·s)") print(f"所需面积: {A_need:.3f} m2,设计面积: 2.56 m2")

从校核结果看,2.56 m² 的床层面积远大于所需面积,余量在 8 倍以上。这个结果说明再生工况不是床层直径的控制条件,吸附塔直径主要由吸附周期的处理气量决定。反过来想,这个富余系数也意味着再生流量在正常设计值附近还有较大调整空间,可以用于缩短加热时间或应对原料气含水量短期升高。

提示:C 值取 0.167 时校核通过,不代表可以把流量无限放大。流速过高会导致床层跳动、分子筛粉化,尤其是自下而上再生的塔,压差波动大时要优先检查再生气进口分布器和床层压降,而不是先怀疑分子筛质量问题。

5. 用116~120℃脱附平台判断再生终点:再生气温差与DCS报警设置

再生过程结束得恰到好处,是分子筛长寿命运行的前提。设计里的温度曲线分四段:A 段是轻烃脱附,B 段是吸附水大量脱附,C 段是重烃脱附,D 段是冷吹降温。其中 B 段有个非常明显的工程特征:床层出口温度在 116~120 ℃区间会出现一个平台,因为水脱附吸热,热再生气提供的热量被大量消耗,升温速率明显变缓。平台的存在和持续长度,直接反映床层内主体水还剩多少。

实际操作中,我一般把出口温度的一阶差分作为判断依据:

# 每1分钟采集一次再生气出塔温度,动态判断再生阶段 t_out = 119.0 # 当前出口温度 t_out_prev = 118.6 # 1分钟前温度 dT_dt = t_out - t_out_prev if 110 < t_out < 125 and dT_dt < 0.8: print("水脱附平台确认,主体水正在脱除") elif t_out > 150 and dT_dt < 0.3: print("重烃脱附阶段,接近加热终点") elif t_out >= 195: print("再生加热完成,可切换冷吹")

这段逻辑的关键在两个阈值:升温速率 0.8 ℃/min 以内识别平台,出口温度 195 ℃以上且升温速率降到 0.3 ℃/min 以内判断加热结束。把这两个条件写进 DCS 报警,比单纯定时切换可靠得多。设计原文给出的热再生气入口 260 ℃、出口 200 ℃,温差 60 ℃偏大;一般推荐入口温度比最终出口温度高 19~55 ℃,典型取 38 ℃。温差过大的代价是床层局部过热,分子筛骨架铝流失加快,使用寿命缩短。若想兼顾加热时间,可以把入口温度降到 238 ℃左右,加热时间可能从 4.5 小时延长到 5 小时以上,这时要重新校核 8 小时周期内是否还有足够的冷吹和切换时间。

冷吹终点同样建议以温度为准而不是以时间为准。床层出口降到 50 ℃左右即可切换吸附;如果冷吹气不是真正的干气,应让冷却气自上而下流过床层,让顶部分子筛先把冷吹气中的水蒸气吸附掉,这样在下一轮吸附周期里,干气出口露点不会因为床层上部预饱和而提前击穿。

本文还有配套的精品资源,点击获取

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询