不同季銨鹽作用下的CO2水合物相平衡

引 言

通過資源化利用的方式來消減CO2氣體排放，是應對全球氣候變化、實現(xiàn)碳中和目標的重要途徑之一[1-2]。CO2氣體水合物是一種由水和CO2在特定環(huán)境下生成的類冰固體，具有獨特的熱力學性質(zhì)和較高的儲冷能力，在海水淡化、沼氣純化、空調(diào)蓄冷、CO2捕集和海底封存等領(lǐng)域都表現(xiàn)出了巨大潛力[3-5]。但自然狀態(tài)下CO2氣體水合物的生成條件較為苛刻（一般高于2.5 MPa，小于10℃），生成速率十分緩慢，嚴重阻礙了水合物技術(shù)的應用和推廣。熱力學促進劑可以通過其與水分子相互作用生成水合物籠的方式來增大水合物的生成推動力，進而改善水合物相平衡條件[6-7]。

目前在國內(nèi)外研究較多的熱力學促進劑中有一類是季銨鹽類添加劑，因其易與水分子結(jié)合，促進水合成籠，從而使得CO2水合物在較為緩和的溫度、壓力范圍內(nèi)快速生成。Sun等[8]測定了四丁基氯化銨（TBAC） + 二氧化碳 + 水混合物的水合物平衡條件，發(fā)現(xiàn)TBAC的存在降低了二氧化碳雙水合物的形成壓力。Mohammadi等[9]探究了TBAC對氣體水合物相平衡條件的影響，結(jié)果表明在TBAC質(zhì)量分數(shù)為5%、15%和22%時能顯著緩和水合物的形成和解離條件，并且從5%增加到22%時，促進作用增強。Wang等[10]研究了CO2和四丁基溴化磷（TBPB）水合物的相平衡行為，在10%、20%、30%、37%（質(zhì)量）的TBPB下生成水合物，可有效降低所需的壓力，但水合物的形成也表現(xiàn)為過冷程度較大和誘導延遲時間較長。Mayoufi等[11]測定了TBPB + CO2 + 水體系的相平衡數(shù)據(jù)，結(jié)果表明，在低濃度(摩爾分數(shù)0.0058)的水中加入TBPB后，281.6 K下CO2水合物的生成壓力降至0.5 MPa，而在相同溫度下，在沒有促進劑的情況下，生成壓力降至3.5 MPa。Wang等[12]通過對四正丁基溴化銨（TBAB）、TBPB和TBAB + TBPB與CO2作用下半水合物的生成動力學研究，表明在一定過冷條件下，TBAB + TBPB溶液的CO2吸收大于TBAB和TBPB溶液，TBAB+TBPB溶液能提高水合物生成動力學和CO2儲存能力。Ye等[13]測定了在壓力0.3~4.4 MPa和溫度279~292 K條件下CO2 + TBAB水合物的平衡數(shù)據(jù)，發(fā)現(xiàn)質(zhì)量分數(shù)在0.05~0.32時，平衡溫度在一定壓力下升高。Fan等[14]采用水合物分離法在281.3 K條件下對67.00%（mol）CH4/CO2混合氣體與TBAB溶液進行CO2捕集，發(fā)現(xiàn)TBAB的加入增大了CO2的捕集效率。Chima-Maceda等[15]采用恒容壓力搜索法測定了1-丙醇/TBAB水溶液中二氧化碳在正己烷存在下的水合物平衡條件，結(jié)果表明，1-丙醇對CO2 + 正己烷水合物具有熱力學抑制作用，而TBAB被證明是一種強有力的熱力學促進劑。

近年來，國內(nèi)外關(guān)于水合物熱力學促進劑的研究[16-18]主要集中在以下兩個方面：一是通過實驗測定不同水合物的生成條件及其所對應的三相平衡點；二是通過對比不同添加劑的促進效果，嘗試尋找一種能夠降低水合物生成能耗的新型、高效、環(huán)保、經(jīng)濟的水合物促進劑。實際上，在實施熱力學強化水合物生成過程中，不同氣體、不同促進劑所生成氣體水合物的相平衡條件都有所不同[19]，因此，必須對不同促進劑作用下的CO2水合物相平衡特征進行深入研究?？紤]到物質(zhì)組成、分子結(jié)構(gòu)及理化性質(zhì)的差異，本文以四丁基氟化銨（TBAF）、四丁基氯化銨（TBAC）、四丁基溴化銨（TBAB）、芐基三乙基氯化銨（TEBAC）為研究對象，探究了其對CO2氣體水合物相平衡條件和穩(wěn)定性的影響。

