1. 引言
随着我国对能源需求量不断增大,石油和煤炭等资源的储量岌岌可危,其开采利用会对环境产生负面影响。而费托合成技术能缓解对环境和能源的担忧,通过将产物继续加工处理后能产生清洁燃料等多种高附加值的化工产品,可以替代石油产品[1]。
费托合成(Fischer-Tropsch)是在温度为200℃~350℃,压力为0.5~5 MPa将合成气(CO和H2)转化成烃类[2] [3]。通常情况下合成气中可能会含有HCl、H2S、PH3、CS2等气态杂质,会导致催化剂中毒进而影响催化剂的活性[4] [5]。目前,对于气态杂质引起的钴基催化剂失活主要集中在硫化物和氮化物,并且微量硫或微量氮就能导致催化剂失活[6]。Visconti等[7]研究了Co/Al2O3催化剂的硫中毒,实验证明,硫含量与催化剂活性呈负相关,当硫含量达到2000 µg/g时会导致催化剂完全失活。Saib等[8]研究了不同氮化物对钴基催化剂的中毒作用,结果表明,钴基催化剂的失活取决于氮化物的种类和含量,其失活是可逆的。除此之外,专门针对合成原料气氛下HCl中毒的热力学研究报道有限。
费托合成催化剂最主要的两个部分是活性成分和载体。铁、钴、钌和镍能作为费托合成催化剂的活性成分,但钌基催化剂成本高,镍基催化剂难以进行催化,故铁和钴被广泛应用于工业[5] [9] [10]。铁基催化剂活性低,适合于生产低碳烯烃[11],而钴基催化剂增长F-T链的能力强、活性高、反应稳定且不易中毒,被认为是一种良好的催化剂[12]。因此,本文主要对钴基催化剂进行研究。此外,良好的催化剂载体可以提供高比表面积[13]。费托合成催化剂的载体主要是SiO2、Al2O3、TiO2、ZrO2等[14]。其中,SiO2的比表面积较高,有利于后续金属的分散[15]。本文以Co/SiO2催化剂为研究对象,并对Co-HCl、Co-HCl-CO、Co-HCl-H2、SiO2-HCl、SiO2-HCl-CO、SiO2-HCl-H2反应体系进行热力学模拟计算,并讨论了可能产生的中毒产物,为延长催化剂寿命和降低成本提供理论参考。
2. 热力学计算与分析
2.1. 热力学计算基础
根据纯物质的物理性质和热化学属性,按照吉布斯自由能最小原理,采用HSC Chemistry 5.0软件对钴基催化剂HCl中毒热力学进行模拟计算。计算前通过前查阅资料对可能产生的化学产物进行了分析[16]。计算时输入的压力为101.325 KPa,温度范围是200℃~350℃,升温步长为10℃。
2.2. 热力学分析
2.2.1. Co-HCl体系
体系中HCl单独存在时,与催化剂中的活性组分Co可能发生的反应见式(1)~(5)。图1是通过HSC Chemistry 5.0软件计算得到的反应式(1)~(5)吉布斯自由能(ΔG)随温度变化的关系。由图1可知,在200℃~350℃时,该体系内仅反应(1)的ΔG < 0,其余反应的ΔG > 0。由此可见,反应(1)能自发向右进行,HCl单独存在时会导致钴基催化剂失活,CoCl2是该体系产生的中毒产物。此外,该反应的ΔG随温度的升高而升高,表明升高温度能缓解钴基催化剂的失活。
Co + 2HCl(g) = CoCl2 + H2(g) (1)
2Co + 2HCl(g) = 2CoH(g) + Cl2(g) (2)
Co + 3HCl(g) = CoH(g) + 1.5Cl2(g) + H2(g) (3)
2Co + 2HCl(g) = 2CoCl(g) + H2(g) (4)
Co + 4HCl(g) = CoCl2 + 2H2(g) + Cl2(g) (5)
Figure 1. Relationship between Gibbs energy and temperature of equation (1)~(5)
图1. 式(1)~(5)的吉布斯自由能与温度的关系
2.2.2. Co-HCl-CO体系
体系中HCl与CO共存时,与催化剂中的活性组分Co可能发生的反应见式(6)~(12)。反应式(6)~(12)的吉布斯自由能与温度的关系见图2。由图2可知,在200℃~350℃时,该体系内仅反应(6)的ΔG < 0,其余反应的ΔG > 0。因此,反应式(6)能自发向右进行。C和CoCl2是该体系产生的中毒产物。
Figure 2. Relationship between Gibbs energy and temperature of equation (6)~(12)
图2. 式(6)~(12)的吉布斯自由能与温度的关系
Co + 2HCl(g) + CO(g) = CoCl2 + H2O(g) + C (6)
2Co + 2HCl(g) + 2CO(g) = 2CoO + 2CCl(g) + H2(g) (7)
3Co + 4HCl(g) + 4CO(g) = Co3O4 + 4CCl(g) + 2H2(g) (8)
2Co + 4HCl(g) + 2CO(g) = 2CoCl2 + O2(g) + 2H2(g) + 2C (9)
Co + 2HCl(g) + 2CO(g) = CoO + H2O(g) + 2CCl(g) (10)
3Co + 8HCl(g) + 8CO(g) = Co3O4 + 4H2O(g) + 8CCl(g) (11)
2Co + 2HCl(g) + CO(g) = C + H2O(g) + 2CoCl(g) (12)
2.2.3. Co-HCl-H2体系
体系中HCl与H2共存时,与催化剂中的活性组分Co可能发生的反应见式(13)~(14)。图3是反应式(13)~(14)的吉布斯自由能与温度的关系。由图3可知,两个反应的ΔG与温度呈负相关,并且ΔG > 0,表明该系列反应不能自发向右进行,不会造成钴基催化剂中毒。
4Co+2HCl(g)+H2 (g)=4CoH(g)+Cl2(g) (13)
4Co + 2HCl(g) + H2(g) = 4CoH(g) + 2Cl(g) (14)
Figure 3. Relationship between Gibbs energy and temperature of equation (13)~(14)
图3. 式(13)~(14)的吉布斯自由能与温度的关系
2.3. SiO2载体与HCl、CO、H2共存催化剂反应体系
2.3.1. SiO2-HCl体系
当SiO2作为催化剂的载体时,体系中仅存在HCl时,其可能发生的反应见式(15)~(36)。图4是反应式(15)~(36)的吉布斯自由能与温度的关系。在图4中可以观察到在200℃~350℃范围内时,所有反应式的ΔG > 0,且ΔG与温度的变化不明显。由此表明,HCl不会对载体SiO2产生影响。
SiO2 + 4HCl(g) = SiCl(g) + 2H2O(g) + 3Cl(g) (15)
SiO2 + 4HCl(g) = SiCl2(g) + 2H2O(g) + 2Cl(g) (16)
SiO2 + 4HCl(g) = Si + 2H2O(g) + 2Cl2(g) (17)
SiO2 + 4HCl(g) = SiCl2(g) + 2H2O(g) + Cl2(g) (18)
2SiO2 + 8HCl(g) = 2SiCl3(g) + 4H2O(g) + Cl2(g) (19)
SiO2 + 4HCl(g) = SiCl4(g) + 2H2O(g) (20)
SiO2 + 4HCl(g) = SiCl3(g) + 2H2O(g) + Cl(g) (21)
2SiO2 + 2HCl(g) = 2SiCl(g) + H2(g) + 2O2(g) (22)
2SiO2 + 2HCl(g) = 2SiCl(g) + H2O(g) + 1.5O2(g) (23)
2SiO2 + 4HCl(g) = 2SiCl2(g) + 2H2O(g) + O2(g) (24)
4SiO2 + 12HCl(g) = 4SiCl3(g) + 6H2O(g) + O2(g) (25)
SiO2 + 4HCl(g) = SiCl4(g) + 2H2(g) + O2(g) (26)
SiO2 + 4HCl(g) = SiCl4(g) + H2O2(g) + H2(g) (27)
SiO2 + 5HCl(g) = SiH(g) + 2H2O(g) + 2.5Cl2(g) (28)
SiO2 + 8HCl(g) = SiH4(g) + 2H2O(g) + 4Cl2(g) (29)
SiO2 + 2HCl(g) = SiO(g) + H2O(g) + Cl2(g) (30)
SiO2 + 4HCl(g) = Si(g) + 2H2O(g) + 2Cl2(g) (31)
SiO2 + 4HCl(g) = Si + 2H2O(g) + 2Cl2(g) (32)
2SiO2 + 8HCl(g) = Si2(g) + 4H2O(g) + 4Cl2(g) (33)
3SiO2 + 12HCl(g) = Si3(g) + 6H2O(g) + 6Cl2(g) (34)
SiO2 + 4HCl(g) = SiCl4(g) + 2OH(g) + H2(g) (35)
SiO2 + 2HCl(g) = SiCl2(g) + 2OH(g) (36)
Figure 4. Relationship between Gibbs energy and temperature of equation (15)~(36)
图4. 式(15)~(36)的吉布斯自由能与温度的关系
2.3.2. SiO2-HCl-CO体系
SiO2载体与HCl和CO共存时,三者可能发生的反应见式(37)~(78)。图5是反应式(37)~(78)的吉布斯自由能与温度的关系。由图5可知,在200℃~350℃时,所有反应式的ΔG > 0。因此,该体系内的反应不会自发向右进行,合成气中的燃料气体CO和杂质气体HCl共存时也不会对载体SiO2产生影响。
2SiO2 + 2HCl(g) + 3CO(g) = 2SiCl(g) + H2O(g) + 3CO2(g) (37)
SiO2 + 2HCl(g) + CO(g) = SiCl2(g) + H2O(g) + CO2(g) (38)
2SiO2 + 6HCl(g) + CO(g) = 2SiCl3(g) + 3H2O(g) + CO2(g) (39)
SiO2 + 4HCl(g) + CO(g) = SiCl4(g) + H2O(g) + CO2(g) + H2(g) (40)
2SiO2 + 6HCl(g) + 3CO(g) = 2SiCl3(g) + H2O(g) + 3CO2(g) + 2H2(g) (41)
SiO2 + 2HCl(g) + 2CO(g) = SiCl2(g) + 2CO2(g) + H2(g) (42)
2SiO2 + 2HCl(g) + 4CO(g) = 2SiCl(g) + 4CO2(g) + H2(g) (43)
SiO2 + 4HCl(g) + 2CO(g) = SiH4(g) + 2Cl2(g) + 2CO2(g) (44)
2SiO2 + 6HCl(g) + 4CO(g) = Si2H6(g) + 3Cl2(g) + 4CO2(g) (45)
SiO2 + 2HCl(g) + CO(g) = SiO(g) + H2(g) + CO2(g) + Cl2(g) (46)
2SiO2 + 2HCl(g) + 4CO(g) = 2SiH(g) + Cl2(g) + 4CO2(g) (47)
SiO2 + 3HCl(g) + CO(g) = SiH(g) + 1.5Cl2(g) + CO2(g) + H2O(g) (48)
SiO2 + 3HCl(g) + 2CO(g) = SiH(g) + 1.5Cl2(g) + 2CO2(g) + H2(g) (49)
SiO2 + 6HCl(g) + 2CO(g) = SiH4(g) + 3Cl2(g) + 2CO2(g) + H2(g) (50)
2SiO2 + 10HCl(g) + 4CO(g) = Si2H6(g) + 5Cl2(g) + 4CO2(g) + 2H2(g) (51)
3SiO2 + 2HCl(g) + 6CO(g) = 2SiH(g) + SiCl2(g) + 6CO2(g) (52)
3SiO2 + 4HCl(g) + 6CO(g) = SiH4(g) + 2SiCl2(g) + 6CO2(g) (53)
5SiO2 + 6HCl(g) + 10CO(g) = Si2H6(g) + 3SiCl2(g) + 10CO2(g) (54)
2SiO2 + HCl(g) + 4CO(g) = SiH(g) + SiCl(g) + 4CO2(g) (55)
5SiO2 + 4HCl(g) + 10CO(g) = SiH4(g) + 4SiCl(g) + 10CO2(g) (56)
8SiO2 + 6HCl(g) + 16CO(g) = Si2H6(g) + 6SiCl(g) + 16CO2(g) (57)
3.5SiO2 + 6HCl(g) + 7CO(g) = Si2H6(g) + 1.5SiCl4(g) + 7CO2(g) (58)
2SiO2 + 4HCl(g) + 4CO(g) = SiH4(g) + SiCl4(g) + 4CO2(g) (59)
5SiO2 + 4HCl(g) + 10CO(g) = 4SiH(g) + SiCl4(g) + 10CO2(g) (60)
4SiO2 + 3HCl(g) + 8CO(g) = 3SiH(g) + SiCl3(g) + 8CO2(g) (61)
4SiO2 + 6HCl(g) + 8CO(g) = Si2H6(g) + 2SiCl3(g) + 8CO2(g) (62)
2SiO2 + 6HCl(g) + 7CO(g) = Si2H6(g) + 3COCl2(g) + 4CO2(g) (63)
SiO2 + 4HCl(g) + 4CO(g) = SiH4(g) + 2COCl2(g) + 2CO2(g) (64)
2SiO2 + 2HCl(g) + 5CO(g) = 2SiH(g) + COCl2(g) + 4CO2(g) (65)
2SiO2 + 2HCl(g) + 6CO(g) = 2SiH(g) + CCl2(g) + 5CO2(g) (66)
SiO2 + 3HCl(g) + 4CO(g) = SiH(g) + 1.5CCl2(g) + 2.5CO2(g) + H2O(g) (67)
4SiO2 + 4HCl(g) + 10CO(g) = 4SiH(g) + CCl4(g) + 9CO2(g) (68)
SiO2 + 4HCl(g) + 4CO(g) = SiH4(g) + CCl4(g) + 3CO2(g) (69)
SiO2 + 16HCl(g) + 4CO(g) = SiH4(g) + 4CCl4(g) + 6H2O(g) (70)
2SiO2 + 28HCl(g) + 7CO(g) = Si2H6(g) + 7CCl4(g) + 11H2O(g) (71)
SiO2 + 10HCl(g) + 2.5CO(g) = SiH(g) + 2.5CCl4(g) + 4.5H2O(g) (72)
SiO2 + 25HCl(g) + 10CO(g) = SiH(g) + 5C2Cl5(g) + 12H2O(g) (73)
5SiO2 + 5HCl(g) + 2CO(g) = 5SiH(g) + C2Cl5(g) + 6O2(g) (74)
SiO2 + 6HCl(g) + CO(g) = SiC(g) + 3Cl2(g) + 3H2O(g) (75)
SiO2 + 6HCl(g) + CO(g) = SiC + 3Cl2(g) + 3H2O(g) (76)
SiO2 + 2HCl(g) + CO(g) = Si + Cl2(g) + H2O(g) + CO2(g) (77)
2SiO2 + 4HCl(g) + 2CO(g) = Si2(g) + 2Cl2(g) + 2H2O(g) + 2CO2(g) (78)
Figure 5. Relationship between Gibbs energy and temperature of equation (37)~(78)
图5. 式(37)~(78)的吉布斯自由能与温度的关系
2.3.3. SiO2-HCl-H2体系
SiO2载体与HCl和H2共存时可能发生的反应见式(79)~(88)。图6是反应式(79)~(88)的吉布斯自由能与温度的关系。由图6可知,在200℃~350℃时,所有反应式的ΔG > 0,并且大部分反应的ΔG值集中在500~1000 KJ/mol。由此表明,该体系内的反应不能自发向右进行。
SiO2 + HCl(g) + 1.5H2(g) = SiCl(g) + 2H2O(g) (79)
SiO2 + 2HCl(g) + H2(g) = SiCl2(g) + 2H2O(g) (80)
2SiO2 + 6HCl(g) + H2(g) = 2SiCl3(g) + 4H2O(g) (81)
2SiO2 + 8HCl(g) + H2(g) = 2SiCl4(g) + 4H2O(g) + 2H(g) (82)
SiO2 + 2HCl(g) + 1.5H2(g) = SiH(g) + 2H2O(g) + Cl2(g) (83)
SiO2 + 2HCl(g) + 3H2(g) = SiH4(g) + 2H2O(g) + Cl2(g) (84)
2SiO2 + 4HCl(g) + 5H2(g) = Si2H6(g) + 4H2O(g) + 2Cl2(g) (85)
SiO2 + 2HCl(g) + H2(g) = Si(g) + 2H2O(g) + Cl2(g) (86)
2SiO2 + 4HCl(g) + 2H2(g) = Si2(g) + 4H2O(g) + 2Cl2(g) (87)
3SiO2 + 6HCl(g) + 3H2(g) = Si3(g) + 6H2O(g) + 3Cl2(g) (88)
Figure 6. Relationship between Gibbs energy and temperature of equation (79)~(88)
图6. 式(79)~(88)的吉布斯自由能与温度的关系
4. 结论
(1) 在200℃~350℃操作温度范围内,体系中仅存在HCl时,与催化剂中的活性组分Co会产生中毒产物CoCl2。
(2) 体系中共存CO时,Co对HCl的抗毒性较差,会产生中毒产物C和CoCl2。
(3) 在200℃~350℃操作温度范围内,HCl无论是单独存在还是与合成气(CO、H2)共存时,均不能与SiO2发生反应。表明费托合成钴基催化剂中,HCl不会对SiO2产生影响。因此,SiO2适合作为钴基催化剂的载体。
NOTES
*通讯作者。