CN110021370B - 一种基于分子模拟的共聚物抗盐性能评价方法及其系统 - Google Patents

一种基于分子模拟的共聚物抗盐性能评价方法及其系统 Download PDF

Info

Publication number
CN110021370B
CN110021370B CN201711182823.0A CN201711182823A CN110021370B CN 110021370 B CN110021370 B CN 110021370B CN 201711182823 A CN201711182823 A CN 201711182823A CN 110021370 B CN110021370 B CN 110021370B
Authority
CN
China
Prior art keywords
polymer
molecules
saline solution
density
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711182823.0A
Other languages
English (en)
Other versions
CN110021370A (zh
Inventor
杨帆
杨小华
王海波
林永学
金军斌
孔勇
王琳
褚奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering filed Critical China Petroleum and Chemical Corp
Priority to CN201711182823.0A priority Critical patent/CN110021370B/zh
Publication of CN110021370A publication Critical patent/CN110021370A/zh
Application granted granted Critical
Publication of CN110021370B publication Critical patent/CN110021370B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Addition Polymer Or Copolymer, Post-Treatments, Or Chemical Modifications (AREA)

Abstract

本发明提供一种基于分子模拟的共聚物抗盐性能评价方法及其系统，该方法包括：构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；构建以聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；建立包含水分子、无机盐分子、以及聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对聚合物盐水溶液模型盒子的密度、温度进行设定；对聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；根据得到动力学平衡状态溶液体系结构模型计算特性黏数。本发明不受实验条件限制，效率高；不受人为因素影响，返回性强，结果可靠。

Description

一种基于分子模拟的共聚物抗盐性能评价方法及其系统
技术领域
本发明涉及油田化学聚合物抗盐性能评价的技术领域，尤其涉及一种基于分子模拟的共聚物抗盐性能评价方法及其系统。
背景技术
抗盐性是油田化学聚合物的重要性能之一，无机盐对聚合物的侵害主要表现为降粘作用，粘度降低导致钻井液流变性变差，从而影响聚合物溶液整体性能，因此，评价聚合物抗盐性能时着重关注聚合物溶液的粘度参数。传统实验方法是通过六速粘度计测试聚合物溶液表观粘度来评价其抗盐性能，属于宏观评价方法，其测试评价受到样品纯度、来源、温度条件等诸多限制，且认为因素影响较大，重复性不佳。
因此，如何更准确地对共聚物抗盐性能进行评价仍是本领域需要攻克的难题之一。
发明内容
本发明的特征和优点在下文的描述中部分地陈述，或者可从该描述显而易见，或者可通过实践本发明而学习。
聚合物是由小分子单体聚合而成的长链大分子，分子中的单键能够内旋转，以致分子具有许多不同的构象，构象不同对应分子尺寸不同，高分子物理中用回转半径R(g)表征分子尺寸，常用特性粘数η表征聚合物溶液粘度，η与R(g)直接通过一定公式相联系。聚合物受到无机盐侵害时，盐离子与聚合物分子间的相互作用会导致聚合物分子尺寸发生改变，宏观表现为粘度的变化。R(g)可通过分子模拟计算得到，继而通过公式可得η，即可表征聚合物的抗盐性能。本发明针对传统评价方法存在的问题以及分子模拟方法的特殊优势和应用前景，提供了一种基于分子模拟的共聚物抗盐性能评价方法及其系统，利用计算机模拟快速计算聚合物尺寸，通过尺寸与特性粘数的关系得到特性粘数这一宏观性能参数。
根据本发明的一个方面，提供一种基于分子模拟的共聚物抗盐性能评价方法，包括：
构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；
构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；
建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定；
对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；
根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
在一个可能的设计中，所述对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型包括：
使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构；
将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综下进行动力学计算，获取盒子密度；
将所述盒子密度与所述预设密度进行对比，获取偏差值；
判断所述偏差值是否大于预设阈值；
若否，则将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算。
在一个可能的设计中，若所述偏差值大于预设阈值，则将所述盒子密度设为预设密度；并返回所述建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定的步骤。
在一个可能的设计中，所述根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数包括：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R；
通过以下公式计算特性黏数η：
Figure BDA0001479530990000031
Figure BDA0001479530990000032
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
根据本发明的另一个方面，提供一种基于分子模拟的共聚物抗盐性能评价系统，包括：
构建模块，用于构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；
合成模块，用于构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；
盒子建立模块，用于建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定；
平衡模块，用于对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；
计算模块，用于根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
在一个可能的设计中，所述平衡模块还用于：
使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构；
将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综下进行动力学计算，获取盒子密度；
将所述盒子密度与所述预设密度进行对比，获取偏差值；
判断所述偏差值是否大于预设阈值；
若否，则将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算。
在一个可能的设计中，所述平衡模块还用于：
若所述偏差值大于预设阈值，则将所述盒子密度设为预设密度；并返回所述建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定的步骤。
在一个可能的设计中，所述计算模块用于：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R；
通过以下公式计算特性黏数η：
Figure BDA0001479530990000041
Figure BDA0001479530990000042
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
本发明提供了一种基于分子模拟的共聚物抗盐性能评价方法及其系统，利用计算机“虚拟实验”进行评价，不受实验条件限制，效率高；不受人为因素影响，返回性强，结果可靠；通过本方法得到的特性粘数，可以在一定程度上比较不同聚合物抗盐能力，对钻井液乃至油田化学用抗盐聚合物的设计和研发具有很好的借鉴意义。
通过阅读说明书，本领域普通技术人员将更好地了解这些技术方案的特征和内容。
附图说明
下面通过参考附图并结合实例具体地描述本发明，本发明的优点和实现方式将会更加明显，其中附图所示内容仅用于对本发明的解释说明，而不构成对本发明的任何意义上的限制，在附图中：
图1为本发明实施例的基于分子模拟的共聚物抗盐性能评价方法的流程示意图。
图2为本发明实施例的得到动力学平衡状态溶液体系结构模型的流程示意图。
图3为本发明实施例的基于分子模拟的共聚物抗盐性能评价方法的流程示意图。
具体实施方式
如图1所示，本发明实施例提供一种基于分子模拟的共聚物抗盐性能评价方法，包括：
101、开始。
102、构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化。
具体实施时，可以利用材料性能模拟软件Materials Studio的3D建模草图软件Skech工具绘制水分子、无机盐分子的三维分子结构模型，通过分子动力学计算模块Forcite进行结构优化，得到最小能量构象。
同样可以利用材料性能模拟软件Materials Studio的3D建模草图软件Skech工具绘制聚合单体的分子模型，通过分子动力学计算模块Forcite进行结构优化，获取聚合单体最小能量构象结构。
103、构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构。
更具体地，可以利用Build工具构建以该聚合物分子最小能量构象结构组成的聚合物分子模型，并进行结构优化；获取聚合物分子最小能量构象结构。
104、建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定。
可以采用Amorphous Cell模块建立包含一定数量水分子、无机盐分子、聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，一般地，聚合物分子最小能量构象结构的个数为1个，水分子的个数为500至3000，无机盐分子的个数可以根据盐水浓度设定，也可以设置在10个以内。本发明并不对上述水分子、无机盐分子、聚合物分子最小能量构象结构的个数做限定。
上述预设密度例如为0.7-1.1g/cm3，预设温度例如为298-493K。需要说明的是上述预设密度、预设温度还可以根据水分子、无机盐分子、聚合物分子的数目自行设定。
105、对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型。
106、根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
107、结束。
在图1对应的实施例的基础上，本发明实施例提供一种基于分子模拟的共聚物抗盐性能评价方法，如图2所示，所述对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型即步骤105包括：
201、开始。
202、使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构。
其中，正则系综是由组成系综的系统与一温度为T的很大的热源相接触并达到热平衡后形成的。选择退火功能后，可以进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过动力学功能进行动力学平衡计算。
203、将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综(NPT)下进行动力学计算，获取盒子密度。
更具体地，是处于等温等压系综(NPT)下进行100-1000ps、100000步数的动力学计算。
204、将所述盒子密度与所述预设密度进行对比，获取偏差值。
205、判断所述偏差值是否大于预设阈值。若否，则进入步骤206；若是，则进入步骤207；
上述预设阈值为5％至10％，例如是7％。
206、将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算。
更具体地，是使体系处于正则系综(NVT)下进行100-1000ps动力学计算。
207、将所述盒子密度设为预设密度；并返回步骤104；
即返回步骤104，并在步骤104中，将在步骤203中计算出的盒子密度作为预设密度，用于对所述聚合物盐水溶液模型盒子的密度进行设定。
208、结束。
在图1或图2对应的任一实施例的基础上，本发明实施例提供一种基于分子模拟的共聚物抗盐性能评价方法，所述根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数即步骤106包括：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R；并通过以下公式计算特性黏数η：
Figure BDA0001479530990000081
Figure BDA0001479530990000082
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
在具体实施时，可以运用Forcite模块的分析功能计算盐水溶液中聚合物回转半径R(g)；更具体地，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到R(g)分布图，对其进行统计平均计算得到R(g)数值。
下面辅以具体的实例进一步说明本发明提供的基于分子模拟的共聚物抗盐性能评价方法：
实例一：评价共聚物P(AM-AA)在1％NaCl离子浓度下抗盐性能
在步骤102中，通过Materials Studio软件绘图工具，绘制水分子、NaCl分子、聚合物单体聚丙烯酰胺(AM)、聚丙烯酸(AA)分子模型，通过Forcite模块进行结构优化，采用COMPASSII力场。
在步骤103中，通过Biuld工具构建AM与AA的无规共聚物记为P(AM-AA)，选择已建好的AM和AA单体模型，输入链长10，输入单体摩尔比为0.5:0.5，建立P(AM-AA)，通过Forcite进行结构优化。
在步骤104中，通过Amorphous Cell模块建立聚合物溶液模型，设置体系密度为1.0g/cm3，选择P(AM-AA)，设置分子个数为1，选择NaCl分子，设置分子个数为3，选择水分子，设置分子个数为1000，温度设为298K，力场选择COMPASSII。得到边长为
Figure BDA0001479530990000091
埃的聚合物盐水溶液体系模型的盒子。
在步骤105中，在Forcite模块中选择退火功能，进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过该模块中的动力学功能进行动力学平衡计算。首先选择NPT系综、298K条件，进行100ps、100000步数的动力学计算，控温方式选择Nose，控压方式选择Berendsen，计算结果中实际密度为1.021g/cm3，小于偏差范围，进行下一步NVT系综动力学计算，参数选择同上，得到动力学平衡状态溶液体系结构模型。
在步骤106中，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到R(g)分布图，对其进行统计平均计算得到R(g)数值，带入公式(1)计算得到P(AM-AA)的特性粘数η。
实例二：评价共聚物P(AM-AA-AMPS)在1％NaCl离子浓度下抗盐性能
在步骤102中，通过Materials Studio软件绘图工具，绘制水分子、NaCl分子、聚合物单体AM、AA、2-丙烯酰胺-2-甲基丙磺酸(AMPS)分子模型，通过Forcite模块进行结构优化，采用COMPASSII力场。
在步骤103中，通过Biuld工具构建AM、AA、AMPS的无规共聚物记为P(AM-AA-AMPS)，选择已建好的AM、AA、AMPS单体模型，输入链长10，输入单体摩尔比为0.5:0.25:0.25，建立P(AM-AA-AMPS)，通过Forcite进行结构优化。
在步骤104中，通过Amorphous Cell模块建立聚合物溶液模型，设置体系密度为1.0g/cm3，选择P(AM-AA-AMPS)，设置分子个数为1，选择NaCl分子，设置分子个数为6，选择水分子，设置分子个数为1000，温度设为298K，力场选择COMPASSII。得到边长为
Figure BDA0001479530990000102
的聚合物盐水溶液体系模型的盒子。
在步骤105中，在Forcite模块中选择退火功能，进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过该模块中的动力学功能进行动力学平衡计算。首先选择NPT系综、298K条件，进行100ps、100000步数的动力学计算，控温方式选择Nose，控压方式选择Berendsen，计算结果中实际密度为1.023g/cm3，小于偏差范围，进行下一步NVT系综动力学计算，参数选择同上，得到动力学平衡状态溶液体系结构模型。
在步骤106中，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到Rg分布图，对其进行统计平均计算得到Rg数值，带入公式(1)计算得到P(AM-AA)的特性粘数η。
实例三：评价共聚物P(AM-AA-AMPS)在1％CaCl2离子浓度下抗盐性能
在步骤102中，通过Materials Studio软件绘图工具，绘制水分子、CaCl2分子、聚合物单体AM、AA、AMPS分子模型，通过Forcite模块进行结构优化，采用COMPASSII力场。
在步骤103中，通过Biuld工具构建AM、AA、AMPS的无规共聚物记为P(AM-AA-AMPS)，选择已建好的AM、AA、AMPS单体模型，输入链长10，输入单体摩尔比为0.5:0.25:0.25，建立P(AM-AA-AMPS)，通过Forcite进行结构优化。
在步骤104中，通过Amorphous Cell模块建立聚合物溶液模型，设置体系密度为1.0g/cm3，选择P(AM-AA-AMPS)，设置分子个数为1，选择CaCl2分子，设置分子个数为3，选择水分子，设置分子个数为2000，温度设为298K，力场选择COMPASSII。得到边长为
Figure BDA0001479530990000101
的聚合物盐水溶液体系模型的盒子。
在步骤105中，在Forcite模块中选择退火功能，进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过该模块中的动力学功能进行动力学平衡计算。首先选择NPT系综、298K条件，进行100ps、100000步数的动力学计算，控温方式选择Nose，控压方式选择Berendsen，计算结果中实际密度为1.018g/cm3，小于偏差范围，进行下一步NVT系综动力学计算，参数选择同上，得到动力学平衡状态溶液体系结构模型。
在步骤106中，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到Rg分布图，对其进行统计平均计算得到Rg数值，带入公式(1)计算得到P(AM-AA)的特性粘数η。
实例四：评价共聚物P(AM-AA-AMPS)在4％CaCl2离子浓度下抗盐性能
在步骤102中，通过Materials Studio软件绘图工具，绘制水分子、CaCl2分子、聚合物单体AM、AA、AMPS分子模型，通过Forcite模块进行结构优化，采用COMPASSII力场。
在步骤103中，通过Biuld工具构建AM、AA、AMPS的无规共聚物记为P(AM-AA-AMPS)，选择已建好的AM、AA、AMPS单体模型，输入链长10，输入单体摩尔比为0.5:0.25:0.25，建立P(AM-AA-AMPS)，通过Forcite进行结构优化。
在步骤104中，通过Amorphous Cell模块建立聚合物溶液模型，设置体系密度为1.02g/cm3，选择P(AM-AA-AMPS)，设置分子个数为1，选择CaCl2分子，设置分子个数为12，选择水分子，设置分子个数为2000，温度设为298K，力场选择COMPASSII。得到边长为
Figure BDA0001479530990000111
的聚合物盐水溶液体系模型的盒子。
在步骤105中，在Forcite模块中选择退火功能，进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过该模块中的动力学功能进行动力学平衡计算。首先选择NPT系综、298K条件，进行100ps、100000步数的动力学计算，控温方式选择Nose，控压方式选择Berendsen，计算结果中实际密度为1.037g/cm3，小于偏差范围，进行下一步NVT系综动力学计算，参数选择同上，得到动力学平衡状态溶液体系结构模型。
在步骤106中，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到Rg分布图，对其进行统计平均计算得到Rg数值，带入公式(1)计算得到P(AM-AA)的特性粘数η。
如图3所示，本发明还提供一种基于分子模拟的共聚物抗盐性能评价系统，包括构建模块301、合成模块302、盒子建立模块303、平衡模块304、计算模块305。其中：
构建模块301，用于构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；具体实施时，可以利用材料性能模拟软件Materials Studio的3D建模草图软件Skech工具绘制水分子、无机盐分子的三维分子结构模型，通过分子动力学计算模块Forcite进行结构优化，得到最小能量构象。
同样可以利用材料性能模拟软件Materials Studio的3D建模草图软件Skech工具绘制聚合单体的分子模型，通过分子动力学计算模块Forcite进行结构优化，获取聚合单体最小能量构象结构。
合成模块302，与所述构建模块301相连，用于构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；
更具体地，可以利用Build工具构建以该聚合物分子最小能量构象结构组成的聚合物分子模型，并进行结构优化；获取聚合物分子最小能量构象结构。
盒子建立模块303，与所述构建模块301以及合成模块302相连；用于建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定；
可以采用Amorphous Cell模块建立包含一定数量水分子、无机盐分子、聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，一般地，聚合物分子最小能量构象结构的个数为1个，水分子的个数为500至15000，无机盐分子的个数可以根据盐水浓度设定，也可以设置在10个以内。本发明并不对上述水分子、无机盐分子、聚合物分子最小能量构象结构的个数做限定。
上述预设密度例如为0.7-1.1g/cm3，预设温度例如为298-493K。需要说明的是上述预设密度、预设温度还可以根据水分子、无机盐分子、聚合物分子自行设定。
平衡模块304，与所述盒子建立模块303相连；用于对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；
计算模块305，与所述平衡模块304相连；用于根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
在本发明的一个实施例中，所述平衡模块具体用于：
使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构；其中，正则系综是组成系综的系统是由N个粒子组成的，同温度为T的很大的热源相接触并达到热平衡。选择退火功能后，可以进行300K到500K的5次循环，得到一系列结构，选择能量最小结构，通过动力学功能进行动力学平衡计算。
将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综下进行动力学计算，获取盒子密度；更具体地，是处于等温等压系综(NPT)下进行100-1000ps、100000步数的动力学计算。
将所述盒子密度与所述预设密度进行对比，获取偏差值；
判断所述偏差值是否大于预设阈值；上述预设阈值为5％至10％，例如是7％。
若否，则将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算。
若是，则将所述盒子密度设为预设密度；并返回所述建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定的步骤。
在本发明的一个实施例中，所述计算模块用于：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R；
通过以下公式计算特性黏数η：
Figure BDA0001479530990000141
Figure BDA0001479530990000142
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
在具体实施时，可以运用Forcite模块的分析功能计算盐水溶液中聚合物回转半径R(g)；更具体地，通过Forcite模块的分析功能，对NVT计算后动力学平衡状态的一系列数据进行分析，得到R(g)分布图，对其进行统计平均计算得到R(g)数值。
本发明提供了一种基于分子模拟的共聚物抗盐性能评价方法及其系统，利用分子模拟软件建立聚合物盐水溶液模型，再利用分子动力学计算系统达到平衡后的分子回转半径，通过回转半径计算得到特性粘数，通过比较特性粘数的大小可以评价不同聚合物抗盐能力。本发明无需通过传统实验，利用计算机“虚拟实验”进行评价，不受实验条件限制，效率高，重复性强，结果可靠，有利于油田化学用抗盐聚合物的快速设计和研发。
以上参照附图说明了本发明的优选实施例，本领域技术人员不脱离本发明的范围和实质，可以有多种变型方案实现本发明。举例而言，作为一个实施例的部分示出或描述的特征可用于另一实施例以得到又一实施例。以上仅为本发明较佳可行的实施例而已，并非因此局限本发明的权利范围，凡运用本发明说明书及附图内容所作的等效变化，均包含于本发明的权利范围之内。

Claims (8)

1.一种基于分子模拟的共聚物抗盐性能评价方法，其特征在于，包括：
构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；
构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；
建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定；
对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；
根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
2.根据权利要求1所述的基于分子模拟的共聚物抗盐性能评价方法，其特征在于，所述对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型包括：
使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构；
将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综下进行动力学计算，获取盒子密度；
将所述盒子密度与所述预设密度进行对比，获取偏差值；
判断所述偏差值是否大于预设阈值；
若否，则将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算；
其中，正则系综是由组成系综的系统与一温度为T的很大的热源相接触并达到热平衡后形成的。
3.根据权利要求2所述的基于分子模拟的共聚物抗盐性能评价方法，其特征在于：
若所述偏差值大于预设阈值，则将所述盒子密度设为预设密度；并返回所述建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定的步骤。
4.根据权利要求1所述的基于分子模拟的共聚物抗盐性能评价方法，其特征在于，所述根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数包括：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R：
通过以下公式计算特性黏数η：
Figure FDA0002725071020000011
Figure FDA0002725071020000012
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
5.一种基于分子模拟的共聚物抗盐性能评价系统，其特征在于，包括：
构建模块，用于构建水分子、无机盐分子以及聚合物单体分子的结构模型，并进行结构优化；
合成模块，用于构建以所述聚合物单体分子组成的聚合物分子的结构模型，并进行结构优化，获取聚合物分子最小能量构象结构；
盒子建立模块，用于建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定；
平衡模块，用于对所述聚合物盐水溶液模型盒子进行分子动力学计算，得到动力学平衡状态溶液体系结构模型；
计算模块，用于根据所述得到动力学平衡状态溶液体系结构模型计算特性黏数。
6.根据权利要求5所述的基于分子模拟的共聚物抗盐性能评价系统，其特征在于，所述平衡模块还用于：
使所述聚合物盐水溶液模型盒子处于正则系综下进行退火，获取聚合物盐水溶液模型盒子能量最小结构；
将所述聚合物盐水溶液模型盒子能量最小结构处于等温等压系综下进行动力学计算，获取盒子密度；
将所述盒子密度与所述预设密度进行对比，获取偏差值；
判断所述偏差值是否大于预设阈值；
若否，则将所述聚合物盐水溶液模型盒子能量最小结构处于正则系综下进行动力学计算；
其中，正则系综是由组成系综的系统与一温度为T的很大的热源相接触并达到热平衡后形成的。
7.根据权利要求6所述的基于分子模拟的共聚物抗盐性能评价系统，其特征在于，所述平衡模块还用于：
若所述偏差值大于预设阈值，则将所述盒子密度设为预设密度；并返回所述建立包含所述水分子、无机盐分子、以及所述聚合物分子最小能量构象结构的聚合物盐水溶液模型盒子，并按预设密度、预设温度对所述聚合物盐水溶液模型盒子的密度、温度进行设定的步骤。
8.根据权利要求5所述的基于分子模拟的共聚物抗盐性能评价系统，其特征在于，所述计算模块用于：
计算所述动力学平衡状态溶液体系结构模型中聚合物分子的回转半径R；
通过以下公式计算特性黏数η：
Figure FDA0002725071020000031
Figure FDA0002725071020000032
其中，NA为阿伏伽德罗常数，M为聚合物分子的分子量。
CN201711182823.0A 2017-11-23 2017-11-23 一种基于分子模拟的共聚物抗盐性能评价方法及其系统 Active CN110021370B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711182823.0A CN110021370B (zh) 2017-11-23 2017-11-23 一种基于分子模拟的共聚物抗盐性能评价方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711182823.0A CN110021370B (zh) 2017-11-23 2017-11-23 一种基于分子模拟的共聚物抗盐性能评价方法及其系统

Publications (2)

Publication Number Publication Date
CN110021370A CN110021370A (zh) 2019-07-16
CN110021370B true CN110021370B (zh) 2020-12-11

Family

ID=67186827

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711182823.0A Active CN110021370B (zh) 2017-11-23 2017-11-23 一种基于分子模拟的共聚物抗盐性能评价方法及其系统

Country Status (1)

Country Link
CN (1) CN110021370B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110556168B (zh) * 2019-09-11 2022-07-05 江西省科学院能源研究所 一种分析水溶性单体在凝胶膜中扩散性质的模拟方法
CN111028891A (zh) * 2019-12-12 2020-04-17 南京科技职业学院 一种基于分子动力学模拟烷基苯磺酸盐表面活性剂油水界面特性的评价方法
CN112802558B (zh) * 2021-01-29 2024-01-26 南京中科格特康科技有限公司 一种基于分子模拟预测用于盐湖提锂的共价有机框架性能的方法
CN113223635A (zh) * 2021-04-09 2021-08-06 中国石油大学(北京) 一种流体剪切黏度确定方法、装置和电子设备
CN113987897B (zh) * 2021-10-12 2023-04-14 广西大学 预测变压器用植物绝缘油物理化学性质的复合模型及应用

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105440211A (zh) * 2014-08-28 2016-03-30 中国石油化工股份有限公司 分子印迹膜的功能单体筛选方法及制备方法
CN106469251A (zh) * 2015-08-17 2017-03-01 杨立威 评估抗菌肽抗菌力的系统及其使用方法
CN106832241A (zh) * 2016-11-07 2017-06-13 北京服装学院 一种利用复合催化剂制备聚对苯二甲酸乙二醇酯的方法
CN107133458A (zh) * 2017-04-24 2017-09-05 西南石油大学 基于分子动力学模拟气体在钻井液中溶解特性的评价方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080040165A1 (en) * 2006-08-08 2008-02-14 Monsanto Technology Llc Transgenic crop financial systems and methods
US10672045B2 (en) * 2015-09-11 2020-06-02 Caterpillar Inc. Systems and methods for processing orders for structural designs

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105440211A (zh) * 2014-08-28 2016-03-30 中国石油化工股份有限公司 分子印迹膜的功能单体筛选方法及制备方法
CN106469251A (zh) * 2015-08-17 2017-03-01 杨立威 评估抗菌肽抗菌力的系统及其使用方法
CN106832241A (zh) * 2016-11-07 2017-06-13 北京服装学院 一种利用复合催化剂制备聚对苯二甲酸乙二醇酯的方法
CN107133458A (zh) * 2017-04-24 2017-09-05 西南石油大学 基于分子动力学模拟气体在钻井液中溶解特性的评价方法

Also Published As

Publication number Publication date
CN110021370A (zh) 2019-07-16

Similar Documents

Publication Publication Date Title
CN110021370B (zh) 一种基于分子模拟的共聚物抗盐性能评价方法及其系统
Wang et al. Progress report on phase separation in polymer solutions
Jang et al. Nanophase-segregation and transport in Nafion 117 from molecular dynamics simulations: effect of monomeric sequence
Carrillo et al. Detailed molecular dynamics simulations of a model NaPSS in water
Chen et al. Promoting effect of common marine cations on hydrate dissociation and structural evolution under a static electric field
Bain et al. Progress in Computer Simulation of Bulk, Confined, and Surface‐initiated Polymerizations
Dalla Bernardina et al. Mechanism of ionization, hydration, and intermolecular H-bonding in proton conducting nanostructured ionomers
Jeon Enhancing electrode wettability in lithium-ion battery via particle-size ratio control
Glagolev et al. Influence of cross-linking rate on the structure of hypercrosslinked networks: Multiscale computer simulation
Latinwo et al. Nonequilibrium work relations for polymer dynamics in dilute solutions
Zhang et al. Salt resistance study and molecular dynamics simulation of hydrophobic-association polymer with internal salt structure
Tang et al. Behaviors of water molecules in polyvinyl alcohol gel amid stretch and temperature changes: A molecular dynamics study
Ma et al. Artificial neural network prediction of metastable zone widths in reactive crystallization of lithium carbonate
Wu et al. Polymer networks by molecular dynamics simulation: Formation, thermal, structural and mechanical properties
Pereyra et al. Modeling polymeric gels: The role of chain flexibility on the structure of physical gels
CN116453603A (zh) 一种构建树枝状聚合物分子结构和耐盐性能定量构效关系的方法、应用
Fyta et al. Multiscale simulation of nanobiological flows
Kwon et al. The breakdown of the local thermal equilibrium approximation for a polymer chain during packaging
CN108984852B (zh) 外力作用下高分子链在表面吸附性质的蒙特卡罗模拟方法
Priezjev et al. Structural relaxation of porous glasses due to internal stresses and deformation under tensile loading at constant pressure
Glagolev et al. Reverse mapping algorithm for multi-scale numerical simulation of polylactic acid
Oh et al. Behaviour of water molecules in Nafion 117 for polymer electrolyte membrane fuel cell by molecular dynamics simulation
Jiang et al. Self-organization of multivalent counterions in polyelectrolyte brushes
Zhang et al. Dynamics simulation on the associative properties of amphiphilic functional monomer modified polyacrylamide copolymers
Zhang et al. Defining the optimal criterion for separating gases using polymeric membranes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant