不同季铵盐作用下的CO2水合物相平衡

1、PAGE PAGE - 18 -不同季铵盐作用下的CO2水合物相平衡引言通过资源化利用的方式来消减CO2气体排放，是应对全球气候变化、实现碳中和目标的重要途径之一1-2。CO2气体水合物是一种由水和CO2在特定环境下生成的类冰固体，具有独特的热力学性质和较高的储冷能力，在海水淡化、沼气纯化、空调蓄冷、CO2捕集和海底封存等领域都表现出了巨大潜力3-5。但自然状态下CO2气体水合物的生成条件较为苛刻（一般高于2.5MPa，小于10），生成速率十分缓慢，严重阻碍了水合物技术的应用和推广。热力学促进剂可以通过其与水分子相互作用生成水合物笼的方式来增大水合物的生成推动力，进而改善水合物相平衡条件6-7

2、。目前在国内外研究较多的热力学促进剂中有一类是季铵盐类添加剂，因其易与水分子结合，促进水合成笼，从而使得CO2水合物在较为缓和的温度、压力范围内快速生成。Sun等8测定了四丁基氯化铵（TBAC）+二氧化碳+水混合物的水合物平衡条件，发现TBAC的存在降低了二氧化碳双水合物的形成压力。Mohammadi等9探究了TBAC对气体水合物相平衡条件的影响，结果表明在TBAC质量分数为5%、15%和22%时能显著缓和水合物的形成和解离条件，并且从5%增加到22%时，促进作用增强。Wang等10研究了CO2和四丁基溴化磷（TBPB）水合物的相平衡行为，在10%、20%、30%、37%（质量）的TBPB下生

3、成水合物，可有效降低所需的压力，但水合物的形成也表现为过冷程度较大和诱导延迟时间较长。Mayoufi等11测定了TBPB+CO2+水体系的相平衡数据，结果表明，在低浓度(摩尔分数0.0058)的水中加入TBPB后，281.6K下CO2水合物的生成压力降至0.5MPa，而在相同温度下，在没有促进剂的情况下，生成压力降至3.5MPa。Wang等12通过对四正丁基溴化铵（TBAB）、TBPB和TBAB+TBPB与CO2作用下半水合物的生成动力学研究，表明在一定过冷条件下，TBAB+TBPB溶液的CO2吸收大于TBAB和TBPB溶液，TBAB+TBPB溶液能提高水合物生成动力学和CO2储存能力。Ye等

4、13测定了在压力0.34.4MPa和温度279292K条件下CO2+TBAB水合物的平衡数据，发现质量分数在0.050.32时，平衡温度在一定压力下升高。Fan等14采用水合物分离法在281.3K条件下对67.00%（mol）CH4/CO2混合气体与TBAB溶液进行CO2捕集，发现TBAB的加入增大了CO2的捕集效率。Chima-Maceda等15采用恒容压力搜索法测定了1-丙醇/TBAB水溶液中二氧化碳在正己烷存在下的水合物平衡条件，结果表明，1-丙醇对CO2+正己烷水合物具有热力学抑制作用，而TBAB被证明是一种强有力的热力学促进剂。近年来，国内外关于水合物热力学促进剂的研究16-18主要

5、集中在以下两个方面：一是通过实验测定不同水合物的生成条件及其所对应的三相平衡点；二是通过对比不同添加剂的促进效果，尝试寻找一种能够降低水合物生成能耗的新型、高效、环保、经济的水合物促进剂。实际上，在实施热力学强化水合物生成过程中，不同气体、不同促进剂所生成气体水合物的相平衡条件都有所不同19，因此，必须对不同促进剂作用下的CO2水合物相平衡特征进行深入研究。考虑到物质组成、分子结构及理化性质的差异，本文以四丁基氟化铵（TBAF）、四丁基氯化铵（TBAC）、四丁基溴化铵（TBAB）、苄基三乙基氯化铵（TEBAC）为研究对象，探究了其对CO2气体水合物相平衡条件和稳定性的影响。1实验部分1.1主要

