Synthesis and Molecular Dynamics Simulation of Amphiphilic Low Molecular Weight Polymer Viscosity Reducer for Heavy Oil Cold Recovery

In order to reduce the viscosity of heavy oil, the performance of emulsifying viscosity reducers is limited. In this study, a new kind of amphiphilic low molecular weight viscosity reducer was prepared by emulsion copolymerization of acrylamide (AM), acrylic acid (AA), 2-acrylamido2-methylpropanesulfonic acid (AMPS), and Butene benzene (PB). The synthesis feasibility and viscosity reduction mechanism of viscosity reducer in heavy oil were explored using Materials Studio software from the perspective of molecular dynamics. The results of the molecular dynamics simulation revealed that the addition of viscosity reducer into heavy oil varied the potential energy, non-potential energy, density and hydrogen bond distribution of heavy oil. Benefiting from its structure, the benzene ring in PB was well embedded in the interlayer structure of asphaltene, contributing to weaken the network structure of the heavy oil. Moreover, the two strong polar groups (COO− and SO3) of AA and AMPS, which constituted the branched chains of the viscosity reducer’s molecular structure, gradually disassembled the network structure from the ‘inward’ to the ‘outward’ of the heavy oil network structure, thereby driving heavy oil viscosity reduction (as clarified by molecular dynamics). Owing to its temperature resistance, this kind of new amphiphilic low molecular copolymer could be an effective viscosity reducer for heavy oil cold recovery at elevated temperatures.


Introduction
China is one of the world's major heavy oil producing countries with abundant geological reserves. More than 70 heavy oil fields have been discovered in the Junggar Basin, Songliao Basin, Nanxiang Basin, and Erlian Basin, and another 12 large and medium-sized oil basins and depressions have been discovered with proven oil-in-place of about 4 billion tons [1]. With depleting conventional light oil reserves, unconventional resources are gaining relevancy for balancing the growing demands of energy. Heavy oil resources, as one of the most relevant unconventional resources, are trapped within underground rocks [2]. Nevertheless, the heavy oil is rich of asphaltene, resin, and paraffin. They aggregate together into macromolecules, which results in high viscosity and poor fluidity as a consequence of compromising the flow and production process [3][4][5]. The heavy oil viscosity typically ranges from 50 to 50,000 mPa.s at reservoir temperature and pressure [6][7][8][9][10][11]. Therefore, reducing the heavy oil viscosity is a prerequisite for improved ward' of the heavy oil network structure as illustrated in Figure 2, driving an viscosity reduction effect. It was worth noting that this molecular structure was from the author's previous study who found that the amphiphilic copolymer wh lecular structure was constituted of AM, and n-butylstyrene could favorably adso emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Material 5.0 molecular simulation software was used to quantitatively monitor the feasi synthesis of viscosity reducer from the perspective of molecular dynamics and to the interaction between heavy oil and viscosity reducer, and to explore the visc duction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic lecular weight copolymer was synthesized as a viscosity reducer in heavy oil.  Carboxyl group/amide group acid group; Water molecule; Heavy oil.

