胡明剛　陶冶　崔哲　田文德

摘 要 目前我國石油化工領(lǐng)域生產(chǎn)過程存在較為嚴(yán)重的油氣蒸發(fā)損耗問題，不僅導(dǎo)致油品資源的損耗，也嚴(yán)重影響了大氣質(zhì)量，常通過吸附方式進(jìn)行油氣回收。針對實(shí)際工業(yè)中油氣吸附效率較低的問題，致力于尋求提高活性炭吸附性能的方法?；谟?jì)算流體力學(xué)軟件中的組分輸運(yùn)及多孔介質(zhì)模型，對油氣吸附過程進(jìn)行建模及動態(tài)模擬，分析活性炭床層內(nèi)部的濃度場分布情況，并采用控制變量法對比空塔氣速對吸附性能的影響。不同氣速下穿透曲線的比對結(jié)果表明：氣速越大穿透曲線越陡峭，相同時間間隔內(nèi)的吸附負(fù)荷越大穿透時間越短。

關(guān)鍵詞 參數(shù)優(yōu)化 CFD模擬 油氣吸附 氣速 穿透曲線 空塔氣速 操作壓力 操作溫度

中圖分類號 TP274? ?文獻(xiàn)標(biāo)識碼 A? ?文章編號 1000-3932（2023）04-0545-05

在石油化工領(lǐng)域，油品組分中存在大量揮發(fā)性較強(qiáng)的輕烴組分，在它們的開采、提煉、儲運(yùn)、應(yīng)用等過程中，由于受技術(shù)、工藝、設(shè)備條件等的限制，部分液態(tài)輕烴組分不可避免地會蒸發(fā)而逸入大氣，油品損耗嚴(yán)重，導(dǎo)致能源浪費(fèi)和經(jīng)濟(jì)損失，還會嚴(yán)重影響大氣質(zhì)量，危害人體健康，并給企業(yè)帶來諸多安全問題［1，2］。因此，有必要對油氣進(jìn)行回收利用，以提高能源利用率、減少經(jīng)濟(jì)損失、保護(hù)生態(tài)環(huán)境和消除安全隱患。

吸附法是目前國內(nèi)較為成熟的、節(jié)能環(huán)保型的油氣回收技術(shù)［3］，但在實(shí)際應(yīng)用中還存在油氣凈化率不夠高、吸附效率有待進(jìn)一步提高的問題。傳統(tǒng)的實(shí)驗(yàn)研究方法大多受設(shè)備規(guī)模、檢測精度等外部條件的制約，存在操作周期長、成本高、人力物力消耗大的缺點(diǎn)。理論研究方法雖然沒有上述缺點(diǎn)，但必須對計(jì)算對象進(jìn)行簡化，在非線性情況下只有少數(shù)流動才能產(chǎn)生分析結(jié)果。

計(jì)算流體力學(xué)（Computational Fluid Dynamics，CFD）是一種可靠的模擬技術(shù)，可以克服實(shí)驗(yàn)方法和理論分析方法的缺點(diǎn)，用于模擬吸附過程中的流體流動狀況和組分輸運(yùn)現(xiàn)象，進(jìn)而優(yōu)化工藝設(shè)備及操作參數(shù)［4，5］。對于以活性炭為吸附劑的吸附過程來說，空塔氣速、操作壓力、操作溫度等參數(shù)都會影響吸附床層的利用效率，許多學(xué)者對該過程進(jìn)行了深入的研究。梁穎做了吸附罐內(nèi)氣體空間濃度分布的CFD數(shù)值模擬，基于有限元法研究了氣體在空間內(nèi)連續(xù)擴(kuò)散的問題，發(fā)現(xiàn)罐內(nèi)的氣體濃度并不是均勻分布的，在高度方向上濃度存在明顯的拐點(diǎn)［6］，這為確定吸附罐內(nèi)活性炭的安裝位置提供了幫助。武華賓設(shè)計(jì)了一套三級冷凝-吸附集成的VOCs處理工藝，并采用Fluent對設(shè)備內(nèi)的流場分布情況進(jìn)行模擬分析，深入探討了設(shè)備選型優(yōu)化和集成工藝參數(shù)設(shè)計(jì)［7］。文獻(xiàn)［8］采用多孔介質(zhì)模型和用戶自定義方程（UDF）相結(jié)合的方法，模擬苯在活性炭上的吸附過程和多孔介質(zhì)中的流場分布，采用層流單相流方程結(jié)合傳質(zhì)方程，建立三維非定常氣體流動模型，吸附速率由用戶自定義的源項(xiàng)添加到模型中。事實(shí)上，CFD模擬方法已成為研究VOCs凈化裝置內(nèi)流體分布狀況及流動規(guī)律的重要手段［9～11］。