6、原料实验过程中所使用的CO2气体是由南京特种气体有限公司提供，纯度为高纯（体积分数为99.99%），主要试剂如表1所示。表1实验试剂Table1Experimentalreagents试剂分子式纯度供应商四丁基氟化铵（TBAF）CHFN0.98上海阿拉丁试剂晶纯实业有限公司四丁基氯化铵（TBAC）CHClN0.99上海阿拉丁试剂晶纯实业有限公司四丁基溴化铵（TBAB）CHBrN0.99上海阿拉丁试剂晶纯实业有限公司苄基三乙基氯化铵（TEBAC）C13H22ClN0.99上海阿拉丁试剂晶纯实业有限公司去离子水H2O实验室自制1.2实验设备与流程鉴于CO2气体水合物生成过程中相态转变的复杂性，采用

7、恒容温度搜索法20-21，通过低温高压装置，分别测定了不同季铵盐促进剂作用下CO2气体水合物的相平衡数据。该实验装置主要由CO2气体钢瓶、恒温水浴槽、搅拌设备、可视高压平衡釜以及数据采集记录系统五部分组成，如图1所示。其中，可视高压相平衡釜的主体材质为316L不锈钢（设计压力为25MPa，体积为220ml），内置机械搅拌装置（最高转速为1440r/min），整体浸没在恒温水浴槽中。平衡釜配有一个0.5%扩散硅芯片压力传感器（1kPa）和两个Pt-100铂金电阻温度传感器（0.1K），主要用于采集釜内压力、液相温度、气相温度等数据，并通过传感器将模拟信号转换为数字信号，传输给计算机实时记录保存实

8、验数据。图1图1实验装置图1CO2气体瓶;2恒温水浴槽;3搅拌装置;4可视高压平衡釜;5数据采集记录系统Fig.1Diagramofexperimentaldevice水合物生成过程：设定初始温度，开启实验数据采集及记录系统，实时记录压力、温度的变化；将季铵盐水溶液置于平衡釜中，缓慢充入CO2气体至实验压力，打开搅拌，开始降温，待釜内压力迅速降低时，说明此时水合物大量生成，停止降温。水合物分解过程：水合物生成约4h后，逐步升温并实时观察水合物形态和数量变化，直至水合物完全分解，体系恢复为气液两相。此时，水合物生成/分解图中的交点温度即为其所对应的压力下的水合物相平衡温度。重复上述步骤三次，以平

9、均值来确定水合物相平衡点。重新设定压力，重复上述过程，即可得到其他相平衡点数据。1.3水合物相平衡点的确定考虑到水合物生成过程中需要较高的过冷度和较长的诱导时间22-23，为准确测定水合物相平衡数据，需要通过水合物完全分解点来确定气-液-固三相相平衡点24-25。本文采用恒容温度搜索法，测定了水合物生成和分解过程中的温度压力（P-T）变化关系，如图2所示。可以看出，水合物生成初始阶段，随着温度的降低，平衡釜内压力连续降低，压力随温度变化的曲线无明显拐点（曲线AB）。随着水合物的不断生成，平衡釜内压力持续下降，出现拐点B后，系统压力迅速降低，说明此时水合物大量生成（曲线BC）。在分解过程中，水合

10、物需要吸收大量热量26，因此初始升温时，釜内压力变化不大（曲线CD）。但随着温度的不断升高，压力不断增大（曲线DE），直至出现P-T相图的交点E，即为水合物的气-液-固三相相平衡点。透过平衡釜可视窗口，可以清晰观察到水合物在生成与分解过程中的状态和数量变化，如图3所示。其中，图3（a）实验初始阶段，气液界面明显，未见有固体生成；图3（b）随着温度的进一步降低，类冰状水合物开始生成，系统处于固液共存状态；图3（c）保持低温状态不变，水合物不断生成，基本观察不到液相的存在；图3（d）随着温度的不断升高，固体水合物随之分解；图3（e）水合物完全分解，体系恢复为气液两相。图2图2CO2气体水合物生成及

