CN109173584A - 一种基于ldf模型的psa优化策略 - Google Patents
一种基于ldf模型的psa优化策略 Download PDFInfo
- Publication number
- CN109173584A CN109173584A CN201811180037.1A CN201811180037A CN109173584A CN 109173584 A CN109173584 A CN 109173584A CN 201811180037 A CN201811180037 A CN 201811180037A CN 109173584 A CN109173584 A CN 109173584A
- Authority
- CN
- China
- Prior art keywords
- oxygen
- bed
- component
- rate
- aperture
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/02—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by adsorption, e.g. preparative gas chromatography
- B01D53/04—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by adsorption, e.g. preparative gas chromatography with stationary adsorbents
- B01D53/047—Pressure swing adsorption
-
- C—CHEMISTRY; METALLURGY
- C01—INORGANIC CHEMISTRY
- C01B—NON-METALLIC ELEMENTS; COMPOUNDS THEREOF; METALLOIDS OR COMPOUNDS THEREOF NOT COVERED BY SUBCLASS C01C
- C01B13/00—Oxygen; Ozone; Oxides or hydroxides in general
- C01B13/02—Preparation of oxygen
- C01B13/0229—Purification or separation processes
- C01B13/0248—Physical processing only
- C01B13/0259—Physical processing only by adsorption on solids
- C01B13/0262—Physical processing only by adsorption on solids characterised by the adsorbent
- C01B13/0274—Other molecular sieve materials
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2253/00—Adsorbents used in seperation treatment of gases and vapours
- B01D2253/10—Inorganic adsorbents
- B01D2253/116—Molecular sieves other than zeolites
Abstract
本发明公开了一种基于LDF模型的PSA优化策略,通过在两个吸附床之间安装可变孔径流量调节阀,利用遗传算法求解下一次吸附时的最佳孔径值,一方面可使氧气浓度在合格的同时,提高氧气回收率,降低能耗;另一方面,避免冲击孔径过小,导致再生吸附效果较差,影响下一次吸附的能力,从而导致生成氧气浓度的降低。
Description
技术领域
本发明涉及变压吸附制氧机领域,具体一种基于LDF模型的PSA优化策略。
背景技术
变压吸附制氧机(简称PSA制氧机),即利用PSA专用分子筛选择性吸附空气中的氮气,排出纯度很高的氧气。吸附常常是在压力环境下进行的,变压吸附提出了加压和减压相结合的方法,它通常是由加压吸附、减压再组成的吸附-解吸系统。在等温的情况下,利用加压吸附和减压解吸组合成吸附操作循环过程。吸附剂对吸附质的吸附量随着压力的升高而增加,并随着压力的降低而减少,同时在减压过程中,放出被吸附的气体,使吸附剂再生,外界不需要供给热量便可进行吸附剂的再生。在定孔状态下,当吸附压力一定时,冲洗流量基本维持不变,氧气浓度及回收率变化幅度较小,当吸氧量增大时无法在短时间内增加氧气浓度;若将定孔孔径增大,虽然氧气浓度会提高,但氧气回收率会大幅下降,造成能耗的提高,且频繁更换定孔不具备实际意义。
发明内容.
本发明要解决的技术问题是克服现有的缺陷,提供一种基于LDF模型的PSA优化策略,为实现上述目的,本发明提供如下技术方案:一种基于LDF模型的PSA优化策略,包括以下步骤:
Step.1利用扩展Langmuir等温线方程计算氧气的平衡吸附量与氮气的平衡吸附量 其中a为氧气,b为氮气,n为混合气体中的组分数,i为混合气体中的气体组分,i=a,b,bi为i组分的吸附等温线常数,Ci为i组分的气相摩尔浓度,qis为i组分的饱和吸附量,上述参数均为已知;
Step.2;利用LDF模型的质量传递率方程计算当前时刻的吸附率与 其中ki为i组分LDF模型传质系数,qi为分子筛中i组分的吸附量,上述参数为已知;
Step.3计算下一时刻i组分的吸附量qi(i,j+1), 其中Δt为时间步长,取值区间为[0.005,0.01],j为时间序列,j=1,2,3…,上述参数均为已知;
Step.4双组分系统的总质量平衡方程为将理想气体状态方程带入上述总质量平衡方程,得到其中ε为床层的空隙率,z为吸附床轴向坐标,C为混合组分总的摩尔浓度,以上参数均为已知;
Step.5根据能量守恒关系,得到有关床温变化的能量平衡方程, 其中λ为床壁的导热系数,cg为空气比热,ρg为空气的密度,cs为吸附相的摩尔浓度,ρb为氮气的密度,H为吸附热,h为管道的热交换系数,Tw为管壁温度,以上参数均为已知;
Step.6设置温度边界条件,令TZ=0=T入口,L为床体长度,联立Step.4和Step.5的方程可计算出当前时间序列i组分的气体流动速度u(i,j)和当前时间序列i组分的气体温度T(i,j);
Step.7通过冲洗孔孔眼的气体流量为通过孔眼的气体流量受临界压比v比的影响而变化,当时,当时,其中k为氧气的传质系数,P上和P下分别为孔眼的上下游压力,T上为上游气体温度;Step.8将理想气体状态方程带入流体质量平衡方程 得到计算氧气浓度的质量平衡方程 其中T为当前时刻的气体温度,ya为氧气浓度,和分别为氧气和氮气的平均吸附量,DL为轴向扩散系数;
Step.9当床体处于吸附状态时,床体的实际出口流量为Q出=Q产-Q冲,氧气回收率为由于采用双床循环吸附,因此床1的Q冲即为床2解吸的入口流量;
Step.8设实际出口流量的氧气浓度下限为α,氧气回收率下限为β,设当前床1处于吸附状态,床2处于解吸状态时,当冲击孔的孔径增加时,冲洗流量增大,对床2的反向清洗更彻底,再生效果好,可提高床2下一次的氧气浓度,但会使当前床1的出口流量降低,氧气回收率降低,能耗升高;当吸氧量增高时,出口流量与氧气浓度均会下降,缩小冲击孔的孔径从而增大出口流量,氧气回收率同样上升,能耗比下降,但冲洗效果下降会影响下一次生成的氧气浓度;通过安装可变孔径流量调节阀,利用遗传算法对下一次吸附的孔径进行实施调节,在满足氧气实际出口流量及浓度的同时,提高氧气回收率,降低能耗;
Step.9设可变孔径流量调节阀的孔径变化范围为[1mm,10mm],精度为0.0001mm,随机在孔径变化范围内取N个不同的孔径值,N的范围为[500,1000],
Fi最小时的孔径值i,令计算个体的生存概率其中i为第i个孔径值;
Step.10按照个体生存概率,选择N/2个孔径值进行复制,将孔径值转化为二进制,以概率pc进行两两交叉操作,交叉位置为二进制孔径值的第2、3、7位,交叉后的孔径值满足[1mm,10mm];设变异概率为pm,变异时将二进制孔径值的第4、5、7位取反;若则令否则不变;继续迭代100次后取Fi最小时的可求出下一次使为最大时的孔径值i。
与现有技术相比,本发明的有益效果是:通过在两个吸附床之间安装可变孔径流量调节阀,利用遗传算法求解下一次吸附时的最佳孔径值,一方面可使氧气浓度在合格的同时,提高氧气回收率,降低能耗;另一方面,避免冲击孔径过小,导致再生吸附效果较差,影响下一次吸附的能力,从而导致生成氧气浓度的降低。
附图说明
图1为本发明的变压吸附双塔循环制氧机的结构示意图;
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。本发明提供一种技术方案:包括以下步骤:
Step.1利用扩展Langmuir等温线方程计算氧气的平衡吸附量与氮气的平衡吸附量 其中a为氧气,b为氮气,n为混合气体中的组分数,i为混合气体中的气体组分,i=a,b,bi为i组分的吸附等温线常数,Ci为i组分的气相摩尔浓度,qis为i组分的饱和吸附量,上述参数均为已知;
Step.2;利用LDF模型的质量传递率方程计算当前时刻的吸附率与 其中ki为i组分LDF模型传质系数,qi为分子筛中i组分的吸附量,上述参数为已知;
Step.3计算下一时刻i组分的吸附量qi(i,j+1), 其中Δt为时间步长,取值区间为[0.005,0.01],j为时间序列,j=1,2,3…,上述参数均为已知;
Step.4双组分系统的总质量平衡方程为将理想气体状态方程带入上述总质量平衡方程,得到其中ε为床层的空隙率,z为吸附床轴向坐标,C为混合组分总的摩尔浓度,以上参数均为已知;
Step.5根据能量守恒关系,得到有关床温变化的能量平衡方程, 其中λ为床壁的导热系数,cg为空气比热,ρg为空气的密度,cs为吸附相的摩尔浓度,ρb为氮气的密度,H为吸附热,h为管道的热交换系数,Tw为管壁温度,以上参数均为已知;
Step.6设置温度边界条件,令TZ=0=T入口,L为床体长度,联立Step.4和Step.5的方程可计算出当前时间序列i组分的气体流动速度u(i,j)和当前时间序列i组分的气体温度T(i,j);
Step.7通过冲洗孔孔眼的气体流量为通过孔眼的气体流量受临界压比v比的影响而变化,当时,当时,其中k为氧气的传质系数,P上和P下分别为孔眼的上下游压力,T上为上游气体温度;Step.8将理想气体状态方程带入流体质量平衡方程 得到计算氧气浓度的质量平衡方程 其中T为当前时刻的气体温度,ya为氧气浓度,和分别为氧气和氮气的平均吸附量,DL为轴向扩散系数;
Step.9当床体处于吸附状态时,床体的实际出口流量为Q出=Q产-Q冲,氧气回收率为由于采用双床循环吸附,因此床1的Q冲即为床2解吸的入口流量;
Step.8设实际出口流量的氧气浓度下限为α,氧气回收率下限为β,设当前床1处于吸附状态,床2处于解吸状态时,当冲击孔的孔径增加时,冲洗流量增大,对床2的反向清洗更彻底,再生效果好,可提高床2下一次的氧气浓度,但会使当前床1的出口流量降低,氧气回收率降低,能耗升高;当吸氧量增高时,出口流量与氧气浓度均会下降,缩小冲击孔的孔径从而增大出口流量,氧气回收率同样上升,能耗比下降,但冲洗效果下降会影响下一次生成的氧气浓度;通过安装可变孔径流量调节阀,利用遗传算法对孔径进行实施调节,在满足氧气实际出口流量及浓度的同时,提高氧气回收率,降低能耗;
Step.9设可变孔径流量调节阀的孔径变化范围为[1mm,10mm],精度为0.0001mm,随机在孔径变化范围内取N个不同的孔径值,N的范围为[500,1000],设适应度函数为且yai≥a,Fi最小时的孔径值i,令计算个体的生存概率其中i为第i个孔径值;
Step.10按照个体生存概率,选择N/2个孔径值进行复制,将孔径值转化为二进制,以概率pc进行两两交叉操作,交叉位置为二进制孔径值的第2、3、7位,交叉后的孔径值满足[1mm,10mm];设变异概率为pm,变异时将二进制孔径值的第4、5、7位取反;若则令否则不变;继续迭代100次后取Fi最小时的可求出使为最大时的孔径值i。具体实施例
设吸附床数量为2,床长154mm,床体内径122mm,原料气管道9mm,剩余管道8mm,DL=2.1cm2/s,λL=12.5w/m.K,ka=30.5/s-1,kb=11.5/s-1s,h=12.5w/m2.K,ρg=12.5kg/m3,cg=12.5kJ/kg.K,Ha=13.188kJ/mol,Hb=18.942kJ/mol,q=6130mol/m3,Ka=4.7,Kb=14.8,孔径变化范围为[1mm,10mm],精度为0.0001mm,随机在孔径变化范围内取800个不同的孔径值,设适应度函数为且yai≥a,pc=0.94,pm=0.02,迭代次数为100,模拟环境为变压吸附循环30次后,下表为重复10次所取得的最佳孔径值及预测氧气浓度和氧气回收率。
从上表可以看出,当孔径值取7.5时预测氧气浓度及氧气回收率达到最高值,由于遗传算法存在易收敛于局部最优解的缺陷,因此可在计算时增加运算次数,取其中最高值即可。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (1)
1.一种基于LDF模型的PSA优化策略,包括以下步骤:
Step.1针对变压吸附双塔循环制氧机,利用扩展Langmuir等温线方程计算氧气的平衡吸附量与氮气的平衡吸附量 其中a为氧气,b为氮气,n为混合气体中的组分数,i为混合气体中的气体组分,i=a,b,bi为i组分的吸附等温线常数,Ci为i组分的气相摩尔浓度,qis为i组分的饱和吸附量,上述参数均为已知;
Step.2;利用LDF模型的质量传递率方程计算当前时刻的吸附率与 其中ki为i组分LDF模型传质系数,qi为分子筛中i组分的吸附量,上述参数为已知;
Step.3计算下一时刻i组分的吸附量qi(i,j+1), 其中Δt为时间步长,取值区间为[0.005,0.01],j为时间序列,j=1,2,3…,上述参数均为已知;
Step.4双组分系统的总质量平衡方程为将理想气体状态方程带入上述总质量平衡方程,得到其中ε为床层的空隙率,z为吸附床轴向坐标,C为混合组分总的摩尔浓度,以上参数均为已知;
Step.5根据能量守恒关系,得到有关床温变化的能量平衡方程, 其中λ为床壁的导热系数,cg为空气比热,ρg为空气的密度,cs为吸附相的摩尔浓度,ρb为氮气的密度,H为吸附热,h为管道的热交换系数,Tw为管壁温度,以上参数均为已知;
Step.6设置温度边界条件,令TZ=0=T入口,L为床体长度,联立Step.4和Step.5的方程可计算出当前时间序列i组分的气体流动速度u(i,j)和当前时间序列i组分的气体温度T(i,j);
Step.7通过冲洗孔孔眼的气体流量为通过孔眼的气体流量受临界压比v比的影响而变化,当时,当时,其中k为氧气的传质系数,P上和P下分别为孔眼的上下游压力,T上为上游气体温度;
Step.8将理想气体状态方程带入流体质量平衡方程 得到计算氧气浓度的质量平衡方程 其中T为当前时刻的气体温度,ya为氧气浓度,和分别为氧气和氮气的平均吸附量,DL为轴向扩散系数;
Step.9当床体处于吸附状态时,床体的实际出口流量为Q出=Q产-Q冲,氧气回收率为由于采用双床循环吸附,因此床1的Q冲即为床2解吸的入口流量;
Step.8设实际出口流量的氧气浓度下限为α,氧气回收率下限为β,设当前床1处于吸附状态,床2处于解吸状态时,当冲击孔的孔径增加时,冲洗流量增大,对床2的反向清洗更彻底,再生效果好,可提高床2下一次的氧气浓度,但会使当前床1的出口流量降低,氧气回收率降低,能耗升高;当吸氧量增高时,出口流量与氧气浓度均会下降,缩小冲击孔的孔径从而增大出口流量,氧气回收率同样上升,能耗比下降,但冲洗效果下降会影响下一次生成的氧气浓度;通过安装可变孔径流量调节阀,利用遗传算法对孔径进行实施调节,在满足氧气实际出口流量及浓度的同时,提高氧气回收率,降低能耗;
Step.9设可变孔径流量调节阀的孔径变化范围为[1mm,10mm],精度为0.0001mm,随机在孔径变化范围内取N个不同的孔径值,N的范围为[500,1000],设适应度函数为且yai≥a,当Fi最小时,令为计算个体的生存概率其中i为第i个孔径值;
Step.10按照个体生存概率,选择N/2个孔径值进行复制,将孔径值转化为二进制,以概率pc进行两两交叉操作,交叉位置为二进制孔径值的第2、3、7位,交叉后的孔径值满足[1mm,10mm];设变异概率为pm,变异时将二进制孔径值的第4、5、7位取反;若则令为否则不变;继续迭代100次后取Fi最小时的可求出使为最大时的孔径值i。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811180037.1A CN109173584B (zh) | 2018-10-10 | 2018-10-10 | 一种基于ldf模型的psa优化策略 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811180037.1A CN109173584B (zh) | 2018-10-10 | 2018-10-10 | 一种基于ldf模型的psa优化策略 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109173584A true CN109173584A (zh) | 2019-01-11 |
CN109173584B CN109173584B (zh) | 2021-08-24 |
Family
ID=64947929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811180037.1A Active CN109173584B (zh) | 2018-10-10 | 2018-10-10 | 一种基于ldf模型的psa优化策略 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109173584B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112520704A (zh) * | 2020-11-26 | 2021-03-19 | 河南省日立信股份有限公司 | 一种六氟化硫气体回收装置及回收率测量方法 |
CN113691893A (zh) * | 2021-08-09 | 2021-11-23 | 武汉理工大学 | 一种压力机能耗远程智能监控平台及方法 |
CN112520704B (zh) * | 2020-11-26 | 2024-04-26 | 河南省日立信股份有限公司 | 一种六氟化硫气体回收装置及回收率测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1522788A (zh) * | 2003-09-12 | 2004-08-25 | 北京科技大学 | 一种变压吸附分子筛制氧装置 |
CN104331636A (zh) * | 2014-11-26 | 2015-02-04 | 华东理工大学 | 一种c8芳烃模拟移动床分离过程的建模和优化方法 |
-
2018
- 2018-10-10 CN CN201811180037.1A patent/CN109173584B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1522788A (zh) * | 2003-09-12 | 2004-08-25 | 北京科技大学 | 一种变压吸附分子筛制氧装置 |
CN104331636A (zh) * | 2014-11-26 | 2015-02-04 | 华东理工大学 | 一种c8芳烃模拟移动床分离过程的建模和优化方法 |
Non-Patent Citations (5)
Title |
---|
A.NALAPARAJU等: "CO2 captureincation-exchangedmetal–organic frameworks:Holistic", 《CHEMICALENGINEERINGSCIENCE》 * |
S. FAROOQ等: ""Numerical Simulation of a Pressure Swing Adsorption Oxygen Unit", 《CHEMICAL ENGINEERING SCIENCE》 * |
崔红社等: "小型变压吸附制氧过程的数值模拟", 《深冷技术》 * |
戴先知等: "变压吸附空分制氧过程非等温模拟", 《低温工程》 * |
王瑞平主编: "《工程热力学》", 31 August 2009, 西北工业大学出版社 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112520704A (zh) * | 2020-11-26 | 2021-03-19 | 河南省日立信股份有限公司 | 一种六氟化硫气体回收装置及回收率测量方法 |
CN112520704B (zh) * | 2020-11-26 | 2024-04-26 | 河南省日立信股份有限公司 | 一种六氟化硫气体回收装置及回收率测量方法 |
CN113691893A (zh) * | 2021-08-09 | 2021-11-23 | 武汉理工大学 | 一种压力机能耗远程智能监控平台及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109173584B (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Simultaneous biogas purification and CO2 capture by vacuum swing adsorption using zeolite NaUSY | |
Kim et al. | Study on a numerical model and PSA (pressure swing adsorption) process experiment for CH4/CO2 separation from biogas | |
Tian et al. | Experiment and simulation study of a dual-reflux pressure swing adsorption process for separating N2/O2 | |
Kumar et al. | A versatile process simulator for adsorptive separations | |
TWI484994B (zh) | 具有多容器床的變壓吸附方法及系統 | |
TW200817080A (en) | Performance stability in shallow beds in pressure swing adsorption systems | |
TW200817079A (en) | Performance stability in rapid cycle pressure swing adsorption systems | |
CN107349746A (zh) | 一种循环操作中压力变化联动控制的变压吸附方法 | |
JP2002001113A (ja) | 圧力変動吸着分離用の吸着剤及び吸着筒並びに装置 | |
CN109985487B (zh) | 一种基于Fluent的变压吸附碳捕集的优化方法 | |
EP2958655B1 (en) | Separation process of gaseous compounds from natural gas with low exergy losses | |
Chou et al. | Simulation of a four-bed pressure swing adsorption process for oxygen enrichment | |
CN107456845B (zh) | 一种变压吸附装置及其控制方法 | |
CN109173584A (zh) | 一种基于ldf模型的psa优化策略 | |
Huang et al. | Commercial adsorbents as benchmark materials for separation of carbon dioxide and nitrogen by vacuum swing adsorption process | |
CN110240129B (zh) | 一种变压吸附制氮与干燥一体式节能装置和方法 | |
Choi et al. | Incorporation of a valve equation into the simulation of a pressure swing adsorption process | |
Li et al. | Optimization design research of air flow distribution in vertical radial flow adsorbers | |
CN204111314U (zh) | 小型三塔结构分子筛制氧机 | |
CN112004774A (zh) | 一种采用变压吸附技术生产氧气的方法 | |
Chai et al. | Efficiency of Nitrogen Desorption from LiX Zeolite by Rapid Oxygen Purge in a Pancake Adsorber. | |
CN109529534B (zh) | 一种净化有机硅分馏尾气并回收有机硅单体的方法及装置 | |
Beeyani et al. | Parametric studies and simulation of PSA process for oxygen production from air | |
Tańczyk et al. | Multicomponent pressure swing adsorption. Part II. Experimental verification of the model | |
CA1153316A (en) | Rapid cycle gas separation system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |