An improved modeling approach for asphaltene deposition in oil wells including particles size distribution

There are several approaches to model Asphaltene deposition process in the wellbore. There are different assumptions to simplify the problem in the previous investigations for specific conditions, limiting the prediction range of the models. In this work, the effect of precipitated asphaltene particles size is included, to extend the available modeling approaches for deposition profile. To do so, two-dimensional partial differential equations based on asphaltene micro aggregates material balance including asphaltene aggregation, diffusion and deposition are numerically discretized and solved to find asphaltene deposition profile, in radial and vertical directions of vertical oil wells. The modeling results are verified with the results of the well-known ADEPT (asphaltene deposition tool in flow lines) model of Kurup et al. (2011). The size dependent diffusion coefficients of Escobedo and Mansoori (2010) are used to extend the base model. In addition, the population balance method (PBM ) was included to improve the aggregation process description with size distribution of asphaltene particles. Based on the developed model a parametric study is performed to study the effect of asphaltene particles average size, flow rate, wellbore radius and fluid viscosity. The model evaluation shows the importance of asphaltene particle size in the deposition profile. In addition, the evaluation results show that as the average asphaltene particle size increases for a given distribution, the amount of deposition in the wellbore decreases.


Introduction
Asphaltene deposition is an important challenge in the petroleum industry, and proper understanding of its nature and deposition behavior is still under investigation.The general operational definition of asphaltene is according to its solubility in different solvents; Asphaltenes are polar and non-volatile components of the crude oil, which are insoluble in low molecular weight alkanes (e.g.n-heptane).These molecules are slightly soluble in the aromatic solvents (e.g.toluene).Asphaltene deposition in the oil industry processes like production, transportation and processing is one of the important flow assurance problem that can make different challenges.One of the best methods to deal with this problem is the modification of production processes and prevention of the precipitation and deposition by the economical, technical and environmental point of view.In fact, other methods such as chemical and mechanical methods are used when the production process modification is impossible.It should be mentioned that different chemicals at different temperature and pressure conditions have different effects for different kind of crude oils that it makes treatment difficult (Rezakazemi et al., 2017).Therefore, there are many researches about the type of chemicals that are compatible with various crude oils to decrease asphaltene precipitation and then deposition (Hutin et al., 2016;Kashefi et al., 2018).
The fluid flow in porous media mostly occurs in the horizontal directions with a limited change in thermodynamic conditions.In contrast, the petroleum fluid experiences considerable changes in pressure, temperature and composition making challenges to deposition modeling in the vertical wellbores.With upward oil flow in the wellbore, the thermodynamic conditions change leading to the asphaltene separation known as precipitation for which there are many thermodynamic modeling researches in the literature.Deposition is the process of asphaltene aggregates sticking (after precipitation) to the wellbore wall.A review of previous literature reveals that there are only a few qualitative and quantitative research studies on this phenomenon in both laboratory and field scale for the deposition process.The information about the asphaltene deposition position and thickness in the production wells as a function of time are limited.Therefore, model verification is difficult.High-pressure conditions and steep slopes in some deviated wells prevent the use of the tools for monitoring of asphaltene deposition.In some cases, the detector itself leads to pressure changes.Therefore, it promotes deposition process due to system instability leading to blockage and sometimes losing the tool.
One of the earliest studies on asphaltene deposition is provided by Ramirez-Jaramillo et al. (2006).The modeling assumptions were in a way that the asphaltene sticks on the walls due to the molecular diffusion in radial direction and the concentration gradient is due to temperature gradient along the wellbore.Based on the results from this model, the deposition thickness cannot exceed 15% of inner effective diameter of the wellbore.This model fails to match with the recent laboratory data, which shows that, the temperature gradient around the wall does not critically affect the deposition rate.Soltani Soulgani et al. (2009Soulgani et al. ( , 2010) ) utilized a simple empirical correlation, and assumed that the deposition rate on the wellbore was controlled by chemical interaction mechanism.Despite the acceptable results obtained from this model, there is no reasonable explanation for the mechanism of the deposition process.The next model was developed by Vargas et al. (2010).This model included some sub-models explaining precipitation, aggregation, deposition and transportation.In fact, according to this model, a deposition simulator was prepared based on a balance equation for turbulent flow associated with thermodynamic precipitation modeling using Perturbed Chain Statistical Associating Fluid Theory (PC-SAFT).Utilizing a simple mathematical model, it predicts the deposition of asphaltene along the flow line in different flow situations (Vargas et al., 2010).The major drawback of this approach was that it assumed a constant rate of precipitation and aggregation along the wellbore.Eskin et al. (2011) developed a model for predicting asphaltene deposition in the wellbore for turbulent flow based on particle mass flux transfer.They studied the deposition problem from fluid mechanic's point of view.In this method, the description of particle transfer into the wall was a function of particle size, density and flow parameters.They used the couette device, consisting of two cylinders, for the verification of their results.The inner cylinder rotates with various speed to apply different flow and particles deposit in the outer immobile wall to model the asphaltene deposition in the wellbore using different particle size distributions with method of particle population balance.This is a strong model from the physical point of view but it has some weaknesses in mathematical simulation.It should be mentioned that in the case of particle size Marques et al. (2008) obtained notable results from laboratory work.In addition, they showed that different sizes of asphaltene aggregates have different aromaticity, aliphatic composition, alkyl chains and metal concentration as well as percent of vanadium to nickel.These results may conclude to the difference in asphaltene deposition rates as well Kurup et al. (2011); Kurup et al. (2012) extended the model proposed by Vargas et al. (2010) through the 1D steady and pseudo steady precipitation simulator (ADEPT).In this model, the precipitation was obtained as a profile for the whole desired length of the well, which is an advantage over Vargas et al. (2010) model, however, aggregation rate was still constant.In addition, they have claimed that the 1D modeling along the wellbore is reasonable although, the radial diffusion has a noticeable effect.
This work describes a comprehensive asphaltene deposition simulator with macroscopic simulation approach.Despite its simplicity, it includes all processes that, the asphaltene particles may encounter in the wellbore.It considers amount, depth, and profile of asphaltene deposition along the wellbore by taking into account the changes in pressure, temperature and particle size.Like the earlier works of Vargas et al. (2010) and Kurup et al. (2011), three processes of precipitation, aggregation and deposition are included in this model.In addition, the proposed model in this work includes the role of aggregates' size in the asphaltene deposition.For this purpose, the amount of precipitation, asphaltene particles aggregation rate, and asphaltene deposition rate are investigated separately in the wellbore (these parameters are assumed constant in the previous works).Asphaltene precipitation and aggregation are explained by thermodynamics phase behavior and kinetic studies respectively; however, the deposition of asphaltene is related to a series of mechanisms related to surface characteristics, phase behavior, flow regime and rheological properties of the flow.During the precipitation process, a portion of asphaltene particles become unstable which leads to instability and phase separation.Increasing the tendency of particles to aggregate and further aggregation by collisions to each other form the aggregates which is considered by Population Balance Method (PBM), in this study, in addition to the work of Vargas et al. (2010) and Kurup et al. (2011).In this work, using the proposed modeling approach, a sensitivity analysis for different parameters is presented for particles average size, flow rate, viscosity, and well radius.Laboratory data and other required inputs are gathered from Tavakkoli (2013) and the results were compared with results obtained by ADEPT.

