劉振民，李 斌，薛永兵，王瀟瀟，王遠洋

(太原科技大學化學與生物工程學院，山西 太原 030024)

在生物柴油生產(chǎn)過程中大約會副產(chǎn)質(zhì)量分數(shù)10%的甘油， 快速發(fā)展的生物柴油工業(yè)已導(dǎo)致越來越嚴重的甘油過剩局面[1]。 近年來，以甘油為平臺化合物生產(chǎn)高附加值的產(chǎn)物引起了國內(nèi)外研究者的關(guān)注[2-3]，其中，碳酸甘油酯（GC）作為一種高附加值的甘油下游產(chǎn)物，具有高沸點、低凝固點、低揮發(fā)性、低可燃性、強極性、無毒無味、無腐蝕、無污染、與水互溶、生物可降解等特性，可廣泛用于涂料、膠黏劑、化妝品的潤飾劑以及藥物制備中的溶劑和氣體分離膜的合成等[4]。

以甘油為原料合成GC的方法主要有CO2法、尿素法和酯交換法等。 CO2法[5]主要以分子篩、離子交換樹脂或二丁基氧化錫為催化劑，CO2與甘油在超臨界的條件下可直接合成出GC，但是該路線反應(yīng)條件苛刻，需要使用助溶劑，且受熱力學的限制，GC產(chǎn)率較低。 尿素法[6]以路易斯酸為催化劑，需要負壓抽出反應(yīng)生成的氨氣來促進甘油的轉(zhuǎn)化，同時該路線存在氨氣回收的問題。 酯交換法[7]是由碳酸乙烯酯/碳酸烷基酯與甘油反應(yīng)合成GC的方法，該路線的甘油轉(zhuǎn)化率高、GC選擇性好， 是一條最具有工業(yè)化應(yīng)用前景的工藝路線。 亨斯曼公司以甘油和碳酸乙烯酯為原料，于2009年實現(xiàn)工業(yè)化生產(chǎn)，但產(chǎn)物GC和乙二醇的沸點相差較小，分離成本較高。 而甘油與碳酸二甲酯反應(yīng)制備GC因反應(yīng)條件溫和、操作設(shè)備簡單，且產(chǎn)物GC和甲醇的沸點相差較小而容易分離提純等優(yōu)點而廣受青睞。

目前, 國內(nèi)外對甘油/碳酸二甲酯（DMC）交換合成GC尚處在研究探索階段， 且關(guān)于GC和DMC的基礎(chǔ)熱力學數(shù)據(jù)尚無文獻報道。 Jiabo Li[8]等采用Benson和R-D基團貢獻法結(jié)合熱力學相關(guān)公式，只給出甘油/DMC 酯交換合成GC 反應(yīng)在298.15K、313.15K、333.15K和353.15K下的平衡常數(shù)?？蝶惥闧1]等采用Benson和ABWY基團貢獻法對酯交換反應(yīng)進行了熱力學計算， 但計算過程較簡略。 本文采用Benson基團貢獻法、Trouton規(guī)則、 對應(yīng)狀態(tài)法理論以及M-P法對甘油/DMC酯交換合成GC的反應(yīng)進行了詳細的熱力學計算和分析， 并采用Aspen 對該反應(yīng)進行模擬驗證，以期對后續(xù)實驗研究和可能的工業(yè)化生產(chǎn)提供一定的理論依據(jù)。

1 熱力學計算流程

1.1 反應(yīng)方程式

甘油/DMC酯交換反應(yīng)合成GC的反應(yīng)如圖式1。

圖式1 甘油/DMC酯交換合成GC的反應(yīng)方程式Scheme 1 Reaction formula of glycerol carbonate synthesis by transesterification of glycerol with DMC

1.2 計算方法

由于Benson基團貢獻法給出的是在氣態(tài)下化合物的生成焓、熵和等壓熱容，而在反應(yīng)條件下，反應(yīng)物均在液相狀態(tài)下進行化學反應(yīng)，因此，計算液態(tài)下化合物的生成焓和等壓熱容，需要知道蒸發(fā)焓的數(shù)據(jù)。 蒸發(fā)焓也稱汽化熱或蒸發(fā)潛熱，它是指同溫度下飽和蒸汽和飽和液體的焓差，即飽和液體汽化生成飽和蒸汽焓的變化值。

物質(zhì)的蒸發(fā)焓△Hν， 可根據(jù)物質(zhì)的沸點（Tb） 可由Trouton規(guī)則計算得到，其計算公式為：

所以，

式中：Tr為對應(yīng)態(tài)溫度，ω為偏心因子。

1.2.4 反應(yīng)體系△r和lnKθ與溫度T的關(guān)系

已知△r， 可根據(jù)基爾霍夫定律， 計算出△r(T)～T變化的函數(shù)。

再由范特霍夫（Vant’Hoff）方程式，可計算出lnKθ～T變化的函數(shù)。

2 計算過程

2.1 DMC和GC熱力學函數(shù)△f、△f與、的計算