1 實驗部分

1.1 主要原料

實驗過程中所使用的CO2氣體是由南京特種氣體有限公司提供，純度為高純（體積分數(shù)為99.99%），主要試劑如表1所示。

表1   實驗試劑

Table 1  Experimental reagents

1.2 實驗設(shè)備與流程

鑒于CO2氣體水合物生成過程中相態(tài)轉(zhuǎn)變的復雜性，采用恒容溫度搜索法[20-21]，通過低溫高壓裝置，分別測定了不同季銨鹽促進劑作用下CO2氣體水合物的相平衡數(shù)據(jù)。該實驗裝置主要由CO2氣體鋼瓶、恒溫水浴槽、攪拌設(shè)備、可視高壓平衡釜以及數(shù)據(jù)采集記錄系統(tǒng)五部分組成，如圖1所示。其中，可視高壓相平衡釜的主體材質(zhì)為316 L不銹鋼（設(shè)計壓力為25 MPa，體積為220 ml），內(nèi)置機械攪拌裝置（最高轉(zhuǎn)速為1440 r/min），整體浸沒在恒溫水浴槽中。平衡釜配有一個0.5%擴散硅芯片壓力傳感器（±1 kPa）和兩個Pt-100鉑金電阻溫度傳感器（±0.1 K），主要用于采集釜內(nèi)壓力、液相溫度、氣相溫度等數(shù)據(jù)，并通過傳感器將模擬信號轉(zhuǎn)換為數(shù)字信號，傳輸給計算機實時記錄保存實驗數(shù)據(jù)。

圖1

圖1   實驗裝置圖

1—CO2氣體瓶;2—恒溫水浴槽;3—攪拌裝置;4—可視高壓平衡釜;5—數(shù)據(jù)采集記錄系統(tǒng)

Fig.1   Diagram of experimental device

水合物生成過程：設(shè)定初始溫度，開啟實驗數(shù)據(jù)采集及記錄系統(tǒng)，實時記錄壓力、溫度的變化；將季銨鹽水溶液置于平衡釜中，緩慢充入CO2氣體至實驗壓力，打開攪拌，開始降溫，待釜內(nèi)壓力迅速降低時，說明此時水合物大量生成，停止降溫。

水合物分解過程：水合物生成約4 h后，逐步升溫并實時觀察水合物形態(tài)和數(shù)量變化，直至水合物完全分解，體系恢復為氣液兩相。此時，水合物生成/分解圖中的交點溫度即為其所對應的壓力下的水合物相平衡溫度。重復上述步驟三次，以平均值來確定水合物相平衡點。重新設(shè)定壓力，重復上述過程，即可得到其他相平衡點數(shù)據(jù)。

1.3 水合物相平衡點的確定

考慮到水合物生成過程中需要較高的過冷度和較長的誘導時間[22-23]，為準確測定水合物相平衡數(shù)據(jù)，需要通過水合物完全分解點來確定氣-液-固三相相平衡點[24-25]。本文采用恒容溫度搜索法，測定了水合物生成和分解過程中的溫度壓力（P-T）變化關(guān)系，如圖2所示?？梢钥闯?，水合物生成初始階段，隨著溫度的降低，平衡釜內(nèi)壓力連續(xù)降低，壓力隨溫度變化的曲線無明顯拐點（曲線AB）。隨著水合物的不斷生成，平衡釜內(nèi)壓力持續(xù)下降，出現(xiàn)拐點B后，系統(tǒng)壓力迅速降低，說明此時水合物大量生成（曲線BC）。在分解過程中，水合物需要吸收大量熱量[26]，因此初始升溫時，釜內(nèi)壓力變化不大（曲線CD）。但隨著溫度的不斷升高，壓力不斷增大（曲線DE），直至出現(xiàn)P-T相圖的交點E，即為水合物的氣-液-固三相相平衡點。透過平衡釜可視窗口，可以清晰觀察到水合物在生成與分解過程中的狀態(tài)和數(shù)量變化，如圖3所示。其中，圖3（a）實驗初始階段，氣液界面明顯，未見有固體生成；圖3（b）隨著溫度的進一步降低，類冰狀水合物開始生成，系統(tǒng)處于固液共存狀態(tài)；圖3（c）保持低溫狀態(tài)不變，水合物不斷生成，基本觀察不到液相的存在；圖3（d）隨著溫度的不斷升高，固體水合物隨之分解；圖3（e）水合物完全分解，體系恢復為氣液兩相。