Methodology
The focus in this work is to extend a deposition model according to the work of Kurup et al. (2011) including the effective mechanisms in the macroscopic transportation of asphaltene particles through the wellbore (convective flow of micro aggregates, aggregation, diffusion, precipitation and deposition).This model has the ability to predict the probability, location and thickness of asphaltene deposition in the flow lines for different operating conditions.Alongside the upward flow in a wellbore and with respect to the different factors and probabilities, the precipitated particles stick to each other and the aggregates are formed with different sizes.It is also worth noting that the tendency to aggregate between asphaltene masses for different oil samples and with identical precipitant are different.This issue indicates the importance of the oil composition on the precipitation and aggregation rate.Finally, it is likely that these aggregates deposit on the walls of the wellbore due to the contact with the physical surfaces available in the flow lines or previously deposited asphaltenes.In other words, without a good understanding of aggregation process, the results for dynamic deposition cannot trusted.
As shown in Figure 1, in this work similar to ADEPT simulator (Kurup et al., 2011), the simulation consists of three modules, thermodynamics, kinetics and dynamics.To obtain amount of asphaltene precipitation along axial length of the wellbore, the pressure and the temperature profiles are obtained through codes for the velocity profile of the two-phase oil and gas flow (velocity calculation for one phase and slippage velocity calculation for the other phase according to flow pattern).The same approach for the calculation of the temperature, pressure and velocity was used in the work of Kurup et al. (2011).Afterwards, for the calculation of soluble asphaltene concentration in equilibrium conditions (i.e., C eq ) in the thermodynamic module, PC-SAFT equation is used.It is in the form of dimensionless concentration with respect to the concentration of dissolved asphaltene at inlet conditions.This equation of state has been successfully used by David Ting et al. (2003), Gonzalez et al. (2007), and Vargas et al. (2009) for asphaltene precipitation.This EOS is effective in the study of phase behavior of asphaltene in oil, because of its ability in the modeling of complex molecules like asphaltene in the high pressure and temperature conditions.In PC-SAFT equations, characterization of crude oil based on SARA analysis was used to aid in computations reduction.Finally, soluble asphaltene concentration was used as an input for the dynamic module.Before any extension, the approach of Kurup et al. (2011) was developed into 2D using our code.In 2D material balance equation used in this study, different terms were included for precipitation, aggregation, and deposition for each control volume of the wellbore discretized domain (Fig. 2), as in Vargas et al. (2010).Despite its simplicity, this dimensionless equation presents the main processes that asphaltene particles experience: This PDE is subjected to different boundary and initial conditions (Kurup et al., 2011): In these relations, C is the dimensionless concentration (c/c 0 ) of precipitated asphaltene particles.c 0 is inlet mass concentration of asphaltene solubilized in the oil phase.i is the size of asphaltene particles.C f is dimensionless concentration of asphaltene dissolved in oil phase.C eq is the dimensionless maximum concentration of asphaltene in oil at equilibrium condition.u is dimensionless time.h is dimensionless radius (r=R).r is the local radius inside the wellbore and R is the radius of wellbore.j is dimensionless length (z=L).z is local distance alongside the wellbore and L is axial length of the wellbore.Da P is precipitation Damkohler number (Da P is calculated without considering size of aggregates), Da D is deposition Damkohler number, and K diss is kinetic of dissolution factor.The dimensionless parameters are described as: Da p ¼ rate of production of microaggregate by precipitation rate of transport of microaggregate by advection ; Da D ¼ rate of consumption of microaggregate by deposition rate of transport of microaggregate by diffusion ; where Pe is Peclet number.v z is average axial velocity of the oil flowing through wellbore, k p is precipitation kinetic constant, k D is deposition kinetic constant, k B is Boltzmann constant, d p is particle diameter, and v is kinematic viscosity.
One of the differences of this study to the work of Kurup et al. (2011) is related to the correlation used for the diffusion term.In the wellbore tubing, there are different flow layers including, turbulent flow core, buffer layer and laminar (wall) sub layer each with different diffusivity coefficients.The asphaltene existing near the tubing wall (i.e., the laminar boundary layer) deposits at a constant rate with the diffusivity coefficient proportional to the distance from the wall.This is another clue for the importance of the 2D modeling.To include it in this research, the correlations presented by Escobedo & Mansoori (2010) are used: In this relation, D B is the Brownian diffusivity and D t is eddy diffusivity.Eddy diffusivity is a function of radial distance, and, is calculated for each of the three main regions in the wellbore, i.e. laminar sub layer, buffer layer, and turbulent core, with different correlations as described in following. where: On the other hand, after applying different diffusivity coefficients for each layer in the wellbore, the velocity in the circular pipe is calculated using turbulent flow distribution function: In Equation ( 2), the term on the left-hand side shows changes in the absorbed asphaltene on the wall with time.In this study, it is assumed that all of the absorbed asphaltene will deposit similar to Kurup et al. (2011).The first term on the right-hand side shows the asphaltene particles advection considered only in axial direction, as diffusion is dominant in radial direction.The second term shows asphaltene transport in radial direction under diffusion mechanism.Diffusion in axial direction is not included as the advection is dominant in the axial direction.The third and fourth terms determine the kinetics of aggregation and precipitation, respectively.The deposition process with constant rate is included into the calculations at the third boundary condition.The parameters description in the Equation ( 2) and the mathematical models used in this work are consistent with the work of Vargas et al. (2010) and Kurup et al. (2011).In the dynamic module, obtaining the particles aggregation kinetic is a prerequisite.However, in this work, this parameter is not constant.Aggregation rate is calculated from PBM for each control volume in the wellbore discretized domain, because, different particle sizes lead to different amounts of absorption.In this work, the geometric PBM form Maqbool (2011) is used.The r A for each class of size is evaluated as following.The mass balance equation was simultaneously solved with equation set describing the size distribution evolution as (Maqbool, 2011): where: In these relations, K i,i is dimensionless collision kernel between two ith aggregates (K i,i is the collision kernel between two ith aggregates).R 0 is geometric scaling factor and is equal to 2.0.R g is universal gas constant and T is absolute temperature.u is viscosity of the medium.d i is the diameter of colliding aggregates ith, b is the efficiency of collision, e is packing factor for random close packing of spheres, Vm is nanoaggregate volume, Mw A is molecular weight of asphaltene, N A is Avogadro constant, N agg is the number of asphaltene molecules in a nanoaggregate, and r 0 is oil density.In Equation ( 19), the term in the left-hand side represents the net rate of generation of ith aggregate obtained by the summation of terms for different generation and depletion mechanism.The first and second terms on the right-hand side show the rate of the generation of ith aggregate, as a result of the reactions of two or more particles with total volume equal to ith aggregate.The third and fourth terms on the right-hand side show the rate of depletion of ith aggregate due to collision with one another or other aggregates.This equation, neglects the breakage and fragmentation terms, otherwise, two adjusting parameters would have been required, making the modeling process more difficult.In addition, the source of precipitation enters the equation continuously and the amount of concentration reduces due to the deposition.The kinetic parameter, ki,j represents the collision kernel between two aggregates with different sizes (m 3 kmol À1 s À1 ).Since it is assumed that nanoaggregates act like primary units for asphaltene aggregations, the Brownian aggregation kernel is used in the calculations.When the pressure changes along the wellbore, the nanoaggregates are destabilized and aggregation produces larger particles.Therefore, the precipitation onset for asphaltene is actually the onset of detection of asphaltene aggregates as they grow to form larger sizes.Besides, with a simplifying assumption, that asphaltene aggregates are spherical and have a packing factor (e) which is close to 0.637 for spheres aggregating by accident.Therefore, the only unknown parameter in modeling of PBM is the efficiency of collision (b) which depends on the interparticle interactions.The parameter b, can be estimated from the regression through laboratory data.However, the model is not sensitive to it.Therefore, its value was taken from Maqbool (2011).Finally, after finding C, dimensionless concentration of precipitated (phase separated) asphaltene particles, in the wellbore, dimensionless deposition flux, j A , is calculated as: The steady state approach is used to solve the equations.In other words, deposition rate remains constant with the time.Therefore, the pressure and temperature profile and the amount of soluble asphaltene at equilibrium conditions for the entire well length are calculated once.The deposition of asphaltene in the well tubing decreases the available area for the flow and increases the velocity of the flow which is included in the calculation.This leads to pressure drop in tubing, pressure changes and other consequences.However, the mentioned changes can be ignored for more calculation speed without losing the solution accuracy.The wellbore pressure drop, relating to the friction, comprises 10% of the whole pressure drop along the vertical well.In addition, the asphaltene deposition rate is slow (in the order of 10 À7 ), and therefore it doesn't create an intense change in pressure and temperature.According to Kurup et al. (2012), there is not a significant difference between steady state and pseudo steady state.However, these differences increase with time at the end of forming depositions.Therefore, considering the calculation time, using a steady state method is efficient for the sake of simple and fast modeling.Hence, in this work, it is assumed that the asphaltene deposition on the tubing's wall doesn't change the pressure and temperature profile in the well.It should be noted that, the effect of diameter reduction because of asphaltene deposition is included on reduction of available area for deposition for the next time step, with subtracting asphaltene deposition thickness from the well radius.In this study, after each time step, radius of the wellbore is corrected as: Where d is the height of deposited asphaltene.Its value is obtained as follows: In addition, this quantity can be obtained from dimensionless parameters according to: In this work, the well was discretized in radial and axial dimensions and finite volume method was used to solve dimensionless partial differential Equation (Eq.2).In this model, sensitivity analysis for different parameters were investigated for the particles average size, flow rate, viscosity and well radius.

Results and discussion
The results description is divided into three main sections.At first, the integrated basic thermodynamic precipitation, and dynamic deposition models are compared with ADEPT results of Kurup et al. (2012) (1D) to evaluate the implementation of the modifications on the data of Tavakkoli ( 2013).After a verification step, the aggregation model of PBM is included.Finally, a sensitivity analysis is performed on the asphaltene particles average size, flow rate, viscosity and well radius.
In the stage of model extension in 2D for a single particle size, a mesh independency study was performed and shown in Figure 3.It should be noted to present the sensitivity of the selected mesh to the timestep length, the mesh independency study results are presented for different timesteps.The fine mesh dimensions are Nz = 200, Nr = 50 and Dt = 3 days.The relative differences of deposited asphaltene of several mesh dimensions form the reference mesh are showed in Figure 3.The mesh with smallest dimensions which does not produce a large deviation was selected for further study.Its dimensions are Nz = 40, Nr = 15 Dt = 30 days.

Precipitation and deposition modeling verification
The experimental data of Tavakkoli (2013) were used and the modeling results were compared to his work using ADEPT of Kurup et al. (2012).The wellbore fluid composition is presented in Table 1.The SARA analysis of the wellbore fluid is in Table 2.The wellbore structure information with all of the used parameters are shown in Table 3.The characterization data for PC-SAFT are calculated as shown in Table 4 according to Tavakkoli (2013).The profile of asphaltene content in the equilibrium condition (Ceq) along the wellbore was obtained using vertical flow pattern of the reservoir fluid and compared with the work of Tavakkoli (2013), as shown in Figure 4 (TOP).Density of the reservoir fluid reduces along the tubing with decrease in pressure and temperature.The effect of expansion of the light components of the reservoir fluid is the same as light components injection into the crude oil.Both of them cause asphaltene to precipitate.From precipitation onset, asphaltene starts to precipitate and the precipitation amount increases up to the bubble point.At pressure below the bubble point, light components vaporize gradually out of the oil.This event causes the solubility of asphaltene in the crude oil to increase until the lower onset point in which asphaltene dissolves again completely.The calculated solubility parameter of the oil supports this speculation (Bahrami et al., 2015).In fact, for each reservoir fluid, there is a minimum solubility point for asphaltene at fluid bubble point as mentioned by Tavakkoli (2013).It should be noted that, according to the figure, the difference(which is not considerable) in the equilibrium concentration between this work and the work of Tavakkoli ( 2013) is due to the fact that both works made separate regression, of PVT data.In fact, for PVT data in this work, we made the regression process using our code.As the regression is a multivariable optimization, it does not produce unique answer.The next step, is the deposition calculation.Since all of the deposition model improvements, in this work, are based on the work of Kurup et al. (2011), the correct implementation of the method is verified at first using the data of Tavakkoli (2013).The deposition flux profiles are plotted in Figure 4 (BOTTOM) in comparison to Tavakkoli (2013).The first figure is from 1D (ADEPT simulator) of Kurup et al. (2011) in Tavakkoli (2013), the second is prepared with this work, 2D approach (constant diffusion coefficient) and C eq from thermodynamic precipitation approach of this work.Two figures are similar.This means, asphaltene deposition profile of 2D becomes very close to 1D (ADEPT simulator) of Kurup et al. (2011) in Tavakkoli (2013).The sensitivity analysis shows that the constant diffusion coefficient is so large that, it doesn't show the radial distribution effects.It is one of the reasons of the similarity between 2D and 1D (ADEPT simulator) of Kurup et al. (2011) in Tavakkoli (2013) In the next step, diffusion coefficients of Escobedo & Mansoori (2010) were utilized (Eq.14).The size effect results in different penetration rates for each size and each layer of the wellbore were calculated in 2D model.Therefore, the particles absorb with different rates compared to those come from the center to the layer near the wellbore.According to Figure 5, the deposition profiles along the wellbore are obtained by choosing different sizes for asphaltene particles, which results in different deposition flux.The deposition flux decreases for larger size distribution of asphaltene particles, which shows the importance of applying the size effect into the absorption equation.Nanometer is the order of magnitude in which asphaltene begins to precipitate (Eskin et al. (2011) showed that this size could be between 1.6-5 nm).In addition, Eskin et al. (2011) declared that in simulating asphaltene deposition inside Couette device, particles larger than 10 micrometers do not settle down.This is the reason for choosing this range for the current study.The differences between 2D and 1D simulation, for each particle size, are provided in Figure 6.This figure shows the asphaltene concentration in the layer close to the wellbore (it was obtained when we discretized the well) proportional to the average concentration along the radius versus particle size.Their values were averaged over the whole length of the wellbore.The asphaltene concentration in the layer close to the wellbore and the average concentration along the radius are equivalent to 2D and 1D models, respectively.In fact, the differences between 2D and 1D simulations become larger for the bigger sizes.When the size of the particles is small, the deposition flux on the wall is equal to the average volume flux in that depth of well, and 1D simulation is similar to 2D.The deviation becomes greater with size.For the aforementioned reasons, there is no doubt that particle size has critical effect on the diffusion coefficient as well as deposition flux.

