Modeling the Wax Deposition Process and Diffusion Kinetics of Wax Crystallization

Due to the depletion of reserves of medium and light oils, the oil, and gas industry has to pay more and more attention to the industrial development of heavy, highly viscous hardto-recover oil deposits. In the world, light oils are extracted no more than 50%, heavy highviscosity oils in the range from 10 to 30%, depending on the characteristics of the oil [1] and every year the share of production of heavy oils enriched with high-melting paraffinic hydrocarbons and Resinous-Asphaltene Substances (RAS) increases. In the production of high-viscosity oils, especially in cold climates, serious problems arise associated with the deposition of wax, which reduces the efficiency of field development. Oil is subdivided into low paraffinic (less than 1.5% by weight), paraffinic (from 1.5 to 6%), high paraffinic (more than 6%). This is an important characteristic influencing the technologies used in oil production. Paraffinic oils have the ability to form deposits on the walls of the well, as a result of which there is a decrease in debt and, in the future, blockage of the well. The practice of oil production in the fields shows that the main areas of accumulation of AsphalteneResin-Paraffin Deposits (ARPD) are well pumps, lifting strings in wells, flow lines from wells, reservoirs of field gathering points [2].

The main factor affecting the intensity of the paraffins crystallization is a decrease in temperature, which, in turn, is determined by the degassing process (heat absorption due to degassing and a decrease in temperature due to the Joule-Thomson effect). Thus, the main "driving effect" of reducing well productivity is the degassing process [5].
It is known that the dissolving capacity of oil in relation to paraffins decreases with decreasing temperature and degassing of oil. In this case, the temperature factor prevails [2]. The rate of heat transfer depends on the temperature difference between the fluid and the surrounding rocks at a certain depth, and the thermal conductivity of the annular space between the lifting pipes and the production column.
With an increase in the rate of oil movement, the intensity of deposits initially increases, which is explained by an increase in flow turbulence and, consequently, an increase in the frequency of formation and separation of bubbles from the pipe surface, floating suspended particles of paraffin and asphalt-resinous substances [2]. In addition, the fluid stream tears off part of the deposits from the pipe walls, which can explain the sharp decrease in deposits in the interval 0-50m from the wellhead. The fluid stream, having high flow velocities is more resistant to cooling, which also slows down the formation of ARPD.
The pipe wall roughness and the presence of solid impurities in the system also contribute to the release of paraffin from oil into the solid phase [2]. In addition to these main factors, the intensity of waxing of pipelines during the transportation of water cuttings of good production can be influenced by the water cut of the product [6] and the pH value of formation waters [2]. Studies [4,7] have shown that the most probable mechanism of paraffin formation is the crystallization mechanism, i.e., in conditions when the gas-liquid flow can transport heavy components, the formation, and further build-up of paraffin deposits occur due to the growth of crystals directly on the surface of underground equipment. Wax formation begins where oil comes in contact with the cold walls of underground equipment.
The necessary conditions for the formation of paraffin deposits are [8]: the presence of high-molecular hydrocarbon compounds in oil-primarily paraffins.

B.
reduction of reservoir pressure to bubble point pressure.
C. decrease in the flow temperature to values at which the solid phase is separated from oil. D. the presence of a substrate with a low temperature, on which high-molecular hydrocarbons crystallize with sufficiently strong adhesion to the surface, excluding the possibility of tearing off the deposits by the flow of a gas-liquid mixture or oil at a given technological mode.
There are many other factors that contribute to or prevent the intensive formation of paraffin deposits. The most significant of them can be attributed to the oil flow rate, the process of gas release when the pressure in the pipeline decreases, the presence of mechanical impurities, which are active crystallization centers, the condition of the equipment surface, the water cut of the well production, etc. Various experimental techniques such as flow profile, cold finger and rotating cylinder is used to study the wax deposition phenomenon [9].

