爱情的文章
第38卷 第1期2018年01月
西安科技大学学报
JOURNALOFXI’ANUNIVERSITYOFSCIENCEANDTECHNOLOGY
Vol.38 No 1Jan 2018
DOI:10.13800/j.cnki.xakjdxxb.2018.0103文章编号:1672-9315(2018)01-0017-09
少量乙烯对甲烷爆炸危险性及反应历程的影响
收稿日期:2017-09-20 责任编辑:高 佳
苹果丢失基金项目:国家自然科学基金(51674193);国家重点研发计划(2017YFC0804T);陕西省自然科学基础研究(2017JM5068)通信作者:罗振敏(1976-),女,山东兖州人,教授,博士生导师,E mail:luozm@xust.edu.cn
罗振敏1,2,3,李 逵1,2,3,毛文龙1,2,3,张 江1,4,程方明1,2,3,王 涛1,2,3
,傅 文1
(1.西安科技大学安全科学与工程学院,陕西西安710054;2.西安科技大学西部矿井开采及灾害防治教育部重点实验室,
陕西西安7
10054;3.西安科技大学陕西省煤火灾害防治重点实验室,陕西西安710054;4.西安科技大学期刊中心,陕西西安710054)
摘 要:工业生产中可燃气体爆炸往往是多元可燃气体在空气中的爆炸,而乙烯作为活性非常高的可燃性气体,对甲烷爆炸的危险性有很大影响,为了进一步考察乙烯影响CH4爆炸的规律,采用标准的可燃气体爆炸极限测试装置,分析了含有不同浓度C2H4时,CH4在空气中的爆炸极限;并基于密度泛函理论(DFT),在B3LYP/6-31G水平下对相关基元反应进行定量分析,计算出相应热力学数据,然后用反应内禀坐标法验证反应路径,结合反应动力学来分析C2H4对甲烷爆炸反应历程的影响。结果表明:少量C2H4可使甲烷爆炸上限和下限均有所下降,且爆炸下限下降更明显,爆炸危险度F值增大;C2H4裂解过程中的·C2H3与O2
反应生成爆炸链反应的中间产物CH2O,增加了CH4爆炸链分支反应;C2H4的存在增加了链反应所需的中间产物,提高了整个反应体系活化中心浓度,促进了甲烷爆炸;计算结果从微观角度很好地解释了实验结果,得出了C2H4/CH4混合体系16步简化反应机理。关键词:安全科学与工程;C2H4/CH4;爆炸极限;爆炸危险度;密度泛函理论;热力学中图分类号:TD712.7 文献标志码:A
Influenceofminorethyleneonmethaneexplosion
hazardandreactioncourse
LUOZhen min1,2,3,LIKui1,2,3,MAOWen long1,2,3,ZHANGJiang1,4
,
CHENGFang ming1,2,3,WANGTao1,2,3,FUWen
1
(1.CollegeofSafetyScienceandEngineering,Xi’an
UniversityofScienceandTechnology,Xi’an710054,China;
2.KeyLaboratoryofWesternMineExploitationandHazardPrevention,MinistryofEducation,
Xi’anUniversityofScienceandTechnology,Xi’an710054,China;3.ShaanxiKeyLaboratoryofPreventionandControlofCoalFire,Xi’anUniversityofScienceandTechnology,Xi’an710054,China;
4.JournalPublishingCenter,Xi’anUniversityofScienceandTechnology,Xi’an710054,China)
Abstract:Inindustrialproduction,theflammablegaseshavegreateffectsontheexplosioncharacteris ticsofnaturalgascompnents.Ethyleneastheflammablegashasaveryhighactivity.Inordertofurtherinvestigatetheinfluenceofethyleneontheexpl
osionofmethane,theexplosionlimitsofmethaneinairwithethyleneatvaryingconcentrationsweretestedandanalyzedbasedonastandardtestdevicefortheexplosionlimitsofcombustiblegases,andrelatedelementaryreactionswasquantitativelyanalyzedbyapplyingtheB3LYP/6-31Gmethodsofdensityfunctionaltheory
,meanwhile,calculatedthecorre 博看网 . All Rights Rerved.
spondingthermodynamicdataandverifiedthereactionpathswithintheintrinsicreactioncoordinatemethod,andanalyzethemicroscopicmechanismofethyleneonmethaneexplosionbindingkinetics.Re suhsshowthatupperexplosionlimitandlowerexplosionlimitofmethanedecreaseandtheexplosivehazarddegreeFincreaseswhenmixingwithasmallamountofethylene,whilelowerlimitofexplosiondropssig
nificantly.C2H3inethylenecrackingprocessreactswithoxygentoproduceanintermediateproductformaldehydeoftheexplosionchain,whichincreasethemethaneexplosivechainbranchingre action.Increasingtheintermediateproductisrequiredforchainreactionintheexistenceofethylene,thus,improvetheconcentrationoftheactivationcenterofentirereactionsystemandpromotethemeth aneexplosion.Theexperimentalresultsarewellexplainedfromamicroperspective,obtainingthe16stepssimplifiedreactionmechanismofC2H4/CH4m
ixture.Keywords:safetyscienceandengineering;methane/ethylene;explosionlimit;explosivehazardde gree;densityfunctionaltheory;thermodynamics
0 引 言
以甲烷为主要成分的天然气作为一种清洁能源日益广泛地应用在工业和生活中,一直以来,由
含甲烷的可燃气体引发的爆炸事故时有发生,造成严重的财产损失和人员伤亡。乙烯作为可燃气体,具有爆炸极限小、压力峰值高和危险度大等特性,若与甲烷混合会增加混合气体的活性,使得甲烷爆炸极限范围明显增大,加剧甲烷爆炸,而造成更大的危害
[1]
。因此,开展C2H4/CH4
混合气体爆炸的研究对安全防爆具有一定指导意义。
目前,国内外学者主要对初始条件、点火能量、障碍物存在、惰性及部分可燃性气体等影响甲
烷爆炸的因素进行研究[2-7]
,在爆炸极限方面,陈
晓坤等采用近球型气体爆炸反应罐测试了在空气中C2H6对甲烷爆炸极限参数的影响[8]
;周利华对
矿井火区可燃性混合气体爆炸三角形判别法及其
爆炸危险性进行了研究
[9]
;VandenSchoor和Ver
plaetsen研究了在温度和压力升高时,空气中C2H6
和C3H8
爆炸极限的变化规律[10-11]
;Bolk等研究了流动状态下C2H6—Air—N2混合气体的爆炸极
限[
12]
。而涉及到乙烯对甲烷爆炸特性的研究较少,仅部分学者开展了考察,研究表明,少量的乙烯气体可使空气中甲烷爆炸的界限加宽,危险度
增大[1]
;C2H4
与CO之间存在着协同氧化诱导效应以及C2H4/CH4混合气体的爆炸极限临界氧浓度
降低,使甲烷爆炸危险性增大[13-14],这些研究也
只是在宏观上进行了定性分析,并未从微观角度深入研究。随着链式反应理论的提出和化学动力学的发展,发现可燃气体爆炸并不是单一化学反应,而是一种复杂的热爆炸和链式反应,自由基的
生成和消失直接影响反应的进行[
15]
。因此,如果从微观分子角度来定量地分析C2H4/CH4混合爆炸反应的机理,那么将会深入了解C2H4影响甲烷爆炸的规律。为此,文中基于气体爆炸极限测试实验,研究乙烯对甲烷爆炸极限的影响,同时应用Gaussian09程序,采用密度泛函理论的B3LYP/6-31G方法优化分子构型,并对相关基元反应作定量分析,从热力学和动力学角度,研究C2H4对甲烷链式爆炸的微观作用机理,得出C2H4/CH4混合体系步简化反应机理,从微观角度分析乙烯对甲烷爆炸的影响规律。
1 实验与计算方法
1.1 实验系统
实验系统采用标准的可燃气体爆炸极限测试装置,如图1所示,主要由爆炸反应管、点火器、循环泵、真空压力计、泄压口等组成。爆炸反应管长为14
00±50mm,管内径60±5mm,底部装置通径为25mm的泄压阀;配气、循环搅拌、点火等步骤均在计算机上自动控制完成,配气精确度为0 1%,采用电压为300VA电压互感器为点火电源,电极间距设置3~4mm,放电时间0.1~1s可调;实验在常压下进行,初始温度18~26℃,环境湿度56~64%RH,分别测试加入不同体积分数(0%,0 4%,0.8%,1.2%,1.6%,2.0%)的C2H4气体后CH4在空气中的爆炸极限,循环搅拌时间为5min,点火放电时间为500ms.1.2 计算方法
运用Gaussian09量子化学程序包,采用密度
泛函理论的B3LYP/6-31G方法[16]
,优化反应物
(
Re)、中间体(IM)、过渡态(TS)和产物(P)的几81 西安科技大学学报 2018年
博看网 . All Rights Rerved.
何构型,计算出相关基元反应的热力学数据如活化能、反应焓变、吉布斯自由能变,并用反应内禀坐标法验证反应路径[17]
,来考察C2H4
对甲烷链式爆炸的微观机理。
2 实验结果及计算分析
故人西辞黄鹤楼烟花三月下扬州2.1 C2H4对空气中CH4
爆炸极限的影响实验测得CH4/air中加入不同体积分数的C2H4气体时CH4的爆炸极限,且计算甲烷的爆炸危险度及爆炸极限范围如图2和图3所示,来研究C2H4对甲烷爆的影响规律。其中,可燃气体的爆炸危险度F是爆炸上下限的平均值和下限值的差
再除以平均值[18],其表达式F=(槡UL-L)/
槡
UL,U为爆炸上限,%;L为爆炸下限,%,F值范围为0~1,F值越接近于1则表明CH4爆炸危险性越大。
由图2可知,随着C2H4
浓度的增加,爆炸极限范围增大,爆炸下限下降比较明显。图3表明C2H4浓度增加,空气中CH4爆炸危险度F值逐渐增大,2%的C2H4使得CH4爆炸危险度增大约59 1%.表明,C2H4在一定程度上增大了CH4在空气中的爆炸危险性。
2.2 C2H4/CH4
链式爆炸反应过程计算及分析对于甲烷燃爆反应机理已有了很多相关报道,若某物质的存在能生成主要中间产物或优先于某些链反应发生,从而生成更多中间产物,那么就提高了反应体系活化中心浓度,促进了甲烷爆
炸[19-23]
。基于此,现借鉴敏感性分析甲烷燃爆反
应的简化机理[24]
,结合C2H4/CH4
混合体系中可能发生的反应,选取其中11个基元反应,从反应发生的难易程度和概率出发,进行反应机理的简化,来分析C2H4
对甲烷爆炸影响的微观作用机理。C2H→ 4·C2H3+·H(1)CH4+O→ 2·CH3+·HO2(2)C2H4+O→ 2·C2H3+·HO2(3)·C2H3+O→ 2CH2O+·HCO(4)CH2O+O→ 2·HO2+·HCO(5)·H+CH2→ O·O+·CH3(6)CH4+·→ H·CH3+H2(7)CH4+·→ O·CH3+·OH(8)C2H4+·→ O·C2H3+·OH(9)CH4+·→ OH·CH3+H2O(10)C2H4+·→ OH·C2H3+H2
O(11)
2.2.1 反应历程分析
图4为各反应中间体、过渡态及生成物构型,红色虚线表示原子间空间距离,并未实际成键。由图4可以看出,IM1的碳碳双键直接裂解到生成物P1,并不存在过渡态TS1,IM2中的一个氧原子与·OH达到成键距离后形成·HO2
东芝笔记本怎么样
;IM3中一个氢第1期罗振敏等:少量乙烯对甲烷爆炸危险性及反应历程的影响
9
1
图1 可燃气体爆炸极限测试装置Fig.1 Testdeviceforexplosionlimitsofcombustiblegases
*+(%"*$,-./0
*+(+"*$(),-./0
*#+(,-,-./0
大拇指戴戒指图2 C2H4对CH4
爆炸极限的影响Fig.2 Influenceofethyleneonexplosionlimitsofmethane
!"#
$%&$%!
#%'#%(
#%#
)*(
)!*( +,-./孕妇体重
0"10
0"(00"(20%200%220%&00%&2)*(
图3 C2H4对CH4爆炸危险度的影响Fig.3 Influenceofethyleneonexplosive
hazarddegreeofmethane
博看网 . All Rights Rerved.
原子发生断裂,与O2结合生成P3;IM4中碳碳双键断裂后与O2的2个氧原子分别结合形成P4,IM5碳氧双键发生断裂后一个氧原子与羟基结合形成P5,IM6中碳氧键发生断裂后氧原子被氢原子取代,形成P6,IM7中的一个碳氢键发生断裂后,2个游离的氢原子结合,形成P7;IM8中游离的氧原子攻击碳氢键,并与断裂后的氢原子结合形成P
8,IM9中同样是游离的氧原子攻击碳氢键,并与断裂后的氢原子结合形成P9,IM10中碳氢键发生断裂后的氢原子与游离的羟基结合形成P10,IM11中游离羟基攻击碳氢键,碳氢单键发生断裂后氢原子与羟基结合生成P11.通过图4所示的计算得知,
各中间体和生成物均存在。
图4 B3LYP/6-31G水平下各反应中间体、过渡态、及生成物的构型(键长单位/10-10
m,键角单位/(°))
Fig.4 Geometriesofintermediate,transitionstatesandproductsattheB3LYP/6-31Glevel
(Bondlengths/10-10
m,bondangles/(°))
2.2.2 热力学与动力学分析
11个基元反应的反应物、过渡态和产物以及正逆反应活化能、反应焓变和吉布斯自由能见表
1~4.其中,Ea为正反应活化能;Eb
为逆反应活化能;ΔHθ为反应焓变;Δ
Gθ
为反应吉布斯自由能变。在初始的CH4/C2H4
混合体系中,首先可能发02 西安科技大学学报 2018年
博看网 . All Rights Rerved.
第1期罗振敏等:少量乙烯对甲烷爆炸危险性及反应历程的影响生的反应只有反应1,2,3,从反应热力学角度出
发,由表1可知,3个反应的ΔGθ均大于0,因此常温常压下均不能自发进行。对比反应(1)C2H→ 4
·C2H3+·H和反应(3)C2H4+O→ 2·C2H3+
·HO2
的热力学数据可看出,2个反应的ΔHθ
都小于0属于放热反应,反应(3)的正反应活化能Ea为141.331kJ/mol,逆反应活
化能Eb为328.608kJ/mol,此反应Ea<Eb,反应有利于正向进行。反应(1)ΔE为28.959kJ/mol,不存在过渡态,虽然反应(3)有利于反应正向进行,但对比化学反应平衡,对于从反应物开始反应时的次序几乎无影响,而反应(
3)的反应活化能小于反应(1)的反应活化能,同样条件下,在此为了方便对比分析,简化描
述为率先发生是反应(1)。反应(2)CH4+O→ 2
·CH3+HO2
的焓变ΔHθ
大于0为吸热反应,正反应活化能Ea为289.645kJ/mol,逆反应活化能Eb为36.494kJ/mol,此反应Ea>Eb,反应有利于逆向进行。对比反应(2)和反应(3)的热力学数据可看出,反应(3)有利于反应正向进行,也就是说更有利于C2H4与O2结合;反应(3)为放热反应,同时反应(3)的正反应活化能Ea小于反应(2),反应(3)发生所需跨越的能垒较小,因此,同样条件下,反应(3)比反应(2)优先发生。故此处可以说反应(1)首先引发了C2H4/CH4混合体系爆炸链锁反应,为后续爆炸链反应提供热量以及关键自由基·H和·C2H3
.表1 反应(1)~(4)的热力学数据
Tab. Thermodynamicdataofreaction(1)~(4)
构型E/HartreeH/HartreeG/HartreeΔE/(kJ·mol-1)Δ
Hθ/(kJ·mol-1
)
Δ
Gθ
/(kJ·mol1
)
IM1
-78.53911
-78.53817
-78.56435
---(1)C2H→ 4·C2H3+
·HTS1---28.959-375.420
44.292P1-78.52808-78.68116-78.54748-
-
-
构型
E/Hartree
H/Hartree
G/Hartree
ΔEa/(kJ·mol-1)ΔEb/(kJ·mol-1
)
半天ΔHθ/(kJ·mol-1)Δ
Gθ
/(kJ·mol-1)IM2-190.73959-190.73411-190.76548
----(2)CH4+O→ 2·CH3+·HO2
TS2-190.62927-190.62291-190.65737289.645
36.494259.268
210.670
P2
-190.64317-190.63536-190.68524
----IM3-228.79956-228.79789-228.83691
----(3)C2H4+O→ 2·C2H3+·HO2
TS3-228.74573-228.74041-228.77437141.331
328.608
-212.797
197.543
P3
-228.87089-228.87894-228.91215
----IM4-228.22651-228.22557-228.25668
---
-
(4)·C2H3+O→ 2CH2
O+·HCOTS4-228.21808-228.21844-228.2420522.133199.853
-177.720-184.573
P4
-228.29420-228.29326-228.32698警察叔叔图片
-
-
-
-
由于反应1比反应3的优先发生,可以使体系中的C2H4大部分反应生成·C2H3,因此在后续的链式反应中,体系中与O2的反应主要体现为反应2,4.由表1可知,链反应(4)·C2H3+
O→ 2·CH2O+HCO正反应活化能Ea为22.133kJ/mol远小于反应2的289.645kJ/mol,逆反应活化能Eb为
199.853kJ/mol,焓变ΔHθ
为-177.720kJ/mol,为放热反应,Δ
Gθ<0,此反应在常温常压下能自发进行,同时,此反应Ea<Eb,反应有利于正向进行。比较反应(2)的热力学数据可以发现,反应(4)比反应(2)更容易发生,也就是说在反应(1)首先引发初始混合体系链锁反应后,·C2H3
可以优先或代替CH4去与O2发生反应,同时为体系提供热量和链分支反应(5)的反应物CH2O以及爆炸过程关键自由基·HO2的生成,促进了CH4爆炸的链引发和链分支反应。因此,C2H4的加入使得CH4在空气中更易发生爆炸反应,爆炸危险性将增大,这也解释了C2H4的加入使得空气中CH4爆炸下限出现明显下降和爆炸危险度F值的增大。
由表2可知,反应(5)CH2O+O→ 2·HO2+·HCO正反应活化能Ea为193.184kJ/mol远小于反应(2)的289.645kJ/mol,逆反应活化能Eb为
50 515kJ/mol,ΔHθ
为141.646kJ/mol,是吸热反应,Δ
Gθ
>0,此反应在常温常压下不能自发进行,1
2 博看网 . All Rights Rerved.