筆者基于CFD軟件中的多孔介質(zhì)模型、組分輸運(yùn)模型和自主編寫的UDF代碼進(jìn)行油氣吸附過程的數(shù)值模擬，多孔介質(zhì)模塊包含一整套物理場接口，可以近似地模擬活性炭中的實(shí)際孔隙。針對單相的多孔介質(zhì)模型不能反映固體吸附劑與氣體吸附質(zhì)之間的質(zhì)量、熱量傳遞問題，通過編寫UDF還原了吸附分離過程中氣固兩相間的傳質(zhì)、傳熱現(xiàn)象，將單相模型改進(jìn)為氣固兩相耦合模型，探索了吸附過程內(nèi)在機(jī)理，優(yōu)化了過程操作參數(shù)。

1 數(shù)值求解方法

由一組偏微分代數(shù)方程給出的模型方程較難求解。當(dāng)氣流軸向擴(kuò)散不明顯時，濃度曲線可能會比較陡峭，這類曲線需要精細(xì)化的網(wǎng)格來進(jìn)行處理，這就需要大量的計(jì)算工作。因此，在優(yōu)化吸附過程時減少計(jì)算工作量就顯得尤為關(guān)鍵，所以應(yīng)該將數(shù)值進(jìn)行離散化處理［12］。筆者采用有限差分法將計(jì)算域進(jìn)行離散化，將一組偏微分代數(shù)方程轉(zhuǎn)換為一組微分代數(shù)方程。在求解過程中，數(shù)值求解方法選擇較為穩(wěn)定的迎風(fēng)差分法［13］，該方法模擬時間較短且模擬結(jié)果合理準(zhǔn)確，可以通過增加節(jié)點(diǎn)數(shù)量來提高精度。

1.1 質(zhì)量平衡方程

1.2 動量平衡方程

1.3 吸附動力學(xué)

2 案例應(yīng)用

將實(shí)際工業(yè)中油氣吸附罐設(shè)備的規(guī)格尺寸等比縮小100倍，在建模軟件中建立幾何模型，然后通過網(wǎng)格劃分軟件對構(gòu)建的幾何體進(jìn)行網(wǎng)格劃分，多孔介質(zhì)區(qū)全部采用結(jié)構(gòu)化的六面體網(wǎng)格并添加邊界層，以更好地模擬真實(shí)流動。采用Mesh軟件提供的正交質(zhì)量和單元畸變度兩種評價指標(biāo)來評價網(wǎng)格質(zhì)量，評價結(jié)果表明所劃分的網(wǎng)格可以滿足計(jì)算精度要求。將劃分完網(wǎng)格的模型導(dǎo)入CFD求解器中進(jìn)行瞬態(tài)求解，求解器設(shè)置見表1。

實(shí)際的油氣組成較為復(fù)雜，為降低模擬的復(fù)雜性，用丁烷（C4H10）、戊烷（C5H12）和空氣的混合物代替油氣，各組分的質(zhì)量分?jǐn)?shù)分別為丁烷0.3、戊烷0.3、空氣0.4。由于C4H10和C5H12間存在競爭吸附現(xiàn)象，且活性碳層對分子量較大的C5H12吸附能力強(qiáng)，優(yōu)先吸附C5H12，所以C4H10比C5H12先穿透，因此將C4H10穿透時間作為衡量吸附性能的指標(biāo)，具體操作方法為監(jiān)控吸附罐出口處C4H10的質(zhì)量分?jǐn)?shù)。

