CN113223635A - 一种流体剪切黏度确定方法、装置和电子设备 - Google Patents
一种流体剪切黏度确定方法、装置和电子设备 Download PDFInfo
- Publication number
- CN113223635A CN113223635A CN202110381562.5A CN202110381562A CN113223635A CN 113223635 A CN113223635 A CN 113223635A CN 202110381562 A CN202110381562 A CN 202110381562A CN 113223635 A CN113223635 A CN 113223635A
- Authority
- CN
- China
- Prior art keywords
- fluid
- data
- model
- shear viscosity
- carrying
- 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.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000000547 structure data Methods 0.000 claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 37
- 230000008878 coupling Effects 0.000 claims description 29
- 238000010168 coupling process Methods 0.000 claims description 29
- 238000005859 coupling reaction Methods 0.000 claims description 29
- 238000004088 simulation Methods 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 15
- 239000010779 crude oil Substances 0.000 claims description 14
- 239000000203 mixture Substances 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 238000012887 quadratic function Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 238000009472 formulation Methods 0.000 claims description 3
- 239000012071 phase Substances 0.000 description 31
- 230000003993 interaction Effects 0.000 description 19
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 14
- 238000006073 displacement reaction Methods 0.000 description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 10
- 230000006872 improvement Effects 0.000 description 8
- 229910002092 carbon dioxide Inorganic materials 0.000 description 7
- 239000001569 carbon dioxide Substances 0.000 description 7
- 239000003921 oil Substances 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 5
- 238000000329 molecular dynamics simulation Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000009878 intermolecular interaction Effects 0.000 description 4
- 238000011084 recovery Methods 0.000 description 4
- 239000008346 aqueous phase Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Business, Economics & Management (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Marketing (AREA)
- Chemical & Material Sciences (AREA)
- Animal Husbandry (AREA)
- Economics (AREA)
- Agronomy & Crop Science (AREA)
- Human Resources & Organizations (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Crystallography & Structural Chemistry (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本说明书实施例公开了一种流体剪切黏度确定方法、装置和电子设备。所述方法包括：获取至少一种组分的单相结构模型；根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；对所述体相结构模型进行结构优化，得到优化结构数据；对所述优化结构数据进行构型平衡处理，得到平衡结构数据；对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。本说明书实施例可以快速高效地确定流体的剪切黏度数据。
Description
技术领域
本说明书实施例涉及石油勘探开发技术领域，特别涉及一种流体剪切黏度确定方法、装置和电子设备。
背景技术
21世纪以来我国进入社会经济快速发展阶段，对能源需求量日益增加。近年我国石油能源自给率不足30％，石油供给安全形势严峻。保持传统油藏石油产量持续稳定和开采非常规油气藏对保障我国能源安全至关重要。利用化学驱技术，将外部流体(例如水，超临界CO2等)注入油藏储层孔隙中能够大幅度提高原油采收率，这项技术涉及到复杂孔隙流体剪切黏度的计算。因此，流体剪切黏度的计算，对于选择驱替流体的种类，评价孔隙中原油驱替效果，和设计原油开采工艺路线具有重大的理论和实践意义。
在相关技术中，可以通过室内实验测量流体的剪切黏度。但是通过实验测量流体黏度的过程，受操作条件(如温度和压强控制，人为误差等)、检测技术和经济成本等因素影响，难以快速高效地实现。
发明内容
本说明书实施例提供一种剪切黏度确定方法、装置和电子设备，以快速高效地确定流体的剪切黏度数据。本说明书实施例的技术方案如下。
本说明书实施例的第一方面，提供了一种剪切黏度确定方法，包括：
获取至少一种组分的单相结构模型；
根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；
对所述体相结构模型进行结构优化，得到优化结构数据；
对所述优化结构数据进行构型平衡处理，得到平衡结构数据；
对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；
对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
本说明书实施例的第二方面，提供了一种剪切黏度确定装置，包括：
获取单元，用于获取至少一种组分的单相结构模型；
搭建单元，用于根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；
第一处理单元，用于对所述体相结构模型进行结构优化，得到优化结构数据；
第二处理单元，用于对所述优化结构数据进行构型平衡处理，得到平衡结构数据；
第三处理单元，用于对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；
分析单元，用于对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
本说明书实施例的第三方面，提供了一种计算设备，包括：
至少一个处理器；
存储有程序指令的存储器，其中，所述程序指令被配置为适于由所述至少一个处理器执行，所述程序指令包括用于执行如第一方面所述方法的指令。
本说明书实施例提供的技术方案，可以采用分子动力学模拟方法确定流体的剪切黏度数据。相较于通过实验测量流体黏度，本说明书实施例可以将复杂繁琐的操作流程简单化，极大地减少了操作流程中的人为因素，能够大幅度节约实验成本，从而为原油开采选择驱替相流体种类，提高原油采收率和设计原油开采工艺路线提供理论和实践指导。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施例中剪切黏度确定方法的流程示意图；
图2为本说明书实施例中利用Berendsen耦合算法得到的水相剪切黏度与K值的一种拟合曲线图；
图3为本说明书实施例中利用Nose-Hoover耦合算法得到的水相剪切黏度与K值的另一种拟合曲线图；
图4为本说明书实施例中剪切黏度确定装置的结构示意图；
图5为本说明书实施例中计算设备的结构示意图。
具体实施方式
下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。
剪切黏度(shear viscosity)，是指稳流状态下剪切应力与剪切速率之比。剪切黏度是液体分子内摩擦的量度，也是物体黏流性质的一项具体反映。促使原油膨胀降低剪切黏度是驱油的重要机理。因此计算剪切黏度对筛选驱替剂，预测油藏开发时的原油驱替效率具有重要意义。
分子动力学模拟方法，是指利用牛顿运动定律，从微观角度研究分子之间相互作用过程。利用分子动力学模拟方法，能准确控制流体的温度、压强等条件，快速计算复杂孔隙流体的剪切黏度。为此，本说明书实施例提供一种流体剪切黏度确定方法。所述剪切黏度确定方法可以基于分子动力学模拟方法实现。所述剪切黏度确定方法可以在台式机、笔记本电脑、平板电脑、移动智能电话、服务器、或者服务器集群等计算设备上实现。
请参阅图1，所述剪切黏度确定方法可以包括以下步骤。
步骤S101：获取至少一种组分的单相结构模型。
在一些实施例中，所述组分可以为流体中的成分。所述流体包括储层孔隙中用于驱替原油的流体。当然，所述流体还可以包括其它的流体，例如化学聚合物等。所述组分可以包括二氧化碳分子、水分子等。所述单相结构模型可以用于表示组分的结构。
在一些实施例中，可以根据流体中单一组分的重复单元，构建至少一种组分的单相结构模型。例如，可以利用Materials Studio软件构建二氧化碳分子、水分子等组分的单相结构模型；可以将构建的单相结构模型转化为分子模拟软件所需要的输入文件。其中，所述分子模拟软件可以包括Gromacs软件等。
步骤S103：根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型。
在一些实施例中，所述配方可以包括所需流体中各种组分的比例。所述体相结构模型可以用于表示流体的结构。
在一些实施例中，可以构建设定大小的模型盒子；可以根据模型盒子的大小和流体中各种组分的比例，确定各种组分的分子数；可以根据组分的分子数，向模型盒子中置入由组分的单向结构模型表示的分子。例如，可以利用Gromacs软件构建大小为7×7×20nm的模型盒子；可以向模型盒子中置入10000个水分子、3000个二氧化碳分子，得到流体的体相结构模型。
步骤S105：对所述体相结构模型进行结构优化，得到优化结构数据。
在一些实施例中，可以在设定的第一模拟参数下，对所述体相结构模型进行能量最小化处理，得到优化结构数据。所述第一模拟参数包括以下至少之一：分子力场、能量限度。在实际中可以采用最速下降算法将所述体相结构模型的能量最小化，得到优化结构数据。
例如，可以利用Gromacs软件，设置分子力场为Charmm36-jul2017.ff，设置能量限度为800kJ mol-1nm-1；采用最速下降算法将所述体相结构模型的能量最小化，得到优化结构数据。
步骤S107：对所述优化结构数据进行构型平衡处理，得到平衡结构数据。
在一些实施例中，可以在设定的第二模拟参数下，对所述优化结构数据进行等温等压系综(NPT系综)模拟，以进行充分的构型和密度充分平衡，得到平衡结构数据。所述第二模拟参数包括以下至少之一：分子间的相互作用、温度、温度耦合算法、压力、压力耦合算法、截断半径、时间步长、模拟时长。其中，分子间的相互作用可以包括相同分子间的相互作用和不同分子间的相互作用。
例如，可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Berendsen算法，设置压力为20Mpa，设置压力耦合算法为Berendsen算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为5ns；对所述优化结构数据进行NPT系综模拟，以进行充分的构型和密度充分平衡，得到平衡结构数据。其中，分子间的库伦相互作用可以采用PME(Particle-MeshEwald)方法计算。
步骤S109：对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据。
在一些实施例中，可以在设定的第三模拟参数下，对平衡结构数据进行正则系综(体积和温度不变，NVT系综)模拟，以进行充分的动力学平衡，得到动力学轨迹数据。所述第三模型参数包括以下至少之一：分子间的相互作用、温度、温度耦合算法、截断半径、时间步长、模拟时长。其中，分子间的相互作用可以包括相同分子间的相互作用和不同分子间的相互作用。在实际应用中，可以采用一种温度耦合算法，对所述平衡结构数据进行动力学平衡处理，得到一种动力学轨迹数据。例如，可以采用Berendsen温度耦合算法或者Nose-Hoover温度耦合算法，对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据。或者，还可以采用两种温度耦合算法，对所述平衡结构数据进行动力学平衡处理，得到两种动力学轨迹数据。所述两种动力学轨迹数据用于确定流体的剪切黏度的上限值和下限值。例如，可以分别采用Berendsen温度耦合算法和Nose-Hoover温度耦合算法，对所述平衡结构数据进行动力学平衡处理，得到两种动力学轨迹数据。其中，通过Berendsen温度耦合算法获得的动力学轨迹数据，可以用于确定流体的剪切黏度上限值。通过Nose-Hoover温度耦合算法获得的动力学轨迹数据，可以用于确定流体的剪切黏度下限值。
例如，可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Nose-Hoover算法或者Berendsen算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为15ns；对所述平衡结构数据进行NVT系综模拟，以进行充分的动力学平衡，得到最后5ns的动力学轨迹数据。其中，分子间的库伦相互作用可以采用PME(Particle-Mesh Ewald)方法计算。
步骤S111：对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
在一些实施例中，所述流体的剪切黏度数据可以包括流体的剪切黏度上限值和/或流体的剪切黏度下限值。具体地，可以对动力学轨迹数据进行横向电流自相关函数分析，得到流体的剪切黏度随k的变化曲线；可以对所述变化曲线进行二次函数拟合，得到流体的剪切黏度数据。其中，所述横向电流可以指横向扰动波计算方法。所述横向电流自相关函数分析，可以是指通过横向扰动波计算方法进行自相关函数分析。所述k可以表示横向动量的波矢量。
值得说明的是，在对动力学轨迹进行分析的过程中，若最终体系混相形成单相，则可以计算体系内所有分子的横向电流自相关函数。若最终体系分相形成多个相，则可以计算指定体相范围内所有分子的横向电流自相关函数。
例如，可以采用Gromacs软件，处理动力学轨迹数据，计算密度分布曲线。所述密度分布曲线可以用于确定所求流体相的区域位置。因此，根据所述密度分布曲线，可以构建所要计算的水相位置范围区域内所有分子的索引组。可以采用Gromacs软件，处理索引组分子平衡过程轨迹数据，可以计算横向电流自相关函数，得到黏度与k的数据文件。可以采用origin软件对所述数据文件进行二次函数拟合，拟合曲线的y轴截距即为剪切黏度值。
在一些实施例中，通过步骤S109，可以获得一个动力学轨迹数据。通过对动力学轨迹数据进行自相关函数分析，可以得到流体的剪切黏度数据。所述流体的剪切黏度数据可以包括流体的剪切黏度上限值或者下限值。或者，通过步骤S109，可以获得两个动力学轨迹数据。通过对所述两个动力学轨迹数据进行自相关函数分析，可以得到流体的剪切黏度数据。所述流体的剪切黏度数据可以包括流体的剪切黏度上限值和下限值。
本说明书实施例的剪切黏度确定方法，可以采用分子动力学模拟方法确定流体的剪切黏度数据。相较于通过实验测量流体黏度，本说明书实施例可以将复杂繁琐的操作流程简单化，极大地减少了操作流程中的人为因素，能够大幅度节约实验成本，从而为原油开采选择驱替流体种类，提高原油采收率和设计原油开采工艺路线提供理论和实践指导。
具体地，本说明书实施例可以具有以下有益效果：
1、通过引入模拟仿真减少实验次数，加快研发进度；
2、仿真建模的便捷性使本说明书实施例适合复杂的流体配方；
3、操作流程明确，简单易用，降低操作门槛，能够快速推广；
4、操作流程化，确保仿真结果的可靠性和重复性；
5、模拟仿真控制条件的准确性和任意性使本说明书实施例能应用广泛。
6、有利于优选高效的驱替流体，提高原油采收率。
在一个场景示例中，在求解水相黏度值上限值时，可以利用Materials Studio软件构建二氧化碳分子、水分子等组分的单相结构模型；可以利用Gromacs软件构建大小为7×7×20nm的模型盒子；可以向模型盒子中置入10000个水分子、3000个二氧化碳分子，得到流体的体相结构模型；可以利用Gromacs软件，设置分子力场为Charmm36-jul2017.ff，设置能量限度为800kJ mol-1nm-1；采用最速下降算法将所述体相结构模型的能量最小化，得到优化结构数据；可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Berendsen温度耦合算法，设置压力为20Mpa，设置压力耦合算法为Berendsen压力耦合算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为5ns；对所述优化结构数据进行NPT系综模拟，以进行充分的构型和密度充分平衡，得到平衡结构数据；可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Berendsen温度耦合算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为15ns；对所述平衡结构数据进行NVT系综模拟，以进行充分的动力学平衡，得到最后5ns的动力学轨迹数据；可以对轨迹数据进行横向电流自相关函数分析，得到如图2所示的剪切黏度随k的变化曲线；可以对曲线进行二次函数拟合，得到待测流体的剪切黏度的上限值。
在另一个场景示例中，在求解水相黏度值下限值时，可以利用Materials Studio软件构建二氧化碳分子、水分子等组分的单相结构模型；可以利用Gromacs软件构建大小为7×7×20nm的模型盒子；可以向模型盒子中置入10000个水分子、3000个二氧化碳分子，得到流体的体相结构模型；可以利用Gromacs软件，设置分子力场为Charmm36-jul2017.ff，设置能量限度为800kJ mol-1nm-1；采用最速下降算法将所述体相结构模型的能量最小化，得到优化结构数据；可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Berendsen温度耦合算法，设置压力为20Mpa，设置压力耦合算法为Berendsen压力耦合算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为5ns；对所述优化结构数据进行NPT系综模拟，以进行充分的构型和密度充分平衡，得到平衡结构数据；可以利用Gromacs软件，设置相同分子间的相互作用和不同分子间的相互作用均为L-J参数，设置温度为323.15K，设置温度耦合算法为Nose-Hoover温度耦合算法，设置截断半径为1.1nm，设置时间步长为1fs，设置模拟时长为15ns；对所述平衡结构数据进行NVT系综模拟，以进行充分的动力学平衡，得到最后5ns的动力学轨迹数据；可以对轨迹数据进行横向电流自相关函数分析，得到如图3所示的剪切黏度随k的变化曲线；可以对曲线进行二次函数拟合，得到待测流体的剪切黏度的下限值。
请参阅图4。本说明书实施例提供一种流体剪切黏度确定装置。所述装置可以包括以下单元。
获取单元201，用于获取至少一种组分的单相结构模型；
搭建单元203，用于根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；
第一处理单元205，用于对所述体相结构模型进行结构优化，得到优化结构数据；
第二处理单元207，用于对所述优化结构数据进行构型平衡处理，得到平衡结构数据；
第三处理单元209，用于对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；
分析单元211，用于对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
下面介绍本说明书电子设备的一个实施例。图5是该实施例中电子设备的硬件结构示意图。如图5所示，该电子设备可以包括一个或多个(图中仅示出一个)处理器、存储器和传输模块。当然，本领域普通技术人员可以理解，图5所示的硬件结构仅为示意，其并不对上述电子设备的硬件结构造成限定。在实际中该电子设备还可以包括比图5所示更多或者更少的组件单元；或者，具有与图5所示不同的配置。
所述存储器可以包括高速随机存储器；或者，还可以包括非易失性存储器，例如一个或者多个磁性存储装置、闪存或者其他非易失性固态存储器。当然，所述存储器还可以包括远程设置的网络存储器。所述远程设置的网络存储器可以通过诸如互联网、企业内部网、局域网、移动通信网等网络连接至所述区块链客户端。所述存储器可以用于存储应用软件的程序指令或模块，例如用于实现本说明书图1或者图2所对应实施例的程序指令或模块。
所述处理器可以按任何适当的方式实现。例如，所述处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific IntegratedCircuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。所述处理器可以读取并执行所述存储器中的程序指令或模块。
所述传输模块可以用于经由网络进行数据传输，例如经由诸如互联网、企业内部网、局域网、移动通信网等网络进行数据传输。
本说明书还提供计算机存储介质的一个实施例。所述计算机存储介质包括但不限于随机存取存储器(Random Access Memory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard Disk Drive,HDD)、存储卡(Memory Card)等等。所述计算机存储介质存储有计算机程序指令。在所述计算机程序指令被执行时实现：本说明书图1或者图2所对应的实施例。
需要说明的是，本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其它实施例的不同之处。另外，可以理解的是，本领域技术人员在阅读本说明书文件之后，可以无需创造性劳动想到将本说明书列举的部分或全部实施例进行任意组合，这些组合也在本说明书公开和保护的范围内。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本说明书各个实施例或者实施例的某些部分所述的方法。
本说明书可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本说明书，本领域普通技术人员知道，本说明书有许多变形和变化而不脱离本说明书的精神，希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。
Claims (10)
1.一种流体剪切黏度确定方法，包括：
获取至少一种组分的单相结构模型；
根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；
对所述体相结构模型进行结构优化，得到优化结构数据；
对所述优化结构数据进行构型平衡处理，得到平衡结构数据；
对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；
对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
2.如权利要求1所述的方法，所述配方包括所需流体中各种组分的比例；
所述搭建流体的体相结构模型，包括：
构建设定大小的模型盒子；
根据模型盒子的大小和流体中各种组分的比例，确定各种组分的分子数；
根据组分的分子数，向模型盒子中置入由组分的单相结构模型表示的分子。
3.如权利要求1所述的方法，所述对所述体相结构模型进行结构优化，包括：
在设定的第一模拟参数下，对所述体相结构模型进行能量最小化处理。
4.如权利要求1所述的方法，所述对所述优化结构数据进行构型平衡处理，包括：
在设定的第二模拟参数下，对所述优化结构数据进行NPT系综模拟。
5.如权利要求1所述的方法，所述对所述平衡结构数据进行动力学平衡处理，包括：
在设定的第三模拟参数下，对所述平衡结构数据进行NVT系综模拟。
6.如权利要求1所述的方法，所述对所述平衡结构数据进行动力学平衡处理，包括：
采用两种温度耦合算法，对所述平衡结构数据进行动力学平衡处理，得到两种动力学轨迹数据，所述两种动力学轨迹数据用于确定流体的剪切黏度的上限值和下限值。
7.如权利要求1所述的方法，所述对所述动力学轨迹数据进行自相关函数分析，包括：
对动力学轨迹数据进行横向电流自相关函数分析，得到流体的黏度随k的变化曲线；对所述变化曲线进行二次函数拟合，得到流体的剪切黏度数据；所述k表示横向动量的波矢量。
8.如权利要求1所述的方法，所述流体包括储层孔隙中用于驱替原油的流体。
9.一种流体剪切黏度确定装置，包括：
获取单元，用于获取至少一种组分的单相结构模型；
搭建单元，用于根据所需流体的配方，利用组分的单相结构模型，搭建流体的体相结构模型；
第一处理单元，用于对所述体相结构模型进行结构优化，得到优化结构数据；
第二处理单元，用于对所述优化结构数据进行构型平衡处理，得到平衡结构数据；
第三处理单元，用于对所述平衡结构数据进行动力学平衡处理，得到动力学轨迹数据；
分析单元，用于对所述动力学轨迹数据进行自相关函数分析，得到流体的剪切黏度数据。
10.一种电子设备，包括：
至少一个处理器；
存储有程序指令的存储器，其中，所述程序指令被配置为适于由所述至少一个处理器执行，所述程序指令包括用于执行如权利要求1-8中任一项所述方法的指令。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110381562.5A CN113223635A (zh) | 2021-04-09 | 2021-04-09 | 一种流体剪切黏度确定方法、装置和电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110381562.5A CN113223635A (zh) | 2021-04-09 | 2021-04-09 | 一种流体剪切黏度确定方法、装置和电子设备 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113223635A true CN113223635A (zh) | 2021-08-06 |
Family
ID=77086835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110381562.5A Pending CN113223635A (zh) | 2021-04-09 | 2021-04-09 | 一种流体剪切黏度确定方法、装置和电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113223635A (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014106554A (ja) * | 2012-11-22 | 2014-06-09 | Toyo Tire & Rubber Co Ltd | 輸送係数の算出方法、輸送係数の算出装置、及び輸送係数の算出プログラム |
| CN106596336A (zh) * | 2016-12-13 | 2017-04-26 | 东北石油大学 | 黏弹性聚合物在多孔介质中弹性湍流特性测定装置及方法 |
| CN106934084A (zh) * | 2015-12-30 | 2017-07-07 | 中国石油天然气股份有限公司 | 一种带底油凝析气藏的相态拟合方法 |
| CN107290506A (zh) * | 2017-07-28 | 2017-10-24 | 中国石油大学(北京) | 一种定量评价储层成岩演化过程中孔隙度时空演变的方法 |
| JP2018049583A (ja) * | 2016-09-14 | 2018-03-29 | 住友重機械工業株式会社 | シミュレーション方法、コンピュータプログラム、及びシミュレーション装置 |
| US20180156902A1 (en) * | 2015-06-22 | 2018-06-07 | Board Of Trustees Of Michigan State University | Shear viscosity imaging with acoustic radiation force |
| CN108829919A (zh) * | 2018-05-02 | 2018-11-16 | 嘉兴学院 | 多尺度模拟核壳结构多相电流体动力学的方法 |
| CN108959844A (zh) * | 2018-06-12 | 2018-12-07 | 苏州创腾软件有限公司 | 采用分子模拟方法评价聚合物材料流动性能的方法及系统 |
| CN109883649A (zh) * | 2019-03-14 | 2019-06-14 | 北京石油化工学院 | 通过模拟模型研究纳米流体在纳米通道内流动行为的方法 |
| CN110021370A (zh) * | 2017-11-23 | 2019-07-16 | 中国石油化工股份有限公司 | 一种基于分子模拟的共聚物抗盐性能评价方法及其系统 |
| CN110055046A (zh) * | 2018-01-19 | 2019-07-26 | 中国石油天然气股份有限公司 | 高温高盐油藏注入水配制的高黏弹性聚合物驱油剂及制备方法 |
| US20190375976A1 (en) * | 2018-06-06 | 2019-12-12 | Elementis Specialties, Inc. | Oil based drilling fluids useful having reduced change in high shear viscosity over a wide temperature range |
| CN112214935A (zh) * | 2019-06-24 | 2021-01-12 | 中国石油化工股份有限公司 | 多孔介质中流体相态的计算方法及系统 |
| CN112380790A (zh) * | 2020-11-11 | 2021-02-19 | 桂林电子科技大学 | 一种emd方法计算不同形状的纳米颗粒的纳米流体粘度 |
-
2021
- 2021-04-09 CN CN202110381562.5A patent/CN113223635A/zh active Pending
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014106554A (ja) * | 2012-11-22 | 2014-06-09 | Toyo Tire & Rubber Co Ltd | 輸送係数の算出方法、輸送係数の算出装置、及び輸送係数の算出プログラム |
| US20180156902A1 (en) * | 2015-06-22 | 2018-06-07 | Board Of Trustees Of Michigan State University | Shear viscosity imaging with acoustic radiation force |
| CN106934084A (zh) * | 2015-12-30 | 2017-07-07 | 中国石油天然气股份有限公司 | 一种带底油凝析气藏的相态拟合方法 |
| JP2018049583A (ja) * | 2016-09-14 | 2018-03-29 | 住友重機械工業株式会社 | シミュレーション方法、コンピュータプログラム、及びシミュレーション装置 |
| CN106596336A (zh) * | 2016-12-13 | 2017-04-26 | 东北石油大学 | 黏弹性聚合物在多孔介质中弹性湍流特性测定装置及方法 |
| CN107290506A (zh) * | 2017-07-28 | 2017-10-24 | 中国石油大学(北京) | 一种定量评价储层成岩演化过程中孔隙度时空演变的方法 |
| CN110021370A (zh) * | 2017-11-23 | 2019-07-16 | 中国石油化工股份有限公司 | 一种基于分子模拟的共聚物抗盐性能评价方法及其系统 |
| CN110055046A (zh) * | 2018-01-19 | 2019-07-26 | 中国石油天然气股份有限公司 | 高温高盐油藏注入水配制的高黏弹性聚合物驱油剂及制备方法 |
| CN108829919A (zh) * | 2018-05-02 | 2018-11-16 | 嘉兴学院 | 多尺度模拟核壳结构多相电流体动力学的方法 |
| US20190375976A1 (en) * | 2018-06-06 | 2019-12-12 | Elementis Specialties, Inc. | Oil based drilling fluids useful having reduced change in high shear viscosity over a wide temperature range |
| CN108959844A (zh) * | 2018-06-12 | 2018-12-07 | 苏州创腾软件有限公司 | 采用分子模拟方法评价聚合物材料流动性能的方法及系统 |
| CN109883649A (zh) * | 2019-03-14 | 2019-06-14 | 北京石油化工学院 | 通过模拟模型研究纳米流体在纳米通道内流动行为的方法 |
| CN112214935A (zh) * | 2019-06-24 | 2021-01-12 | 中国石油化工股份有限公司 | 多孔介质中流体相态的计算方法及系统 |
| CN112380790A (zh) * | 2020-11-11 | 2021-02-19 | 桂林电子科技大学 | 一种emd方法计算不同形状的纳米颗粒的纳米流体粘度 |
Non-Patent Citations (5)
| Title |
|---|
| JIHUI JIA 等: "Molecular Dynamics Study on CO2 Foam Films with Sodium Dodecyl Sulfate: Effects of Surfactant Concentration, Temperature, and Pressure on the Interfacial Tension", 《ENERGY FUELS》 * |
| PING XUE 等: "Molecular dynamics simulation of thickening mechanism of supercritical CO2 thickener", 《CHEMICAL PHYSICS LETTERS》 * |
| 倪海欧 等: "NaNO3-KNO3-NaNO2三元混合相变熔盐结构与物性的分子动力学模拟", 《储能科学与技术》 * |
| 薛萍: "分子动力学模拟研究CO2_EOR技术", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
| 赵欣 等: "基于分子动力学模拟的流体剪切力学特性分析", 《北京理工大学学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rongy et al. | Mixing from Fickian diffusion and natural convection in binary non‐equilibrium fluid phases | |
| Orozco et al. | Optimization of intermolecular potential parameters for the CO2/H2O mixture | |
| Horowitz et al. | A concurrent efficient global optimization algorithm applied to polymer injection strategies | |
| Yang et al. | Three-dimensional computational fluid dynamics simulation of valve-induced water hammer | |
| Peng et al. | Study on force schemes in pseudopotential lattice Boltzmann model for two‐phase flows | |
| Lorentzen et al. | Estimation of production rates with transient well-flow modeling and the auxiliary particle filter | |
| Yu et al. | GPU accelerated lattice Boltzmann simulation for rotational turbulence | |
| Yan et al. | Calculation of minimum miscibility pressure using fast slimtube simulation | |
| Hu et al. | Prediction of interfacial structure and tension of binary mixtures containing carbon dioxide | |
| Dias et al. | Strain-injection and crack-path field techniques for 3D crack-propagation modelling in quasi-brittle materials | |
| Yuan et al. | Research on theoretical model of dynamic bulk modulus of gas-containing hydraulic oil | |
| Henderson et al. | Modeling and analysis of the isothermal flash problem and its calculation with the simulated annealing algorithm | |
| Asadi et al. | Central-moments-based lattice Boltzmann for associating fluids: A new integrated approach | |
| Zhang et al. | Modeling and simulation of pigging for a gas pipeline using a bypass pig | |
| CN113223635A (zh) | 一种流体剪切黏度确定方法、装置和电子设备 | |
| Wang et al. | Investigation on friction models for simulation of pipeline filling transients | |
| Wang et al. | Numerical simulation of one-dimensional two-phase flow using a pressure-based algorithm | |
| Mu et al. | Mass-conserved density gradient theory model for nucleation process | |
| Al-Arfaj | Estimation of the normal boiling points of water and 1, 2-dimethoxyethane through the viscosity–temperature dependence in the corresponding binary mixtures | |
| Sharafi et al. | Mathematical model for prediction of oil recovery of core-flood tests in process of viscoelastic polymer flooding | |
| Tsypkin et al. | Linear stability of a filtration flow with gas–oil interface within the brinkman approach | |
| Liang et al. | Modeling low Reynolds number incompressible flows with curved boundaries using SPH | |
| Civan | Including nonequilibrium relaxation in models for rapid multiphase flow in wells | |
| Shiri et al. | The effects of grain geometry on waterflooding and viscous fingering in micro-fractures and porous media from a lattice Boltzmann method study | |
| Khait et al. | A GPU-based integrated simulation framework for modelling of complex subsurface applications |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210806 |