圖2

圖2   CO2氣體水合物生成及分解過程P-T曲線

A→B降溫過程階段;B→C恒溫水合物生成階段;C→D水合物維持階段;D→E水合物分解階段

Fig.2   P-T curve of formation and decomposition process of CO2 hydrate

圖3

圖3   CO2氣體水合物形成過程

Fig.3   CO2 hydrate formation process

2 結(jié)果與討論

2.1 數(shù)據(jù)計算及可靠性驗證

2.1.1 實驗數(shù)據(jù)計算

（1）平均驅(qū)動力ΔTave$?T_{ave}$的計算

水合物生成的平均驅(qū)動力ΔTave$?T_{ave}$是衡量添加劑對CO2氣體水合物促進效果的重要參數(shù)[27]，如式（1）所示：

ΔTave=1n∑i=1nΔTi$?T_{ave}=\frac{1}{n}\sum _{i=1}^n?T_i$(1)

其中，n是數(shù)據(jù)點的數(shù)量，ΔTi$\Delta T_i$是第i$i$個點在相同壓力下由促進劑所引起相平衡溫度的差值。

（2）相變潛熱的計算

水合物相變潛熱是氣體水合物生成或分解過程中放出或吸收的單位熱量，通常采用Clausius-Clapeyron方程來計算[28]，如式（2）所示：

d(lnP)d(1/T)=?ΔHdZR$\frac{d\left(lnP\right)}{d(1/T)}=-\frac{?H_d}{ZR}$(2)

其中，P、T分別為氣體水合物生成時的相平衡壓力和溫度；R為氣體常數(shù)（8.314 J·mol-1·K-1）；Z為氣體壓縮因子，通過PR方程計算；ΔHd$\Delta H_d$為水合物的生成相變潛熱。

該方程假設(shè)0℃以上時，氣體水合物的生成相變潛熱ΔHd$?H_d$僅取決于氣體分子所占據(jù)的水合物空腔結(jié)構(gòu)，與氣體種類、濃度無關(guān)[29]。根據(jù)水合物相平衡數(shù)據(jù)，以lnP對1/T作圖，并進行線性擬合，其所對應的斜率大小不僅能夠用于推測水合物結(jié)構(gòu)，還可用于計算水合物生成相變潛熱，相變潛熱計算值越大，說明該水合物越穩(wěn)定。同時，曲線的擬合度R2可用作判斷實驗數(shù)據(jù)是否滿足一致性檢驗[30-31]。

2.1.2 實驗可靠性驗證

為驗證裝置的可靠性，實驗測定了CO2氣體水合物的相平衡數(shù)據(jù)，并將其與文獻數(shù)據(jù)[27-28, 32-33]進行對比，如圖4所示?？梢钥闯?，實驗值與文獻值吻合良好，相對誤差不超過3%，說明恒容溫度搜索法可以準確測定CO2氣體水合物的相平衡數(shù)據(jù)。

圖4

圖4   CO2水合物相圖

Fig.4   Phase diagram of CO2 gas hydrate

2.2 季銨鹽對CO2水合物相平衡的影響