Aggregation modeling using PBM
To include different particle sizes effect, the asphaltene particles were classified into 12 classes according to size.Including, the particle sizes bigger than the last class in our classification didn't have any effect on the deposition of asphaltene.Afterwards, the absorption equation (Eq.2) and population balance (Eq.19) were simultaneously solved for each size in all of the blocks.As the deposition rate was considered steady and time independent, a number of proper time intervals were selected to reach steady state conditions.In Figure 7, the thickness of deposited asphaltene along the wellbore versus time is shown.As it indicates after one month, the radius of the wellbore reduces about 0.015 cm, which is lower than the model with constant diffusion coefficient and no population balance.Figure 8 shows the contribution of particles with different sizes in the deposition flux.The shape of deposition flux for each size is almost like the general trend.It is shown that, initially, the deposition rate is increased with size and after a specific point, it is inversed.Besides, the results show that the asphaltene particle of class eight has the highest deposition rate.In fact, it is concluded that for this case, size distribution effect is well associated with the effect of size abundance and they cannot be separated.Moreover, the distribution of bigger sizes in the wellbore is higher in this stage of deposition; therefore, increasing in deposition flux is to some extent due to the abundance of those sizes.In Figure 9, the results of the 2D simulation with population balance are compared with the results of Tavakkoli (2013) (1D with constant diffusion coefficient) for 30 days.In this figure, deposition rate and deposition trend in the models are compared.The results clearly indicate that the deposition rate in comparison with ADEPT in Tavakkoli (2013) is smaller.This difference in deposition rate is due to applying the kinetic of aggregation effect (i.e.size distribution) in the simulation.In ADEPT (Kurup et al.,   L. Eyni et al.: Oil & Gas Science and Technology -Rev. IFP Energies nouvelles 73, 25 (2018) 2011), the amount of the kinetic aggregation constant was considered so small that the probability of sticking the particles to one another was almost zero.As the ADEPT simulation considers the single particle size, the selection of such kinetic aggregation constant is reasonable.Nevertheless, by applying size distribution and assigning the probability of particle sticking, the larger aggregates are formed.As it is mentioned before, the diffusion rate for the particles with larger size is smaller; therefore, it has a lower deposition rate.Thus, the size distribution, with larger particle sizes, leads to a lower deposition rate.The trends of diagrams are similar with the minor difference because of kinetic of aggregation and size distribution effect.In fact, as it was shown, each class of considered size distributions has its own trend and the overall effects results in a distinct trend from ADEPT.To show the size effect separately in Figure 10, the deposition flux is plotted for different average sizes for 30 days.The deposition flux decreases with increasing in average size of asphaltene aggregates.It should be mentioned that some of inhibitors, which are used for preventing the deposition of asphaltene, break the bigger asphaltene aggregations and in fact, unlike their injection goals, they lead to asphaltene deposition worsening Vargas et al. (2010).To demonstrate this effect, for average amount of different size distributions, average deposition flux was obtained and plotted in Figure 11.It means that the deposition flux starts from a maximum amount and reduces until it reaches its maximum change at a point and then the slop of change reduces so that, finally, it reaches the minimum amount and after that, there is no change in deposition flux.
The importance of 2D modeling with PBM (this work), in which a range of particle sizes is included, over the 2D modeling without PBM, in which a fixed size equal to average particle size, is demonstrated for 30 days in Figure 12.

Sensitivity analysis
Using the extended deposition model, the effect of different parameters on the deposition can be analyzed.In this work, a sensitivity analysis of three different flow rates was    Comparison of 2D modeling using with fixed particle size (size effect is included on the diffusion coefficient) and 2D modeling using variable particle size (PBM) for 30 days.conducted using two different particle size distributions as showed in Figure 13.According to this figure, with increasing the flow rate within the same time (30 days), two effects on deposition flux can be seen; the first one (nonaligned effect) is convection, in which, as it increases, it carries more particle and prevents deposition.Therefore, the amount of deposition increases with increasing the flow rate.This effect is in part clear at the ascending part of the diagram.In second mechanism (aligned effect), velocity gradient in the radial direction increases with flow rate and it leads to increase the concentration gradient in the radial direction (At the beginning, asphaltene particles are generated by precipitation.This is modeled as a volume source.Therefore, all the points in radial direction have the same chance.Suppose asphaltene is precipitated in the wellbore when the onset point is reached.Because of velocity gradient, particles in the middle of the wellbore flow faster than the particles near the wall.Therefore, concentration distribution is stretched.This results in sharp gradients perpendicular to the flow, which diffusion can then iron out.The use of radial gradient in velocity distribution as well as radial diffusion coefficient in the first equation manages this effect.).Then this gradient leads the particle to move to the wall of wellbore.In fact, the penetration of particles increases toward the wall.The second mechanism prevails the first one in a part of the path.The reason of this dominance is the increase in suspended asphaltene concentration along the well.Therefore, the usage of chokes with smaller opening diameter can help decreasing the deposition.It is worth mentioning that in 1D model of Kurup et al. (2011), only, the effect of first mechanism is investigated.Therefore, it is also another reason for the superiority of the 2D model.The first mechanism in ascending part of figure for the bigger average size of asphaltene aggregates is lower than smaller average size of asphaltene aggregates, which shows that the investigation on the effect of parameters without considering average size of asphaltene aggregates may not be completely correct.Therefore, proposed modeling can enhance our understanding of size effect on deposition.After that, a similar sensitivity analysis was conducted for two different particle size distributions, for three different viscosities, as showed in Figure 14.The viscosity influences two parameters, diffusion and kernel, in which, the former and the latter have favorable and adverse effects on deposition promotion, respectively.As it is mentioned before, diffusion does not show the radial distribution effects.In fact, it reduces the mass transfer resistance in the direction of the well radius.Hence, the amount of deposition increases with increasing diffusion.Besides, the aggregation increases with the kernel.Therefore, the amount of deposition increases with decreasing the kernel.Due to the aforementioned reasons, with increasing the viscosity within the same time (30 days), the deposition flux increases in the smaller average size of asphaltene aggregates.However, with increasing the average size of asphaltene aggregates, the mutual effects of the viscosity and average size of asphaltene aggregates on deposition can be seen.Once again, it shows that the investigation on the effect of parameters without considering average size of asphaltene aggregates may not be completely correct.In the next step, another sensitivity analysis was conducted for three different well radiuses at two different particle size distributions, as showed in Figure 15.The flow rate   increases with the radius in a constant velocity.In addition, as we showed before, the amount of deposition increases with increasing the flow rate.Therefore, the amount of deposition increases with decreasing the well radius.Finally, in all sensitivity analyses, the deposition flux decreases with increasing in average size of asphaltene aggregates for each constant flow rate, viscosity, well radius and temperature, which confirms the result of Figure 10.

