- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents
Advances in Mechanical Engineering
Volume 2012 (2012), Article ID 973652, 10 pages
Numerical Study on the Effect of Wax Deposition on the Restart Process of a Waxy Crude Oil Pipeline
1National Engineering Laboratory for Pipeline Safety/Beijing Key Laboratory of Urban Oil and Gas Distribution Technology, China University of Petroleum, Beijing 102249, China
2PetroChina Pipeline Research and Development Center, Langfang 065000, China
Received 16 October 2012; Accepted 20 November 2012
Academic Editor: Bo Yu
Copyright © 2012 Qing Miao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
A new model for wax deposition is established based on the loop tests data of 9 typical crude oils in conjunction with the concept of the effective wax precipitation which occurs in the heat boundary layer in the pipeline. The thickness of the wax deposition of a waxy oil pipeline is predicted by the established model, and the results of the model agree well with the field data. Subsequently, the models for the shutdown and restart processes of a waxy crude oil pipeline are established. Then the restart process of a wax crude oil pipeline is investigated numerically, the temperature drop during the shutdown process and the transient inlet pressure are presented, and the effect of the wax deposition on the safety of the restart process is clarified.
Wax precipitation in crude oil pipelines at or below the wax appearance temperature (WAT) can lead to gelling that inhibits flow by causing significant non-Newtonian behaviour and increasing effective viscosities as the temperature of a waxy crude oil approaches its pour point . When the temperature of the pipeline wall is below WAT, the deposition of a layer of paraffin molecules would form near the pipe wall, constricting the flow. In the practical pipeline transportation engineering, pigging is often conducted periodically for a hot waxy crude oil pipeline to keep its throughput and operation pressure relatively steady. But for some hot crude pipelines, the throughputs of which decline with the reductions of the outputs of the upstream oil fields; pigging operations are not done regularly or even not done at all for a long period. For these pipelines, the existence of the wax deposition in them could help reduce the increment of heat loss resulting from the reduction of the throughput and even make it possible to keep the pipeline operating safely without increasing the heat capacities of heat stations. But the heat preservation of wax deposition does not imply that it will also keep or improve the safety of a long-term unpigged hot crude pipeline during its shutdown and restart. Since the wax deposition decreases the accumulation heat of a hot crude pipeline in the soil around itself during running, which is a disadvantage to the thermal conditions after the shutdown of the pipeline, the influences of thicker wax deposition on the shutdown and restart of a long-term unpigged hot crude pipeline cannot be deduced directly from normal operation and should be investigated quantitatively in detail for assessing the safety of such a pipeline.
Some researchers, such as Del Carmen García and Urbina  and Del Carmen García , have studied correlations between the properties of crude oils and their flowing properties, including the precipitation and deposition of wax during flow. Models have been developed to predict the onset of wax precipitation and the deposition of wax along pipeline walls. However, accurately modelling deposition in pipelines can be a complex and difficult undertaking, because, while precipitation is mainly a function of thermodynamic variables such as composition, pressure, and temperature, deposition is also dependent on flow hydrodynamics, heat and mass transfer, and solid-solid and surface-solid interactions . Some researches about the wax deposition could be referred to [5–15].
The studies on the calculation of the restart process of oil pipeline are extensive. Some valuable researches can be found in [16–21]. But the impacts of wax deposition on the restart process are rare. In this paper, the impact of the wax deposition formed during the running without pigging on the safety of the restart process of a pipeline is analyzed and some significant conclusions are drawn.
The layout of this paper is as follows. The characteristic of wax deposition is to be investigated first. Then the models for the shutdown and restart processes of a waxy crude oil pipeline will be established. Lastly, the effects of wax deposition on the shutdown and restart processes of a waxy crude oil pipeline will be investigated.
2. Mathematical Model for Predicting the Wax Deposition
When the crude oil is transported in the pipeline, since the ambient temperature at a certain length of the pipeline is usually lower than the crude temperature at the same length, the temperature near the inner pipe wall cools faster than that far from the inner wall and causes a temperature gradient pointing to inner pipe wall. And the temperature at inner pipe wall firstly reduces to WAT, and the wax precipitates near the wall in the first place and then crystallizes and deposits on the wall. The wax deposition near the wall results in the reduction of the liquid wax concentration in this area and causes a solved wax concentration gradient pointing to the inner wall, which impels the liquid wax molecules in the central area of pipe cross-section moving to inner wall and precipitating and depositing there continuously. This temperature gradient induced solved wax concentration is just the impetus of the wax deposition, and thus it is deduced that only the wax precipitation occurred in a region where an oil temperature gradient exists and contributes to the wax deposition, or say, is effective to the wax deposition. From this issue, even if all the oil temperatures along a radius in some pipe cross section have reduced below the WAT so the wax precipitation has occurred anywhere in the cross section; however, only the wax precipitation in the area near the pipe wall in which heat flux points outwards has contribution to the wax deposition on the pipe wall, instead of the whole amount of wax precipitation in the cross section, while the wax precipitations in other areas will flow downstream in wax crystals forms.
According to the above analyses, the concept of effective wax precipitation (EWP) is proposed, which associated with the wax precipitation occurred in the region where an oil temperature gradient exists and contributes to the wax deposition. Since the flow states in most industrial crude pipelines are in turbulent regimes, the temperature gradient only exists in laminar region near the wall. As for the laboratory loop test for wax deposition (which is the focus of our study), because the length of the test tube in a loop is usually about 1.5 meters to 2 meters and the heat boundary layer in the tube would not expand fully according to heat transfer method, the temperature gradient also exists in the heat boundary layer near the wall, shown in Figure 1.
Since the ambient temperature around the test tube changes axially in practical test process, the oil temperature varies along the test tube and the heat boundary layer thickness also changes along the test tube. In order to calculate the EWP in test tube accurately, the length of the test tube would be divided into several segments. Then the thickness of the heat boundary layer of each section could be calculated according to heat transfer method and thus the whole effective mass of crude which flows through the heat boundary layer during a test can be calculated according to the following formulae: where indicates the thickness of heat boundary layer of the th section; and represent the mass flow rate for the th section and the mass flow rate of the whole test tube; indicates the effective inside diameter for the th section; represents the test time for wax deposition; indicates the whole effective mass of crude which flows through the heat boundary layer during a test.
Once is obtained, the EWP for the tested crude during the test time can be determined from it according to the wax precipitation characteristics of the crude.
For studying the relationship between wax precipitation and its deposition, a variable named effective deposition ratio (EDR) should be defined to connect the EWP with the wax deposition, which is defined as where is the whole amount of wax deposition and is for EWP. In the following, effective deposition ratio will be abbreviated as EDR.
In order to detect the hidden relationship between EDR and EWP, a test loop shown in Figure 2 was devised to obtain the necessary data. In the figure, and are two differential pressure gauges used to measure the pressure drop of test tube 1 and 2; cooling circulator 1 and 2 are temperature control equipments used to control the ambient temperatures around test tube 1 and 2; mass flow meter is used to measure the mass flow rate of the loop; oil jar provides sample source and compensates wax concentration; pump provides power to keep sample circulating.
With the experiment data of nine different crude oils, it is discovered that the relationships between effective deposition ratio (EDR) and effective wax precipitation (EWP) for each kind of representative crudes studied and tested in the test loop obey the same rule; that is, the EDR decreases with the increase of the EWP and, when the EWP is less, the EDR changes rapidly; when the EWP increases to some extent, the EDR begins to change slowly. Further, put all the EDR and EWP data from 9 kinds of crude oils together, the tendency concluded above is clear at a glance (see Figure 3). It is further found that the statistic rule on the relationship of EDR and EWP is obvious and could be expressed in power law mathematically as follows: where represents the EWP; and are two constants.
Based on the results above, the significance for introducing the concept of EWP is emphasized. Through the definition of EWP, the universal rule on the relationship between EWP and wax deposition available to many crude oils is found quantitatively. The concept of EWP uncovers the internal connection of the wax deposition with the wax precipitation and provides a method to divide the part contributing to wax deposition from the whole wax precipitation in pipeline or in test tube. More important, the rule common to the wax deposition of many kinds of crude oils in the test loop makes it possible that a common wax deposition model could be created for many kinds of crude oils and a large number of loop tests necessary for wax deposition prediction previously could be left out.
In fact, the outline of the common wax deposition model has been given in formula (3) theoretically. But that form might not be used for practice yet for the reason that the EWP is unknown for a new kind of crude never tested in the test loop. From the definition and calculation of EWP, it is explicit that EWP is connected with the thickness of heat boundary layer and the thickness of heat boundary layer strongly depends on the velocity and temperature distributions in the test tube at different test conditions, so the factors such as the viscosity of crude, shear stress at inner pipe wall, radial temperature gradient at pipe wall, oil temperature in central area of test tube, and the velocity of oil flow which have impacts on velocity and temperature distributions would all influence the determination of EWP. In other words, the EWP is a function of all factors mentioned above. From the above discussion, the practical common wax deposition model should have the following form: where are constants needed to be determined; and are the viscosity of tested crude at test condition and shear stress at inner pipe wall, respectively; and are radial temperature gradient at pipe wall and the amount of wax precipitation from tested crude for unit oil temperature reduction, respectively; and are oil temperature in the central area of test tube and the average velocity of oil flow, respectively.
On the basis of wax deposition loop test results of 9 kinds of crude oils and applying the method of linear multi-variable stepwise regression to check and analyze the formula (4), the final form of the common wax deposition model is determined as follows: Formula (5) is the common wax deposition model which is available to most kinds of crude oils. From now on, an efficient and sample-saving wax deposition model can be created for any crude in the case of no wax deposition loop tests.
3. Mathematical Models for the Shutdown and Restart Processes of a Waxy Crude Pipeline
During pipeline shutdown, heat transfer in crude oil can be divided into three stages as (1) natural-convection-dominated stage, (2) joint control of natural convection and heat conduction, and (3) pure heat conduction . With the following assumptions the temperature drop process is then simplified as a heat conductive process. (a)The variation of the specific heat capacity of the crude oil is applied to describe the latent heat of wax crystallization during shutdown.(b)Stagnation point introduced in the second stage describes the internal pipeline as a liquid region with natural convection and a solid region with thermal conductivity.(c)Natural convection can be solved as conduction after an equivalent thermal conductivity is applied . By this means, we can convert the calculation of natural convection into a conductive heat transfer. The equivalent thermal conductivity is dependent on the natural convection coefficient, the temperature of oil, and the temperature gradient on the interface. Then the heat transfer equation of the oil is then shown as in which, is a temperature gradient at the interface.
Though the equivalent thermal conductivity is not a real heat conduction coefficient, it guaranteed energy conservation and greatly simplified the calculation. This approximation is verified by field data of crude oil pipelines .
Heat transfer equations of wax deposition layer, pipe wall, and corrosion protective coating as well as the heat conduction equation of the soil are the same as those in the normal operation process. Matching conditions are shown in (7) while boundary conditions are the same as those in the normal operation process. Consider where to are shown in Figure 4 and CPC in Figure 4 is the abbreviation for the corrosion protective coating.
During the process of restart, the thermodynamic model is the same with that during normal operation, which can be referred to . In addition, the hydrodynamic model for the restart process is as follows.
Momentum equation: where is the average velocity of the oil flow; is the restart pressure along the line; is the inclination angle of the pipeline; is the density of the crude; is the gravity acceleration; is the effective inside diameter of the pipe; is the shear stress at inner pipe wall, Pa.
The shear stress at inner pipe wall depends on the flow pattern and flow behavior of the crude and can be determined according to the following formulae.(1) When the crude belongs to Newtonian fluid: : : (2) When the crude belongs to non-Newtonian fluid: : : where and are the Reynolds number and Metzner-Reynolds number, respectively; is the volumetric flow rate; and are the viscosity of the crude and friction coefficient, respectively; and are consistency coefficient and flow behavior index of pseudoplastic fluid, respectively.(3) When the crude presents a thixotropic behavior, can be determined by solving the Houska thixotropic model. Initial and boundary conditions: where and are two constants.
4. Numerical Methods
Here a thermal influence zone  is adopted as a computational domain and then the numerical simulation is performed. The pipeline is divided into segments in axial direction, shown in Figure 5. A Delaunay triangulation method is applied in a software package to generating the grids in a Cartesian coordinate system (as shown in Figure 6) specified by the parameters of pipeline burial depth and radius defined by users. The soil domain was divided into nonoverlapped triangles by the method mentioned above. In order to simulate the temperature field more accurately, a mesh that is dense near the pipeline but coarse far away from it is used to improve the computational efficiency while resolving the temperature gradient that is high near the pipeline and low away from it. Grid independence was tested by applying three sets of grids of different sizes, and the numerical simulation results show that the results are independent of mesh size. Due to the temperature drop during the shutdown and restart processes, the temperature is not uniform in the radial direction. The temperatures of oil, wax deposit, and corrosion-inhibiting coating were calculated, respectively, in a polar coordinate system shown in Figure 7. Then the governing equations are discretized by the combination of the finite volume method and the finite difference method, and details can be found in [24, 26, 27].
5. Results and Discussion
5.1. Characteristic of the Wax Deposition of a Waxy Crude Oil Pipeline
A hot crude oil pipeline transporting Daqing crude in the Northeast of China had not been pigged over the past four years due to the reductions on throughputs. The hot crude pipeline is 436.3 km long and its diameter is mm. The maximum allowed inlet oil temperature of each pipeline interval is 70°C and the maximum allowed operation pressure (MAOP) is 4.32 MPa. The design annual throughput of the pipeline is . During the past four years, the throughputs of this pipeline reduced year after year from to , then to and . With the reduction of the throughput and the growth of the wax deposition, the thermal insulation effect of successively increasing wax deposition appeared gradually, which kept the outlet oil temperatures of each pipeline interval almost unchanged without adding more heat capacities to the heat stations. Based on the mathematical wax deposition model created for Daqing crude, numerical simulation for the wax deposition of this long-term unpigged hot crude pipeline is conducted. The results for the first interval of the pipeline are shown in Figure 8.
From Figure 8, it can be found that the wax deposition is not identical along the pipeline. In the middle of the fifth year during no pigging, the maximum wax deposition of the first pipeline interval is about 117 mm and the wax deposition at the outlet of the same interval is near 80 mm. The simulation shows that the position where the maximum wax deposition occurs has moved towards the outlet of the interval at the end of the fourth year and the displacement is about 11 km, which results in the widening of the range of the maximum wax deposition. The phenomenon results from the nonequality of the wax deposition of Daqing crude in this hot pipeline and is the accumulation of the point-by-point movements of the maximum wax deposition during the fourth year. The mechanism lies in when the throughput and the inlet oil temperature of a pipeline interval keep constant, the wax deposition will make the oil temperatures along the line increase slowly. For a hot crude pipeline along which the wax deposition distributes unevenly, the oil temperature at the position of the maximum wax deposition occurring initially would increase and the temperature difference between the oil flow and the inner pipe wall there would reduce. The next downstream position will replace the initial position as a new maximum wax deposition position. Thus the location of maximum wax deposition will move downstream successively, which can bring a span of maximum wax deposition.
In the case of a long-term unpigged hot crude pipeline with reducing throughputs, whether the location of maximum wax deposition will move depends on two effects jointly: one is the reductions of the oil temperatures caused by lower throughput, another is the increments of the oil temperatures due to the wax deposition. Generally, at the beginning of the operation without pigging, the wax deposition is thinner and the heat loss caused by low throughput is stronger than the heat preservation of the wax deposition, which may result in the location of maximum wax deposition moving upstream. With the growth of the wax deposition, the insulation of wax deposition enhances and the heat loss rate arising from the reduction of throughput decreases due to the reduction of the inner pipe diameter, which will make the location of maximum wax deposition moving downstream. When the maximum wax deposition grows to a certain value (depending on pipeline sizes and operation conditions), the two effects balances and the location of maximum wax deposition will stop moving. Referring to Figure 8, it can be reasoned that the pipeline runs nearby the equilibrium point in the former three years under no pigging and over that point in the fourth year. It needs to be noted that the unpigged hot pipeline transporting Daqing crude began to transport a kind of light crude after no pigging of more than 4 years. Some maintenance and repairs were done for the pipeline to satisfy the new demand, so it was possible to measure the real wax deposition in situ. The following Figures 9 and 10 exhibit the maximum wax deposition and the wax deposition at the outlet of the first pipeline interval in situ, respectively, verifying the accuracy of the simulation.
5.2. Temperature Drop Characteristics of the Waxy Crude Oil Pipeline
Based on the above model, numerical simulations for the temperature drops along the pipeline with the shutdown duration of the unpigged hot crude pipeline carrying Daqing crude are completed. The results for the first pipeline interval in December of the 4th year of no pigging are shown in Figure 11.
To investigate the impact of the wax deposition formed through a long-term operation on the temperature drop of the hot crude pipeline, the oil temperatures at the outlet of the first pipeline interval after 20-hour shutdown in December of each year during four years under no pigging are selected from the simulation results. Also the wax depositions at the same position in the same time are selected from the wax deposition simulations. All the data are drawn in the same figure, see Figure 12.
From Table 1, the insulation effect of the increasing wax deposition is obvious in spite of the successive reduction of throughputs during no-pigging operation. When the throughputs decreased slowly and the inlet oil temperatures kept constant (in Decembers of the 1st and 2nd years), the increasing wax deposition ensured the same outlet oil temperatures. While the throughputs declined a lot and the inlet oil temperature was enhanced a few (in December of the 3rd year, especially in December of the 4th year), the wax deposition made the outlet oil temperatures relatively steady.
But from Figure 12, it can be concluded that the wax deposition formed during normal operation will not always be positive to keeping the oil temperature after the shutdown of the pipeline. The existence of a turning point corresponding to the 2nd year proves this issue. It can be explained as follows: on the one hand, the wax deposition slows down the heat transfer between oil and the surrounding soil and thus reduces the oil temperature drop rate during the shutdown. On the other hand, the reduction of the flow rate reduces the outlet oil temperature during the normal operation, leading to a lower oil temperature with the same shutdown duration. So the final impact of the wax deposition on the temperature drop of a hot crude pipeline after its shutdown depends on the two effects simultaneously: the insulation enhancement of the wax deposition and the flow rate reduction in the pipeline, each of which advances oppositely. Also, there is a balance point between the two effects. For the case in Figure 12, the balance point corresponds to the wax deposition of about 30 mm at the outlet of the first pipeline interval.
5.3. The Effect of Wax Deposition on the Restart Process of the Waxy Crude Oil Pipeline
Combining the hydrodynamic model with the thermodynamic model introduced above, numerical simulation for the restart process of the long-term unpigged hot crude pipeline is carried out. The results aimed at the first pipeline interval are selected for analyses, which show how the inlet pressures vary with the duration of restart for different shutdown durations in different years of no pigging (see Figures 13, 14, 15, 16, and 17).
Inspection of Figure 13 reveals that the safe maximum shutdown duration for the first pipeline interval in December of the first year of no pigging is 40 hours, at which the inlet pressure already reaches the MAOP of the pipeline, namely, 4.32 MPa. When the shutdown durations lie between 40 and 44 hours, the retaining times of maximum inlet pressure on MAOP delay with the shutdown time and are 0.5 hour, 1 hour, 2.3 hours, and 3.8 hours, respectively, as shown in Figure 14. The longer the duration of the maximum inlet pressure on MAOP is, the larger the risk of restart failure is. So, the safe maximum shutdown duration for the case in Figure 13 should be determined as 40 hours despite the fact that the pipeline interval may be restarted successfully above 40 hours but below 44 hours. Once the pipeline interval discussed was shut down over 44 hours under operation conditions in December of the first no-pigging year, it would not be restarted normally according to the operation regulations for the pipeline and emergency procedures should be activated.
From Figure 15, it is clear that the safe maximum shutdown duration for the first pipeline interval in December of the second year of no pigging is 42 hours. Between 42 and 45 hours of shutdown durations, the pipeline interval might be restarted but the risk of the restart failure is large. Over 45 hours, the pipeline interval could not be restarted normally.
From Figure 16, the safe maximum shutdown duration for the first pipeline interval in December of the third year of no pigging is 25 hours. Over 26 hours, the pipeline interval could not be restarted normally.
From Figure 17, the safe maximum shutdown duration for the first pipeline interval in December of the fourth year of no pigging is 18 hours. The pipeline interval could not be restarted normally when the shutdown delayed over 19 hours.
The restart safeties of the first pipeline interval during no pigging are derived from the above analysis and listed in Table 2.
In Table 2, SMSD is the abbreviation for safe maximum shutdown duration, IPRF for interval of possible restart failure, and CSDURN for critical shutdown duration unable to restart normally.
From Figure 13 to Figure 17, it could be found that the SMSDs of the first pipeline interval in Decembers during no pigging reduced with the growth of the wax deposition except in the second year. One of the causes is the final impact of wax deposition on the thermodynamics process during shutdown, which is discussed in Section 5.2. Another one arises from the fact that the growth of the wax layer reduces the effective inner diameter of the pipe resulting in the increasing of the restart pressure, even if the reduction of the throughput can counteract partly this tendency. Furthermore, the IPRF almost does not exist since the second year, which implies that there is no buffer between safe restart and the limit state of restart failure after the second year of no pigging, namely, the restart safety margin of the pipeline reduces.
In this paper, the characteristic of the wax deposition of a waxy crude oil pipeline is investigated. A new model for wax deposition is established and validated by the field data. Besides, the effects of wax deposition on shutdown and restart process of the crude oil pipeline are studied and the following conclusions are obtained.
The wax deposition in a long-term unpigged hot crude pipeline decreases the heat loss of the pipeline to the surroundings at any time and the heat accumulation in the soil as well. The insulation of the wax deposition is beneficial to the heat preservation of a hot crude pipeline during its normal running, especially under decreasing throughputs. However, the reduction of the flow rate reduces the outlet oil temperature during the normal operation, leading to a lower oil temperature with the same shutdown duration. Finally, the increasing insulation for the pipeline and the reducing flow rate caused by the growth of wax deposition jointly affect the safety of the shutdown and restart of such a special pipeline.
The study is financially supported by the National Science Foundation of China (no. 51134006).
- K. S. Pedersen and H. P. Rønningsen, “Influence of wax inhibitors on wax appearance temperature, pour point, and viscosity of waxy crude oils,” Energy and Fuels, vol. 17, no. 2, pp. 321–328, 2003.
- M. Del Carmen García and A. Urbina, “Effect of crude oil composition and blending on flowing properties,” Petroleum Science and Technology, vol. 21, no. 5-6, pp. 863–878, 2003.
- M. Del Carmen García, “Paraffin Deposition in Oil Production,” in SPE International Symposium on Oilfield Chemistry, pp. 161–167, SPE, Houston, Tex, USA, February 2001.
- A. Hammami, J. Ratulowski, and J. A. P. Coutinho, “Cloud points: can we measure or model them?” Petroleum Science and Technology, vol. 21, no. 3-4, pp. 345–358, 2003.
- P. Singh, R. Venkatesan, H. S. Fogler, and N. Nagarajan, “Formation and aging of incipient thin film wax-oil gels,” AIChE Journal, vol. 46, no. 5, pp. 1059–1074, 2000.
- R. Venkatesan and H. S. Fogler, “Comments on analogies for correlated heat and mass transfer in turbulent flow,” AIChE Journal, vol. 50, no. 7, pp. 1623–1626, 2004.
- H. S. Lee, Computational and rheological study of wax deposition and gelation in subsea pipelines [Ph.D. thesis], University of Michigan, Ann Arbor, Mich, USA, 2008.
- A. Matzain, Multiphase flow paraffin deposition modeling [Ph.D. dissertation], The University of Tulsa, Tulsa, Okla, USA, 1999.
- H. K. Mirazizi, Investigation of single-phase paraffin deposition characteristics under turbulent flow [M.S. thesis], The University of Tulsa, Tulsa, Okla, USA, 2011.
- M. Lashkarbolooki, A. Seyfaee, F. Esmaeilzadeh, and D. Mowla, “Experimental investigation of wax deposition in Kermanshah crude oil through a monitored flow loop apparatus,” Energy and Fuels, vol. 24, no. 2, pp. 1234–1241, 2010.
- R. Hoffmann and L. Amundsen, “Single-phase wax deposition experiments,” Energy and Fuels, vol. 24, no. 2, pp. 1069–1080, 2010.
- J. F. Tinsley, The effects of polymers and asphaltenes upon wax gelation and deposition [Ph.D. dissertation], Princeton University, Princeton, NJ, USA, 2008.
- H. O. Bidmus and A. K. Mehrotra, “Solids deposition during “cold flow” of wax-solvent mixtures in a flow-loop apparatus with heat transfer,” Energy and Fuels, vol. 23, no. 6, pp. 3184–3194, 2009.
- R. Tiwary and A. K. Mehrotra, “Deposition from wax-solvent mixtures under turbulent flow: effects of shear rate and time on deposit properties,” Energy and Fuels, vol. 23, no. 3, pp. 1299–1310, 2009.
- N. V. Bhat and A. K. Mehrotra, “Modeling the effect of shear stress on the composition and growth of the deposit layer from “waxy” mixtures under laminar flow in a pipeline,” Energy and Fuels, vol. 22, no. 5, pp. 3237–3248, 2008.
- J. Sestak, M. E. Charles, M. G. Cawkwell, and M. Houska, “Start-up of gelled crude oil pipelines,” Journal of Pipelines, vol. 6, no. 1, pp. 15–24, 1987.
- M. G. Cawkwell and M. E. Charles, “An improved model for start-up of pipelines containing gelled crude oil,” Journal of Pipelines, vol. 7, no. 1, pp. 41–52, 1987.
- H. Petter Rønningsen, “Rheological behaviour of gelled, waxy North Sea crude oils,” Journal of Petroleum Science and Engineering, vol. 7, no. 3-4, pp. 177–213, 1992.
- C. Chang, Q. D. Nguyen, and H. P. Rønningsen, “Isothermal start-up of pipeline transporting waxy crude oil,” Journal of Non-Newtonian Fluid Mechanics, vol. 87, no. 2-3, pp. 127–154, 1999.
- M. R. Davidson, Q. D. Nguyen, C. Chang, and H. P. Rønningsen, “A model for restart of a pipeline with compressible gelled waxy crude oil,” Journal of Non-Newtonian Fluid Mechanics, vol. 123, no. 2-3, pp. 269–280, 2004.
- G. Vinay, A. Wachs, and I. Frigaard, “Start-up transients and efficient computation of isothermal waxy crude oil flows,” Journal of Non-Newtonian Fluid Mechanics, vol. 143, no. 2-3, pp. 141–156, 2007.
- W. Li and J. J. Zhang, “Studies of temperature drop of buried waxy crude pipelines in shutdown,” Oil & Gas Storage and Transportation, vol. 23, pp. 4–8, 2004 (Chinese).
- X. H. Yang, Design and Management of Oil Pipelines, China University of Petroleum Press, Dongying, Shandong, China, 2006.
- B. Yu, C. Li, Z. Zhang et al., “Numerical simulation of a buried hot crude oil pipeline under normal operation,” Applied Thermal Engineering, vol. 30, no. 17-18, pp. 2670–2679, 2010.
- X. G. Cui and J. J. Zhang, “Determination of the thermal influence zone of buried hot oil pipeline on steady operation,” Journal of the University of Petroleum, vol. 28, no. 2, pp. 75–78, 2004.
- C. Xu, B. Yu, Z. Zhang, J. Zhang, J. Wei, and S. Sun, “Numerical simulation of a buried hot crude oil pipeline during shutdown,” Petroleum Science, vol. 7, no. 1, pp. 73–82, 2010.
- B. Yu and C. Xu, “Study on restart of crude oils batch pipelined with different outlet temperature,” Oil & Gas Storage and Transportation, vol. 28, pp. 4–16, 2009.