油氣在多孔介質(zhì)內(nèi)吸附傳質(zhì)過程中的丁烷質(zhì)量分?jǐn)?shù)云圖如圖1所示，紅色區(qū)域代表活性炭已處于飽和狀態(tài)，藍(lán)色區(qū)域表示活性炭還未被利用，中間部分代表傳質(zhì)區(qū)。具體數(shù)值如圖1比色柱所示。

吸附罐不同軸向位置的壓力分布情況如圖2所示，可以看出，氣體進(jìn)、出口封頭處壓力大小不變，多孔介質(zhì)區(qū)域壓強(qiáng)均勻下降。適當(dāng)?shù)膲航涤兄跉怏w均勻分布，壓降過大易造成安全隱患。

如果吸附過程處于理想情況下，即不存在傳質(zhì)阻力、吸附速率趨于無窮大、氣固兩相能夠瞬間達(dá)到平衡狀態(tài)，穿透曲線的形狀將呈直角折線。但實(shí)際中由于存在傳質(zhì)阻力、流動速度分布不均等因素的影響，穿透曲線為S型曲線。

筆者對比了空塔氣速分別為0.2、0.3、0.4、0.5 m/s的情況下，對吸附穿透曲線造成的影響，進(jìn)而指導(dǎo)工業(yè)實(shí)際操作參數(shù)的選用，如圖3所示。

認(rèn)為當(dāng)出口濃度C等于0.05倍的入口濃度，即Ct /C0=0.05時，活性炭床層被穿透，達(dá)到“穿透點(diǎn)”，對應(yīng)的流動時間即為穿透時間。出口濃度C等于0.95倍的入口濃度，即Ct /C0=0.95時，活性炭床層飽和，達(dá)到“流干點(diǎn)”。圖3中的模擬結(jié)果表明，當(dāng)其他模擬條件一定時，在相同的時間間隔內(nèi)，氣流空塔氣速的大小對床層吸附的性能影響較大。隨著氣流空塔氣速v增加，穿透曲線斜率變大，吸

附傳質(zhì)速率加快，穿透曲線向左平移，吸附床達(dá)到“穿透點(diǎn)”的時間縮短，達(dá)到“流干點(diǎn)”的時間也減小，床層的吸附負(fù)荷變大。

3 結(jié)束語

當(dāng)氣流空塔氣速從0.2 m/s增加到0.5 m/s，穿透時間從3.6 s縮短到1.0 s。隨著氣速增加，吸附速率加快，穿透曲線斜率也隨之變大，穿透曲線更為陡峭，床層穿透時間和達(dá)到平衡的時間明顯減小，等量吸附劑處理的吸附質(zhì)越多，即吸附負(fù)荷越大。然而，氣速增加后單位時間內(nèi)經(jīng)過床層的氣體量增加，床層的無效層厚度和傳質(zhì)區(qū)長度也延長，床層的利用率降低。綜合考慮氣速對吸附床層的無效層厚度和傳質(zhì)區(qū)長度的影響，是下一步工作需要完善的重要內(nèi)容。

參 考 文 獻(xiàn)

［1］ 王教凱.油氣處理廠VOCs排放特征及綜合管控對策研究［J］.石油石化綠色低碳，2020，5（6）：40-44.

［2］ 王思宇，黃玉虎，胡瑋，等.加油站油氣處理裝置作用及VOCs排放現(xiàn)狀［J］.環(huán)境科學(xué)研究，2021，34（2）：279-285.

［3］ 黃義峰.油氣的危害及回收方法［J］.廣州化工，2018，46（10）：103-104；107.