Conclusion
In this work, the dynamic deposition of asphaltene in the oil wells including the precipitation, aggregation and deposition of asphaltene with changes in the temperature, pressure, and size along the axial length of the wellbore was simulated and analyzed through an improved modeling approach.The proposed modeling approach extends the available models to 2D numerical simulation including different particle size effects through PBM.These three parts were integrated in a simple and comprehensive 2D simulator.The model was verified with a common asphaltene deposition tool to evaluate its basic modeling capabilities to obtain the asphaltene deposition profile along tubing.The analysis showed that larger aggregate sizes, makes significant reduction of the diffusion rate and consequently lower deposition rate.Moreover, the relative importance of the aggregation and deposition phenomena were analyzed through the proposed model including size distribution effect.This result revealed the impracticality or even adverse effect of some asphaltene precipitation inhibitors.Besides, the critical size effect of asphaltene aggregate sizes and mutual effects of flowrate on the deposition process were included in the parametric study.Effects of flowrate were related to the effect of convection and velocity gradient in the radial direction, which the former and the latter have favorable and adverse effects on deposition promotion, respectively.Moreover, effects of the critical size were related to the effect of diffusion and collision kernel, which as we explained in results and discussion section the former and the latter have favorable and adverse effects on deposition promotion, respectively.At the end, a sensitivity analysis of viscosity and radius were conducted using two different particle size distributions.In fact, mutual effects of viscosity and average size of asphaltene aggregates make a different result for various particle size distributions in sensitivity analysis.Therefore, to predict the asphaltene deposition accurately in production wells with pressure and temperature depletion, it is recommended to consider the asphaltene size distribution.