文獻中尚無關(guān)于DMC、GC的焓值和熵值的數(shù)據(jù)，故對酯交換反應(yīng)進行熱力學計算，需要對DMC和GC的焓值和熵值進行估算。

表1 Benson法[10-12]估算DMC的標準生成焓值和熵值Table 1 Standard molar enthalpy of formation and standard molar entropy of dimethyl carbonate estimated by Bensonmethod

利用Benson 基團貢獻法求得的氣態(tài)DMC 的△f為-602.40 kJ/mol，而通過分子力學計算方法，經(jīng)過優(yōu)化DMC分子的構(gòu)型得到DMC的△f為-579.02 kJ/mol，二者較為吻合，說明Benson基團貢獻法的準確性較高[13]。 考慮DMC的蒸發(fā)焓，根據(jù)Trouton規(guī)則，按照公式△H=0.088Tb計算，則得到液態(tài)DMC的△f為634.37 kJ/mol，與文獻值-640.40 kJ/mol基本一致[1]。說明利用Benson基團貢獻法和Trouton規(guī)則求得的DMC的△f準確度高，且計算過程簡便。

表2 Benson法估算GC的標準生成焓和熵值Table 2 Standard molar enthalpy of formation and standard molar entropy of glycerol carbonate estimated by Benson method

通過Benson 基團貢獻法，由表2數(shù)據(jù)可計算出：

根據(jù)Trouton規(guī)則，△Hν=0.088Tb=55.18 kJ/mol，可進一步得到：

2.2 DMC和GC熱力學函數(shù)和的計算

要獲得反應(yīng)的熱力學性質(zhì)隨溫度的變化關(guān)系，就需要知道反應(yīng)體系中各物質(zhì)的熱容隨溫度的變化關(guān)系，甲醇和甘油的熱容值隨溫度的變化可以從數(shù)據(jù)手冊中查詢，但DMC和GC的熱容值則需要通過熱力學估算求得。

Benson 基團貢獻法給出的每種基團分別在300K、400K、500K溫度下的數(shù)值如表3所示。

表3 DMC和GC的Benson 法基團的值Table 3 Heat capacity values of DMC and GC estimated by Benson method

表3 DMC和GC的Benson 法基團的值Table 3 Heat capacity values of DMC and GC estimated by Benson method

?

表4 DMC和GC的～T之間的關(guān)系式Table 4 Relationship between heat capacities of DMC and GC and temperature estimated by the corresponding state method

表4 DMC和GC的～T之間的關(guān)系式Table 4 Relationship between heat capacities of DMC and GC and temperature estimated by the corresponding state method

?

根據(jù)Stermling和Brown提出對應(yīng)狀態(tài)法理論，即和之間的關(guān)系式計算DMC和GC的， 這里可采用M-P法[14]簡便計算DMC和GC的Tc和偏心因子ω，結(jié)果見表5。

表5 DMC和GC的Tc和ω值的計算Table 5 Critical temperature and acentric factor of DMC and GC estimated by M-P method

表6 DMC和GC的Cθp,l～T之間的關(guān)系式Table 6 Relationship between heat capacities of DMC and GC and temperature estimated by the corresponding state method

3 結(jié)果與討論

3.1 標準狀況下反應(yīng)的焓變、 熵變和吉布斯自由能變以及平衡常數(shù)的計算

除了DMC和GC的熱力學性質(zhì)需要估算外，甲醇和甘油的熱力學性質(zhì)可以從化學化工手冊中查詢得到。

表7 反應(yīng)物和產(chǎn)物的熱力學數(shù)據(jù)[14-15]Table 7 Thermodynamic data of reactants and products

根據(jù)標準狀況的數(shù)據(jù)計算， 可得到以下數(shù)據(jù)：酯交換焓變△r=10.42kJ/mol， 反 應(yīng)熵 變△r=124.65 J·mol-1·K-1，反應(yīng)吉布斯自由能變=平衡常數(shù)Kθ=5.05*104。

通過標準狀況下的熱力學數(shù)據(jù)計算表明，在101.325kPa，298.15K下，甘油和DMC酯交換合成GC反應(yīng)的焓變=10.42kJ/mol，表明該反應(yīng)在此條件下是吸熱反應(yīng)。 從反應(yīng)的熱力學方面考慮，溫度升高，平衡向GC生成的方向移動，甘油的轉(zhuǎn)化率增大，故適當提高反應(yīng)溫度可以促進GC的生成；從反應(yīng)的動力學方面考慮，反應(yīng)溫度升高，提高分子的運動動能， 增加甘油和DMC分子的有效碰撞數(shù)，反應(yīng)速率加快，縮短反應(yīng)時間；因此，從熱力學和動力學兩個方面考慮，提高反應(yīng)溫度對甘油的轉(zhuǎn)化率和GC的生成都有促進作用。

