2,4-二硝基甲苯热解自催化特性鉴别及其热解动力学

物理化学学报( Wuli huaxue xuebao)Acta Phys. -Chim. Sin. 2013, 29(3), 479-485479[Article]doi:10.3866PKU.WHXB201212141www.whxb.pku.edu.cn2,4-二硝基甲苯热解自催化特性鉴别及其热解动力学鲍士龙陈网桦陈利平高海素吕家育(南京理工大学化工学院安全工程系,南京210094摘要:为研究2,4-二硝基甲苯(2,4-DNT)的热危险性及其分解反应的特征,利用差示扫描量热仪(DSC)对该物质进行了动态扫描测试,得到其起始分解温度T范围为272.4-303.5°C,分解热ΔH约为2.22kJ·g.在此基础上,采用瑞士安全技术与保障研究所提出的快速鉴别法(瑞士方法)及数值模拟技术,对其分解反应的特性参数进行了推算,结果表明其分解具有自催化性.采用 Malek法分析了该物质分解反应的最概然机理函数并得出了相关动力学参数,表明其分解具有自催化性且符合 Sestak-Berggren双参数自催化模型(SB模型,这与瑞士方法所得结论一致.采用等温DSC测试获得了该物质的‘钟形’热解曲线,从而验证了两种方法的结论关键词:自催化分解反应;差示扫描量热仪;快速鉴别; Malek法;分解反应动力学中图分类号:O642;0643Identification and thermokinetics of autocatalytic ExothermicDecomposition of 2, 4-DinitrotolueneBAO Shi-Long CHEN Wang- Hua CHEN Li-PilGAO Hai-SLU Jia-Yu(Department of Safety Engineering, Institute of Chemical Technology, Nanjing University of Technology and Science,Nanjing 210094, P. R. ChinaAbstract: The thermal stability and autocatalytic decomposition of 2, 4-dinitrotoluene(2, 4-DNT)isinvestigated under dynamic and isothermal conditions using differential scanning calorimetry(DSC).Thetemperature range of initial exothermic temperature(To)is 272.4-3035C, and its decomposition enthalpyAHa)is about 2.22 KJ. g. An identification method based on a numerical simulation technique from theSwiss Institute for the Promotion of Safety and Security(Swiss method) is used to determine thecharacteristic parameters of the decomposition reaction, revealing that the decomposition of 2, 4-DNTpotentially autocatalytic. The Malek method is used to determine the most probable mechanism functionand kinetic parameters of 2, 4-DNT decomposition. The Sestak-Berggren model with two parameters issuitable to describe the autocatalytic decomposition of 2, 4-DNT, which is consistent with the results of theSwiss method and isothermal experimental results. Isothermal DSC experiments confirmed that thedecomposition of 2, 4-DNT is autocatKey Words: Autocatalytic decomposition reaction; Differential scanning calorimeter; Rapididentification; Malek method; Decomposition kinetics1引很长一段时间内都不会出现明显的温升现象,而自催化分解反应是指分解产物在反应过程中旦产生明显温升,则迅速发展为热爆炸.这是导致充当催化剂的反应.与一般的具有n级动力学分解具有自反应性的化学物质(如爆炸性物质、有机过氧反应不同,自催化分解反应在发生热自燃热爆炸前化物、自反应性物质等)易发生热自燃与热爆炸的重Received: September 28, 2012: Revised: December 13, 2012: PublislC YHCnMHGThe project was supported by the National Natural Science Foundation of China (51204099)国家自然科学基金(51204099)资助项目C Editorial office of Acta Physico-Chimica SinicaActa Phys. -Chim. Sin. 2013VoL 29要原因.在如今危险化学品安全生产状况虽然持续类似2,4二硝基甲苯这样的含能化学品分解反好转,但是在形势依然严峻的大背景下,为了正确应的特性描述对其生产工艺的研发及储运是至关识别、评价其潜在危险性,并在此基础上针对性地重要的,而在研究其分解反应方面DSC被广泛应采取防控措施,从而确保此类化学物质在储运及使用.2本文从实验和理论分析两方面对2,4DNT分解用等过程中的安全,必须对其反应机理及特性、鉴反应的特征及动力学参数进行了研究.在此基础别方法及表征等进行深入研究上,考虑到等温测试是对自催化分解反应进行实验国外学者对自催化分解反应进行了大量研究,验证的最可靠方法,对该物质在等温条件下热行为并取得了许多研究成果,提出了三种描述自催化反进行了测试,并将测试结果与理论分析进行了印证应的模型,分别是 Prout- Tompkins模型(简称PT模比较.型), Benito-Perez模型(简称BP模型)和源自柏林学院的柏林模型4虽然近二十年来,国内(含台湾地区)2仪器与测试条件的热分析研究得到了迅猛发展,但总体而言,在物实验样品:2,4-DNT(纯度99%,国药(沃凯)试质物料自催化分解反应的热测试与表征方面的研剂)设备:DSC-1专业型,瑞士 METTLER TOLEDO究开展较少.上世纪90年代,国内学者先后研究讨公司生产论了JH9105炸药与奥克托今炸药(HMX)的自催化实验条件:不锈钢密封坩埚,实验中测试气氛特性,但这些研究所没有得出相关的动力学参数及维持不变,分别以50和100mL·s的N2为反应气和模型362000年以后, Zhang等'利用非等温差示扫描干燥气.动态实验中升温速率分别为1、2、4、8°C量热仪(DSC)研究了硝化纤维素在自催化模型下的min,样品量依次为1.00、1.08、1.48、1.13mg,温度热爆炸临界升温速率,之后胡荣祖与其合作者3范围为25-380°C.等温实验中样品量为1.43mg,应用非等温DSC进行数学分析得出了硝化棉的临等温实验温度为230°C界升温速率和自催化分解反应的动力学方程,赵凤起等"针对3,4-二(4-硝基呋喃氨基)2-氧呋喃进行3实验结果与讨论了研究,得岀了其自催化分解的热爆炸临界升温速3.1动态DSC曲线分析率,高红旭等则分析并计算了N-脒基脲二硝酰胺图1是不同升温速率下2,4DNT分解的DSC曲自催化分解情况下的动力学参数及相关的热稳定线,相关数据列于表1中.图2为不同升温速率下的性参数.文献ˉ对几种含能材料的自催化分解动2,4-DNT转化率α与温度的关系曲线.力学进行了分析并推算了相关的安全参数,但都没从图1可以看出,2,4-DNT仅有一个放热峰.起有给出自催化分解反应的鉴别方法.台湾的欧新荣始分解温度和峰值温度随升温速率的增大而升高.等研究了丙烯酸类的自催化特性及安全储存条其分解热的平均值为22kJ·g,大于800J·g-,根件,采用复杂线性回归计算对自催化反应的动力学据文献的判别条件,该分解反应一旦发生,引发事参数进行拟合,并推算了各种安全参数.台湾的Chen等采用细致的多级模型代替了经典的两级自16B(°cmin)催化模型,并进行了一系列的模拟计算.此外,鲜有文献涉及自催化热分解反应研究.2,4-二硝基甲苯(简称2,4DNT)是一种重要的有机合成中间体,在农药、医药、火炸药等相关领域应用广泛.国内外对其研究很多,但主要集中于其动力学参数的计算(非热分析方法)、生物或金属离子对2,4DNT催化降解作用、2,4DNT在自然界中的污染扩散等方面.5虽然文献∞对其热分解过200220240260280300320340360380400程及稳定性进行了研究,但均按一般n级反应模式中国煤化工进行处理,笔者对此存在争议,这也是本文研究的rig1 Dynam出发点之CNMHG24DNTNo. 3鲍士龙等:24-二硝基甲苯热解自催化特性鉴别及其热解动力学481表12,4-DNT热解反应动态DSC实验数据别.具体来说,动态DSC测试时满足以下条件Table 1 Experimental data of the exothermicition of 2. 4-dnt determined by dscda=Aexp( rt)f(a)g=△H出=△H1合-1)o276.62.25305.7其中,βC·min)为升温速率,a为转化率(含义与图31702中一致),A(s-)为指前因子,R为理想气体常数T: initial exothermic temperature: T:: thermal decomposition peakE(kJ·mol)为表观活化能,qW·kg)是放热速率temperature; AHd: decomposition enthalpy△HJ·g)为反应的分解热1.00.9 B/(Cmin)取DSC曲线刚刚偏离基线时的温度TC)为起始分解温度,φ(W·kg)是起始分解温度时对应的放热速率.在T时a非常小(可以认为接近0,f(a)约等一-8于1,则有式(3qo=△ Haoexp(将式(3代入式(2)得到式(4q=90exp(-E(1-1)f(a)220240260280300320340360380当a较小时,fa)可以用1-a来近似,即转化率很低时图22,4DNT转化率(a)-温度关系曲线自催化热解反应动力学可以用一级反应模型,于是Fig2 Relation curves between fraction conversion(a)可以采用式(5)加以描述and temperature of 2, 4-DNTq= go exp(-(-1)(1-a)故的严重度很大.分解反应的转化率与反应时间是相互对应的,又由于温度与时间之间存在线性关DSC曲线的q值与T值是一一对应的.用式(5)描述系,因此转化率与温度存在一一对应的关系.图2为条DSC曲线的起始阶段时,曲线的T与q为固定转化率-温度关系曲线,四条曲线形状及变化趋势值,若q取该曲线上某一实测点的数值,则T也唯基本一致,表明四次测试的反应历程基本一致确定(相应的a也为固定值.此时对于该实测点T、a32自催化热分解反应的鉴别T、φ与R都为定值,只有E为变量.从式(5)可以看3.21瑞士方法的原理出E与q成正比,q值越大则E值越大(式(5)中T小般说来,物质在动态测试下发生自催化热分于等于T,1/小于等于1/).对式(5)求导得到式(6)解的放热曲线有一定的特征,例如峰形较为对称,dr-yo exp(ErI)(1-aE形状较为尖锐而狭窄.最简单的思路是通过对峰的形状特征如峰高、峰宽及其比值等进行定量描述来从式(6)可以看出E与dqdT也成正比,q值增长越快判断其是否为自催化反应.但该方法仅用到了很少则E值越大反应自催化作用越强,反应开始时相比的数据点统计学意义不大,而且受到测试设备、测于一般的n级反应其q值越大且增长速度越快,在试条件的影响很大.瑞士安全技术与保障研究所提式(5)和式(6)中表现为表观活化能也越大因此可以出了一种更为可靠有效、且有较高统计学意义的方采用一级动力学模型,通过调整表观活化能(反映了法2——基于简单一级反应模型的自催化反应鉴别放热曲线的上升陡度)反过来对动态DSC实验获方法(简称瑞士方法):在动态测试条件下,初期出现取的放热峰的开始阶段(本文取a为0001-0.1范围的放热曲线反映了绝热条件下反应物料的行为且内)的数据进行拟合,所得表观活化能越高,自催化自催化分解反应初期,转化率很低,其反应动力学可能性越大可以用一级反应模型加以描述.3因此可利用一级瑞士方法的模型源干式(5)表达式见反应模型对分解反应初期的DSC曲线进行分析,获式(7)中国煤化工取反应的表观活化能,并以此对自催化反应进行鉴CNMHG482Acta Phys. -Chim. Sin. 2013VoL 29(+△)=qcA-R(7(+△)-7)偏离度(a)的定义为△H(t+△n)Psi-pei(10)△Ha方程中T,RE,△H与前面意义相同△H(+△)表e的单位为mW,式(0)中n为实际测试中a范围为示t+Δt时刻的分解产生的热量,该值无法实际测0.001-0.1时所取数据点的个数(具体个数由具体实得,可通过式(8)进行估算验决定).p和p为同一温度下放热速率q的模拟值与实验值.偏离度的实际含义为其在同一温度下放△H(t+△)=△H(t)+q(△t将式(8)带入式(7)得到式(9)热速率q的模拟值与实验值的平均差值.E拟合计算时,先设定E的范围(本文取50-500q(+△)= go exp(RT(t+△)T0kJ·mol-),控制时间步长为1s,调节E值以0.1kJmol为步长向上递增,对实测曲线起始阶段进行拟(1H合,得到偏离度最小的E值,当其数值大于220kJ模型中用到了两个实测值,一是T时所对应的放热mol时推测认为该物质具有自催化性速率q,另一个为分解热△H式(8)中,△H(n)表示到t32.2瑞士方法结果时刻反应的放热量,其迭代计算式为式(8),(△H(1)+结合动态DSC实验数据,进行数值计算与拟合q(1)△△H表示转化率a7(+△)表示+△时刻的温得到的结果如图3.这里需要特别说明的是,此处拟度.由(9)可以看出该方程中有E、t、T、△H(1)、q(1)五合计算得到的活化能并没有任何实际的理化意义个变量,温度T是时间t的线性函数,△H(1)也由时间只是一个数值指标决定,q(1)为迭代计算中上一次计算的结果即时刻t对动态升温速率为1、2、4、8°C·min的2,4-时的放热速率,因此对于固定的E,计算可以用t来DNT放热曲线进行拟合,所得到的表观活化能E分控制进程别为3585、281.8、447、3574kJ·mol,对应的偏本文通过计算机程序得出了实验曲线与拟合离度为0.020、0.22、0.020、0.040mW,其中升温速曲线之间偏离度最小时所对应的E值率为2°C·min的曲线拟合结果偏差较大,分析后(a)B=1°cminexperimental curvexperimental curvefitted curve0.0r210240270300330360210240270300330360390T°cexperimental curve 16(d)B=8C-min-'experimental curve-· fitted curve· fitted curveE≥108642021024027030033036210240中国煤化工图3不同升温速率下2,4DNT的拟合曲线与实验CNMHGFig 3 Fitted and experimental curves of 2, 4-dNT at different ncauig iNo. 3鲍士龙等:24-二硝基甲苯热解自催化特性鉴别及其热解动力学483认为主要是该实验在开始阶段升温较为剧烈.拟合表2应用 Ozawa法求得2,4-DNT的活化能值所得E均大于220kJ·mol,说明2,4-DNT分解反应Table 2 Values of E for 2, 4-DNT by Ozawa method具有潜在自催化性aE/(kJ·mol)aE(kJ·mol)aE/(kJ·mol)0.10140.70.4014260.7013703.3 Malek法推断最概然机理函数及动力学参数0.15140.70,4514260.75132.6331Maek法介绍152.9137.0为深入了解2,4DNT分解反应机理,结合动态132.6142.6DSC数据应用具有较好普适性的 Malek法对0.35140.7132.64DNT的热分解动力学进行了研究. Malek法是由E: activation energy; E(average)=138.9 kJ. me定义的归一化函数y(a)和Z(a)确定机理函数f(a)和通过所得曲线可以得出特征值aM和a,其中aMG(a)的一种较好的方法.美国马里兰大学 Criado教授、塞尔维亚 Jankovic教授以及国际热分析学会均和a分别为y(ax)-a和Za-a曲线峰值处对应的转化认为Mak法可以明显缩小动力学机理函数的可能率,记动态DSC曲线峰值处转化率为a利用这三个范围,该方法几乎不存在假设和近似条件,可以较特征值就可以推断出最概然机理函数fa)由拟合好地提高机理函数形式研究的准确性.确定机理曲线所得的特征值见表3函数fa)和G(a)的表达式需要用到定义函数y(a)和由表3可知a