《重质油组分及其溶解机理的分子动力学研究毕业论文.doc》由会员分享，可在线阅读，更多相关《重质油组分及其溶解机理的分子动力学研究毕业论文.doc（67页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 重质油组分及其溶解机理的分子动力学研究Molecular Dynamics Study of the Composition and the Dissolution Theory of Heavy Oil 学科专业： 物理学研究方向： 凝聚态物理作者姓名： 指导教师： 二一三年四月Molecular Dynamics Study of the Composition and the Dissolution Theory of Heavy Oil A Thesis Submitted for the Degree of MasterCandidate：Yu Gongqi Supervisor：
2、Prof. Chu JunCollege of ScienceChina University of Petroleum (East China)关于学位论文的独创性声明本人郑重声明：所呈交的论文是本人在指导教师指导下独立进行研究工作所取得的成果，论文中有关资料和数据是实事求是的。尽我所知，除文中已经加以标注和致谢外，本论文不包含其他人已经发表或撰写的研究成果，也不包含本人或他人为获得中国石油大学(华东)或其它教育机构的学位或学历证书而使用过的材料。与我一同工作的同志对研究所做的任何贡献均已在论文中做出了明确的说明。 若有不实之处，本人愿意承担相关法律责任。学位论文作者签名： 日期： 年 月
3、日学位论文使用授权书本人完全同意中国石油大学(华东)有权使用本学位论文(包括但不限于其印刷版和电子版)，使用方式包括但不限于：保留学位论文，按规定向国家有关部门(机构)送交、赠送和交换学位论文，允许学位论文被查阅、借阅和复印，将学位论文的全部或部分内容编入有关数据库进行检索，采用影印、缩印或其他复制手段保存学位论文。保密学位论文在解密后的使用授权同上。学位论文作者签名： 日期： 年 月 日指导教师签名： 日期： 年 月 日摘 要重质油聚沉特性等相关研究已经取得了较大进展，但受实验条件的限制,有关沥青质溶解性的微观研究尚未取得实质性突破，严重制约了重质油资源的利用。本文以分子动力学、弗洛里-哈金
4、斯理论等为基础，对重质油溶解行为进行计算机模拟，并对结构、能量等定量数据进行计算分析，以期完善重质油溶解理论并将其更好的用于重质油有效利用的指导。首先，以塔河重油(THAR)和绥中重油(SZAR)分子结构为基础，在MS软件中构建重质油组分模型。对各个组分进行密度模拟，得到了与实验数值高度吻合的数据，相对误差均在3.0%以下。通过密度模拟，我们确定了模拟方法的可行性和软件力场、物态方程、流程设置的准确性，并对沥青质组分有了整体性认识。其次，利用分子动力学方法对沥青质在不同溶剂中的溶解行为进行模拟，通过对混合能、平均配位数以及最低能量结构图的研究，得到了沥青质溶解的定量数据，利用弗洛里-哈金斯理论
5、对数据进行分析。沥青质在脂肪类溶剂中的混合能均为正值，说明溶解过程是吸热的，此类溶剂不能溶解沥青质聚沉。而胶质、苯、硝基苯与沥青质的混合能都为负值，表明混合过程放热，较易溶解，因此聚沉抑制剂的研发应以芳香类溶剂为主。此外，模拟结果表明升高温度可以有效增强聚合物的溶解，但温度要控制在合理的范围之内，为了更好的进行重油解聚，温度应保持在600K左右。最后，从能量和构型两方面对沥青质溶解性进行分子动力学研究。能量方面，先后对不同条件下沥青质聚集体的相互作用能、溶解度参数、构型能进行对比分析。结果表明，随着沥青质聚集程度的增加，溶解度参数逐渐减小；范德瓦尔势能和静电能是沥青质聚集的主要作用，而加入溶剂
6、和升高温度可以减弱这些作用，从而增强沥青质溶解。构型方面，分析了沥青质分子的径向分布函数、分子弯曲以及分子在不同溶剂中回转半径、分子间距等的变化。研究表明升高温度后沥青质的结构更加舒展，烷烃支链的伸展可以有效降低聚沉的发生。分子弯曲表明沥青质分子间存在强大的-作用，加入芳香溶剂时，由于溶剂破坏了沥青质分子间-的作用，沥青质分子的间距增大，分子间作用减弱。通过这些模拟得到了沥青质溶解的定量数据，完善了重质油溶解机理，为重质油的有效利用提供了微观理论指导。关键词：重质油聚沉，分子动力学，弗洛里-哈金斯理论，溶解机理，计算机模拟Molecular Dynamics Study of the Comp
7、osition and the Dissolution Theory of Heavy Oil Yu Gongqi (Physics)Directed by Prof. Chu JunAbstractThe research about composition and dissolution theory of Heavy oil has made great progress. Because the limitation of experimental conditions that the microscopic studies of solubility have not obtain
8、ed the breakthrough. All that severely restricted the use of Heavy oil. We simulate the component of Heavy oil with the theories of MD and F-H in MS. The paper should complete the study of the microscopic mechanism and make the guidance for the inhibition of coagulation. Firstly, we select the two t
9、ypes of Heavy oil as the research object (THAR and SZAR). We build a major component of heavy oil model in the MS software. In order to determine appropriate calculation parameters, we simulate the THAR and SZAR component density. We obtained the reference matches the height data and the relation er
10、rors are 3.0%. Through the first simulation, we determine the feasibility and accuracy. We build a foundation for the calculation by the study about the unit of Heavy oil.Secondly, we studied the mixture properties of the asphalt in the five different solvents with the extension F-H and the MD. Afte
11、r simulation, we study the binding energy, the interaction parameter, free energy and so on. We know about the asphalt coagulation ability of depolymerization. The binding energy and interaction parameter show that the solvent insoluble asphaltene within the scope of 100-600K. On the contrary, the b
12、inding energy and interaction is negative. The simulation calculation confirmed that the raising temperature has the positive effect for the solution of heavy oil, but the temperature should be controled in the reasonable scope. If we want get a good depolymerization and the better activity of the c
13、atalyst, the temperature should be at about 600 k. At last, we respectively study the energy and the configuration with MS. In energy aspects: we research the asphalt aggregate under the condition in the different temperature, different solvents. After that, we analyzed the interaction energy, solub
14、ility parameters. We know the role of energy and - in coagulation. In configuration aspects, we simulated the molecules of RDF, the bending of the layered structure, the change of molecular spacing and so on. The result shows that higher temperature can make the structure more stretch. The extension
15、 of branched can reduce the happening of coagulation. When we add aromatic solvent, the space between asphaltene molecules will increase. The aromatic solvent destroys the - bond of asphaltene molecules. These simulation experiments have found the quantitative data of the asphalt coagulation. We per
16、fected the dissolving mechanism of heavy oil and offered the theoretical guidance for the effective use of Heavy oil.Key words：Heavy oil aggregate, Molecular dynamics, F-H, Dissolution mechanisms, Computer simulation目 录第一章 前 言11.1 国际石油资源现状11.2 我国石油发展现状11.3 重质石油定义21.4 重质油研究现状31.4.1 重质油研究历程31.4.2几种重要的
17、重质油模型51.5 研究意义及研究内容71.6 本章小结8第二章 MD方法简介及重质油密度模拟92.1 MD方法简介92.1.1 MD方法的概念92.1.2 MD方法研究现状92.1.3 MD方法在石油领域应用102.2 Material Studio 软件简介112.3 重质油组分密度模拟132.3.1 计算参数与计算流程132.3.2常见化合物密度模拟132.3.3重质油组分模型构建152.3.4 SZAR组分的密度模拟162.3.5 THAR组分的密度模拟182.3.6数据分析192.4本章小结20第三章 F-H理论对沥青质溶解性的研究213.1模拟实验原理与计算流程213.1.1研究应
18、用的原理213.1.2模拟溶剂的选择213.1.3模拟流程的设定223.2主要计算公式223.3数据结果和分析233.3.1 SZAR沥青溶解性分析233.3.2 THAR沥青溶解性分析253.3.3沥青质混合结构的分析273.3.4混合能随温度变化的分析283.4本章小结29第四章 能量角度对溶解机理的模拟和分析304.1 内聚能和溶解度参数计算304.1.1构建模型和计算流程304.1.2计算公式与数据分析334.2聚集体构型能的模拟研究344.2.1构建模型与参数设定344.2.2模拟结果与计算分析364.3本章小节40第五章 结构角度对溶解机理的模拟和分析415.1回转半径与弯曲角度的
19、模拟分析415.1.1回转半径的研究分析415.1.2分子弯曲的研究分析435.2 径向分布函数与分子间距的模拟分析445.2.1 径向分布函数的计算445.2.2 分子间距的研究分析465.3 预防重质油聚沉的理论指导485.4本章小结48结论49参考文献51致 谢57v中国石油大学（华东）硕士学位论文第一章 前 言1.1 国际石油资源现状石油是世界上最重要的能源资源之一，它的主要成分为烷烃、环烷烃和芳香烃等。通过对原油的炼制可以获得汽油、煤油、柴油等液体燃料，以及应用于各种机械的润滑剂、气态烃等。通过化工加工后ou后后，可获得合成性纤维、合成性橡胶、肥料、医用药物、塑料制品、农用药物和家用
20、洗涤剂等。它在汽车、运输、石油化工、建筑、家庭日用等各行各业中都有大量应用，因此被人们称为国家的“经济血液”和能源界的“黑色黄金”1。为了争夺石油能源，历史上曾多次发生战争。今后只要没有一种新的燃料能够取代石油，国际上对石油的争夺战就不会停止，由此可见石油资源对于国家发展的重要性。目前，国际石油产量同比增长极小甚至呈现下降趋势,尽管探测到的全球原油储量正在上升（美权威刊物的统计显示，2012年国际石油的储量从2011年的1.65万亿桶增加到了1.79万亿桶），但这些增幅主要来自重质原油、页岩气等各种开采难度很大的非常规油藏。尽管全球原油的产量在未来30年内可能仍呈现慢增长态势，但到2040年前
21、后将逐步转入下降期2。今后随着石油、天然气勘探和开发难度的增加，势必使开发成本进一步上升，同时石油进出口在一定程度上将受国际政治经济和区域政治经济的影响，未来对于石油资源的争夺将会更加激烈。因此，对于非常规石油（例如重质油、页岩气）的研究就显得至关重要。1.2 我国石油发展现状根据最新数据显示我国石油储量位居世界第九，约为60亿吨。但与前几位相比储量上有很大的差距（仅为沙特阿拉伯的1/6），而且大多数储量为非常规石油资源。现今我国正处于高速发展时期，能源需求量大。据海关统计，2012年我国进口原油量高达2.7亿吨，比去年同期增加7%，且都在较高价位，进口均价大约为874.9美元/吨，上涨37%
22、（图1-1，显示了前两年我国石油进口的量价）。石油资源把握着我国工业建设的命脉，如果主要依靠进口石油，经济发展将受到国外制约条件的严重影响，而要想实现能源的自给，对于本国非常规石油资源开采和利用的研究将具备重要的实际意义3。图1-1 我国进口原油数据Fig1-1 The data of China imported oil我国的石油资源多是中等偏重型的原油，重质油年产量占原油总产量的20%-25%，而且将在今后呈逐年增加趋势，预计2030年将达到38%左右。而世界范围内，虽然石油储量的60%未被开采利用，但这些剩余石油存储资源也主要以重质原油以及油砂等非常规储量存在4。为了解决这些问题，很多石
23、油科研工作者已经对非常规石油进行了各种研究，对于未解决的问题，依然需要人们更多的研究。1.3 重质石油定义重质原油定义：按照国际能源界规划的一般标准，将原油分为轻质、中质、重质以及特稠等几大类。重质油一般指常规开发技术很难开采的、具有较大密度和粘度的原油。而关于重质油定义的讨论开始于1980年，当时美国能源部成立了油砂定义专题研讨会，用以引导重质油等非常规资源的规范化处理5。1982年在委内瑞拉举行的第二届国际重质原油和沥青油砂会议上正式公布了重质油的定义：即在15.6及0.1MPa下，密度为934kg/m3-1000kg/m3的原油。而在石油行业中一般的划分为：轻质原油的密度小于852kg/
24、m3；中质原油的密度在852kg/m3到930kg/m3之间；重质原油的密度在931kg/m3到998kg/m3之间；超稠原油的密度一般在998kg/m3以上6。重质原油主要是由含有氮、硫、金属元素等杂原子的芳香中心以及脂肪支链共同组成的结构较为复杂的化合物7。由于各地区石油的形成情况不同，地理环境、气候等因素导致了地下油藏组成有较大差别，所以不同地区的重质油组成也有差异。由上面的定义我们可以看出，重质油的密度是比较大的，这主要是因为重质油的主要成分是芳香性的胶质、沥青质，脂肪性的饱和分等，这些大分子结构致使重质油密度较大。如果按物质的极性给重油的主要组分进行一个排位的话，大体是：沥青质胶质饱
25、和分。其中，重质油主要成分中沥青质是结构最复杂、分子量最大、极性最强的组分，沥青质分子的性质一定程度上决定了重质油的稳定性、流变性、相分离性、浸润性等性质8。沥青质化学组成复杂、结构大，如果对此不进行深入研究，在提炼和加工重质油的过程中都将具有很大的难度。因此对于重质油的研究，主要在于沥青质大分子性质的研究。1.4 重质油研究现状1.4.1 重质油研究历程在20世纪30年代时，Nellensteyn和Pfeiffer提出了原油沥青质的胶体模型。即重质油中的大分子相互结合形成胶束，主要成分都胶溶于石油中9-13，这是首次提出沥青质不是以单独的分子存在于石油中，而是以胶束胶溶于石油中，其中胶质分子
26、是沥青质的胶溶剂。自此，国外科学家运用了各种方法对此进行了研究。从20世纪40年代开始，人们使用了多种不同的方法去研究重质油的组成和性质。例如：利用烷基转移和还原降解、选择氧化等方法的化学方法，它通过分析研究化学降解物可以推出分子热反应性能、主要反应途径、沥青质结构等14-16。另外还有一些物理方法，主要包括红外光谱、核磁共振和质谱检测等方法17。在20世纪50年代末期，科学家首次将核磁共振技术应用于重质油的研究中，在测定物质平均结构的应用中取得了很多成果。它已成为当今研究重质油和沥青质结构与组成最常用和最成功的方法之一。这种技术用来测量原油成分、化工馏分以及从煤炭中提取的燃料高分子结构，都取
27、得了巨大成功17。由NMR测定、区分并定量给出了原油组分中的烯、内烯以及多取代内烯等的比例，并运用此项技术分析估算出了沥青质中多核芳烃的平均芳环数。这些数据对于石油炼制、石油化工等诸多领域都是很有用的。这些研究成果都表明，核磁共振波谱法是研究带有脂肪性侧链的大分子稠环结构最为有力的工具。它已成为研究石油组分等的重要分析手段，在各个研究领域发挥着重要作用。在20世纪60年代中期，质谱方法被用来研究石油分子的结构和性质，其中高压电离质谱法在石油分子的结构表征中占据了很重要的地位，在表征离子对、中间馏分、汽油、柴油等组分的结构时也取得了很大成功。随着研究的需要和技术的进步，LVEI-HRMS（低压碰
28、撞电离技术和高分辨率质谱技术的联合应用）也在石油分子的芳香性分析中得到了极大应用18。在20世纪70年代末，在国外科学家Yen19等人的工作中，对固体粉末的沥青质进行了x射线衍射测定(通过5核物质撞击得到沥青质芳香结构)，他们实验得到的芳香族组成的属性和峰值与金刚石的属性和峰值相当一致，Dickie在实验中也获得了相似的结果。基于这些结果，他们利用宏观结构提出了沥青质等大分子的聚合薄片模型，这种模型是由芳香族叠放而成，这个模型是受到金刚石结构的启发而得到，沥青质以重叠的片状结构出现，这个结构是由一个多核结构组成，以芳环片段为中心，烷基侧链连接的大分子结构。他们提出沥青质聚合体是通过3-6个芳香
29、薄片，可能以-作用聚集在一起。但经过后来的模拟得到的一个重要的结果证实这些聚集仅仅通过一种作用是无法形成的，芳香族的相互作用还有其他的情况，这与Hunter20所做的另外的实验结果相一致。因此，不同方向的作用可能会形成不同构型的沥青质聚集体，这些研究发现与Sheu21曾经在甲苯溶剂中所做实验的结果是一致的（基于对64个沥青质分子的模拟，聚集结构主体在3到11个芳香环范围内浮动）。Sheu也发现了沥青质聚集体不仅形成了面面型聚集而且还形成了其他结构的聚集体，这种聚集是容易形变并且在外观上也是不规则的。我们对于这些科学前辈们的结论进行下面的概述：沥青质聚集胶束的结构主要是面对面型(-)、侧向性（T
30、型及型）19。他们还定义这些单一的结构为晶胞或粒子，实验表明重质油是通过这种胶粒形成胶核的，芳香烃混合物的这种结构和沥青质的聚集结构已经通过经典的分子动力学模拟进行了理论上的研究，研究表明在真空条件和有溶剂下会发生诱导聚集和自发聚集，这些聚集形态在自然界是无定形的22。20世纪80年代初期开始，我国科学家也陆续开展了重质油组分、组成结构的研究工作，提出了改进的Brown-Ladner方法23。并用该法研究了多种重质油组分的平均组成，取得了很多优秀的成果。为了对重质油进行更加全面的研究，我国教育部与各石油生产、教学单位在2003年也专项成立了重质油国家重点实验室，并在重质油化工、重油轻质化工艺研
31、究、重油转化催化剂、催化材料、新型清洁燃料生产、重质油的高附加值应用等诸多方面进行了很多科研工作，取得了丰硕的成果。但是由于起步晚、技术落后等原因，我国在石油组分分析、聚集抑制剂、重质油加工等方面与国外先进技术还是有一定差距。我国的石油科研工作者们也正努力提高技术，逐步缩短差距。通过前人不断的努力和各种物理、化学技术的应用，人们对于重油组分的宏观结构有了清晰的认识。这些科学研究及文献参考数据也是本课题继续深入研究重质油结构的前提条件。1.4.2 几种重要的重质油模型(1) 如图1-2，1-3是科学家Yen 19等人提出的重质油模型，其中沥青质的模型主要结构包括缩合芳环构成的芳香片、烷基侧链官能
32、团、杂原子取代基(含硫、氧、氮)等。沥青质芳香薄片依靠分子间作用力、分子内的各种作用在立体构型上形成三维结构，不稳定时即以聚集薄片为核心形成胶束聚沉。这种模型奠定了沥青质分子结构的研究基础，并且也能解释沥青质的胶溶理论，这一结构模型已经为大量的宏观和微观研究结果所证实。图1-2 重油组分示意图-YenFig1-2 Element of asphaltene by Yen图1-3 重油胶束示意图-YenFig1-3 Aggregation of asphaltene by Yen(2) 如图1-4，1-5是中国科学院张胜飞等进行的粗粒化模型的研究成果8，是现代的沥青质分子的自组织和连续网络等新的
33、理论结构模型的一个典范。图1-4 沥青质组分示意图-粗粒化Fig1-4 Element of asphaltene by Shengfei Zhang 图1-5 重油组分示意图-粗粒化Fig1-5 Element of asphaltene by Shengfei Zhang (3) 如图1-6(a-f)是重质油国家重点实验室24-27的杨朝合教授等利用尺寸实验表征法构造的THAR（塔河重质油）和SZAR（绥中重质油）组分模型图,从a到f依次是THAR沥青质模型、SZAR沥青质模型、THAR芳香分模型、THAR胶质模型、SZAR芳香分模型、SZAR胶质模型。由于实验采用尺寸实验表征法使得重油各
34、组分模型更加精确，并且石油原料采自我国油田，使得研究更具有实践意义，因此本文将主要以此为基础进行模型的构建等分析工作。(a)THAR沥青质 (b)SZAR沥青质(c)THAR芳香分 (d)THAR胶质(e)SZAR芳香分 (f)SZAR胶质图1-6 重油组分示意图-重质油实验室Fig 1-6 Element of asphaltene by HONL1.5 研究意义及研究内容(1) 研究意义：石油聚沉对石油的开采、运输和最后加工利用各个环节都会产生严重影响，而聚沉主要是由于沥青质引起的。沥青质的聚沉会降低石油流动性，导致油藏坏死、堵塞输油管线等诸多严重问题，主要可概括为：沥青质在油藏中会附着于
35、岩层中，使得原油的流动性减弱，造成油藏减产甚至坏死。在石油开采时，沥青质聚沉会造成分离器、提油设备和管道等器具的堵塞，破坏器具导致原油无法提取。在原油的加工时，聚沉会导致催化剂失活，降低工业生产效率和工业产量，造成资源的浪费。目前世界石油资源正逐渐往重质化方向发展，国内外很多油田的产量均趋于递减状态28。由于这些危害，石油重质化已成为当今世界能源开采利用的热点和难点。国内外石油科研工作者对此进行了很多深入的研究，在重质油的组成，胶质、沥青质的平均分子结构以及重质油胶粒的缔合结构等方面都取得了很多成果。但由于实验条件等的限制，很难做出分子水平上聚沉机理的研究29。而利用分子动力学对此进行模拟，可
36、以很好的填补真实实验的空缺，给出更多重质油组分的信息，为重质油资源的更好利用提供更多科学依据。因此，本课题具备很重要的实践意义。(2) 研究内容：第二章将在国内外石油科研工作者结构分析试验的基础之上，构建适用于分子模拟的沥青质、胶质、芳香分和饱和分的分子模型，能够在分子水平上对其进行定量描述。对常见化合物进行动力学模拟，确定合适的力场及边界条件，然后对各种组分进行密度的MD模拟，以体现模拟的可行性和物态方程的准确性。第三章运用分子动力学方法、分子力学、Extended FloryHuggins方法相结合的方式，对SZAR沥青质和THAR沥青质在不同溶液中的平均束缚能、平均配位数和混合能进行计算
37、分析，以验证沥青质在不同溶剂中的溶解性大小。这是F-H理论在重质油研究领域的首次应用。通过本章研究为解聚机理的分析打下基础。第四章构建沥青质单体、二聚体、三聚体的模型（聚集程度不断提高），通过MS模拟获得沥青质的溶解度参数、范德华势能(Evdw)、内聚能、静电能(Ee)等能量，从能量角度对解聚机理进行研究。第五章模拟计算沥青质大分子的回转半径，以确定分子在不同状态下的伸展情况。分析沥青质分子的角度弯曲、RDF（径向分布函数）、分子间距等结构性质在不同情况下的数据，从结构角度对解聚机理进行研究，并给出预防重质油聚沉的理论指导。1.6 本章小结本章对国内外石油资源的应用现状进了调研分析，对重质油聚
38、沉的危害和研究现状进行了简单介绍，并由此明确了重质油研究中存在的问题。针对定量数据不足、重质油溶解机理尚未成型等现有成果的局限性，本文将以分子动力学、弗洛里-哈金斯理论为基础，对重质油的组分、结构、溶解机理等进行计算机模拟研究。57第二章 MD方法简介及重质油密度模拟2.1 MD方法简介2.1.1 MD方法的概念MD (Molecular Dynamics-分子动力学)，是一门结合物理、数学和化学的综合性技术。MD是分子模拟的一种方法，它利用牛顿力学的原理对系统的运动进行模拟，再从运动系统的不同状态中随机抽取样本计算构型积分30，利用构型积分为基础，完成对模拟体系热力学性质、动力学性质和其他宏
39、观性质等的计算研究31。我们所用的MD方法，就是用计算机运算来表示统计力学的方法，主要涉及粒子运动的动力学问题。作为实验的一个辅助手段，MD模拟就是对电子及中心核共同形成的系统，求解其运动方程(例：牛顿方程，哈密顿方程，拉格朗日方程等)，系统中的每个原子核运动时都被受到其它全部原子核及电子的作用,模拟并分析体系中各个粒子受力情况，用经典的运动理论和量子的方法来求解体系中各个粒子在随机某个时刻的位置、速度，以确定粒子当时的运动状态，进而推断系统的结构和性质32-34。2.1.2 MD方法研究现状1957年，Alder等人第一次提出了经典的分子动力学方法，并首先在“硬球”液体模型上应用了该方法35
40、。他们用这种方法验证了由Kirkwood在1939年所预言的“刚性球系统从液相状态到结晶相状态的转变”（统计力学原理），这就是人们所熟知的Alder相变。1963年，Rahman在利用经典MD方法研究了连续性液体模型的性质，得出了与实际试验相似的结果36-39。这种方法逐渐被人们所认可和熟知。1972年，Less等科学家发展了经典的MD方法，并扩展应用到了含有速度梯度的非平衡系统中40-43。1980年，Andersen等人在研究中发现并应用了恒压下的MD方法，这种方法对后来的NVT、NPT的研究产生了深远影响44-45。1983年，Gill等人将这种方法推广应用到了包含温度梯度的非平衡系统，
41、从而形成了后来常被应用的非平衡系统MD方法体系46。1984年初，Nose等人在恒压MD方法的启发下创建了恒温MD方法，首次实现了恒温下物质实验的模拟，取得了极大成功47。1985年，为了解决构建势函数模型难度较大的半导体和金属材料等物质，Car等研究者提出并试验了将电子论和MD方法统一起来进行研究，最终得出了第一性原理的MD方法48，这是目前为止应用最广、准确性最高的一种方法，也是本文研究的一个原理基础。1991年，Cagin等又提出了巨正则系统MD方法，这种方法能够用于处理吸附等问题的模拟49-51。人们在摸索中不断发现并解决问题，逐步发展完善了这种方法，经过几十年的发展到今天我们所应用的
42、MD方法。而从20世纪80年代后期开始到现在，大型计算机技术飞速发展，加上理论上多体势函数的提出和发展，使得分子动力学模拟技术有了进一步的发展。本文在总结前人研究的基础上，根据计算机模拟的条件，选择新的重油组分模型。根据第一性原理确定合适的力场和正确的周期性边界条件，进行动力学模拟。并依据化学物性数据及试验数据改变组分模型，得出可能影响重油组分性质的杂原子等因素，进行新的模拟对比并总结影响重油组分关系的因素。研究了重质油的溶解性、静电相互作用、整体相互作用能等，找到了使重质油聚集的主要作用力，同时对不同溶剂中重质油聚集体的伸展情况进行了对比，逐步找到了重油解聚剂的发展方向。2.1.3 MD方法
43、在石油领域应用虽然MD方法很早就开始在科研中应用，但是分子模拟技术在石油研究中的应用却开始于20世纪的末期。分子模拟技术在油品添加剂开发中的应用，可以帮助科学家们从微观角度分析和了解油品添加剂的作用机理及复合效应，能够实现对添加剂结构、性能等的深入研究。由于在节省时间、节约实验以及完成实验所不能的预测等方面的优势，使得这种开发研究具有非常好的应用前景。特别是近些年计算机技术的飞速发展，使得与此相关的文献报道也呈井喷增加趋势。1990年，国外科学家第一次将MD模拟技术用于油田试剂分子的设计中，虽然实验受当时技术及原理限制并未取得突破性进展，但为后来的工作奠定了一定基础。2003年6月，Timot
44、hy等人用MD模拟技术，模拟了润滑油中改性剂（有机钼）的溶液相反应，根据模拟结果显示二烷基-二硫代氨基甲酸钼和二烷基-二硫代磷酸锌得反应受到的主要影响，来自软硬酸和溶剂之间的交互作用，二烷基-二硫代氨基甲酸钼和二烷基-二硫代磷酸锌在润滑油中的存在形式是离子对52。分子模拟的结果较好的解释了实验中所观察到的现象。2004年，Nicholas等人分别用量子力学法与分子动力学法对几种二烷基-二硫代磷酸锌盐(ZDDP)的抗磨添加剂热分解进行了模拟研究。对于ZDDP单体、ZDDP二聚体、含有不同取代基的交联异构体(LI-ZDDPs)进行了模拟计算53-56，得到了关于其反应的重要数据。很好的解释了试验的
45、不合理现象，并弥补了实际试验的缺陷。伴随着全世界范围石油能源的不断消耗，对非常规石油能源的利用和开发老油藏的开发潜力，已经是全世界各国增加生产的主要举措。进入新世纪后我国的石油科研工作也日趋完善。据推测，我国未来石油类化学发展基本的趋势为：研究从单一平衡逐渐走向非平衡态和平衡态共同研究；从单向性走向多向性；从完全定性逐渐发展为精确定量。而在微观和定量上，MD模拟具有独特的优势，它描述的微观物理体系中，每个分子、原子都各自服从经典的牛顿力学，通过应用微观数学模型描述的物理体系，去建立完整的重油组分分子模型，选择一定的周期性边界条件。然后进行动力学的模拟（软件遵从第一性原理），对有用数据进行提取并
46、分析得出的数据，推出实际理论。因为MD方法不会有任何的随机因素存在，故得出的数据较试验数据更具理论研究意义，可以与实际试验相结合，得到更加完善的理论。我们所应用的主要是依据物质的化学物性数据、试验数据，改变组分模型，去发现可能影响重油组分性质的因素。因此采用MD模拟等方法，运用高速发展的计算机来定量、微观的对石油相关领域的问题进行研究，将成为未来石油研究的一个大趋势。2.2 Material Studio 软件简介Material Studio 5.0软件是Accelrys公司为材料科学领域、研发领域所开发的新一代计算机应用软件28。MS中的Amorphous cell模块能够方便快速的构建分
47、子的3D模型，深入进行有机物质、无机晶体、无定型材料、胶体聚合物等的分析，可以在催化试剂、固体化学、聚合物化学、结晶结构学、石油化学等能源和材料领域进行性质预测和聚合物的建模等工作57。软件开发是基于分子理学、动力学、第一性原理等进行。操作灵活、界面可视、直观具体，并且很好的利用了现今飞速发展的计算机科学。Accelrys是美国一家著名的软件公司，它是世界上提供分子模拟、化学信息学和生物信息学等学科软件解决方案的最大电子供应商，其中MS软件是它最成功的作品之一，量子力学密度泛函理论、分子动力学、半经验的量化计算方法、分子力学方法和介观动力学模拟方法等是其软件平台所遵循的一些基本原理。MS软件可以在任何一台装有Linux、Unix、Windows、Win7等操作系统的计算机上进行运行，软件中常用的模块有：Dmol3、Castep、Vamp（这三个模块主要运用量子力学的规律）