马强 林日亿 韩超杰 刘浩
1.中国石油大学(华东)储运与建筑工程学院 2.中国石油天然气管道工程有限公司 3.潍坊市政工程设计研究院有限公司
注蒸汽热采是开采稠油最常见的技术,通过向油藏中注入高温蒸汽加热原油,降低其黏度,增强原油流动性,达到开采的目的[1-2],同时也会促使原油、地层水以及地层矿物之间发生一系列复杂的化学反应,进而产生CO2、H2S等气体[3-5]。H2S是一种有浓烈臭味的剧毒气体,当H2S的体积分数达到0.002%时,可闻到强烈的臭味;当工作环境中H2S的体积分数在0.1%以上时,会造成人员死亡;而且,H2S具有很强的化学活性,对开发运输储藏器具具有极强的腐蚀作用[6-7]。在注蒸汽过程中,越来越多的油气井产生了H2S,危害工作人员的生命安全,在对设备造成极其严重的腐蚀的同时,会出现H2S沉积堵塞现象,导致地面脱硫设备投资增加问题,同时也对环境造成污染。
研究认为,油气藏中的H2S可能来自于:硫化物的细菌还原(bacterial sulfate reduction,BSR)、有机硫化物的热裂解(thermal decomposition of sulfate,TDS)和硫酸盐热化学还原(thermochemical sulfate reduction,TSR)。H2S的毒性对微生物的影响导致通过BSR方式很难生成大量的H2S,而油藏中含硫有机物的含量通常极低,这也使得通过TDS方式产生的H2S不会很多。因此,目前主要认为H2S由TSR反应产生,为此国内外大量学者进行了TSR生成H2S的研究。
丁康乐等[8]以MgSO4为研究对象,分别在450 ℃、480 ℃、510 ℃时探究了TSR反应机理。实验表明,温度和烷烃链长度是影响TSR进行的两个重要因素。
Zhang等[9]通过实验探究在初始无低价硫存在情况下碳氢化合物的TSR反应动力学参数。结果表明:H2S的生成量与CaSO4的用量成正比;pH值可以显著影响动力学参数,随着pH值的降低,反应速率显著增加。Zhang等[10]在定温恒压条件下对含有不同官能团的模型化合物与MgSO4的TSR反应,确定了烃类型对反应的影响。实验分别对含有不同官能团(无硫)的单个有机化合物在MgSO4的存在下进行了反应,并测定了H2S生成量。实验结果表明:不同的烃类型对反应活化能影响很大,硫酸盐还原过程中所涉及的烃类反应活性的差异主要取决于所呈现的官能团的类型;模型化合物反应顺序为:1-辛烯> 1-辛醇> 1-辛酮 >正辛烷>辛烷酸>正辛基苯>二甲苯。
Xia等[11]通过前人的实验数据以及现场数据,研究了气态烃类对TSR的自催化作用。结果表明,非自催化反应主要是由于缺少必要的中间产物。计算了受TSR影响的气态烃和碳同位素分馏的动力学参数,表明气体烃化物TSR反应需要高温条件。
理论计算认为,TSR温度在25 ℃以上就可以发生,但是并没有实际证明在100 ℃以下能够进行TSR[12]。Goldhaber和Orr[7]通过试验数据及现场数据的分析认为,100~140 ℃为TSR反应的启动温度。Machel[13]对烃类气体参与TSR的不同工况下的研究归纳得到,有机物成分、催化剂以及溶液中硫酸根的浓度等因素对TSR需要的最低温度范围有重要影响。
Orr等[14]提出SO42-和气体烃类在较高的温度下发生TSR,生成H2S,并总结了CaSO4和CH4的反应式。戴金星[15]认为,高温高压条件下硫酸盐与烃类等有机物能发生氧化还原反应,进而生成了大量的H2S,这与Orr研究结果相似。Worden和Smalley[16]通过研究碳酸盐油气藏中生成H2S的化学反应,结合标志元素硫、碳的同位素数据,给出了烃类与CaSO4的3个反应方程。
根据以上学者提出的TSR反应方程,可以总结出TSR的一般形式:
硫酸盐+烃类→CaCO3(CO2)+H2S+H2O+其他产物
(Ⅰ)
通过文献调研可知,国内外学者只针对TSR反应的影响因素进行了研究,未对具体反应路径进行详细阐述。因此,采用不含硫模型化合物正十六烷作为反应物与4种硫酸盐(MgSO4、Na2SO4、CaSO4、Al2(SO4)3)为研究对象,采用气相色谱、傅里叶红外管光谱仪、元素分析仪及水溶液离子分析仪等手段,通过实验探究TSR反应过程的关键步骤、明确反应过程的具体路径。由于硫酸盐是反应中最初存在的唯一的硫源,因此形成的H2S一定是由TSR产生的。这样排除了其他反应生成H2S对TSR产生H2S的干扰,可以更为清晰地探究TSR反应机理。此外,正十六烷是唯一碳源,也由此可以探明不同硫酸盐反应体系对CO2生成的贡献率大小。本实验结果可为稠油注汽热采过程中H2S的预防和治理提供依据。
实验仪器以及所需试剂,分别如表1和表2所列。
表1 主要实验仪器Table 1 Main experimental instruments仪器名称型号误差(最小刻度)产地在线式H2S浓度检测仪PN-2000±3%深圳FID/TCD气相色谱仪GC3800美国耐高温高压抗腐蚀哈氏合金钢反应釜江苏电子天平岛津AUY2200.0001g日本水溶液离子分析仪DIONEX-500离子色谱仪美国DIONEX傅里叶变换红外光谱仪美国NEXUSCHNS/O元素分析仪Elementar Vario EL III德国Elementar安捷伦7890A-SCD气相色谱仪7890A-SCD美国
表2 实验药品试剂Table 2 Experimental reagents药品名称规格厂家正十六烷分析纯上海无水MgSO4分析纯美国无水Na2SO4分析纯美国无水CaSO4分析纯美国无水Al2(SO4)3分析纯美国去离子水分析纯美国NEXUS高纯N299.999%(φ)青岛
实验装置如图1所示。
实验步骤:①称取特定药品,装入石英管。将石英管加入反应釜并反复冲充入1 MPa高纯N2,排尽空气并维持反应所需的高压条件;②当反应按设定的温度、时间完成后,收集气体并依次检测H2S及其他气体成分、收集液相产品检测其物质成分。
温度是影响TSR反应的主要控制因素。地层中硫酸盐的种类比较多,如石膏(CaSO4·2H2O)、明矾石(Al2(SO4)3)、芒硝(Na2SO4·2H2O)及苦土石(MgSO4·H2O)等。通过大量正交实验,确保硫酸盐含量不会对H2S产生影响,确定硫酸盐用量。经过预实验证明,在低于250 ℃时,反应进行程度较小,检测较为困难。因此,实验温度分别设定为250 ℃、275 ℃、300 ℃、325 ℃。将预先设定好的物料加入高温高压反应釜反应24 h,反应完成后,首先利用H2S检测仪检测H2S体积分数,之后对气相进行收集检测,结果见表3。
由表3可知,随着温度的升高,H2S生成量逐渐上升,证明反应随着温度的升高趋于剧烈。这是由于随着反应温度的升高,当温度达到反应物所需活化能时反应得以顺利进行,进而生成更多的H2S。从图2可看出,不同硫酸盐的反应体系生成的H2S量不同,即不同的硫酸盐TSR反应对H2S的贡献率不尽相
表3 正十六烷TSR反应气相产物分析Table 3 Analysis of TSR gas phase of n-hexadecane(mL·g -1(C16H34))名称MgSO4Al2(SO4)3Na2SO4CaSO4250 ℃275 ℃300 ℃325 ℃300 ℃300 ℃300 ℃CH43.925 25.913 6.116 0.808 4.238 4.745 17.995C2H61.614 49.250 67.633 73.485 27.976 48.089 83.386C2H41.243 41.560 39.882 32.359 13.265 21.748 25.546C3H80.000 0.000 0.000 0.000 0.000 4.066 0.000C3H645.919 0.000 0.000 0.000 0.000 0.000 0.000i-C4H107.009 5.478 27.240 18.213 37.324 34.122 33.973n-C4H100.000 2.288 1.167 0.622 7.183 19.256 0.837n-C4H83.078 5.531 2.559 0.000 9.015 8.523 2.641i-C5H121.275 3.640 4.729 6.752 14.214 20.315 1.473n-C5H120.000 0.008 2.761 0.008 3.544 0.003 0.000C+61.942 5.131 1.065 1.556 3.076 4.339 2.612H23.957 0.855 0.819 0.935 0.983 0.845 5.114CO0.000 0.000 3.087 11.993 6.204 0.000 0.389CO212.931 5.983 7.800 12.681 46.566 11.019 2.897烯烃50.241 47.091 42.441 32.359 22.280 30.271 28.187烷烃13.823 86.577 109.646 99.888 94.479 130.595 137.664H2S0.0050.0070.0140.0220.1480.0110.022
同,其顺序为:Al2(SO4)3>CaSO4>MgSO4>Na2SO4,这是由于Al3+所带电荷数较高,使SO42-电子云偏离较大,从而使得Al2(SO4)3反应体系中的活化能更低,生成的H2S量更多。由以上分析可以推测,在TSR反应体系中硫酸盐的自催化作用对H2S的生成起决定性作用。但在不同硫酸盐反应体系中,正十六烷对CO2的贡献率顺序为:Al2(SO4)3>Na2SO4>MgSO4>CaSO4。通过H2S和CO2生成量顺序可以看出,不同的反应体系其反应机理并不相同。
此外还可看出,无论是同一体系不同温度,还是同一温度不同体系,烯烃变化范围较小而烷烃的变化范围较大,证明TSR反应中,烯烃的反应活性比烷烃高得多。随着反应温度的升高,低碳分子的量逐渐增多,说明碳链较长的烃类比碳链较短的活性和还原性强,亦即随着烃类碳数的增加,碳氢化合物的氧化速率增加。由此可以证明,石油中长碳链碳氢化合物的活化能比短碳链的碳氢化合物要低。
采用安捷伦7890A-SCD气相色谱仪,对硫化物定性实验,得到硫化物标样的保留时间[17-20],如表4所列。
表4 标准硫化物的保留时间Table 4 Retention time of standard sulfide硫化物名称保留时间/min硫化物名称保留时间/min甲硫醇3.3092-甲基噻吩12.143乙硫醇3.8833-甲基噻吩12.443异丙硫醇4.613C4/C5噻吩17.200~19.400叔丁基硫醇5.366二乙基二硫醚19.790丙硫醇5.697苯丙噻吩19.987甲基乙基硫醚5.8862-甲基苯丙噻吩21.3712-丁硫醇7.4913-甲基苯丙噻吩22.615噻吩7.550正十六烷23.643异丁硫醇7.9104-甲基苯丙噻吩25.341二乙基硫醚8.8495-甲基苯丙噻吩26.341丁硫醇9.706C2苯丙噻吩30.320~38.170C2/C3噻吩10.870~16.800C3苯丙噻吩38.880~42.580二苯硫醚47.502C4/C5苯丙噻吩43.370~48.520
通过收集正十六烷与MgSO4反应体系的油相产物,采用安捷伦7890A-SCD色谱仪对其硫化物进行定性检测,结果如图3所示。
根据表4和图3可以得出,正十六烷与MgSO4反应体系油相产物中硫化物有甲硫醇、异丙硫醇、叔丁硫醇、噻吩/异丁硫醇、二乙基硫醚、C2噻吩、C3噻吩、C4噻吩、2-甲基苯丙噻吩、正十六烷、4-甲基苯丙噻吩、5-甲基苯丙噻吩、C2苯丙噻吩、C4苯丙噻吩、C5苯丙噻吩及二苯硫醚。从硫化物类型可知,正十六烷与MgSO4反应后生成的硫化物主要为噻吩类、少数的硫醇和硫醚类物质。
(Ⅱ)
(Ⅲ)
(Ⅳ)
MgSO4接触离子对的形成机制主要有3个步骤,具体如式(Ⅴ)所示。
(Ⅴ)
结合正十六烷的结构及反应后产物分析结果,可得到正十六烷与MgSO4的TSR反应体系的反应过程,如图4所示。
通过正十六烷的TSR反应过程可知,正十六烷参与TSR反应过程主要进行了质子化作用、热解反应、硫代硫酸盐向有机硫化物转化、H2S自催化作用及硫化物热解和水解等反应。这与Meshoulam等[24]的研究结论一致。
为了探究时间和温度对TSR生成H2S的影响,对正十六烷和MgSO4的TSR反应动力学进行了研究。从稠油及模型化合物TSR反应生成H2S的规律可以看出,随时间的增加,H2S生成速率整体呈现先增后减的变化规律,所以反应适用于如式(1)的动力学方程[25-26]。
(1)
式中:x为H2S的转化率,%;t为时间,d;n为总包反应的反应级数;k为反应速率,1/d;T为反应温度,K。
(2)
式中:A为指前因子,1/d;R为摩尔气体常量,J/(K·mol);E为表观活化能,J/mol。
对式(1)进行积分求解,可得:
(3)
式中:c为积分常数。
已知H2S生成率和时间,改变反应级数n,求得对应下的线性回归系数,根据最大的回归系数确定该反应体系的反应级数,得到不同反应条件下的反应速率k。将式(2)变换为:
(4)
根据不同的反应速率和反应温度,绘制回归曲线,得到斜率和截距,计算求得各反应体系生成H2S的动力学参数。热模拟实验的具体反应条件及实验结束后检测到的H2S生成含量见表5。
从图5可看出,正十六烷和MgSO4的反应级数为0.8。根据H2S生成率和时间的关系进行线性拟合,并求出其回归直线,如图6所示。
表5 正十六烷与MgSO4反应的H2S生成含量Table 5 Generation of H2S for the reaction between n-hexadecane and magnesium sulfate反应体系反应温度/℃反应时间/hH2S体积分数/%正十六烷+MgSO4+水250240.002 5 300240.005 1 350240.013 1 250480.004 6 300480.009 9 350480.027 8 250720.006 4 300720.013 9 350720.039 7
由拟合直线的斜率和截距求得正十六烷与MgSO4反应生成H2S的指前因子为0.585 7 1/d,表观活化能为61.498 kJ/mol。
(1) 由正十六烷的TSR反应产物研究可知,主要产物为MgO、CO2、H2S、烃类气体以及硫醇类、噻吩类和硫醚类有机硫化物。
(2) 4种金属盐TSR生成H2S量顺序为:Al2(SO4)3>CaSO4>MgSO4>Na2SO4;生成CO2量顺序为:Al2(SO4)3>Na2SO4>MgSO4>CaSO4。通过比较不同体系H2S和CO2生成顺序可知,金属阳离子电荷数越大,自催化作用越强;不同硫酸盐体系反应路径不同。
(3) 推导了正十六烷和MgSO4的TSR反应路径,包括质子化作用、热解反应、硫代硫酸盐向有机硫化物转化、H2S自催化作用及硫化物热解和水解等反应。
(4) 正十六烷和MgSO4反应体系的指前因子为0.585 7 1/d,反应活化能为61.498 kJ/mol。