Fig. 1 .
Fig. 1.Schematic diagram of the all steps of developed simulator, Inputs*: Oil & Gas characterization-Flow rate & Geometry-Precipitation, Aggregation, and Deposition constant rate-Experimental data

Fig. 4 .
Fig. 4. (TOP) Profile of asphaltene content in equilibrium condition (comparison to (Tavakkoli, 2013)) along the wellbore, (BOTTOM) Asphaltene deposition results of 2D model with constant diffusion coefficient and precipitation model parameters in this work and 1D ADEPT in Tavakkoli (2013) profile along the wellbore

Fig. 6 .
Fig. 6.Asphaltene concentration in the layer near to the wellbore proportional to the average concentration along the radius versus to size of particles.m versus d, m=Asphaltene concentration in the layer near the wellbore proportional to the average concentration along the radius, d p = size of particles.Their values were averaged over the whole length of the wellbore.

Fig. 8 .
Fig. 8. Asphaltene deposition profile (proposed 2D modeling approach with Escobedo & Mansoori (2010) size dependent diffusion coefficients and PBM) along the wellbore for each considered particle sizes; contribution of particles with different size in the deposition flux.

Fig. 10 .
Fig. 10.Asphaltene deposition profile (proposed model) along the wellbore for different average sizes for 30 days.

Fig. 11 .
Fig. 11.Average deposition flux versus average amount of different size distributions.

Fig. 12 .
Fig.12.Comparison of 2D modeling using with fixed particle size (size effect is included on the diffusion coefficient) and 2D modeling using variable particle size (PBM) for 30 days.

Fig. 13 .
Fig. 13.Sensitivity analysis to the wellbore flow rates for different particle size distributions after 30 days.

Fig. 14 .
Fig. 14.Sensitivity analysis to the viscosity for different particle size distributions after 30 days.

Fig. 15 .
Fig. 15.Sensitivity analysis for the effect wellbore radius for different particle size distributions, deposition profile after 30 days.

Table 3 .
(Tavakkoli, 2013) properties and parameters(Tavakkoli, 2013) Boltzmann constant, 1.38066 Â 10 23 J/K K i,j dimensionless collision kernel between ith size aggregate and jth size aggregate R radius of the wellbore, cm r the local radius inside the wellbore, cm R ́geometric scaling factor R g the universal gas constant, 8.314 J/mol/K