1. 引言
全球油页岩按油换算后其储存量大约有6902亿吨,达世界石油资源储量的3倍。如图1所示,在美国,页岩油的储量有5369亿吨 [1]。近年来美国在有关深藏页岩油开采技术方面获得了突飞猛进的进步,其关键技术特征为水平井技术和水力压裂技术。2014年其页岩油的产量高达美国本土常规石油产量6亿吨的三分之一2.1亿吨。由此国际油气资源供应的状况被完全改变,媒体将其称之为“页岩革命” [2]。我国的油页岩按油换算后储量为476亿吨,达全国可开发石油储量233亿吨的2倍 [3]。而今,在中国对油页岩的开发应用已经实现了商业化。油页岩主要挖掘埋藏在100~400 m处的浅层油页岩,主要选用在近地表面露天挖掘的方法来开采。截止到2014年,中国通过油页岩炼制得到的页岩油产量已经达到100万吨 [4]。预计到2020年,页岩油的产能规模将突破300万吨 [5]。针对我国发展替换油气措施的新兴产业、减缓油气资源不足的问题,油页岩的开发利用将起到关键性的支撑作用。
Figure 1. Shale oil-equivalent reserves in the world [1]
图1. 世界页岩油资源储量 [1]
油页岩的热解干馏工艺分为地上干馏和地下干馏 [6]。中国只有地上干馏工艺实现了工业化应用,代表技术为抚顺炉干馏技术。然而抚顺炉技术存在油页岩利用率低(80%)、产品收率低(65%)和油中轻质组分(沸点小于350℃)含量低(20%)等缺点 [4]。中国科学院过程工程研究所发明了一种间接加热内构件移动床干馏炉,并且建造了处理规模为1000吨/年的油页岩的中试装置 [7]。相比抚顺炉技术,新技术使得油页岩的利用率提高20%,油收率提高25%,页岩油中轻质组分的含量提高46% [8]。
通过间接加热内构件移动床技术干馏油页岩所需的热量主要来源于干馏气的循环燃烧。然而分离油气分离单元得到的干馏气经燃烧尚不足以提供足够的热量供油页岩干馏,所以需要补充其它燃料气气源,如焦炉气、天然气或者炼厂气等。然而,燃料气的来源和运输长距离是待攻克的难题。由于油页岩开采过程中伴有褐煤产生,所以文章以间接加热内构件移动床油页岩干馏过程为基础,提出煤辅助油页岩高效热解过程,如图2所示,热解过程包括油页岩的干馏单元、煤气化单元、分离油气单元和半焦燃烧供热发电单元。经过提前研磨和筛分的油页岩,颗粒小于10 mm的页岩在干馏炉中进行热解反应。页岩颗粒热解所需的热量一部分来自于热解产生的干馏气全部循环,这部分热量不足,尚需提供额外燃料气。油页岩开采过程伴生煤气化产生的合成气燃烧补充油页岩热解所需额外的热量。热解反应后,在干馏炉的顶部有油气混合物产生,混合物中的粉尘经旋风分离器去除。油气混合物经过脱除粉尘后进入油气分离器,经分离器分离后成为页岩油和干馏气。在燃烧室内,干馏气和煤气化合成气燃烧释放的热量用于页岩颗粒在干馏炉内的热解。另外,热解过程中会有半焦产生,其从干馏炉底部排出,经干燥后再次进入循环的流化床系统,燃烧后生成的高压蒸汽用于蒸汽轮机发电。
Figure 2. Coal-assisted oil shale pyrolysis process
图2. 煤辅助油页岩热解过程
本文第二节首先进行了建模,模拟了煤辅助油页岩高效热解干馏过程的关键单元(干馏单元);第三节分析了煤辅助油页岩高效热解过程的技术和经济性能,将其结果和已有的抚顺式油页岩干馏作对比,指出煤辅助油页岩高效热解过程的优势和不足。
2. 建模与模拟
文章采用Aspen Plus (Version 7.2)流程模拟软件对煤辅助油页岩热解进程进行模拟,重点是油页岩干馏单元的建模,模拟其干馏过程。油页岩和煤的工业分析见表1,其中M、FC、V和A分别表示水分、固定碳、挥发性和灰分,ar表示油页岩的接收基。
Table 1. Proximate analyses of oil shale and coal
表1. 油页岩和褐煤的工业分析
文章采用Aspen Plus (Version 7.2)流程模拟软件对煤辅助油页岩热解进程进行模拟,重点是油页岩干馏单元的建模,模拟其干馏过程。油页岩的热解规模为300万吨/年(375 t/h)。原料油页岩选用吉林桦甸含油率为9%的油页岩,原料煤采用油页岩开采伴生的褐煤,油页岩和褐煤的工业分析数据如表1所示 [4] [8]。
油页岩干馏单元
间接加热内构件移动床干馏炉的模拟计算程序步骤和相应的Aspen Plus模拟流程如图3所示,干馏炉由干燥段和干馏段构成,并分段对其进行建模与模拟。本文选用RYield反应器模块对干燥段进行模拟,设定反应器模块的温度为200℃,压力为0.1 MPa。另外,采用RCSTR反应器模块对干馏段进行模拟,设定反应器的温度为550℃,压力为0.1 MPa。
Figure 3. Aspen simulation calculation procedure
图3. 移动床干馏炉模拟计算程序步骤
干馏炉中干燥段发生的反应是页岩颗粒表面水和内部毛细管水份的蒸发,热解页岩中油母质主要发生在干馏段,热解反应方程如式(1):
Kerogen → H2 + H2O + H2S + NH3 + CO + CO2 + CH4 + C2H4 + C2H6 + C3H8+ C4H10O + Semi-coke(1)
油页岩的热解过程的动力学方程,如式(2)所示 [9]:
(2)
其中x指转化率,t指反应时间,k指反应比速,n指反应级数。根据阿伦尼乌斯方程,速率常数k的表达式,如式(3)所示 [10]:
(3)
其中A代表指前因子,Ea指表观活化能,R指气体常数(R = 8.314 J/mol K),T表示反应温度。依据文献 [4] [11],Ea为158 KJ/mol、A为2.03 × 109、n为1.0。动力学方程可表达为下式(4):
(4)
通过对间接加热内构件移动床技术进行模拟,获得单元油页岩经干馏热解后各产物的产率分布,结果如表2所示。对比模拟数据与实验数据 [8],相对误差小于8%。
Table 2. Comparison simulation and industrial data
表2. 移动床干馏单元页岩热解产物产率分布
3. 技术经济分析
3.1. 技术分析
油页岩干馏单元页岩使其热解所需的热量由干馏气和煤气化合成气经循环燃烧产生热量提供。本文以干馏炉为基础,计算其热力学平衡。页岩带入的热量和页岩热解所需的热量等于页岩热解反应热和页岩热解各产物带出的热量,热平衡方程如(5)所示:
(5)
其中Qshale表示原料页岩带入的热量,Qret表示页岩热解反应所需的热量,Qpyr表示页岩热解产生的热,Qoil、Qgas、Qsemi-coke和Qvap分别表示页岩油、半焦、干馏气以及水分在干馏炉出口处带出的热量。通过如下式(6)计算原料页岩带入的热量Qshale,如下式所示:
(6)
其中mshale表示页岩质量流量,Cp,shale表示页岩比热容。根据有关作者的研究 [8],桦甸页岩的比热容Cp,shale为1.05 kJ/(kg∙℃),通过式(6)可计算得到Qshale为60.8 kW。其次,油页岩热解反应产生的热Qpyr可由下式(7)计算,
(7)
其中
表示蒸发水分的质量流量,
表示单位页岩热解反应产生的热量。根据相关文献可知 [4],
为235 kJ/kg。可代入式(7)计算Qpyr值是263.0 kW。Qoil热解页岩油在干馏炉出口处带出的热量可由下式计算,如式(8)所示:
(8)
其中moil表示干馏热解页岩油的质量流量,Cp,oil表示所需比热容,Tret表示热解温度,roil表示汽化潜热。由相关文献可知 [8],页岩油的比热容Cp,shale值为2.26 kJ/(kg∙℃),roil汽化潜热值为368.46 kJ/kg。Qoil可由式(8)计算获得123.0 kW。Qgas表示热解产物干馏气在出口处带出的热量,其计算如式(9)所示:
(9)
其中mgas表示热解产物干馏气的质量流量,Cp,gas表示干馏气比热容。干馏气在出口处温度可经模拟得到为360℃,Cp,gas表示热容,其大小为6.17 kJ/(kg∙K)。Qgas可由式(9)计算获得,其值为305.3 kW。热解产物半焦带出的热量Qsemi-coke的计算,如式(10)所示:
(10)
其中msemi-coke代表产物半焦的质量流量,Cp,semi-coke代表比热容。根据已有研究可知其比热容Cp,semi-coke为1.8 kJ/(kg∙℃) [8],则Qsemi-coke可由式(10)计算得到,大小为883.3 kW。Qvap表示水份蒸发所带出的热量,计算式如式(11)所示:
(11)
其中msemi-coke为水分蒸发掉的质量流量。经式(11)计算获得Qvap为148.4 kW。代入热平衡方程式(5)可计算得到一台干馏炉热解油页岩所需1662.2 kW的热量。包括90台干馏炉的干馏过程,因此整个干馏过程热解页岩油所需Qret为149.6 MW的热量。页岩热解所需热量由干馏气和煤气化合成气燃烧提供。文章利用RStoic反应器模型来模拟气体燃炉,设定反应温度750℃,压力为0.1 MPa。Qcomb表示干馏气燃烧所放出热量,其计算式如式(12)所示:
(12)
其中wk为干馏气中各组分的占比,ΔhL表示低位热值。其中在干馏过程中油气分离所得到的干馏气经全部循环燃烧释放的热量可由式(12)计算得到,其热量为89.6 MW。由于页岩热解干馏需要149.6 MW的热量,因此需要通过煤气化合成气来补充额外的热量Qsupl,其计算式如式(13)所示:
(13)
通过式(13)计算可知,需要额外补充的热量为60 MW,才能足够满足整个油页岩干馏过程所需的热量。需要的煤气化合成气的量的计算,如式(14)所示:
(14)
其中Vsyn表示合成气的体积流量,
和
分别表示各组分的体积占比和低位热值。所需合成气的量Vsyn利用(14)式计算获得,需5.3 × 104 Nm3/h合成气,耗损煤炭的量为16.9 t/h。
3.2. 经济性能分析
经济性能分析指标是生产成本,生产成本(Production cost, PC)主要涉及原料购买所需费用、公用工程费用、折旧、生产管理、以及设备维修所需费用,总生产费用通过如下式(15)所示估算:
PC = CR + CU + CO&M + CD + CPOC + CAC + CDSC (15)
其中CR为原料成本,CU为公共工程所需金额,CO&M为设备维修和操作耗费,CD和CPOC分别为折旧费和工厂管理费,CAC、CDSC为行政费和销售费用。其中依据模拟可计算出CR原料费和CU公用工程用度。操作人工费用的计算参考作者之前的工作 [3],根据生产所需成本的10%来计算。根据20年的直线折旧来计算折旧金额,4%的残值 [12]。
Figure 4. Production costs of the IHMB and Fushun refinery processes
图4. 煤辅助油页岩热解(a)和抚顺式(b)油页岩热解过程的生产成本概率分布
通过计算得到煤辅助油页岩热解过程和传统的抚顺热解过程的生产成本,如图4所示为单位页岩油生产所需成本的概率分布图,从图中可以看出:1) 两个过程的平均生产成本分别为2636 CNY/t和3024 CNY/t,新过程的生产成本相比相同处理规模的抚顺式炼制过程的生产成本低12%,主要由于间接加热内构件移动床炼制过程的原料利用率高和油收率高。2) 具有经济竞争力的原油价格新过程价格约为51 $/bbl,抚顺式干馏原油售价58 $/bbl。全球油价目前稳定在50 $/bbl附近,在这种情况下,间接加热内构件移动床炼制过程具有一定的经济竞争力,相比传统的抚顺式炼制间接加热内构件移动床技术炼制具有绝对的经济竞争力。然而当国际油价恢复到达100 $/bbl的高位时,抚顺式和新的油页岩炼制技术过程均拥有较强的经济竞争力。
4. 结论
为解决传统抚顺炉技术瓶颈,中国科学院过程工程研究所研发了新的间接加热内构件移动床技术。新的移动床技术热解油页岩所需的热量一部分来自于热解产生的干馏气全部循环,这部分热量不足,尚需要提供额外燃料气,这部分燃料气由油页岩开采过程伴生煤气化合成气提供。本研究经过对煤辅助油页岩热解反应建模和模拟分析,得到以下结论:1) 热力学分析表明热解375 t/h的油页岩需要149.6 MW的热量,其中干馏气全部循环可释放89.6 MW的热量,剩余60 MW的热量由煤气化合成气供给。2) 经济分析表明采用新过程生产单位页岩油需成本2636元/t,相比抚顺式,其成本降低12%。归因于间接加热内构件移动床技术炼制时原料油页岩利用率和油收率分别比抚顺炉炼制高出20%和25%。
NOTES
*第一作者。