11、分解过程P-T曲线AB降温过程阶段;BC恒温水合物生成阶段;CD水合物维持阶段;DE水合物分解阶段Fig.2P-TcurveofformationanddecompositionprocessofCO2hydrate图3图3CO2气体水合物形成过程Fig.3CO2hydrateformationprocess2结果与讨论2.1数据计算及可靠性验证2.1.1实验数据计算（1）平均驱动力TaveTave的计算水合物生成的平均驱动力TaveTave是衡量添加剂对CO2气体水合物促进效果的重要参数27，如式（1）所示：Tave=1ni=1nTiTave=1ni=1nTi(1)其中，n是数据点的数量，T

12、iTi是第ii个点在相同压力下由促进剂所引起相平衡温度的差值。（2）相变潜热的计算水合物相变潜热是气体水合物生成或分解过程中放出或吸收的单位热量，通常采用Clausius-Clapeyron方程来计算28，如式（2）所示：d(lnP)d(1/T)=HdZRd(lnP)d(1/T)=-HdZR(2)其中，P、T分别为气体水合物生成时的相平衡压力和温度；R为气体常数（8.314Jmol-1K-1）；Z为气体压缩因子，通过PR方程计算；HdHd为水合物的生成相变潜热。该方程假设0以上时，气体水合物的生成相变潜热HdHd仅取决于气体分子所占据的水合物空腔结构，与气体种类、浓度无关29。根据水合物相平衡

13、数据，以lnP对1/T作图，并进行线性拟合，其所对应的斜率大小不仅能够用于推测水合物结构，还可用于计算水合物生成相变潜热，相变潜热计算值越大，说明该水合物越稳定。同时，曲线的拟合度R2可用作判断实验数据是否满足一致性检验30-31。2.1.2实验可靠性验证为验证装置的可靠性，实验测定了CO2气体水合物的相平衡数据，并将其与文献数据27-28,32-33进行对比，如图4所示。可以看出，实验值与文献值吻合良好，相对误差不超过3%，说明恒容温度搜索法可以准确测定CO2气体水合物的相平衡数据。图4图4CO2水合物相图Fig.4PhasediagramofCO2gashydrate2.2季铵盐对CO2水

14、合物相平衡的影响季铵盐热力学促进剂（TBAF、TBAC、TBAB、TEBAC）作用下的CO2水合物相平衡数据列于表2、图5、图6。其中，Ted、Ped分别为气-液-固三相平衡时的温度和压力，TaveTave表示以温度为驱动力的平均驱动力。可以看出，TBAF、TBAC、TBAB均可以明显改变CO2水合物相平衡条件，且不同促进剂的影响效果也不尽相同。这是因为季铵盐分子结构中的阳离子易于填充到水合物笼中，进而占据CO2水合物空腔结构中较大的晶穴空腔，使得水合物结构变得更加稳定有序；同时，季铵盐分子中的阴离子通过取代CO2水合物结构中水分子的方式，形成半笼型结构34，增大固体水合物对CO2的吸附能力，

15、进而使CO2气体水合物生成条件更加温和，相平衡条件向更高的温度和更低的压力区域偏移。表2促进剂作用下CO2水合物相平衡条件Table2PhaseequilibriumconditionofCO2hydrateunderadditiveaction促进剂5%（质量）7%（质量）10%（质量）12%（质量）15%（质量）Ted/KPed/MPaTave/KTed/KPed/MPaTave/KTed/KPed/MPaTave/KTed/KPed/MPaTave/KTed/KPed/MPaTave/KTBAF284.250.4411.00285.350.6511.65286.350.8612.9528

16、6.850.813.40287.350.7413.65287.551.34287.951.38288.251.41289.051.38289.751.35288.751.95289.051.92289.351.9290.351.98291.352.06289.452.46289.852.46290.252.46291.252.56292.252.63289.852.76290.552.92291.253.03292.853.35294.353.68290.753.45291.453.51292.153.57293.353.67TBAC273.350.782.20277.850.924.3528

17、1.751.076.95281.350.717.85280.950.358.15276.651.29280.051.38282.851.46282.951.06282.950.65278.651.75281.451.8283.651.86283.851.43284.051.01280.552.12283.052.29284.852.46285.152.12285.451.74281.752.48284.152.87285.653.25286.152.88286.652.51282.752.91285.053.41286.753.91287.153.58287.453.25283.653.362

18、86.153.97288.054.28288.654.19TBAB279.450.426.35281.050.577.45282.050.369.60283.550.439.85283.050.4910.30281.450.71282.450.81284.850.84284.650.66284.450.86282.751.31284.451.31285.351.03285.550.98285.751.16284.051.73285.651.62286.451.52286.551.42286.651.52284.451.98286.552.36287.052.03287.452.07287.75

19、2.12285.452.45287.152.67287.652.61288.152.62288.352.53286.052.91287.553.08288.353.25288.753.23289.053.12286.753.57287.953.54288.853.87289.153.7289.453.53287.054.07288.253.82289.454.32289.554.22289.854.12287.254.23288.754.36TEBAC272.751.130.25273.851.230.55274.951.330.75274.251.220.95273.451.071.1527

20、6.251.75276.451.74276.651.73276.151.65275.651.57278.352.24278.452.14278.452.04278.251.98278.051.92279.452.75279.952.68280.352.62280.252.48279.652.34280.653.09281.053.02281.352.95281.352.91281.452.86281.753.65281.753.46281.753.26282.253.3282.753.34282.854.03283.253.93283.553.83283.753.79283.853.74285

21、.354.26图5图5季铵盐作用下CO2水合物相平衡条件Fig.5PhaseequilibriumconditionforCO2hydrateunderquaternaryammoniumsaltaction图6图6不同质量分数促进剂作用下CO2水合物相平衡条件Fig.6PhaseequilibriumconditionforCO2hydrateunderdifferentaccelerantmassfraction图5（a）（d）分别给出了不同质量分数条件下TBAF、TBAC、TBAB和TEBAC对CO2水合物的促进效果。可以看出，TBAF、TBAC、TBAB和TEBAC四种季铵盐促进剂都可

22、以扩大水合物的稳定区。相同条件下，随着季铵盐配比浓度的增加，改善CO2水合物相平衡条件的效果呈现出先增后减的趋势，因此热力学促进剂添加量存在最佳值。当季铵盐水溶液质量分数为5%10%时，随着质量分数的增加，CO2水合物相平衡温度显著升高，其中质量分数为10%TBAF溶液对CO2水合物促进效果最好，平均驱动力Tave为12.95K。当季铵盐水溶液质量分数升高至10%15%后，随着质量分数的进一步增加，尤其是使用量大于其与水完全生成水合物的理论配比时，过量的四丁基卤化铵将会产生抑制效果，导致促进效果减弱。图6（a）（e）分别给出了相同质量分数下TBAF、TBAC、TBAB和TEBAC对CO2水合物

23、的促进效果。可以看出，在相同条件下，CO2气体水合物生成的相平衡温度由高到低分别为：TBAFTBABTBACTEBAC，其中TBAF、TBAC、TBAB的促进效果明显高于TEBAC。从微观角度来看，TEBAC分子结构中的苯环占据了较大空间，导致水分子无法在其周围生成活泼水分子簇，进而削弱了CO2气体分子和游离水的迁移程度，因此TEBAC促进效果最差。TBAF分子结构中阴离子F-与CO2分子之间范德华作用力更强，使得CO2更容易进入到不稳定的水合物分子簇中，进而生成稳定的气体水合物，因此TBAF促进效果最为明显。但过量TBAF能够与水生成固体季铵盐水合物，增大CO2进入液相的阻力，进而抑制水合物

24、生成，导致其促进效果减弱。综合考虑促进剂的质量分数、水合物特性以及促进效果，压力范围在0.354.50MPa内，选择质量分数范围在10%12%内的TBAB作为CO2气体水合物的促进剂。2.3季铵盐作用下CO2水合物相变潜热的计算TBAF、TBAC、TBAB和TEBAC作用下的CO2水合物相变潜热计算结果见表3。可以看出，相同浓度下，CO2水合物相变潜热由大到小分别为：TBABTBAFTBACTEBAC，这与由相平衡数据得到季铵盐促进效果的优越性顺序基本一致。其中，TBAB和TBAF作用下CO2水合物相变潜热相接近且明显高于其他季铵盐，说明其促进效果最好，所对应的水合物生成条件也最为温和。同时，

25、TBAF、TBAC和TBAB三种促进剂的斜率计算值在-28000-20000之间，进一步说明四丁基季铵盐与水结合能够生成半笼型气体水合物，这与Sun等8研究结果相符。TEBAC的斜率计算值大于-10000，说明TEBAC的作用下生成的CO2水合物结构更接近于S型，与四丁基季铵盐的半笼型水合物存在较大差异。拟合度R2值均大于0.99，说明所测相平衡数据具有很好的线性拟合关系，满足热力学一致性检验原则。此外，TBAC计算偏差较大，为27.73%，这可能是因为TBAC在低浓度时所生成半笼型水合物易受Cl-的影响，使TBAC对CO2气体分子的吸附能力减弱，进而导致低浓度时TBAC作用下的水合物相变潜热

26、较小。表3季铵盐作用下CO2水合物的相变潜热Table3DissociationenthalpiesofCO2salthydrateunderquaternaryammoniumaction促进剂质量分数/%AslopeHd/(kJmol-1)R2TBAF5-23641163.510.99377-23457167.520.997510-28785205.490.994412-23947171.120.992715-24775175.540.9972TBAC5-1092574.750.99307-20674126.890.994810-21348147.980.994612-22486156.4

27、00.995915-27154194.800.9951TBAB5-26843198.500.99857-25296181.370.992310-24526178.820.996112-23530174.140.998215-23850172.880.9910TEBAC5-8882.859.050.99427-9704.564.490.997310-9588.263.060.995212-9162.161.960.995615-9904.565.370.99732.4气-液-固三相平衡数据的计算气体水合物的相平衡条件是衡量其能否稳定存在的重要标准，本文利用Chen-Guo模型35，结合PR状态方

28、程36和改进Joshi经验活度模型37，分别对TBAF、TBAB、TBAC和TEBAC作用下CO2气体水合物热力学相平衡数据进行了计算。2.4.1Chen-Guo热力学模型气体水合物相平衡条件预测模型多是基于范德瓦尔-普朗特38统计热力学基本模型发展而来的，但计算结果始终存在较大误差。不同于范德单吸附机理的水合物形成/分解机理，Chen-Guo模型认为水合物成核过程中同时进行着以下两种动力学过程：准化学反应动力学过程和吸附动力学过程。因此，在对各类气体水合物体系的相平衡条件进行预测时，Chen-Guo模型的计算结果与实验数值符合得更好。方程形式如下：f=f0(1)f=f01-(3)=Cf1+C

29、f=Cf1+Cf(4)C=Xexp(YTZ)C=XexpYT-Z(5)f0=f0(T)f0(P)f0(aw)f0=f0(T)f0(P)f0(aw)(6)f0(T)=Aexp(BTC)f0T=AexpBT-C(7)f0(P)=exp(PT)f0P=expPT(8)f0(aw)=a1/2wf0aw=aw-1/2(9)式中，值仅与水合物结构相关，对于S型结构为1/3，对于S型结构为2；CO2气体计算常数X、Y、Z分别为1.64641011Pa-1、2799.66、15.9；A、B、C数值分别为963.72/1015Pa、-6444.5、36.67；=V2R=V2R，可作为常数处理，对于S型水合物结构

30、=0.4242K/105Pa=0.4242K/105Pa，对于S型水合物结构=1.0224K/105Pa=1.0224K/105Pa；对于纯水体系液相活度awaw为1，对于促进剂体系，需要通过活度模型计算。2.4.2PR状态方程气相逸度的计算结果会直接影响气体水合物相平衡数据的预测精度。现有研究表明39，PR状态方程更适合用于预测CO2气体水合物相平衡条件。方程形式如下：P=RTVba(T)V2+2bVb2P=RTV-b-a(T)V2+2bV-b2(10)a=aR2T2.5cPca=aR2Tc2.5Pc(11)b=bRTcPcb=bRTcPc(12)2.4.3改进Joshi经验活度模型Greenspan等40认为液相活度与促进剂的浓度、系统的温度、压力相关。促进剂浓度和系统温度对水合物活度影响比较大，而压力变化对活度的影响比较小。因此，本文采用Joshi等37提出的液相活度awaw经验关联式，并参照Zhao等41的改进方法，忽略压力对液相活度awaw的影响，得到新的活度经验方程，如式（13）所示：aw=(1xm