Approaches to mathematical modeling and mechanisms of the wax deposition process
The mathematical modeling approach is widely used to predict and monitor wax deposition, either directly numerically or using software that has an internal mathematical model. Wax deposition models are designed based on assumptions and options.
Four wax deposition mechanisms need to be considered when modeling wax deposition growth in pipelines, namely molecular diffusion, shear dispersion, gravity deposition, and Brownian diffusion [10]. Molecular diffusion, which can be described as mass transfer and energy balance, is universally recognized as dominant. At the same time other mentioned mechanisms have been identified as having little impact on this phenomenon.
Some researchers have proposed a shear dispersion mechanism to analyze the wax deposition behavior, especially in laminar flow. However, other researchers reported based on several experiments that this mechanism plays an insignificant role in the deposition of wax [11]. The gravity settling mechanism assumes that wax crystals will settle to the bottom of pipelines because the crystals are denser than oil. However, this effect appears insignificant since experimental studies comparing vertical and horizontal liquid flows did not reveal differences in the amount of deposited wax [12]. Once the wax crystals are deposited and suspended in the oil, the wax crystals will behave according to the Brownian motion rules. Since the effect of motion is likely to transport crystals to areas of lower wax concentration, this can influence the wax deposition behavior. However, its effect is generally considered minimal and is usually neglected during modeling of the overall deposition mechanism. Consequently, the last three mechanisms are not widely accepted, especially for computer simulation procedures [13][14][15].
Various studies have shown that molecular diffusion is considered the dominant mechanism in the wax deposition process [16,17], and shear stress is significant only when the oil temperature is much lower than WAT and the wax concentration in the oil is high [18].
In turn, the process of wax deposition due to molecular diffusion can be divided into the following five stages [19]: precipitation of dissolved paraffin molecules and formation of a layer of deposits on the surface of the pipe wall. b.
formation of radial gradient concentration of dissolved paraffin components, which leads to the diffusion of paraffin molecules from the bulk of oil to the wall. c.
the deposition of paraffin components on the surface of the existing sediment, which helps to increase its thickness.
internal diffusion of paraffin molecules inside the sediment, which leads to an increase in the proportion of solid paraffin in the sediment.

e.
counter diffusion of dewaxed oil from the deposit. Calculates the mass transfer rate of wax using the heat transfer analogy. Takes into account the phenomenon of shear and includes the effect of shear deposition.

