CN109378041B - Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof - Google Patents

Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof Download PDF

Info

Publication number
CN109378041B
CN109378041B CN201811467326.XA CN201811467326A CN109378041B CN 109378041 B CN109378041 B CN 109378041B CN 201811467326 A CN201811467326 A CN 201811467326A CN 109378041 B CN109378041 B CN 109378041B
Authority
CN
China
Prior art keywords
beads
petroleum
ring
heteroatom
average molecular
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
CN201811467326.XA
Other languages
Chinese (zh)
Other versions
CN109378041A (en
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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201811467326.XA priority Critical patent/CN109378041B/en
Publication of CN109378041A publication Critical patent/CN109378041A/en
Application granted granted Critical
Publication of CN109378041B publication Critical patent/CN109378041B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a mesoscale simulation petroleum molecule coarse graining model and a construction method and a prediction method thereof. The construction method comprises the following steps: firstly, obtaining average molecular structures of petroleum samples with different component types, sequentially judging whether the structures contain structures such as fused aromatic ring system atoms, heteroatom ring system atoms, naphthenic ring system atoms, side chain atoms and the like, and replacing by adopting equivalent replacement beads to finally construct a mesoscale simulation petroleum molecule coarse graining model. The prediction method comprises the following steps: obtaining solubility parameters of the replacement beads by molecular dynamics simulation, and obtaining conservative force parameters by associating the Flory-Huggins theory with the dissipative particle dynamics theory; a mesoscale simulated petroleum molecule coarse graining model is utilized, and a dissipation particle dynamics method is combined for simulation to obtain a stable equilibrium state molecular dispersion system structure; the microstructure properties of the petroleum sample are predicted through statistical analysis, and the method has a good prediction effect particularly on the structural properties which are difficult to observe.

Description

Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof
Technical Field
The invention belongs to the field of petroleum processing, and particularly relates to a mesoscale simulated petroleum molecule coarse graining model and a construction method and a prediction method thereof.
Background
The colloidal dispersion system on the microscopic scale of petroleum has important significance in the processes of petroleum exploitation, storage and transportation and processing, and is also a key factor influencing the stability of petroleum products. Petroleum is a complex mixture composed of various hydrocarbon compounds and non-hydrocarbon compounds, and the macroscopic stability of the petroleum is the comprehensive expression of molecular aggregation behavior on the molecular level and is closely related to the chemical composition and relative content of the petroleum.
The molecular composition of a petroleum system is complex and difficult to measure, and structural properties on a plurality of micro-molecular levels cannot be obtained through experiments, so that the microstructure of the petroleum system needs to be researched by adopting a molecular simulation means. The traditional molecular simulation method such as molecular dynamics has accurate simulation result but small scale, and can not complete the calculation of molecular motion behavior of a huge petroleum molecular system, so that a dissipative particle dynamics method on mesoscopic scale is needed to simulate the petroleum system. Even so, currently, a suitable method for establishing a medium-scale coarse-grained model of petroleum molecules is still lacked, and cannot be linked with a composition model in molecular management.
Disclosure of Invention
In order to solve the problems, the invention aims to provide a method for constructing a mesoscale simulation petroleum molecule coarse graining model; the invention also aims to provide a mesoscale simulation petroleum molecule coarse graining model constructed by the construction method; it is also an object of the present invention to provide a method for predicting microstructural properties of petroleum molecules.
The purpose of the invention is realized by the following technical means:
on one hand, the invention provides a construction method of a mesoscale simulation petroleum molecule coarse graining model, which comprises the following steps:
acquiring average molecular structures of petroleum samples with different component types;
judging whether the average molecular structure of the petroleum sample contains a six-membered aromatic ring structure, if so, positioning the six-membered aromatic ring structure in the center of the condensed aromatic structure, and replacing the six-membered aromatic ring structure with an equivalent replacement bead; if not, judging whether a heteroatom structure is contained;
judging whether the average molecular structure of the petroleum sample contains a binary aromatic ring structure and a quaternary aromatic ring structure, if so, dividing a four-membered ring structure fragment and a binary ring structure fragment at intervals around a replacement bead of a six-membered aromatic ring structure, replacing the quaternary ring structure and the binary ring structure with an equivalent replacement bead, and then judging whether the average molecular structure of the petroleum sample contains a heteroatom structure; if not, directly judging whether the structure contains a heteroatom structure;
step four, carrying out the step two or the step three, judging whether the average molecular structure of the petroleum sample contains a heteroatom structure, if so, replacing the heteroatom structure by adopting equivalent replacement beads, and judging whether the average molecular structure contains a naphthenic ring structure; if not, directly judging whether the naphthenic ring structure is contained;
step five, carrying out step four, judging whether the average molecular structure of the petroleum sample contains a naphthenic ring structure, if so, replacing the naphthenic ring structure by adopting equivalent replacement beads, and then replacing the side chain structure by adopting equivalent replacement beads with side chain structures; if not, directly adopting equivalent replacement beads of the side chain structure to replace the side chain structure, and obtaining the medium-scale simulation petroleum molecule coarse graining model.
In the above construction method, preferably, the replacement beads include aromatic ring beads, cycloalkane ring beads, heteroatom beads, and side chain beads;
the aromatic ring beads include six-membered aromatic ring beads (A6 beads), four-membered aromatic ring beads (A4 beads), and two-membered aromatic ring beads (A2 beads);
the cycloalkane ring beads include six-membered cycloalkane ring beads (N6 beads), five-membered cycloalkane ring beads (N5 beads), quaternary cycloalkane ring beads (N4 beads), ternary cycloalkane ring beads (N3 beads), and binary cycloalkane ring beads (N2 beads);
the heteroatom beads include sulfur heteroatom-containing beads (S beads), oxygen heteroatom-containing beads (O beads), nitrogen heteroatom six-membered ring beads (NI6 beads), and nitrogen heteroatom five-membered ring beads (NI5 beads);
the side chain beads include side chain two consecutive carbon atom beads (R2 beads) and side chain three consecutive carbon atom beads (R3 beads).
In the method of the present invention, the assignment and replacement processes of replacement beads are respectively in the following order: replacing a fused aromatic ring system atom, replacing a heteroatom ring system atom, replacing a cycloalkane ring system atom, and replacing a side chain atom;
replacement of the fused aromatic ring system atoms is locking the a6 beads at the center of the fused aromatic ring system, and then replacing the a4 and a2 beads at intervals therearound;
the replacement of the heteroatom ring system atom is to replace the heteroatom to be a heteroatom bead on the basis of the constructed coarse grained structure of the condensed aromatic ring system;
the replacement of the ring system atoms of the naphthenic rings is to replace different segments of the naphthenic ring systems containing carbon atoms to corresponding beads of the naphthenic rings on the basis of the constructed coarse grained structure of the condensed aromatic ring systems containing heteroatoms;
the replacement of the side chain atom is to replace different side chain fragments containing carbon atoms to corresponding side chain beads on the basis of the constructed coarse grained structure of the condensed aromatic ring system containing the heteroatom and the naphthenic ring atom.
In the above construction method, preferably, the specific method for obtaining the average molecular structure of the petroleum samples of different types is as follows:
determining the element content in the petroleum sample; determining structural parameters of the petroleum sample; and constructing the average molecular structure of the petroleum sample.
In the above construction method, preferably, the component types of the petroleum sample include one or more of asphaltene, pectin, aromatic components, and saturated components.
In the above construction method, preferably, the element includes one or more of carbon, hydrogen, oxygen, nitrogen, and sulfur.
In the above construction method, preferably, the structural parameter includes one or more of aromatic carbon ratio, unit aromatic ring number, unit aromatic carbon number, unit naphthenic ring number, nitrogen atom number, oxygen atom number, sulfur atom number, unit naphthenic carbon number, and unit branched chain carbon number.
In the above-described construction method, preferably, the element content includes one or more of a carbon mass fraction, a hydrogen mass fraction, a sulfur mass fraction, a nitrogen mass fraction, and an oxygen mass fraction.
In the above construction method, the specific operation of obtaining the average molecular structure of different types of petroleum samples is disclosed in patent application publication No. CN107704720A (a construction method of a petroleum average molecular structure model and a property prediction method), which is incorporated herein by reference in its entirety.
On the other hand, the invention also provides the mesoscale simulated petroleum molecule coarse graining model constructed by the construction method.
In yet another aspect, the present invention also provides a method for predicting microstructural properties of petroleum molecules, comprising the steps of:
step 1, obtaining solubility parameters of replacement beads by adopting molecular dynamics simulation, and associating Flory-Huggins theory and dissipative particle dynamics theory to obtain conservative force parameters;
step 2, combining the mesoscale simulated petroleum molecule coarse graining model constructed by the construction method and the obtained conservative force parameter with a dissipative particle dynamics method, and obtaining a stable equilibrium state molecular dispersion system structure through dissipative particle dynamics simulation; the microstructure properties of the petroleum sample are predicted by statistical analysis.
In the above method, preferably, the microstructure property includes one or more of a number of sub-aggregates, a molecular aggregation rate, and a molecular aggregation morphology.
According to the scheme provided by the invention, various types of structural parameters can be obtained by adopting a nuclear magnetic resonance method which is relatively simple and convenient to operate, the mass fraction of each element is obtained through element analysis, the true molecular structure information is obtained through an atomic force microscope, the average molecular structure of the petroleum sample is calculated by taking the data as original data (see the specific operation process), and a mesoscale simulation petroleum molecule coarse graining model is constructed. And (3) obtaining a bead solubility parameter by adopting molecular dynamics simulation, associating the Flory-Huggins theory with the dissipative particle dynamics theory to obtain a conservative force parameter, and carrying out dissipative particle dynamics simulation on the basis to predict the microstructure property of the petroleum sample. The division concept of guiding the structure in the molecular management field to be lumped is introduced into the mesoscale simulation, so that the direct cutting of molecules is avoided, and the dissipative particle dynamics method is directly butted with the molecular management. The field of molecular management refers to the study of the composition and properties of petroleum at the molecular level; the division concept of structure-oriented lump refers to the use of 22 structure vectors to clearly describe the structure of petroleum molecules, a very efficient method for dividing petroleum radicals and representing petroleum molecules.
Compared with the prior art, the scheme of the invention has the following advantages:
(1) the division concept of guiding the structure in the molecular management field to be lumped is introduced into the mesoscale simulation, so that the direct cutting of molecules is avoided, and the dissipative particle dynamics method is directly butted with the molecular management.
(2) Starting from a mesoscale simulation petroleum molecule coarse graining model and conservative force parameters, some microstructure properties which are difficult to measure, such as molecule aggregation number, molecule aggregation morphology and the like, of a petroleum sample can be predicted.
(3) The scheme of the invention can predict the micro equilibrium state configuration of the petroleum sample and the related solution system thereof. Only a small amount of experimental data are needed in the whole process, structural data which cannot be obtained in the experimental process can be obtained, and the method is small in workload, low in cost and fast.
Drawings
FIG. 1 is a flowchart of a method for constructing a mesoscale simulation petroleum molecular coarse graining model in example 1;
FIG. 2 is an example of replacement beads in example 1;
FIG. 3 is a diagram of the four-component average molecular model and the coarse grained model of the petroleum sample in example 1;
FIG. 4 shows the distribution results of the microstructure of the petroleum sample in example 1 at room temperature (1, asphaltene molecules; 2, colloidal molecules; 3, aromatic molecules; 4, saturated molecules);
FIG. 5 is a graph of the sample average molecular model and the coarse grained model in example 2;
FIG. 6 shows the distribution results of the microstructure of the asphaltene solution in example 2 at room temperature (1, asphaltene molecules; 5, toluene molecules; 1, asphaltene molecules; 2, colloid molecules; 3, aromatic molecules; 6, water molecules);
FIG. 7 is a graph of asphaltene aggregation rate as a function of concentration for example 2;
FIG. 8 is a graph of the sample average molecular model and the coarse grained model in example 3;
FIG. 9 shows the distribution of the microstructure of the mixed system of the petroleum sample and water in example 3 at room temperature.
Detailed Description
The technical solutions of the present invention will be described in detail below in order to clearly understand the technical features, objects, and advantages of the present invention, but the present invention is not limited to the practical scope of the present invention.
In the following examples, sulfur and nitrogen contents were measured in elemental analysis data of each petroleum sample by using an American ANTEK7000 sulfur-nitrogen analyzer; the sulfur content was measured by ultraviolet fluorescence (ASTM 5453); the nitrogen content was determined by chemiluminescence (ASTM 5762); the carbon and hydrogen content is measured by adopting a Flash EA 1112 organic microelement analyzer; gel Permeation Chromatography (GPC) was used for the relative molecular mass distribution. The experimental apparatus was a Waters GPC515-2410System, USA, the mobile phase was Tetrahydrofuran (THF), flow rate: 1mL/min, detector temperature: the standard sample was Polystyrene (PS) at 30 ℃. The four components are separated by liquid-solid adsorption chromatography (LSAC). Density measurement was carried out by the pycnometer method (GB 2540-81).
Example 1
The embodiment provides a method for constructing a mesoscale simulation petroleum molecule coarse graining model, and fig. 1 is a flow chart of the method, and the method comprises the following steps:
acquiring average molecular structures of petroleum samples with different component types; the method specifically comprises the following steps:
the component types of the petroleum sample of this example include four types: asphaltenes, gums, aromatic and saturated fractions. The specific method for obtaining the average molecular structure of the petroleum sample is as follows:
determining the element content in the petroleum sample; determining structural parameters of the petroleum sample; the average molecular structure of different types of petroleum samples was constructed.
Wherein: the elements include carbon, hydrogen, oxygen, nitrogen and sulfur. The structural parameters comprise aromatic carbon rate, unit aromatic ring number, unit aromatic carbon number, unit naphthenic ring number, nitrogen atom number, oxygen atom number, sulfur atom number, unit naphthenic carbon number and unit branched chain carbon number. The mass fraction of the element includes a carbon mass fraction, a hydrogen mass fraction, a sulfur mass fraction, a nitrogen mass fraction, and an oxygen mass fraction.
Specific procedures for obtaining the average molecular structure of different types of petroleum samples are described in patent application publication No. CN107704720A (a method for constructing a model of the average molecular structure of petroleum and a method for predicting properties), which is incorporated herein by reference in its entirety.
And step two, setting the replacement beads by adopting a molecular management method (the replacement beads defined in the embodiment are shown in FIG. 2). Judging whether the average molecular structure of the petroleum sample contains a six-membered aromatic ring structure, if so, positioning the six-membered aromatic ring structure in the center of the molecular condensed aromatic structure, and replacing with an equivalent replacement bead A6 bead; if not, judging whether a heteroatom structure is contained;
judging whether the average molecular structure of the petroleum sample after replacement contains a binary aromatic ring structure and a quaternary aromatic ring structure, if so, dividing a four-membered ring structure fragment and a binary ring structure fragment at intervals around the A6 bead, replacing the quaternary aromatic ring structure and the binary aromatic ring structure by equivalent replacement beads A4 and A2, and then judging whether the average molecular structure contains a heteroatom structure; if not, directly judging whether the structure contains a heteroatom structure;
step four, carrying out the step two or the step three, judging whether the average molecular structure of the petroleum sample contains a heteroatom structure, if so, replacing the average molecular structure with equivalent replacement beads (S, O, NI6, NI5 and the like), and judging whether the average molecular structure contains a naphthenic ring structure; if not, directly judging whether the naphthenic ring structure is contained;
step five, carrying out the step four, judging whether the average molecular structure of the petroleum sample contains a naphthenic ring structure, if so, replacing the naphthenic ring structure by adopting equivalent replacement beads (N6, N5, N4, N3, N2 and the like), and then replacing the side chain structure by adopting equivalent replacement beads (R2, R3 and the like) with side chain structures; if not, directly adopting equivalent replacement beads (R2, R3 and the like) with side chain structures to replace the side chain structures, and respectively obtaining various medium-scale simulation petroleum molecule coarse graining models of different types of petroleum samples. The results are shown in fig. 3, and fig. 3 shows the mesoscale simulated petroleum molecule coarse graining model diagrams corresponding to the molecular structures of the four different types of components of the petroleum sample in the present example.
The present embodiment also provides a method for predicting a property of a microstructure of a petroleum molecule, comprising the steps of:
step 1, obtaining solubility parameters of replacement beads by molecular dynamics simulation, and obtaining conservative force parameters by associating Flory-Huggins theory and dissipative particle dynamics theory as shown in table 1; the resulting conservative force parameters are shown in table 2:
table 1 solubility parameter table
Figure BDA0001890139150000061
Figure BDA0001890139150000071
TABLE 2 conservative force parameters table
A2 A4 A6 N2 N3 N4 N5 N6 NI5 NI6 O R2 R3 S T W
A2 25.0
A4 25.0 25.0
A6 25.2 25.1 25.0
N2 27.6 29.1 31.0 25.0
N3 27.2 28.0 29.5 25.1 25.0
N4 27.5 28.2 29.2 25.2 25.0 25.0
N5 28.9 29.9 31.2 25.0 25.0 25.1 25.0
N6 34.1 35.6 37.5 25.8 26.5 26.8 25.9 25.0
NI5 25.9 25.7 25.4 32.1 31.2 32.2 34.9 42.7 25.0
NI6 27.5 27.0 26.6 36.3 34.5 36.2 39.6 49.1 25.3 25.0
O 28.8 28.0 27.3 40.0 37.7 37.9 41.5 51.7 25.7 25.1 25.0
R2 32.2 35.4 38.9 26.1 27.1 28.1 27.0 25.3 39.4 45.5 51.3 25.0
R3 31.8 33.3 36.4 25.7 26.3 26.9 26.1 25.0 38.0 42.7 47.9 25.1 25.0
S 25.1 25.3 25.7 27.3 26.5 26.4 27.6 31.9 26.9 29.0 30.3 32.5 30.8 25.0
T 40.3 40.9 41.8 39.9 39.4 39.2 25.8 42.5 44.5 48.3 49.8 44.6 42.9 39.6 39.0
W 78.8 99.4 106.6 114.5 99.1 151.8 168.1 202.4 78.7 74.4 79.5 135.7 150.1 112.4 169.6 25.0
Note: t is toluene solvent beads and W is water beads.
Step 2, combining the mesoscale simulation petroleum molecule coarse graining model of different components and the obtained conservative force parameter with a dissipative particle dynamics method, and the steps are as follows: firstly, the simulation system obtains energy optimization through a steepest descent method under the condition of constant temperature and constant volume, bead kinematic calculation is carried out by applying a dissipative particle dynamic theory after the low energy is achieved so as to achieve a final equilibrium state, namely the energy of the system is stabilized at a low value, and a final simulation result is output. Obtaining a stable equilibrium state molecular dispersion system structure through dissipative particle dynamic simulation; the microstructure properties (number of sub-aggregates, molecular aggregation rate and molecular aggregation morphology) of the petroleum sample were predicted by statistical analysis. The density and four component ratios used for the simulations are shown in table 3, for example. The distribution of the microstructure of the petroleum sample at room temperature in the micro-equilibrium state is obtained, as shown in FIG. 4.
TABLE 3 oil sample analytical data
Figure BDA0001890139150000072
Figure BDA0001890139150000081
As can be seen from table 3 and fig. 4: the simulation result shows that four-component molecules in the petroleum sample system are in disperse distribution, no large amount of molecular aggregation occurs, and the simulation result is consistent with the stability of the sample obtained by experimental observation.
Example 2
In this example, the crude granulation model obtained by using a toluene solution of asphaltene as a sample was prepared in the same manner as in example 1, as shown in FIG. 5. The resulting conservative force parameter table was prepared as in example 1 and is shown in Table 1 above.
And (3) carrying out dissipative particle dynamics simulation on the coarse graining model and corresponding conservative force parameters, and establishing 5 different concentration simulation systems according to 20g/L, 60g/L, 100g/L, 150g/L and 200 g/L. The distribution of the microstructure of the asphaltene-containing toluene solution at room temperature was obtained as shown in FIG. 6. The resulting asphaltene aggregation rate versus concentration curve is shown in FIG. 7.
As can be seen from fig. 5, 6 and 7: the number of asphaltene aggregates in toluene and the rate of asphaltene aggregation increased with increasing asphaltene concentration, consistent with experimental results. The aggregate morphology of asphaltenes is similar to that observed in the experiment.
Example 3
In this example, the mixed system of the petroleum sample and water in example 1 was used as a sample, and the coarse grained model obtained in example 1 was subjected to the same method as shown in FIG. 8. The resulting conservative force parameter table is shown in Table 1, following the procedure of example 1.
And (3) carrying out dissipative particle dynamics simulation on the coarse graining model and corresponding conservative force parameters, wherein the dynamic simulation is carried out according to the water-oil ratio of 0.5: 1. 1: 1. 5: 1. 10: 1. 20: 1 and 50: 1 6 different simulation systems were set up. The micro-equilibrium structure distribution of the mixed system of the petroleum sample and the water at normal temperature is obtained, as shown in fig. 9.
As can be seen from fig. 8 and 9: with the increase of the water content, the system configuration is changed from a water-in-oil configuration to an oil-in-water configuration, asphaltene molecules are distributed between oil-water interfaces, good surface activity is shown, and the simulation result is consistent with the experimental result.
As can be seen from the above examples, the solution of the present invention can predict the configuration of the microequilibrium state of the petroleum sample and the related solution system. Only a small amount of experimental data are needed in the whole process, structural data which cannot be obtained in the experimental process can be obtained, and the method is small in workload, low in cost and fast.
The above-mentioned embodiments are only for illustrating the technical ideas and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and to carry out the same, and the invention is not limited to the embodiments, i.e. equivalent variations or modifications made within the spirit of the present invention are still within the scope of the present invention.

Claims (7)

1. A method of predicting microstructural properties of a petroleum molecule, characterized by: the method comprises the following steps:
step 1, obtaining solubility parameters of replacement beads by adopting molecular dynamics simulation, and associating Flory-Huggins theory and dissipative particle dynamics theory to obtain conservative force parameters;
step 2, combining a mesoscale simulated petroleum molecule coarse graining model and the obtained conservative force parameters with a dissipative particle dynamics method, and obtaining a stable equilibrium state molecular dispersion system structure through dissipative particle dynamics simulation; predicting the microstructure property of the petroleum sample through statistical analysis;
the method for constructing the mesoscale simulated petroleum molecule coarse graining model comprises the following steps:
acquiring average molecular structures of petroleum samples with different component types;
judging whether the average molecular structure of the petroleum sample contains a six-membered aromatic ring structure, if so, positioning the six-membered aromatic ring structure in the center of the condensed aromatic structure, and replacing the six-membered aromatic ring structure with an equivalent replacement bead; if not, judging whether a heteroatom structure is contained;
judging whether the average molecular structure of the petroleum sample contains a binary aromatic ring structure and a quaternary aromatic ring structure, if so, dividing a four-membered ring structure fragment and a binary ring structure fragment at intervals around a replacement bead of a six-membered aromatic ring structure, replacing the quaternary ring structure and the binary ring structure with an equivalent replacement bead, and then judging whether the average molecular structure of the petroleum sample contains a heteroatom structure; if not, directly judging whether the structure contains a heteroatom structure;
step four, carrying out the step two or the step three, judging whether the average molecular structure of the petroleum sample contains a heteroatom structure, if so, replacing the heteroatom structure by adopting equivalent replacement beads, and judging whether the average molecular structure contains a naphthenic ring structure; if not, directly judging whether the naphthenic ring structure is contained;
step five, carrying out step four, judging whether the average molecular structure of the petroleum sample contains a naphthenic ring structure, if so, replacing the naphthenic ring structure by adopting equivalent replacement beads, and then replacing the side chain structure by adopting equivalent replacement beads with side chain structures; if not, directly replacing the side chain structure with the equivalent replacement bead of the side chain structure to obtain a mesoscale simulation petroleum molecule coarse graining model;
the replacement beads comprise aromatic ring beads, naphthenic ring beads, heteroatom beads and side chain beads;
the aromatic ring beads comprise six-membered aromatic ring beads, four-membered aromatic ring beads and binary aromatic ring beads;
the cycloalkane ring beads include six-membered cycloalkane ring beads, five-membered cycloalkane ring beads, four-membered cycloalkane ring beads, three-membered cycloalkane ring beads, and two-membered cycloalkane ring beads;
the heteroatom beads comprise sulfur-containing heteroatom beads, oxygen-containing heteroatom beads, nitrogen-containing heteroatom six-membered ring beads and nitrogen-containing heteroatom five-membered ring beads;
the side chain bead comprises a bead with two continuous carbon atoms on the side chain and a bead with three continuous carbon atoms on the side chain.
2. The method of claim 1, wherein the average molecular structure of the different types of petroleum samples is obtained by:
determining the element content in the petroleum sample; determining structural parameters of the petroleum sample; and constructing the average molecular structure of the petroleum sample.
3. The method of claim 2, wherein: the component types of the petroleum sample include one or more of asphaltenes, gums, aromatic components, and saturates.
4. The method of claim 2, wherein: the elements include one or more of carbon, hydrogen, oxygen, nitrogen and sulfur.
5. The method of claim 2, wherein: the structural parameters comprise one or more of aromatic carbon rate, unit aromatic ring number, unit aromatic carbon number, unit naphthenic ring number, nitrogen atom number, oxygen atom number, sulfur atom number, unit naphthenic carbon number and unit branched chain carbon number.
6. The method of claim 2, wherein: the element content includes one or more of a carbonaceous mass fraction, a hydrogen mass fraction, a sulfur mass fraction, a nitrogen mass fraction, and an oxygen mass fraction.
7. The method of claim 1, wherein: the microstructural properties include one or more of a number of sub-aggregates, a rate of molecular aggregates, and a morphology of molecular aggregates.
CN201811467326.XA 2018-12-03 2018-12-03 Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof Active CN109378041B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811467326.XA CN109378041B (en) 2018-12-03 2018-12-03 Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811467326.XA CN109378041B (en) 2018-12-03 2018-12-03 Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof

Publications (2)

Publication Number Publication Date
CN109378041A CN109378041A (en) 2019-02-22
CN109378041B true CN109378041B (en) 2021-01-08

Family

ID=65375268

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811467326.XA Active CN109378041B (en) 2018-12-03 2018-12-03 Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof

Country Status (1)

Country Link
CN (1) CN109378041B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111899807B (en) * 2020-06-12 2024-05-28 中国石油天然气股份有限公司 Molecular structure generation method, system, equipment and storage medium
CN111892938B (en) * 2020-06-12 2022-08-30 中国石油天然气股份有限公司 Method, system, device and storage medium for determining molecular composition of crude oil
CN111739589B (en) * 2020-06-29 2022-07-05 青岛科技大学 Mesoscale-based simulation method for dissolving lignin by eutectic solvent
CN112216349A (en) * 2020-10-13 2021-01-12 中国民航大学 Method for constructing rigid asphaltene molecular structure in dissipative particle dynamics simulation
CN113223624B (en) * 2021-02-05 2023-04-28 中南大学 Cross-scale simulation method for predicting microstructure evolution in colloid shearing motion process
CN113223626B (en) * 2021-03-22 2023-06-09 中国石油大学(北京) Method and device for determining molecular scale reactor model

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104285227A (en) * 2012-05-30 2015-01-14 埃克森美孚研究工程公司 Method for reconciling model of composition of material from petroleum stream
CN106934137A (en) * 2017-03-05 2017-07-07 北京工业大学 A kind of coarse Molecular Dynamics method for analyzing Graphene assembly
CN108509724A (en) * 2018-04-03 2018-09-07 嘉兴学院 A kind of method of multi-scale Simulation nano particle heterogeneous fluid characteristic

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104285227A (en) * 2012-05-30 2015-01-14 埃克森美孚研究工程公司 Method for reconciling model of composition of material from petroleum stream
CN106934137A (en) * 2017-03-05 2017-07-07 北京工业大学 A kind of coarse Molecular Dynamics method for analyzing Graphene assembly
CN108509724A (en) * 2018-04-03 2018-09-07 嘉兴学院 A kind of method of multi-scale Simulation nano particle heterogeneous fluid characteristic

Also Published As

Publication number Publication date
CN109378041A (en) 2019-02-22

Similar Documents

Publication Publication Date Title
CN109378041B (en) Measure simulation petroleum molecule coarse graining model and construction method and prediction method thereof
Fang et al. Study on the asphaltene precipitation in CO2 flooding: a perspective from molecular dynamics simulation
Santos et al. Molecular dynamics simulation of n-alkanes and CO2 confined by calcite nanopores
Mullins et al. Asphaltenes explained for the Nonchemist
Bousige et al. Realistic molecular model of kerogen’s nanostructure
Chong et al. Excavating hidden adsorption sites in metal-organic frameworks using rational defect engineering
Wang et al. Influence of pressure on spontaneous imbibition in tight sandstone reservoirs
Zhang et al. A new diminishing interface method for determining the minimum miscibility pressures of light oil–CO2 systems in bulk phase and nanopores
Abutaqiya et al. Accurate modeling of asphaltene onset pressure in crude oils under gas injection using Peng–Robinson equation of state
Zhu et al. Correlating asphaltene dimerization with its molecular structure by potential of mean force calculation and data mining
Forsythe et al. A geological model for the origin of fluid compositional gradients in a large Saudi Arabian oilfield: An investigation by two-dimensional gas chromatography (GC× GC) and asphaltene chemistry
Costa et al. Aggregation of model asphaltenes: a molecular dynamics study
Liu et al. Thermodynamic modeling of the equilibrium partitioning of hydrocarbons in nanoporous kerogen particles
Zhang et al. Molecular insight into the oil charging mechanism in tight reservoirs
Okello et al. Towards estimation and mechanism of Co2 adsorption on zeolite adsorbents using molecular simulations and machine learning
Eftekhari et al. A molecular dynamics approach to investigate effect of pressure on asphaltene self-aggregation
Hu et al. Impact of maturity on kerogen pore wettability: A modeling study
Li et al. Conformational properties of a polymer tethered to an interacting flat surface
Desgranges et al. Coarse-Grained Model and Boiling Point Prediction for Asphaltene Model Compounds via HMC-WL Simulations
Xiong et al. Molecular dynamics simulation and experimental studies of the wettability behaviors of shales
Oulebsir et al. Diffusion of supercritical fluids through single-layer nanoporous solids: Theory and molecular simulations
Steinmetz et al. Experimental and mesoscopic modeling study of water/crude oil interfacial tension
Jensen et al. The impact of electrostatics in bulk Linde Type A zeolites
Sun et al. A permeation model of shale gas in cylindrical-like kerogen pores at geological conditions
Jin et al. Hydrocarbon phase behavior under nanoconfinement and Shale-PVT model development

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