季銨鹽熱力學促進劑（TBAF、TBAC、TBAB、TEBAC）作用下的CO2水合物相平衡數(shù)據(jù)列于表2、圖5、圖6。其中，Ted、Ped分別為氣-液-固三相平衡時的溫度和壓力，ΔTave$\Delta T_{ave}$表示以溫度為驅(qū)動力的平均驅(qū)動力。可以看出，TBAF、TBAC、TBAB均可以明顯改變CO2水合物相平衡條件，且不同促進劑的影響效果也不盡相同。這是因為季銨鹽分子結(jié)構(gòu)中的陽離子易于填充到水合物籠中，進而占據(jù)CO2水合物空腔結(jié)構(gòu)中較大的晶穴空腔，使得水合物結(jié)構(gòu)變得更加穩(wěn)定有序；同時，季銨鹽分子中的陰離子通過取代CO2水合物結(jié)構(gòu)中水分子的方式，形成半籠型結(jié)構(gòu)[34]，增大固體水合物對CO2的吸附能力，進而使CO2氣體水合物生成條件更加溫和，相平衡條件向更高的溫度和更低的壓力區(qū)域偏移。

表2   促進劑作用下CO2水合物相平衡條件

Table 2  Phase equilibrium condition of CO2 hydrate under additive action

圖5

圖5   季銨鹽作用下CO2水合物相平衡條件

Fig.5   Phase equilibrium condition for CO2 hydrate under quaternary ammonium salt action

圖6

圖6   不同質(zhì)量分數(shù)促進劑作用下CO2水合物相平衡條件

Fig.6   Phase equilibrium condition for CO2 hydrate under different accelerant mass fraction

圖5（a）~（d）分別給出了不同質(zhì)量分數(shù)條件下TBAF、TBAC、TBAB和TEBAC對CO2水合物的促進效果?？梢钥闯觯琓BAF、TBAC、TBAB和TEBAC四種季銨鹽促進劑都可以擴大水合物的穩(wěn)定區(qū)。相同條件下，隨著季銨鹽配比濃度的增加，改善CO2水合物相平衡條件的效果呈現(xiàn)出先增后減的趨勢，因此熱力學促進劑添加量存在最佳值。當季銨鹽水溶液質(zhì)量分數(shù)為5%~10%時，隨著質(zhì)量分數(shù)的增加，CO2水合物相平衡溫度顯著升高，其中質(zhì)量分數(shù)為10% TBAF溶液對CO2水合物促進效果最好，平均驅(qū)動力?Tave為12.95 K。當季銨鹽水溶液質(zhì)量分數(shù)升高至10%~15%后，隨著質(zhì)量分數(shù)的進一步增加，尤其是使用量大于其與水完全生成水合物的理論配比時，過量的四丁基鹵化銨將會產(chǎn)生抑制效果，導致促進效果減弱。

圖6（a）~（e）分別給出了相同質(zhì)量分數(shù)下TBAF、TBAC、TBAB和TEBAC對CO2水合物的促進效果?？梢钥闯?，在相同條件下，CO2氣體水合物生成的相平衡溫度由高到低分別為：TBAF > TBAB > TBAC > TEBAC，其中TBAF、TBAC、TBAB的促進效果明顯高于TEBAC。從微觀角度來看，TEBAC分子結(jié)構(gòu)中的苯環(huán)占據(jù)了較大空間，導致水分子無法在其周圍生成活潑水分子簇，進而削弱了CO2氣體分子和游離水的遷移程度，因此TEBAC促進效果最差。TBAF分子結(jié)構(gòu)中陰離子F-與CO2分子之間范德華作用力更強，使得CO2更容易進入到不穩(wěn)定的水合物分子簇中，進而生成穩(wěn)定的氣體水合物，因此TBAF促進效果最為明顯。但過量TBAF能夠與水生成固體季銨鹽水合物，增大CO2進入液相的阻力，進而抑制水合物生成，導致其促進效果減弱。綜合考慮促進劑的質(zhì)量分數(shù)、水合物特性以及促進效果，壓力范圍在0.35~4.50 MPa內(nèi)，選擇質(zhì)量分數(shù)范圍在10%~12%內(nèi)的TBAB作為CO2氣體水合物的促進劑。

2.3 季銨鹽作用下CO2水合物相變潛熱的計算