University of Michiqan Model
LedaFlow from Kongsberg Simulates wax crystallization and wax deposition on the pipe wall. It has two main limitations: it does not take into account the effect of pressure on WAT, which determines the onset of wax settling and does not simulate the restriction of flow due to the formation of wax deposits, thus failing to reproduce the inlet pressure increment. [13,19] FloWax FloWax from KBC Complete Composite Wax Deposition Model which takes into account thermodynamic wax deposition, wax diffusion based on the heat and mass transfer analogy and the shear effect. [26] Various commercial instruments have been developed predict correctly and describe the wax deposition phenomenon [13,[20][21][22][23][24][25][26]. All wax deposition models follow the same formula for estimating wax deposition thickness. What changes in each model is the calculation of the parameters, which depends on the experimental data. A brief description of the model used by the software is shown in table 1.

Diffusion kinetics of paraffin crystallization
Wax crystallization begins with nucleation, where in oil, solid asphaltene particles with a size of 1nm or a solid phase can be used as nuclei. In [27], an empirical model is used to describe the processes of crystallization and wax deposition, in contrast to which a diffusion model corresponding to the diffusion kinetics of the formation and deposition of waxes from oil is proposed below. The crystallization process is a diffusion process, and the condition is satisfied on the surface of the nuclei [28] ( ) Assuming that dr Vdt = the equation equation of kinetics crystallization of paraffins can be represented as is the crystallization rate constant, which also depends on temperature T, C-is the concentration of paraffin in the volume of the boundary layer, C * -is the equilibrium concentration of paraffin in the liquid and solid phases, depending on temperature, β L -is the mass transfer coefficient, D -is coefficient of molecular diffusion . The solution to this equation can be represented as

RDMS.000857. 15(2).2021
If the maximum thickness of the deposited paraffin layer is equal to the pipe radius, then there is a complete blockage of the pipe section. Using the experimental data [27], it is possible to estimate the coefficients entering the equation (4)  (where ∆ T is the absolute temperature difference between the surface temperature and the crystallization temperature). Practical calculations show that in the indicated limit of temperature change, the exponent changes insignificantly 0.65 n = .
The results of experimental studies [27] and calculated values by formula (4) for different temperatures are shown in Figure 1. As follows from Figure 1, the concentration of paraffins eventually tends to the established state.

Model for the formatıon of a dense layer of partıcles on a surface
The formation of a dense layer of paraffin particles on the pipe surface is determined by their transfer to the surface and the nature of the surface itself (roughness). The layer of particles is compacted under the influence of external stresses (pressure) and coagulation phenomena.
The mass of the deposited layer of paraffins on the surface of the cylindrical pipe is determined based on their geometric considerations in the form In the following description, the index s will apply to all parameters of the deposited layer of paraffin particles. R 0 -is the radius of the clean pipe (no deposits), s δ -is thickness of the deposited paraffin layer, s ρ -is layer density, and L-is pipe length.
Using the expression (5), the change in mass of the deposited layer of paraffins is defined as The solution of (7) under the initial conditions (for a clean pipe) is defined in the form Unlike paraffin particles that form in the volume of the boundary layer, asphalt-resinous substances are contained in oil in a free state. Experimental studies of the sedimentation of asphaltresinous substances [29] in the tubes of heat exchangers showed that the thickness of the deposits can be described by expression . Here b m and k are empirical constants determined from experimental data.
As shown in figure 2, with an increase in the flow rate, the limiting thickness of the deposits decreases, and with an increase in temperature, it increases.

Temperature profile in deposited pipe
The deposition of paraffin particles and asphalt-resinous substances with low coefficient of thermal conductivity contributes to reduce of the oil flow temperature over time with increasing thickness, which has a positive effect on the rate of wax crystallization. Since the formation and deposition process is slow, heat transfer in the deposited pipe can be described by the following equation: where is ( ) where is ( ) r ϕ the initial temperature distribution over the pipe section, depending on the initial temperature of the flow in the pipe. Using the method of separation of variables, solution (18) can be represented as Here ( ) 0 J r is the zero-order Bessel function of the real argument, , Thus, the series (19) quickly converges and it is enough to limit itself to the first term, i.e. the temperature distribution in the pipe with deposition can be represented as In [30], of flow temperature change with paraffin deposition are given. Since we are limited only to the first member of the series,  , c T is the oil flow temperature in the pipe. Figure 3 shows the temperature distribution in the pipe for different deposit thicknesses.

Results and Discussion
The presented diffusion models of crystallization kinetics (4) and deposition of paraffins on the surface of pipes (7,8) are migration-diffusion in nature, where both diffusion of paraffin from the liquid phase to the solid phase and diffusion transfer of paraffin particles to the surface plays an important role. The formation of a dense layer of particles on the surface of pipes is determined by many factors, among which it is important to note hydrodynamic, mass exchange and thermodynamic conditions, rheological properties of the dispersed medium, adhesion compatibility of particles with a streamlined surface, physical and chemical transformations in the boundary layer, particle sizes and concentration, the orientation of the streamlined surface, etc. Analysis of such processes showed that deposition on the inner surface of the pipes to form a dense layer of particles with a low thermal conductivity degrades heat exchange with the external medium, resulting in a decrease in heat transfer and heat transfer coefficients. Introducing a dimensionless thickness of deposits in the form, (in these expressions the subscript "0" refers to a clean pipe). As follows from these formulas, if, in the case of laminar flow, an increase in the thickness of the deposits in the pipes leads to hydrodynamic instability, then for a turbulent flow, the same factor contributes to an increase in energy dissipation, a decrease in the scale of turbulence and, associated with this, attenuation in the intensity of turbulence. From this, it follows that crystallization, precipitation, and formation of a dense layer on the surface of the pipes have a reverse effect on the flow hydrodynamics and heat -mass transfer in the pipes. If a decrease in the temperature in the core of the flow is observed over time, then the formed paraffin particles together with asphaltenes particles have a significant influence on the structure formation in the flow and thereby change the rheological properties of the oil [31]. Formation of a dense layer of particles of asphalt, and paraffin substances on the surface of pipes significantly affects the change in temperature profile in time, length and radius (13). A decrease in the temperature of the flow in the pipe results in a shift of the crystallization front to the center of the pipe, as a result of which the process of crystallization of paraffins is carried out in the volume of the pipe. The formed paraffin particles in the flow volume can also be transferred diffusely to the pipe surface, confirming the migration-diffusion deposition mechanism. At a high concentration of paraffin particles in the flow volume, it is possible to form coagulation structures that pass into aggregates, clusters of aggregates and eventually form together with particles of asphaltene a frame, leading to an increase in the rheological properties and viscosity of the oil. In the case of coagulation of paraffin particles in the flow volume, their free settling is possible, characterized by a migration-gravitational or mixed mechanism.
Computational modeling can be a powerful tool for predicting and interpreting wax deposition phenomena after validation with available experimental data. Most laboratory experiments are based on a single-phase flow, which often gives deviations from multiphase flow cases. Simulation results are not always applicable to real field cases. For a more reliable interpretation, the scaling effect must be taken into account. Current research needs to focus on more complex environments. One way to overcome this problem is to validate these models with an additional field dataset.