在101.325kPa，298.15K下，甘油和DMC酯交換合成GC反應(yīng)的吉布斯自由能變=-26.85kJ/mol，為負值，表明該反應(yīng)在標準狀況下可自發(fā)進行。

3.2 不同溫度下反應(yīng)的、lnKθ及計算

由于標準摩爾生成焓數(shù)據(jù)來源的限制，常常只能直接計算出298.15K時的標準摩爾反應(yīng)焓（298.15K）。 從表5和表6中的～T關(guān)系式可以計算得到總反應(yīng)的：

把上式代入到基爾霍夫公式（8）式并進行積分得到：

按照范特霍夫（Vant’Hoff）方程式：

總反應(yīng)的標準平衡常數(shù)可由積分求得，即

圖1 反應(yīng)焓變隨溫度變化的關(guān)系Fig. 1 Enthalpy change of transesterification reaction vs temperature

圖2 反應(yīng)的吉布斯自由能變隨溫度變化的關(guān)系Fig. 2 Gibbs free energy change of transesterification reaction vs temperature

圖3 平衡常數(shù)隨溫度變化的關(guān)系Fig. 3 Equilibrium constant of transesterification reaction vs temperature

從圖3中l(wèi)nK～T的關(guān)系圖中也可以看出， 溫度在300～390K范圍內(nèi)反應(yīng)的平衡常數(shù)隨溫度的增大而增大，也說明升高反應(yīng)溫度有利于反應(yīng)的進行。

4 Aspen模擬驗證

由Asp en 熱力學數(shù)據(jù)模擬得到在298.15K、101.325kPa的條件下，該反應(yīng)的焓變=13.72kJ/mol；反應(yīng)吉布斯自由能變= =-30.52kJ/mol，與上述結(jié)果相一致，驗證了估算的正確性。 進一步還模擬了在101.325kPa、300～400K條件下該反應(yīng)的焓變和吉布斯自由能變，見表8。

表8 Aspen 模擬該反應(yīng)的焓變和吉布斯自由能變Table 8 Aspen simulates the reaction in enthalpy and gibbs free energy

由表8可以看出在101.325kPa、300～400K條件下， 該反應(yīng)的焓變, 反應(yīng)吉布斯自由能變，說明在此溫度范圍內(nèi)，反應(yīng)為可自發(fā)進行的吸熱反應(yīng)，與上述計算結(jié)果相一致。

5 結(jié)論

本文采用Benson基團貢獻法計算DMC和GC的標準摩爾生成焓、 標準摩爾熵和摩爾等壓熱容，采用Trouton規(guī)則計算DMC和GC的蒸發(fā)焓。 采用熱力學基本公式計算得到標準狀態(tài)下甘油與DMC酯交換合成GC的焓變=10.42kJ/mol、 熵變=124.65J/(mol·K)、 吉布斯自由能變=-26.85kJ/mol和平衡常數(shù)lnKθ=10.83。

在300～390K范圍內(nèi)， 進一步分析該反應(yīng)的焓變、 吉布斯自由能變和平衡常數(shù)隨溫度的變化關(guān)系，分析結(jié)果表明，反應(yīng)的焓變隨溫度的升高而降低，且其值總是大于零，說明在此溫度范圍內(nèi)，該反應(yīng)為吸熱反應(yīng)，升高溫度有利于反應(yīng)的進行；反應(yīng)的吉布斯自由能變隨溫度的升高而降低，且其值總為負值， 進一步表明升高溫度有利于反應(yīng)的進行；反應(yīng)的平衡常數(shù)隨溫度的增大而增大，也說明升高反應(yīng)溫度有利于反應(yīng)的進行。 同時，采用Aspen對該反應(yīng)的熱力學數(shù)據(jù)進行模擬，得到該酯交換反應(yīng)的焓變?yōu)檎?、吉布斯自由能變?yōu)樨撝?，表明在標準狀態(tài)下， 酯交換反應(yīng)為吸熱反應(yīng)且可自發(fā)進行，與計算結(jié)果相一致。

符號說明

△Hi—基團i的焓值，kJ/mol；△Si—基團i的 熵值，J·mol-1·K-1；—基團i的熱容值，J·mol-1·K-1；—反應(yīng)的焓變， kJ·mol-1；—反應(yīng)的吉布斯自由能變，kJ/mol；—理想氣體的標準摩爾生成焓，kJ/mol；—理想氣體的標準摩爾體熵，J·mol-1·K-1；—理想氣體的標 準摩爾熱容，J·mol-1·K-1；—液體的標準摩爾生成焓，kJ/mol；—液體的標準摩爾熵，J·mol-1·K-1；—液體的標準摩爾熱容，J·mol-1·K-1；△Hν—化合物的蒸發(fā)焓，kJ/mol；Tb—化合物的沸點，K；Tc—化合物的臨界溫度，K；ω—偏心因子；Kθ—平衡常數(shù)；δ—對稱數(shù)；η—對映(光學)異構(gòu)體；Scorrect—化合物熵值的環(huán)式修正值，J·mol-1·K-1。