TBAF、TBAC、TBAB和TEBAC作用下的CO2水合物相變潛熱計算結(jié)果見表3?？梢钥闯?，相同濃度下，CO2水合物相變潛熱由大到小分別為：TBAB > TBAF > TBAC > TEBAC，這與由相平衡數(shù)據(jù)得到季銨鹽促進效果的優(yōu)越性順序基本一致。其中，TBAB和TBAF作用下CO2水合物相變潛熱相接近且明顯高于其他季銨鹽，說明其促進效果最好，所對應的水合物生成條件也最為溫和。同時，TBAF、TBAC和TBAB三種促進劑的斜率計算值在 -28000~-20000之間，進一步說明四丁基季銨鹽與水結(jié)合能夠生成半籠型氣體水合物，這與Sun等[8]研究結(jié)果相符。TEBAC的斜率計算值大于-10000，說明TEBAC的作用下生成的CO2水合物結(jié)構(gòu)更接近于SⅠ型，與四丁基季銨鹽的半籠型水合物存在較大差異。擬合度R2值均大于0.99，說明所測相平衡數(shù)據(jù)具有很好的線性擬合關(guān)系，滿足熱力學一致性檢驗原則。此外，TBAC計算偏差較大，為27.73%，這可能是因為TBAC在低濃度時所生成半籠型水合物易受Cl-的影響，使TBAC對CO2氣體分子的吸附能力減弱，進而導致低濃度時TBAC作用下的水合物相變潛熱較小。

表3   季銨鹽作用下CO2水合物的相變潛熱

Table 3  Dissociation enthalpies of CO2 salt hydrate under quaternary ammonium action

2.4 氣-液-固三相平衡數(shù)據(jù)的計算

氣體水合物的相平衡條件是衡量其能否穩(wěn)定存在的重要標準，本文利用Chen-Guo模型[35]，結(jié)合PR狀態(tài)方程[36]和改進Joshi經(jīng)驗活度模型[37]，分別對TBAF、TBAB、TBAC和TEBAC作用下CO2氣體水合物熱力學相平衡數(shù)據(jù)進行了計算。

2.4.1 Chen-Guo熱力學模型

氣體水合物相平衡條件預測模型多是基于范德瓦爾-普朗特[38]統(tǒng)計熱力學基本模型發(fā)展而來的，但計算結(jié)果始終存在較大誤差。不同于范德單吸附機理的水合物形成/分解機理，Chen-Guo模型認為水合物成核過程中同時進行著以下兩種動力學過程：準化學反應動力學過程和吸附動力學過程。因此，在對各類氣體水合物體系的相平衡條件進行預測時，Chen-Guo模型的計算結(jié)果與實驗數(shù)值符合得更好。方程形式如下：

f=f0(1?θ)α$f=f^0{\left(1-\theta \right)}^{\alpha }$(3)θ=Cf1+Cf$\theta =\frac{Cf}{1+Cf}$(4)C=Xexp(YT?Z)$C=Xexp\left(\frac{Y}{T-Z}\right)$(5)f0=f0(T)f0(P)f0(aw)$f^0=f^0\left(T\right)f^0\left(P\right)f^0\left(a_w\right)$(6)f0(T)=Aexp(BT?C)$f^0\left(T\right)=Aexp\left(\frac{B}{T-C}\right)$(7)f0(P)=exp(βPT)$f^0\left(P\right)=exp\left(\frac{\beta P}{T}\right)$(8)f0(aw)=a?1/λ2w$f^0\left(a_w\right)=a_w^{-1/{\lambda }_2}$(9)

式中，α$\alpha$值僅與水合物結(jié)構(gòu)相關(guān)，對于SⅠ型結(jié)構(gòu)α$\alpha$為1/3，對于SⅡ型結(jié)構(gòu)α$\alpha$為2；CO2氣體計算常數(shù)X、Y、Z分別為1.6464×1011 Pa-1、2799.66、15.9；A、B、C數(shù)值分別為963.72/1015 Pa、-6444.5、36.67；β=ΔVλ2R$\beta =\frac{\Delta V}{{\lambda }_{2}R}$，可作為常數(shù)處理，對于SⅠ型水合物結(jié)構(gòu)β=0.4242 K/105 Pa$\beta =0.4242K/{10}^{5}Pa$，對于SⅡ型水合物結(jié)構(gòu)β=1.0224 K/105 Pa$\beta =1.0224K/{10}^{5}Pa$；對于純水體系液相活度aw$a_w$為1，對于促進劑體系，需要通過活度模型計算。