Materials and Instruments
Acrylamide (AM, 99.9%), acrylic acid (AA, 99.5%), 2-acrylamido-2-methylp sulfonic acid (AMPS, 99%) were obtained from Tianjin Beilian Fine Chemical D ment Co., Ltd, Tianjin, China. Butene benzene (PB, 98.9%) was from Beijing Ba Technology Co., Ltd, Beijing, China. Ammonium persulfate [(NH4)2SO4, 99.5%], sulfite (Na2SO3, 99.5%), sodium hydroxide (NaOH, 96%), tween-80, petroleum eth ing point 30-69 °C), dichoromethane (CH2Cl2, 99.5%), chloroform (CHCl3, 99.5%), (99.5%), anhydrous ethanol (99.7%), mercaptan (99.5%), and isopropanol (99.5 procured from Tianjin Fuchen chemical reagent factory, Tianjin, China. Kerosene alumina (100-200 mesh, specific surface area >150 m 2 /g, pore volume:0.23 cm cm 2 /g), glass indicator adsorption column, muffle furnace were purchased from ward' of the heavy oil network structure as illustrated in Figure 2, driving an viscosity reduction effect. It was worth noting that this molecular structure was from the author's previous study who found that the amphiphilic copolymer wh lecular structure was constituted of AM, and n-butylstyrene could favorably adso emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Material 5.0 molecular simulation software was used to quantitatively monitor the feasi synthesis of viscosity reducer from the perspective of molecular dynamics and to the interaction between heavy oil and viscosity reducer, and to explore the visc duction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic lecular weight copolymer was synthesized as a viscosity reducer in heavy oil.  Carboxyl group/amide group acid group; Water molecule; Heavy oil.   O/W structure, but also to disassemble the network structure from the 'inward' to 'outward' of the heavy oil network structure as illustrated in Figure 2, driving an effective viscosity reduction effect. It was worth noting that this molecular structure was inspired from the author's previous study who found that the amphiphilic copolymer whose molecular structure was constituted of AM, and n-butylstyrene could favorably adsorbed the emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Materials Studio 5.0 molecular simulation software was used to quantitatively monitor the feasibility of synthesis of viscosity reducer from the perspective of molecular dynamics and to simulate the interaction between heavy oil and viscosity reducer, and to explore the viscosity reduction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic low molecular weight copolymer was synthesized as a viscosity reducer in heavy oil.  Network structure of heavy oil and mechanism of viscosity reduction of amphiphilic polymer. Hydrophobic group containing benzene ring; Carboxyl group/amide group/sulfonic acid group; Water molecule; Heavy oil. O/W structure, but also to disassemble the network structure from the 'inward' to 'outward' of the heavy oil network structure as illustrated in Figure 2, driving an effective viscosity reduction effect. It was worth noting that this molecular structure was inspired from the author's previous study who found that the amphiphilic copolymer whose molecular structure was constituted of AM, and n-butylstyrene could favorably adsorbed the emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Materials Studio 5.0 molecular simulation software was used to quantitatively monitor the feasibility of synthesis of viscosity reducer from the perspective of molecular dynamics and to simulate the interaction between heavy oil and viscosity reducer, and to explore the viscosity reduction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic low molecular weight copolymer was synthesized as a viscosity reducer in heavy oil.   O/W structure, but also to disassemble the network structure from the 'inward' to 'outward' of the heavy oil network structure as illustrated in Figure 2, driving an effective viscosity reduction effect. It was worth noting that this molecular structure was inspired from the author's previous study who found that the amphiphilic copolymer whose molecular structure was constituted of AM, and n-butylstyrene could favorably adsorbed the emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Materials Studio 5.0 molecular simulation software was used to quantitatively monitor the feasibility of synthesis of viscosity reducer from the perspective of molecular dynamics and to simulate the interaction between heavy oil and viscosity reducer, and to explore the viscosity reduction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic low molecular weight copolymer was synthesized as a viscosity reducer in heavy oil.   O/W structure, but also to disassemble the network structure from the 'inward' to 'outward' of the heavy oil network structure as illustrated in Figure 2, driving an effective viscosity reduction effect. It was worth noting that this molecular structure was inspired from the author's previous study who found that the amphiphilic copolymer whose molecular structure was constituted of AM, and n-butylstyrene could favorably adsorbed the emulsified oil droplets in heavy oil wastewater [27][28][29]. In addition, the Materials Studio 5.0 molecular simulation software was used to quantitatively monitor the feasibility of synthesis of viscosity reducer from the perspective of molecular dynamics and to simulate the interaction between heavy oil and viscosity reducer, and to explore the viscosity reduction mechanism in heavy oil. Based on molecular dynamics, an amphiphilic low molecular weight copolymer was synthesized as a viscosity reducer in heavy oil.
Secondly, Materials Studio 5.0 software was used to obtain a total of 16 representative molecular structures of 4 distinctive heavy oil samples in the 3D Atomistic document [31], and the molecular number proportion to scale modeling was adjusted in Amorphous cell. The temperature, number of cycles, and initial density were fixed at 298 K, 10, and 0.6 g/cm 3 , respectively. The structural modeling of heavy oil was optimized from three distinctive aspects including, energy, geometry, and dynamic as illustrated in Figure 3.

Establishment of Heavy Oil Structure Model
Firstly, the composition of SZ36-1 was examined according to SH/T0 termination of petroleum asphaltene components" and SY/T 6282-1997" [30 ical approach of fluid physical property for heavy-oil reservoirs-measurem flow rheological characterization for crude oil" [24,25]. In detail, the analy gents were re-distilled in a round-bottomed flack equipped with a fraction refluxing condenser and collector, while heated around the boiling point. the neutral alumina was activated in a porcelain evaporating dishes and h for 6 h. Afterwards, 50 mg of dehydrated heavy oil was dissolved with 3 tilled petroleum ether, then a precipitate of asphaltene was removed by dried. The asphaltene was adsorbed with activated alumina and placed int tion column. The saturated hydrocarbons, aromatic hydrocarbons and were washed with 40~60 mL re-distilled petroleum ether, mixture of dichlo petroleum ether (v/v, 1/2), and mixture of chloroform and anhydrous eth respectively. The solvents were removed using rotary vacuum dryer, an were dried in a vacuum oven under constant temperature.
Secondly, Materials Studio 5.0 software was used to obtain a total of 16 molecular structures of 4 distinctive heavy oil samples in the 3D Atomistic and the molecular number proportion to scale modeling was adjusted in A The temperature, number of cycles, and initial density were fixed at 298 g/cm 3 , respectively. The structural modeling of heavy oil was optimized f tinctive aspects including, energy, geometry, and dynamic as illustrated in

Kinetic Simulation of Viscosity Reducer Synthesis Process
In this study, the molecular models of AM, AA, AMPS, and PB were 3D Atomistic document. The structure and performance of the obtained co optimized from three distinctive aspects including, energies, geometries, an the Forcite calculation module to ensure the consistency in molecular stru formance of copolymer with the desired one as depicted in Figure 4. After

Kinetic Simulation of Viscosity Reducer Synthesis Process
In this study, the molecular models of AM, AA, AMPS, and PB were established in 3D Atomistic document. The structure and performance of the obtained copolymers were optimized from three distinctive aspects including, energies, geometries, and dynamics in the Forcite calculation module to ensure the consistency in molecular structure and performance of copolymer with the desired one as depicted in obtained. Subsequently, the Amorphous cell module (temperature of 333 K and number of cycles of 10) was employed to establish the viscosity reducer model and the calculation was carried out using Forcite analysis. The dynamic simulation of viscosity reduction mechanism of copolymer in heavy oil was therefore established. The simulation was performed in the Compass force field using the ratio of heavy oil to viscosity reducer of 1:1, the temperature was controlled at 333 K, with cycle number of 10, and the optimized spatial structure of the interaction between viscosity reducer and heavy oil was obtained as represented in Figure 5.
optimizing the copolymers molecular structure in the polymer builder too tio of AM:AA:AMPS:PB = 30:2:3:3, chain length of 10, and chain segment obtained. Subsequently, the Amorphous cell module (temperature of 333 of cycles of 10) was employed to establish the viscosity reducer model and was carried out using Forcite analysis. The dynamic simulation of visco mechanism of copolymer in heavy oil was therefore established. The simu formed in the Compass force field using the ratio of heavy oil to viscosity the temperature was controlled at 333 K, with cycle number of 10, and the tial structure of the interaction between viscosity reducer and heavy oil w represented in Figure 5.   optimizing the copolymers molecular structure in the polymer builder too tio of AM:AA:AMPS:PB = 30:2:3:3, chain length of 10, and chain segment obtained. Subsequently, the Amorphous cell module (temperature of 333 of cycles of 10) was employed to establish the viscosity reducer model and was carried out using Forcite analysis. The dynamic simulation of visco mechanism of copolymer in heavy oil was therefore established. The simu formed in the Compass force field using the ratio of heavy oil to viscosity the temperature was controlled at 333 K, with cycle number of 10, and the tial structure of the interaction between viscosity reducer and heavy oil w represented in Figure 5.   complete dispersion of monomers with the pH of the mixture maintained between 9-10. Once the components were uniformly dispersed, 0.01-0.05% (based on total monomers mass) of chain transfer agents (mercaptan and isopropyl alcohol) were added, followed by nitrogen bubbling for 30 min. Subsequently, 1-1.5% of KPS (based on total monomers mass) was added to the mixture and the reaction was conducted at 60 • C under continuous nitrogen atmosphere purging for 8 h. Once the reaction terminated, the product was washed 3-5 times with a mixture of acetone/ethanol (2/1, v/v), dried at 80 • C, and crushed into powder. The optimal mole ratio of AM:AA:AMPS:PB = 60:5:6:9 was used to prepare the viscosity reducer (P 2 ).

Characterization and Viscosity Reduction Effect
FT-IR spectra was taken on Nicolet 560 FT-IR analyzer in the frequency range of 4000-400 cm −1 . The thermogravimetry analysis (TGA) was performed Dupont SDTQ600 thermogravimetry analyzer with heating rate of 20 K/min from room temperature to 600 • C under N 2 atmosphere. The molecular weight was measured on Agilent gel chromatography analysis system at 30 • C, the aqueous solution 0.1 mol/L was used as the mobile phase at a flow rate of 1 mL/min and the injection volume was 50 µL.
A certain amount of P 2 was added into the heavy oil, stirred for 20 min at room temperature, and the apparent viscosity of mixture was examined by Brookfield DV3T rheometer (Brookfield Engineering Lab, Bolfe, MA, USA) at a desired temperature.

Establishment of Heavy Oil Structure Model
Heavy oil is composed primarily of asphaltenes, and colloidal macromolecule compounds, but the specific composition is still unclear. At present, two mainstream approaches are used to establish the microscopic model of heavy oil. The first one refers to the construction of model based on the average molecular structure [31]. In this way, the B-L method is employed to draw the average molecular structure of the four components in heavy oil regarding the data from hydrogen spectrum and X-ray diffraction, and the complex molecular structure can be transformed into four distinctive structures. However, this does not picture the specific molecular structure of heavy oil. Hence, the complex interaction in-between heavy oil molecules cannot be simulated. The second one lies on the construction of physical properties model which is identical to that of heavy oil by exploiting the known four representative molecules [32]. After analyzing the composition of heavy oil, four main components including, saturated hydrocarbons, aromatic hydrocarbons, resins and asphaltene constituted the heavy oil as listed in Table 1. Referring to [33], a total of 16 molecular structures representing the 4 components of heavy oil were obtained. Based on the actual composition of heavy oil, the number of molecules of each component was adjusted so that the composition of the constructed heavy oil model was consistent with that of SZ36-1 as shown in Table 2. The heavy oil model was therefore established using the method described in Section 2.2.1.

Dynamic Simulation of Viscosity Reducer Synthesis Process
As earlier described in Section 2.2.2, both temperature variation and total energy variation were simulated during the synthesis process of viscosity reducer, and a total of 50 dynamic models were acquired as indicated in Figures 6 and 7. From the curve of temperature variation represented in Figure 6, it can be observed that the temperature fluctuated in small range close to the set synthesis temperature (333 K), which revealed that the set temperature of the periodic system had little influence on the simulated results. The curve of variation of the total energy was also presented in Figure 7. It was found that, at a fixed reaction temperature of 333 K, the total energy of the system also slightly fluctuated within a small range. We could not only observe that this fluctuation occurred around 670 kcal/mol, but also the energy of C-C double bond was 145.7 kcal/mol, inferring that under this reaction condition, the rupture of C-C double bond could be envisioned, and the polymerization reaction was therefore theoretically feasible.

Dynamic Simulation of Viscosity Reduction Mechanism
As it can be observed in Figure 8, compared to the corresponding curve of heavy oil, that of the viscosity reducer with heavy oil significantly increased. This sharp variation in potential energy could be due to the attachment of asphaltene and resin molecules to the viscosity reducer which resulted in the forming of a large spatial steric hindrance which was constituted of an abundant number of branched-chain structures of the viscosity reducer molecules. It can be inferred that the viscosity reducer molecules which was constituted of a sulfonic group, carboxyl group, and hydrophobic structure with benzene ring mutually induced strong adsorption of benzene ring structure in asphaltene and colloid molecules. From the same Figure, we could also observe an obvious increase in negativity of non-potential energy which might be attributed to the two strong polar monomers (AA and AMPS) as a consequence of intensifying the interaction between the viscosity reducer and inward of heavy oil. As a matter of fact, the strong hydrophilicity of the sulfonic group in AMPS was beneficial to maintain the O/W emulsion via the formation of more stable hydrogen bonding. Subsequently, the network structure was gradually disassembled from the "inward" to "outward" of the heavy oil structure. Under such circumstances, the

Dynamic Simulation of Viscosity Reduction Mechanism
As it can be observed in Figure 8, compared to the corresponding curve of heavy oil, that of the viscosity reducer with heavy oil significantly increased. This sharp variation in potential energy could be due to the attachment of asphaltene and resin molecules to the viscosity reducer which resulted in the forming of a large spatial steric hindrance which was constituted of an abundant number of branched-chain structures of the viscosity reducer molecules. It can be inferred that the viscosity reducer molecules which was constituted of a sulfonic group, carboxyl group, and hydrophobic structure with benzene ring mutually induced strong adsorption of benzene ring structure in asphaltene and colloid molecules. From the same Figure, we could also observe an obvious increase in negativity of nonpotential energy which might be attributed to the two strong polar monomers (AA and AMPS) as a consequence of intensifying the interaction between the viscosity reducer and inward of heavy oil. As a matter of fact, the strong hydrophilicity of the sulfonic group in AMPS was beneficial to maintain the O/W emulsion via the formation of more stable hydrogen bonding. Subsequently, the network structure was gradually disassembled from the "inward" to "outward" of the heavy oil structure. Under such circumstances, the long carbon chain in AMPS, asphaltene, and colloid molecules were longitudinally extended resulting in the formation of planar and easily combed linear structure, thus preventing the occurrence of three-dimensional assembly structure and eventually the viscosity reduction. From Figure 9, it was noticed that the overall density of both heavy oil and viscosity reducer with heavy oil fluctuated within a small range. However, compared to the overall density of heavy oil (1.00 g/cm 3 ), that of the viscosity reducer with heavy oil (0.96 g/cm 3 ) was slightly lower, translating that the asphaltene and colloid compounds in heavy oil were removed. Thus, this phenomenon contributed not only to lowering the density of heavy oil, but also the declining of heavy oil viscosity could be envisioned (which will be investigated later). Total energy (Kcal/ mol) Frame Potential energy of heavy oil system with viscosity reducer added Non-potential energy of heavy oil system with viscosity reducer added Potential energy of heavy oil system Non-potential energy of heavy oil system It is generally believed that asphaltenes molecules are shaped planar and have dense aromatic rings, side-chains, and a large number of heteroatoms which facilitate the susceptibility to undergo hydrogen bonding and self-aggregation. Moreover, the basic structure of resin molecules is based on the aromatic system composed of multiple rings and there are surrounded by naphthenic rings. Both aromatic rings and naphthenic rings are imbued with several normal or isoalkyl side chains of distinctive lengths, and their molecular weight, polarity, and susceptibility to undergo hydrogen bonding are all second to asphaltene [35]. The length of hydrogen bond distribution of heavy oil and in the presence of viscosity reducer were investigated and indicated in Figure 10. As it can be observed, the hydrogen bond of heavy oil was mainly distributed between 1.11-1.14 E (A), which accounted for about 43%, followed by 1.53-1.55 E (A), which represented about 25.7%. These findings suggested that there was not only a large number, but also a stable hy-drogen bonds in the heavy oil. From the corresponding curve of the length of hydrogen bond distribution of heavy oil with viscosity reducer shown in the same figure, we could determine that the hydrogen bonds were mainly distributed between 1.1-1.11 E (A), which accounted for almost 30%, followed by 1.53-1.57 E (A) with a rate of 6.9% compared to the length of hydrogen bond distribution of heavy oil, which in the presence of viscosity reducer exhibited multiple distribution peaks. This could be ascribed to the amphiphilic function of the viscosity reducer, which contributed to break down the "planar" structure in heavy oil as a result to weaken the interaction between asphaltene and resin. It is worth noting that, compared to the hydrogen bond distribution in the heavy oil, in the presence of viscosity reducer was slightly wider, mainly attributed to a significant decrease in the proportion, indicating that the viscosity reducer can significantly reduce the distribution proportion of hydrogen bonds in the heavy oil. After collapsing the rigid structure, in addition to the viscosity reducer, various components were involved in the formation of hydrogen bonding, contributing to effectively reduce the viscosity of heavy oil.   [36]. The broad absorption peak noticed at 622.54 cm −1 was due to the absorption characteristic of benzene ring. In the light of these results, the four monomers were involved in the copolymerization reaction. acted monomers. The initial decomposition temperature was recorded as 250 °C, while that of polyacrylamide was 210 °C [36], suggesting that the viscosity reducer had superior temperature resistance. This was mainly due to introduction of benzene ring structure of PB and sulfonic group which contributed to improve the thermal stability of viscosity reducer, thus indicating that the viscosity reducer could find use in high temperature heavy oil formation.

TGA Analysis of Viscosity Reducer
The TGA analysis curve of viscosity reducer and polyacrylamide taken from room temperature to 600 • C in N 2 atmosphere was displayed in Figure 12. The first region be-low100 • C with a large slope was ascribed to elimination of water molecules and unreacted monomers. The initial decomposition temperature was recorded as 250 • C, while that of polyacrylamide was 210 • C [36], suggesting that the viscosity reducer had superior temperature resistance. This was mainly due to introduction of benzene ring structure of PB and sulfonic group which contributed to improve the thermal stability of viscosity reducer, thus indicating that the viscosity reducer could find use in high temperature heavy oil formation. polymerization reaction.

TGA Analysis of Viscosity Reducer
The TGA analysis curve of viscosity reducer and polyacrylamide taken from room temperature to 600 °C in N2 atmosphere was displayed in Figure 12. The first region be low100 °C with a large slope was ascribed to elimination of water molecules and unre acted monomers. The initial decomposition temperature was recorded as 250 °C, while that of polyacrylamide was 210 °C [36], suggesting that the viscosity reducer had superior temperature resistance. This was mainly due to introduction of benzene ring structure o PB and sulfonic group which contributed to improve the thermal stability of viscosity re ducer, thus indicating that the viscosity reducer could find use in high temperature heavy oil formation.

Molecular Weight Test of Viscosity Reducer
As determined by gel chromatography method, the results of molecular weight of viscosity reducer were listed in Table 3. It can be observed that, the weight-average molecular weight (Mw) of viscosity reducer varied between 8-100,000. The low molecular weight enables the dispersion of viscosity reducer in heavy oil, which favor the facile migration into the 'inward' of heavy oil as a result to better weaken, and eventually destroy, the network structure of the heavy oil.

Effect of Viscosity Reducer on Heavy Oil Viscosity
In order to ensure the fluidity of heavy oil, the experimental temperature was maintained at 70 • C, distinctive concentration of P 2 solutions was added into the heavy oil and the final viscosity was taken. Figure 13 revealed the contribution of viscosity reducer concentration on the viscosity of heavy oil. As it can be observed, with increasing P 2 concentration from 0 to 0.5 wt.%, the viscosity of heavy oil dropped by almost 86.14%, as to reflect its favorable viscosity reduction effect. As a matter of fact, the viscosity reducer (P 2 ) had a low molecular weight and sufficiently long molecular chain. Moreover, the molecular structure of P2 was composed of strong polar groups COO − (AM, AA) and SO 3 − (AMPS), which aided the migration of P 2 into the 'inward' of heavy oil as a consequence to disassemble, but also to effectively oppose the water molecules from entering the internal structure. Thus, the O/W structure predominated over the W/O one as a consequence to favor the fluidity of heavy oil. In this way, the benzene ring structure which constituted the P 2 molecular structure was well embedded in the asphaltene interlayer [37]. From the same plot, we could notice that increasing the P 2 concentration contributed to gradually increase the apparent viscosity of heavy oil, which might be due to entanglement of the long chain of viscosity reducer in the O/W emulsion. the fluidity of heavy oil. In this way, the benzene ring structure which constituted the P2 molecular structure was well embedded in the asphaltene interlayer [37]. From the same plot, we could notice that increasing the P2 concentration contributed to gradually increase the apparent viscosity of heavy oil, which might be due to entanglement of the long chain of viscosity reducer in the O/W emulsion.

Conclusions
In this study, considering SZ36-1 as the research object, a combination of experimental method and molecular dynamic simulation method were employed to establish the structural model of heavy oil, and dynamic simulation of viscosity reduction mecha-

Conclusions
In this study, considering SZ36-1 as the research object, a combination of experimental method and molecular dynamic simulation method were employed to establish the structural model of heavy oil, and dynamic simulation of viscosity reduction mechanism. The feasibility of the synthesis method of low molecular weight amphiphilic polymer viscosity reducer was confirmed by molecular dynamics and the molecular structure was characterized. By introducing the hydrophobic benzene ring and strong hydrophilic groups into the polymer molecular chain, the amphiphilic low viscosity reducer could interact with the 'inward' resulting to imbue the 'outward' with various anisotropies, which weakened and destroyed the dense network structure of heavy oil and eventually drive the viscosity reduction of heavy oil. Thus, with increasing low molecular weight concentrations from 0 to 0.5%, the heavy oil viscosity reduced by about 86.14% at 70 • C. The experimental results also proved that low molecular weight amphiphilic structure viscosity reducer exhibited good viscosity reduction effect. For heavy oil of SZ36-1 (formation temperature at 50-65 • C), the polymer can meet the need for heavy Oil Cold Recovery, which has a good reference for the flowable heavy oil.
Author Contributions: Conceptualization, L.X. and W.R.; data curation, Y.C. and W.G.; methodology C.M. and M.Z.; investigation, H.Z. and X.L. All authors have read and agreed to the published version of the manuscript.