［4］ 魏新利，李慧，張軍.計(jì)算流體動力學(xué)（CFD）在化工領(lǐng)域的應(yīng)用［J］.化工時刊，2006（2）：63-65.

［5］ 舒長青，王友欣.計(jì)算流體力學(xué)在化學(xué)工程中的應(yīng)用［J］.化工管理，2014（6）：29.

［6］ 梁穎.基于有限元法的罐內(nèi)氣體空間濃度分布模擬［J］.當(dāng)代化工，2020，49（4）：732-737.

［7］ 武華賓.基于Fluent的冷凝-吸附集成VOCs治理技術(shù)仿真模擬及設(shè)計(jì)［D］.廣州：華南理工大學(xué)，2020.

［8］ HASHEMIPOUR M H.Experimental study and artificial neural network simulation of methane adsorption on activated carbon［J］.Korean Journal of Chemical Engineering，2012，29：601-605.

［9］ 卜令兵，張劍鋒，楊云.變壓吸附空分制氧應(yīng)用進(jìn)展［J］.氣體分離，2007，46（3）：12-14.

［10］ 姜曉明，聶宇宏，金云學(xué)，等.基于組分輸運(yùn)模型的SCR煙氣脫硝系統(tǒng)模擬及優(yōu)化［J］.熱力發(fā)電，2016，45（11）：106-113.

［11］ 楊亞超.真空變壓吸附捕集二氧化碳建模及動態(tài)特性研究［D］.上海：華東理工大學(xué)，2016.

［12］ 吳琦，杜強(qiáng)，張金華，等.多孔介質(zhì)輕烴吸附數(shù)值研究［J］.當(dāng)代化工，2013，42（4）：513-515；519.

［13］ 張志娟，袁益讓.非飽和水流問題的迎風(fēng)差分法及其數(shù)值模擬［J］.應(yīng)用數(shù)學(xué)學(xué)報，2009，32（3）：385-399.

［14］ 陳勇.氣固兩相流變壓吸附制氧的CFD模擬［D］.大連：大連理工大學(xué)，2013.

［15］ 解晨.甲烷在活性炭上的吸附平衡及充放氣研究［D］.廈門：集美大學(xué)，2012.

［16］ 彭育志.氫氣純化變壓吸附的穿透曲線及循環(huán)熱效應(yīng)［D］.武漢：武漢理工大學(xué)，2016.

［17］ 趙毅，沈艷梅，白鷺.固定床CO2吸附的穿透曲線模型及傳質(zhì)動力學(xué)［J］.華北電力大學(xué)學(xué)報（自然科學(xué)版），2012，39（4）：107-112.

（收稿日期：2022-12-06，修回日期：2023-01-07）

Key Parameters Optimization in Oil & Gas Adsorption Process

Based on CFD Simulation

HU Ming-gang TAO Ye2， CUI Zhe TIAN Wen-de

（1. Qingdao Nuocheng Chemicals Safety Technology Co.， Ltd.；

2. College of Chemical Engineering， Qingdao University of Science & Technology）

Abstract? ?The oil & gas evaporation loss bothers petrochemical production at home， which leads to the loss of oil resources and seriously influences the air quality and the way of adsorption is employed generally for oil & gas recovery. In this paper， considering the poor efficiency of oil & gas adsorption in industry， theways to improve adsorption performance of activated carbon was proposed. In which， having component transport and porous media models in computational fluid dynamics software based to model and simulate the oil & gas adsorption process was implemented， including analyzing the concentration field distribution within activated carbon bed and having control variable method adopted to compare the effect of superficial gas flow rate on the adsorption performance. The comparison of the penetration curves at different gas velocities shows that， the higher gas velocity can bring about steeper penetration curve and in the same time interval， the higher adsorption load results in shorter penetration time.

Key words? ? parameter optimization， CFD simulation， oil & gas adsorption， gas velocity， penetration curve，? superficial gas flow rate， operating pressure， operating temperature