2.4.2 PR狀態(tài)方程

氣相逸度的計算結(jié)果會直接影響氣體水合物相平衡數(shù)據(jù)的預測精度?，F(xiàn)有研究表明[39]，PR狀態(tài)方程更適合用于預測CO2氣體水合物相平衡條件。方程形式如下：

P=RTV?b?a(T)V2+2bV?b2$P=\frac{RT}{V-b}-\frac{a\left(T\right)}{V^2+2bV-b^2}$(10)a=φaR2T2.5cPc$a={\phi }_a\frac{R^2{T}_c^{2.5}}{P_c}$(11)b=φbRTcPc$b={\phi }_b\frac{R{T}_c}{P_c}$(12)

2.4.3 改進Joshi經(jīng)驗活度模型

Greenspan等[40]認為液相活度與促進劑的濃度、系統(tǒng)的溫度、壓力相關(guān)。促進劑濃度和系統(tǒng)溫度對水合物活度影響比較大，而壓力變化對活度的影響比較小。因此，本文采用Joshi等[37]提出的液相活度aw$a_w$經(jīng)驗關(guān)聯(lián)式，并參照Zhao等[41]的改進方法，忽略壓力對液相活度aw$a_w$的影響，得到新的活度經(jīng)驗方程，如式（13）所示：

aw=(1?xm)exp(K1?K2T)$a_w=\left(1-x_m\right)exp\left(K_1-\frac{K_2}{T}\right)$(13)

式中，xm$x_m$是促進劑的摩爾分數(shù)；K1$K_1$和K2$K_2$為擬合參數(shù)，結(jié)果如表4所示。

表4   擬合參數(shù)K1和K2

Table 4  Values of fitting parameters K1 and K2

圖7（a）~（d）給出了不同溫度條件下，TBAF、TBAB、TBAC和TEBAC作用下CO2氣體水合物的相平衡壓力?？梢钥闯觯嬎阒蹬c實驗值的一致性較高，相平衡壓力計算值和實驗值之間最大平均相對誤差為7.50%，說明Chen-Guo模型結(jié)合改進Joshi經(jīng)驗活度模型能夠很好地預測季銨鹽熱力學促進劑作用下的CO2水合物相平衡數(shù)據(jù)。

圖7

圖7   季銨鹽作用下CO2水合物相平衡計算結(jié)果

Fig.7   Phase equilibrium calculation results for CO2 hydrate under quaternary ammonium salt action

3 結(jié) 論

（1）設(shè)計并搭建一套高壓低溫水合物相平衡的實驗裝置，采用恒容溫度搜索法，探討了在壓力0.35~4.50 MPa和溫度272.75~294.35 K范圍內(nèi)，四種季銨鹽促進劑對CO2氣體水合物相平衡特征的影響。結(jié)果表明，隨著促進劑濃度增加，CO2水合物的相平衡條件向高溫、低壓區(qū)域轉(zhuǎn)移。相同條件下，促進劑作用下CO2水合物的相平衡溫度由高到低分別是：TBAF > TBAB > TBAC > TEBAC。

（2）基于Clausius-Clapeyron方程，計算了不同體系的相變潛熱，并探討了其對水合物穩(wěn)定性的影響?？梢钥闯?，相同條件下，對于季銨鹽促進劑，穩(wěn)定性由高到低：TBAB > TBAF > TBAC >TEBAC。其中TBAB、TBAF的計算相變潛熱相接近，說明其穩(wěn)定性相似。

（3）利用Chen-Guo模型，結(jié)合 PR 狀態(tài)方程和改進Joshi 經(jīng)驗活度模型，分別計算了TBAF、TBAB、TBAC和TEBAC作用下CO2水合物熱力學相平衡數(shù)據(jù)，計算結(jié)果與實驗結(jié)果吻合良好，最大平均相對誤差為7.50%。