多环芳烃在超临界环己烷中的溶解：分子动力学模拟研究.pdf

第41卷第2 期2024年4月J.At.Mol.Phys.,2024,41:021003(9pp)原子与分子物理学报JOURNAL OF ATOMIC AND MOLECULAR PHYSICS多环芳烃在超临界环已烷中的溶解：分子动力学模拟研究Vol.41 No.2Apr.2024赵向勃，付峰，陈权，高晓明，闫挺1（1.延安大学化学与化工学院，延安7 16 0 0 0；2.厦门海关技术中心，厦门36 10 2 6）摘要：重油中大量存在的多环芳烃会经过缩合反应会形成稠环芳烃，而溶解扩散的限制会导致稠环芳烃之间的聚合从形成残焦，使得对重油的利用率大大降低，所以多环芳烃热解过程的溶解扩散是重油改质的关键因素之一因此，基于环己烷对烃类良好的溶解性，本文采用分子动力学模拟的方法研究了多环芳烃及其混合物在超临界环己烷中的溶解行为，结果表明在不同温度与密度下超临界环已烷对于具有NAP均具有良好的溶解性，而温度对于Bghip溶解的影响较小，在多环芳烃混合物的油滴溶解过程中,NAP优先溶解到环己烷相中，而Bghip则集中在粒径减小的油滴中：通过计算体系中多环芳烃的径向分配函数与溶剂化自由能发现环己烷与多环芳烃间之间能相互吸引，环烷在多环芳烃分子周围形成的溶剂壳能抑制多环芳烃分子间的聚合：温度升高以及密度的将降低会导致多环芳烃的内聚能密度降低，增加了多环芳烃与超临界环己烷之间的相互作用：环己烷密度越小以及温度的升高可以促进多环芳烃或多环芳烃混合物在超临界环己烷中的溶解。关键词：超临界环己烷；多环芳烃；溶解；分子动力学中图分类号：TQ021.4Dissolution of polycyclic aromatic hydrocarbons in super critical cyclohexane:文献标识码：AA molecular dynamics simulation studyD0I:10.19855/j.1000-0364.2024.021003ZHAO Xiang-Bo,FU Feng,CHEN Quan,GAO Xiao-Ming,YAN Ting1.College of Chemistry and Chemical Engineering,Yan an University,Yanan 716000,China;2.Techndogical Centre,Xiamen Customs District,Xiamen 361026,ChinaAbstract:A large number of polycyclic aromatic hydrocarbons in heavy oil will undergo a condensation reactionto form polycyclic aromatic hydrocarbons(PAHs),and the limitation of dissolution and diffusion will cause thepolymerization between polycyclic aromatic hydrocarbons to form residual coke,which greatly reduces the utiliza-tion rate of heavy oil.Therefore,dissolution and diffusion in the pyrolysis process of PAHS is one of the key fac-tors for the improvement of heavy oil.Based on the good solubility of cyclohexane to hydrocarbons,the moleculardynamics simulation method was used to study the dissolution behavior of polycyclic aromatic hydrocarbons andtheir mixtures in supercritical cyclohexane.Critical cyclohexane has good solubility for NAP,while temperaturehas little effect on the dissolution of Bghip.During the dissolution of oil droplets in the PAH mixture,NAP waspreferentially dissolved into the cyclohexane phase,while Bghip was concentrated in the oil droplets with reducedparticle size.By calculating the radial distribution function and solvation free energy of PAHs in the system,it isfound that cyclohexane and PAHs can attract each other,and the solvent shell formed by cyclohexane aroundPAH molecules can inhibit of intermolecular polymerization PAHs.The increase in temperature and the decrease收稿日期：2 0 2 2-0 7-30基金项目：国家自然科学基金（2 2 0 6 8 0 37)；陕西省自然科学基金（2 0 2 1JQ-615）作者简介：赵向勃（1999一），男，硕士研究生，E-mail：Zh a o XB0 996 16 3.c o m通讯作者：闫挺.E-mail：9992 7 7 7 16 3.c o m021003-1第41卷in density will lead to a decrease in the cohesive energy density of PAHs,increasing the interaction betweenPAHs and supercritical cyclohexane.The lower density of cyclohexane and the increase in temperature can pro-mote the dissolution of polycyclic aromatic hydrocarbons or mixtures of polycyclic aromatic hydrocarbons in super-critical cyclohexane.Key words:Supercritical cyclohexane;Polycyclic aromatic hydrocarbons;Dissolution;Molecular dynamic1 引 言近年来，尽管太阳能、氢能、风能等新能源发展迅速，但传统化工原料仍然在人类的社会发展中发挥着重要的作用：原油经过长期的过度开采及使用，轻质油含量在逐渐减小而重油和超重油储量已占据我国现有石油储量的50%以上1重油具有粘度高、杂原子含量高的特点，传统的加氢脱碳工艺在重油加工中存在诸多局限性2-4脱碳过程会产生大量残焦，造成资源浪费5而重油中的S、N等杂原子在加氢过程中会使催化剂失活超临界流体(SCFs)技术是一种新型的加氢脱碳工艺，可以在不使用催化剂的情况下使得重油改质，抑制焦炭的生成6 超临界流体已经广泛用于重质油轻质化、减黏等方面的研究7-11超临界水(SCW）具有低的介电常数和良好的扩散性能，使得重油在SCW中有着良好的溶解性能，再结合潜在的供氢能力，重油在超临界流体中的轻质化与减黏主要集中在SCW领域中12-14然而，近年来研究发现SCW在重油热解的过程中的供氢效果并不明显，SCW供氢能力有待深人了解15。同时，Zhu等16,17 基于量子力学的理论计算否定了SCW在涉及碳氢自由基的反应中给予氢自由基的可能性，所以越来越多的人认为SCW在重油热解过程中主要作为溶剂发挥了物理作用，然而，水和烃类的混合体系表现出复杂的相行为对其复杂相行为的检测和预测一直是研究的重点18 目前，水和众多低分子量碳氢化合物的相行为已经有了较大的研究19然而，重烃混合物的研究较少，因为由于混合物的不透明、乳化现象均可能干扰相边界测量此外，烃类高温高压下的热解过程采样困难，对这类混合物的相行为测量具有挑战性在之前的研究中，根据Van Konynenburg 和 Scott 分类方案2 0，Athabasca沥青+水的伪二元混合物表现出b 型相行为2 1.2，并测量了沥青在水中的溶解度结果表明在含水的沥青混合物中加入溶剂时，相对于水原子与分子物理学报和沥青，溶剂通常具有较低的密度和较低的临界温度，相对于沥青具有较高的水溶性：预测富烃液相中的混合相行为和水含量会发生显著变化.随着对重油一超临界流体体系相结构、反应机理和反应动力学的研究，研究者认识到不同的相结构及其形成的扩散环境对超临界流体中重油的非催化热解起着至关重要的作用：Morimoto等2 3基于汉森溶解度参数（HSP）提出了SCW的静态介电常数大于2.2，HSP氢键组分小于10MPas是SCW与重油良好混相的必要条件He等2 4,2 5利用基团贡献法，分析了2 7 36 44 K、26 MPa压力下原油和水混合物的相平衡，在压力为2 5MPa、温度为6 40 6 6 0 K的条件下，研究了软沥青、沥青质和水的混合物的相平衡，揭示了单相和两相区域的存在而 Guo 等2 6 通过混相指数解释了SCW中重油的改质，发现混相指数的降低导致软沥青产量增加的同时沥青质和残焦的收率减少。在恒定的混相指数下，温度增加可以导致软沥青产率的增加但是残焦的产率也会有所提升，即在不结焦的高温条件下最大限度的降低混相指数可以提高软沥青的产率。所以重油一超临界流体体系的混相可以揭示重油超临界流体体系中的相结构和溶解行为分子动力学模拟（M D)可以研究超临界状态下的复杂系统2 7.2 8。Jian等2 9通过MD模拟研究了SCW的结构和扩散性能从室温到超临界温度的升高对水中氢键的影响比沿着超临界等温线密度的变化要明显得多.计算证实，随着温度的升高和密度的降低，每个分子的平均氢键数和团簇大小都在减少，Xin等30 通过分子模拟研究了沥青质在超临界流体中的溶剂化反应，结果表明排斥Van der Waals 相互作用对沥青质与水分子相互作用的性质和程度有主要影响：SCW是重油改质中常用的流体，但由于多环芳烃与 SCW 的相互作用表现为斥力，不利于多环芳烃的溶解而超临界环已烷与多环芳烃之间的相互作用为引力因此，本文利用分子动力学方法研究了不同环数的多环芳烃以及多环芳烃混合021003-2第2 期第41卷物在超临界环已烷中的溶解过程，揭示了多环芳烃在超临界环已烷中的溶解机理：首先采用MD模拟方法研究了纯多环芳烃或多环芳烃混合物在超临界环已烷中的溶解行为，考察了多环芳烃的尺度、多环芳烃混合物的组成对溶解过程的影响.通过对多环芳烃在超临界环已烷中的短程溶剂结构和溶剂化自由能的研究，表征了多环芳烃与超临界环已烷环境的相互作用：然后，通过计算多环芳烃间的内聚能密度来揭示了多环芳烃本身之间的相互作用本文将重油在超临界流体中改质中超临界流体的筛选，操作范围的确定提供理论指导.2模拟与理论方法首先MD仿真的可靠性取决于所用的力场，COMPASS力场在烃类/SCW体系中的有效性已被证实31,32，因此本文所用力场采用COMPASS力场应用于多环芳烃与超临界环已烷体系，所采用的计算均基于Materials Studio2019软件的Forcite模块，由于重油当中所含多环芳烃的结构变化较多，因此本文只考虑没有杂环掺人和烷基取代基存在的多环芳烃在本次研究中，选择2 环的萘（NA P)为轻多环芳烃，选择6 环的苯并芘（Bghip）为重多环芳烃如图1所示这些多环芳烃的详细信息如表1所示。图1多环芳烃分子结构Fig.1Molecular structures of different PAHs.在模拟中，首先构建充满超临界环已烷分子的立方体溶剂盒来模拟超临界环已烷环境，如图2所示，环境温度为6 346 7 4K，密度为0.1g/cm0.3g/c m，当模拟多环芳烃的溶解时，将含有纯多环芳烃或多环芳烃混合物的油滴放在溶剂盒的中心在研究多环芳烃的近程溶剂结构和溶剂化自由能时，将一个多环芳烃分子放人随机溶剂盒中，溶剂盒的剩余空间装满了特定状态下的环已烷分子。多环芳烃在超临界环已烷中溶解的典型动态赵向勃，等：多环芳烃在超临界环已烷中的溶解：分子动力学模拟研究分子量茶NAP苯并花Bghip3丝结果与讨论3.1多环芳烃在超临界环已烷中的溶解3.1.1溶剂密度对多环芳烃溶解度的影响在6 7 4K的温度下，采用MD模拟了NAP和Bghip的多环芳烃液滴在0.10=0.30 g/cm三种超临界环已烷密度下的溶解情况模拟时间在500ps时的结果如表2 所示.温度6 7 4K下不同密度的压力在0.0 0 5GPa0.0 0 8 C Pa 之间，根据观察，在不同密度的环已烷中，NAP油滴分子能够迅速地均匀分散在溶剂盒中随着环已烷密度的增大，多环芳烃与环已烷的体系的混相指数逐渐减小：同时，环已烷密度的变化对NAP和Bghip的溶解有重要影响021003-3第2 期女图2 超临界环已烷溶剂盒Fig.2 Schematic diagram ofSC-cyclohexanecube box.模拟过程如下所述：首先，构建纯多环芳烃的本体结构，对溶剂盒进行结构优化并通过NPT系综的退火操作获得了纯多环芳烃的稳定本体结构。从多环芳烃本体结构中取出一个半径为15A的油滴分子，将其放置在一个立方体溶剂盒的中心，立方体溶剂盒边长为8 0 A.此外，溶剂盒中充满了环已烷分子，其密度由环已烷分子的数量来调节环已烷密度从0.1g/cm到0.3g/cm不等.然后，基于NVT系综进行持续时间为50 0 ps的分子动力学模拟，以保证动力学平衡在模拟过程中，温度由NHL恒温器调节静电法和范德华求和法分别设置为vanderWaalsbased和atombased.表1多环芳烃模型性质Table 1Model properties of PAHs名称简写熔点K128.1705353276.3307553沸点K490第41卷当环已烷密度为0.10 g/cm时，NAP分子快速溶解于环烷相中，球形油滴形态严重偏离最初的形态增加环已烷密度至0.30 g/cm时，油滴扩散的效果较好，但是相比于0.10 g/cm的超临界环已烷体系仍然存在NAP分子的聚集与此同时，对于Bghip，在0.10 g/cm时，Bghip分子迅速溶解于超临界环已烷中。但是相比于NAP分子溶解性稍差，在0.30 g/cm模拟结束时，溶剂盒中只发现少量解离的Bghip分子。从图中可以看出NAP油滴可在超临界环已烷中完全分离，但仔细观察表2 的图像可以发现，在溶剂箱中仍存在NAP分子的区域富集，表明NAP发生了团聚由于富集区域的位置在整个模拟过程中随机变化，因此在超临界环已烷中形成的NAP团聚体是不稳定的而当重多环芳烃溶解到环已烷相时，Bghip在溶剂盒中出现稳定的团聚，其主要以二聚体、三聚体等类似于沥青质 一T堆积的形式存在。表2 多环芳烃在不同密度超临界环已烷中的溶解情况Table 2Dissolutions of PAH droplets in SC-cyclohexanewith different densities at temperature of 674 K0.1g/cm30.3g/cm3NAP梦车Bghip3.1.2温度对多环芳烃溶解度的影响在温度为6 346 7 4K、超临界环已烷密度为0.10-0.30 g/cm时，采用MD模拟方法研究了NAP和Bghip液滴在超临界环已烷中的溶解情况，模拟时间50 0 ps情况下的结果如表3和表4所示当环已烷密度为0.30 g/cm时，温度升高20/40K对NAP液滴溶解的影响有轻微的提高.NAP分子从油滴中溶解在环已烷相中的数量基本相同在环已烷密度为0.10 g/cm时，在6 34和654K温度下也能观察到类似的现象，但随着温度的进一步升高，情况发生了较大的变化在674K温度下，大量的NAP分子会以单体的形式溶解到环已烷相中。当环已烷密度为0.30 g/cm时，温度升高原子与分子物理学报表3萘在不同温度下的溶解行为Table 3Dissolutions of NAP droplets in SC-cyclohexane atdifferent temperatures634 K0.1g/cm30.2g/cm30.3g/cm320/40K对Bghip液滴的溶解与NAP分子同样具有明显的影响大量的Bghip分子从油滴中解离并溶解在环已烷相中，在环已烷密度为0.10 g/cm时，在6 34K和6 54K温度下也能观察到类似的现象，但随着温度的进一步升高，发生了的变化在0.2g/cm3674K温度下，大量的Bghip分子以单体、二聚体或更大规模的聚集体的形式溶解到环已烷相中，此时Bghip液滴的剩余部分可以看作是一个大尺度的团簇，类似于重油热解过程中的结焦前驱体.从表3和表4的图像可以看出，较低环已烷密度与适宜的温度有利于 NAP、Bg h i p 在超临界环已烷中的溶解。表4苯并在不同温度下的溶解行为Table 4 Dissolutions of Bghip droplets in SC-cyclohexane atdifferent temperatures;634 K0.1g/cm30.2g/cm30.3g/cm3021003-4第2 期654 K674 K654 K674 K好2第41卷3.1.3多环芳烃混合物在超临界环己烷中的溶解上文对纯NAP分子和Bghip分子的溶解行为进行了研究，而多环芳烃具有不同尺度，组成差异较大的特点，采用MD模拟方法研究了多环芳烃混合物即BghipsoNAP2o（下标代表油滴中苯并芘与萘的比例）混合油滴在超临界环已烷中的溶解行为，如表5所示。表5BghipsoNAP2o油滴在不同密度超临界环已烷中的溶解行为Table 5 Dissolutions of Bghipso NAP2o droplets in SC-cyclo-hexane at different temperatures634K0.1g/cm30.2g/cm30.3g/cm3相比较于纯Bghip 与 NAP体系，NAP与Bghip混合物体系在超临界环已烷环境下的溶解有重要影响模拟50 0 ps后，BghipsoNAP2o的油滴在溶剂箱中有明显的聚集成团簇状在不同温度下，温度越高BghipsoNAP2o的含量相较于Ops时的含量有着明显的降低，说明温度升高会促进多环芳烃混合物在环已烷中的溶解，BghipsoNAP20油滴大部分完全分离，只有小部分的聚集Bghip初始含量有了明显的降低，促进了多环芳烃混合物在环已烷相中的溶解，油滴发生完全解离环已烷相中以NAP为主要成分的多环芳烃分子数量增加，从而减小了油滴的尺寸，导致NAP和Bghip在溶剂箱中均以单体或聚集体的形式扩散.通过对比表456 的图像，轻PAHs比重PAHs在超临界环已烷中的溶解性更好。3.1.4多环芳烃在超临界环烷中的溶解过程从表6 中可以看出，环已烷密度的增加对多环芳烃混合物在超临界环已烷中的溶解有不利的影响当环已烷密度为0.10 g/cm时，溶剂箱中仅存多环芳烃单体和团聚体，而当环已烷密度为赵向勃，等：多环芳烃在超临界环己烷中的溶解：分子动力学模拟研究654 K674 K第2 期0.20g/cm或0.30 g/cm时，溶剂箱中虽然也存在多环芳烃单体和团聚体，但是相比于0.10 g/cm来说，溶解性似乎更差，并且在低环已烷密度下还是在高环已烷密度下进行模拟，都能观察到类似的现象，即轻多环芳烃优先溶解到环已烷相中。以环已烷密度为0.10 g/cm时的模拟为例仅经过12 5ps的模拟，就有相当数量的BghipsoNAP20分子优先溶解到环已烷相中而0.2 0 g/cm和0.30g/cm在12 5ps明显可以看出BghipsoNAP20分子存在明显的聚集现象。表6 6 7 4K下BghipsoNAP2o在超临界环已烷中的溶解过程Table 6 Dissolution procedures of Bghipso NAP2o oil droplets inSC-cyclohexane with varied densities at temperatureof 674 K.0 ps0.1g/cm30.2g/cm30.3g/cm33.2多环芳烃与超临界环已烷相互作用3.2.1径向分配函数多环芳烃与超临界环已烷的混溶是由溶质一溶剂与溶质一溶质相互作用决定的超临界环已烷中多环芳烃的近程溶剂结构和溶剂化自由能反映了体系中溶质一溶剂间的相互作用为了解释多环芳烃的近程溶剂结构，计算了6 7 4K温度下、0.100.30 g/c m环已烷密度中gBghip-H（r）和gNAP-n(r)的径向分布函数（RDF）如图3和图4所示.在三种环已烷密度下，0.10 g/cm的gBclip-H(r)与其他两者分布不同在Bghip分子的油滴半径内，gBghip-H(r)的值保持为0 在2.1A以上半径处，gBgip-H(r)的值逐渐增大，而在0.2 0 g/cm和0.30 g/cm时在10.0 A附近开始接近1，而后持续保持而Fujisawa和Kajimoto等34根据光谱表征，得出如果溶质和溶剂之间存在强吸引的非键相互作用，则容易在溶质周围形成溶剂化壳.在0.10 g/cm时候gBglip-H(r)发现了对应溶剂壳021003-5125 ps375 ps500 ps第41卷存在的峰，则可以认为多环芳烃与超临界环已烷间具有较强的相互作用，0.10 g/cm的环已烷可作为催化剂促进Bghip的溶解，其详细的机理值得进一步研究而0.2 0 g/cm和0.30 g/cm则没有产生溶剂峰，说明高密度下重多环芳烃与超临界环己烷相互作用较弱而在三种密度下的gNAP-H(r)的RDF分布一致，在2.3A处开始增加，而后在10.0 A附近开始趋向于1，继而保持说明轻多环芳烃与超临界环已烷的相互作用较弱，其结果与SCW有着相似的结果.2.52.01.5+(a)H-dia1.0-0.50.0-0.5-图3674K下0.10-0.30 g/cm密度范围中超临界环已烷中Bghip径向分配函数Fig,3 gi-n(r)in SC-cyclohexane at tempera-ture of 674 K and cyclohexane densities from0.10 to 0.30 g/cm32.52.0-1.5-(a)HdNa1.0-0.5-0.0-0.5-图46 7 4K下0.10 0.30 g/cm密度范围中超临界环已烷中NAP径向分配函数Fig.4 g Ap-n(r)in SC-cyclohexane at temperatureof 674 K and cyclohexane densities from 0.10to 0.30 g/cm3.2.2溶剂化自由能计算了NAP和Bghip的溶剂化自由能，如图原子与分子物理学报5所示6 34K时，NAP溶解于超临界环已烷相中，其溶剂化自由能绝对值为4.7 9 8 9 kJ/mol.Bghip分子的溶剂化自由能绝对值为2 1.0 46 3kJ/mol，其绝对值明显大于NAP分子当环已烷密度为0.30 g/cm时，NAP的溶剂化自由能绝对值为7.17 2 9kJ/mol，高浓度的环已烷密度对多环芳烃分子在超临界环已烷环境中的溶解度有消极影响将Bghip和NAP在三种不同环已烷密度的超临界环已烷中溶解，结果表明，环已烷密度的变化影响多环芳烃的溶解度，随着环已烷密度的增加，溶剂化自由能的绝对值显著增加，在0.30 g/0.1 g/cm0.2g/cm0.3g/cm050r(Angstrom)第2 期cm的环已烷密度下，溶剂化自由能绝对值约为7.729kJ/mol和16.8 16 1kJ/mol.另外，在相同温度不同环已烷密度下，Bghip中的溶剂化自由能绝对值相对于超临界环已烷环境中的轻多环芳烃的溶剂化自由能绝对值较大，而溶剂化自由能绝对值越小说明溶解性越好，所以在低密度下和高温下超临界环已烷对NAP和Bghip有着良好的溶解性.尽管多环芳烃的尺度和环已烷的热力学状态发生了变化，但溶剂化自由能的计算值始终为负1015r(Angstrom)10200.1 g/cm0.2g/cm0.3g/cm1520值Bghip分子的溶剂化自由能的绝对值总是比NAP分子更大，证实了在超临界状态下重质多环芳烃比轻多环芳烃更加稳定，同样也意味着NAP在超临界环已烷中的溶解度应该比Bghip效果好.3.2.3内聚能密度由前文的MD模拟结果可知，多环芳烃与超临界环已烷环境下的混相指数同时由另一个因素决定，即多环芳烃分子间的相互作用：为了了解多环芳烃与超临界环已烷环境之间的相互作用，本文计算了NAP和Bghip在超临界环已烷环境下对应温度和压力下的内聚能密度（Cohesive energydensity，CED），计算得到的CED如图6 所示。超临界环已烷环境下中多环芳烃的CED越高，代表着多环芳烃的溶解越困难在相同环己烷密度下6 环的Bghip的CED和2 环NAP的CED相差不大，并且伴随温度的升高两者的CED在逐渐减小，说明在超临界环已烷环境下对NAP 和Bghip都有着不错的溶解能力，并且温度的升高对于NAP和Bghip的溶解性有所提高，这与其前文的MD结果相对应而对于不同密度来说，0.30g/cm的Bghip分子和NAP分子是0.10 g/cm的Bghip分子和NAP分子的8 倍左右这也证实Bghip分子比NAP分子更加稳定且不容易被溶解，021003-6第41卷51015-20-25(a)Solvation free energy at different temperatures赵向勃，等：多环芳烃在超临界环己烷中的溶解：分子动力学模拟研究NAPBghipl634654Temperature(K)（a）不同温度下的溶剂化自由能第2 期3.5E+07NPBghipl3.0E+072.5E+07-2.0E+071.5E+071.0E107-5.0E+060.0F+00674(a)不同密度下多环芳烃内聚能密度(a)CED of PAHs at different densitiesNAPBghipl0.10.2Density(g/cm)0.3NAPBghipl5101015-15202025-250.1(b)不同密度下溶剂化自由能(b)Solvation free energy at different densities图5多环芳烃在超临界环己烷中的溶剂化自由能Fig,5 Solvation of polycyclic aromatic hydrocarbons insupercritical cyclohexane,随着环已烷温度的增加，Bghip在超临界环已烷环境之间的吸引作用增强相比之下，Bghip的内聚能密度对环已烷密度变化更为敏感。在674K温度下，环已烷密度从0.10 g/cm增加到0.30g/cm，导致Bghip的CED大大增加了2.64E7j/cm，当环已烷密度为0.30 g/cm时，随着温度从6 34K升高到6 7 4K，Bg h i p 的CED只降低了2.2 E5 j/cm，表明多环芳烃之间的相互作用增加并且随着环已烷密度的增大超临界环已烷对NAP和Bghip的溶解能力减弱，而温度的增加，会对超临界环已烷环境下NAP和Bghip的溶解有着积极的影响。4结论本文采用分子动力学模拟方法研究了多环芳0.10.20.3Density(g/cm)0.2Density(g/cm)（b）不同温度下多环芳烃内聚能密度(b)CED of PAHs at different temperatures图6 多环芳烃内聚能密度Fig.6 CEDs of PAHs under different temperaturesand densities.烃在超临界环已烷下的溶解度行为，有助于了解重油在超临界环已烷中的溶解行为：结果表明在不同温度与密度下超临界环已烷对于具有2 个芳香环的NAP(轻多环芳烃)均具有良好的溶解性.而温度对于含有6 个芳香环的BghiP的油滴（重多环芳烃)在超临界环已烷中的溶解影响较小：在由多环芳烃混合物组成的油滴溶解过程中，轻环芳烃优先溶解到环已烷相，重环芳烃则集中在粒径减小的油滴中同时，温度升高以及密度的将降低会导致多环芳烃的内聚能密度降低，增加了多环芳烃与超临界环已烷之间的相互作用因此，环已烷密度越小以及温度的升高可以促进多环芳烃或多环芳烃混合物在超临界环已烷中的溶解。参考文献：1Aleklett K,Jakobsson K,Lardelli M,et al.The peak021003-70.3第41卷of the oil age:declining world oil production will halteconomic growth-sydney ideas-the university of syd-ney J.Publications Office,2010,38:1398.2Hosseinpour M,Fatemi S,Ahmadi S J,et al.Thesynergistic effect between supercritical water and redoxproperties of iron oxide nanoparticles during in-situcatalytic upgrading of heavy oil with formic acid.Iso-topic study J.Appl.Catal.B:Environ.,2018,230:91.3 Hao Q,Gong J H,Tan X C,et al.Dissolution of pol-ycyclic aromatic hydrocarbons in subcritical and super-critical water:a molecular dynamics simulation studyJ.Chem.Eng.Sci.,2018,195:958.4Ancheyta J.Modelling of processes and reactors for up-grading of heay petroleum,J.Focus.Catal.,2013,11:8.5Bello S S,Wang C,Zhang M,et al.A review on thereaction mechanism of hydrodesulfurization and hydro-denitrogenation in heavy oil upgrading J.Energ.Fu-el.,35,14:10998.6Fatemi S,Ahmadi S J.Catalytic cracking of petroleumvacuum residue in supercritical water media:impact of-Fe,O,in the form of free nanoparticles and silica-supported granulesJ.Fuel.,2015,159:538.7Jedlovszky P,Vincze A,Horvai G.Full description ofthe orientational statistics of molecules near to inter-faces.water at the interface with CCl J.Phys.Chem.Chem.Phys.,2004,6:1874.8Limmen N.ReaxFF-molecular dynamics simulationsof non-oxidative and non-catalyzed thermal decom-position of methane at high temperatures J.Phys.Chem.Chem.Phys.,2010,12:7883.9Agrawalla S,Duin A.Development and application ofa ReaxFF reactive force field for hydrogen combustionJ.J.Phys.Chem.A,2011,115:960.10 Han Y,Jiang D D,Zhang J L,et al.Development ap-plications and challenges of ReaxFF reactive force fieldin molecular simulationsJ.Front.Chem.Sci.Eng.,2015,10:16.11 Zhang Y M,Li J L,Wang J P,et al.ReaxFF MDSs-based studies on gasification of glucose in supercriti-cal water under microwave heating J.Int.J.Hydro-gen.Energ.,2016,41:13390.12 Akiya N,Savage P E.Roles of water for chemical re-actions in high-temperature water J.J.Chem In-form.,2010,33:293.13 Bai B,Liu Y G,Meng X H,et al.Experimental in-vestigation on gasification characteristics of polycarbon-原子与分子物理学报ate(PC)microplastics in supercritical waterJ.J.Energy.Inst.,2020,93:624.14 Chen J,Fan Y,Zhao X,et al.Experimental investi-gation on gasification characteristic of food waste usingsupercritical water for combustible gas production:ex-ploring the way to complete gasification J.Fuel.,2020,263:116735.15Hosseinpour M,Fatemi S,Ahmadi S J,et al.Isotopetracing study on hydrogen donating capability of super-critical water assisted by formic acid to upgrade heavyoil:computer simulation vs.experimentJ.Fuel.,2018,225:161.16Amani M J,Gray M R,Shaw J M.The phase behaviorof athabasca bitumen+toluene+water ternary mix-tures,J.Fluid.Phase.Equilib