Wax Deposition PDF
Wax Deposition PDF
Wax Deposition PDF
Experimental Characterizations,
Theoretical Modeling, and Field Practices
Emerging Trends and Technologies in
Petroleum Engineering
Series Editor
Abhijit Y. Dandekar
PUBLISHED TITLES:
Wax Deposition: Experimental Characterizations, Theoretical Modeling, and
Field Practices, Zhenyu Huang, Sheng Zheng, H. Scott Fogler
Wax Deposition
Experimental Characterizations,
Theoretical Modeling, and Field Practices
Zhenyu Huang
Sheng Zheng
H. Scott Fogler
This book contains information obtained from authentic and highly regarded sources. Reasonable
efforts have been made to publish reliable data and information, but the author and publisher cannot
assume responsibility for the validity of all materials or the consequences of their use. The authors and
publishers have attempted to trace the copyright holders of all material reproduced in this publication
and apologize to copyright holders if permission to publish in this form has not been obtained. If any
copyright material has not been acknowledged please write and let us know so we may rectify in any
future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information stor-
age or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copy-
right.com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222
Rosewood Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that pro-
vides licenses and registration for a variety of users. For organizations that have been granted a photo-
copy license by the CCC, a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are
used only for identification and explanation without intent to infringe.
Visit the Taylor & Francis Web site at
http://www.taylorandfrancis.com
and the CRC Press Web site at
http://www.crcpress.com
Contents
Series Preface...........................................................................................................ix
Preface.................................................................................................................... xiii
Authors....................................................................................................................xv
1. Introduction......................................................................................................1
1.1 Background of Wax Deposition...........................................................1
1.2 Overview of Wax Testing, Modeling, and Management................. 5
v
vi Contents
6.4.3
Effect of the Operating Temperatures............................... 122
6.4.3.1 Analysis on the North Sea Condensate............. 122
6.4.3.2 Analysis on Another Oil...................................... 123
6.4.3.3 Importance of the Wax Precipitation Curve...... 125
6.4.3.4 Carbon Number Distribution of the Oil............ 126
6.5 Applying Wax Deposition Models to Investigate Deposit
Compositions...................................................................................... 127
6.6 Summary............................................................................................. 129
Bibliography......................................................................................................... 149
Appendix A: Nomenclature.............................................................................. 159
Index...................................................................................................................... 163
Series Preface
ix
x Series Preface
needing a different and unique approach from the drilling, production, and
reservoir engineering perspectives.
The petroleum engineering academic and industry community is also
aggressively pursuing nanotechnology with the hope of identifying innova-
tive solutions for problems faced in the overall process of oil and gas recovery.
In particular, a big spurt in this area in the last decade or so is evident from
the significant activities in terms of research publications, meetings, formation
of different consortia, workshops, and dedicated sessions in petroleum engi-
neering conferences. A simple literature search for a keyword nanotechnology
on http://www.onepetro.org, managed by the Society of Petroleum Engineers
(SPE), returns over 250 publications dating from 2001 onward with the bulk
of them in the last 5 or 6 years. Since 2008, SPE also organized three differ-
ent applied technology workshops specifically focused on nanotechnology in
the E&P industry. An Advanced Energy Consortium with sponsorships from
some major operators and service companies was also formed in 2007 with the
mission of facilitating research in micro and nanotechnology materials and
sensors having the potential to create a positive and disruptive change in the
recovery of petroleum and gas from new and existing reservoirs. Companies
such as Saudi Aramco have taken the lead in taking the first strides in evalu-
ating the potential of employing nanotechnology in the E&P industry. Their
trademarked ResbotsTM are designed for deployment with the injection flu-
ids for in situ reservoir sensing (temperature, pressure, and fluid type) and
intervention, eventually leading to more accurate reservoir characterization
once fully developed. Following successful laboratory core flood tests, they
conducted the industrys first field trial of reservoir nanoagents.7
The foregoing is clearly a statement of the new wave in the petroleum
engineering field, which is being created by emerging trends in unconven-
tional resources and new technologies. The publisher and its series editor are
fully aware of the rapidly evolving nature of these key areas and their long-
lasting influence on the current state and future of the petroleum industry.
The series is envisioned to have a very broad scope that includes but is not
limited to analytical, experimental, and numerical studies and methods and
field cases, thus delivering readers in both academia and industry an author-
itative information source of trends and technologies that have shaped and
will continue to impact the petroleum industry.
References
1. Retrieved from http://www.eia.gov/analysis/studies/worldshalegas/ (accessed
date June 10, 2013).
Series Preface xi
Abhijit Dandekar
University of Alaska Fairbanks
Preface
Wax deposition has become one of the most common flow assurance prob-
lems in the petroleum industry. As petroleum resources shift from onshore
reservoirs toward offshore subsea production, the industry is currently fac-
ing unprecedented challenges to maintain flow assurance for petroleum
production, in which the strategy to prevent or mitigate wax deposition has
become increasingly costly and complicated.
The goal to manage the issue of wax deposition involves answering the
following three questions:
They are typical questions to be answered not only for wax-related issues
but also for a variety of many other common production chemistry problems
in general flow assurance practices. The answer to the first question gener-
ally involves only a few fluid testing procedures that are relatively simple,
while much more understandings on production chemistry and fluid flow
are required to address the second question. The answer to the third ques-
tion requires not only knowledge that are shown in chemical engineering
textbooks but also significant operational experience, and the decision mak-
ers have to fully understand the production capability of the field and the
effectiveness as well as implication of any mitigation/remediation methods.
While there are several books that provide general knowledge on flow
assurance, a book that specifically addresses the issue of wax deposition is
still not yet available. This book is the first one that covers the entire spec-
trum of knowledge on wax deposition phenomena. It provides a detailed
description of the thermodynamic and transport theories for wax deposi-
tion modeling and a comprehensive review of the laboratory testing to help
establish appropriate control strategies for the field. It provides a progressive
introduction to help flow assurance engineers to understand the process of
wax deposition, to be familiar with the various methods to identify its sever-
ity, and to eventually control this problem. For engineering students, practic-
ing engineers, and researchers in the field of flow assurance, this book serves
as an in-depth discussion of how fundamental principles of thermodynam-
ics, heat, and mass transfer can be applied to solve a problem common to the
petroleum industry.
Going back to the three key questions that were raised earlier in this
Preface, we hope to provide valuable information in this book that could help
the readers to address these questions. Chapter 1 presents the background
xiii
xiv Preface
xv
xvi Authors
1
2 Wax Deposition
the world. A schematic of this shift is shown in Figure 1.1 (Huang, Senra,
Kapoor, & Fogler, 2011).
Taking the United States as an example, large offshore reservoirs that are
mainly by the coastlines of Louisiana, Texas, California, and Alaska have
quickly become one of the most crucial elements to the United States stra-
tegic development of energy resources (Economic analysis methodology for the
5-year OCS Oil and Gas Leasing Program for 20122017, 2011). While 20 million
bbl of oils were produced from offshore in the Gulf of Mexico in 1995, this
number has risen to 1400 million bbl in 2007 (Bai & Bai, 2012). The offshore
petroleum fluids are usually transported in long-distance pipelines, which
range from tens to hundreds of kilometers before they eventually reach
onshore processing facilities (Golczynski & Kempton, 2006). The oil typi-
cally comes out of the reservoir at a temperature around 160F and is cooled
significantly as it is transported through the pipes on the ocean floor, where
the water temperature is around 40F. This temperature difference between
the oil in the pipeline and the surrounding water on the ocean floor (160F
to 40F) becomes the driving force that causes the oil in the pipeline to cool
down. As the oil temperature decreases, the waxy components can precipi-
tate out of the oil and form deposits on the pipe wall. The problem of wax
deposition in the subsea pipeline has caused a series of problems for the flow
assurance industry, including increased pressure drop needed for oil trans-
portation and potential blockage of the pipeline. An example of a plugged
Onshore Offshore
FIGURE 1.1
A schematic of the change from onshore to offshore in petroleum production in the late twen-
tieth century. (From Huang, Z. et al., AIChE J. 57, 841851, 2011.)
Introduction 3
FIGURE 1.2
An example of a pipeline being plugged by wax deposits on the wall. (From Singh, P. et al.,
AIChE J., 46, 10591074, 2000.)
4 Wax Deposition
FIGURE 1.3
Areas reported to have wax deposition problems across the world, including the Gulf of
Mexico (Kleinhans et al., 2000), the North Slope (Ashford et al., 1990), the North Sea (Labes-
Carrier et al.,2002; Rnningsen, 2012), North Africa (Barry, 1971), Northeast Asia (Bokin et al., 2010;
Ding et al., 2006), Southern Asia (Agrawal et al., 1990; Suppiah et al., 2010), and South America
(Garcia, 2001).
25 $24.9
15
$12.4
10
5 $5.8
$2.9
$1.9 $1.0
0
7 14 30 60 90 180
Pigging frequency (days)
FIGURE 1.4
Cost of deferred production due to pigging estimated by Niesen (2002). (Niesen, V.: The real
cost of subsea pigging. E&P Mag. 2002. pp. 9798. Copyright Wiley-VCH Verlag GmbH & Co.
KGaA. Reproduced with permission.)
Wax precipitation
curve
Deposition
model input
Deposition Tuning Wax deposition
testing modeling
Deposition
model output
FIGURE 1.5
General overview of wax testing, modeling, and management procedures.
Wax appearance
temperature (WAT)
Temperature
FIGURE 1.6
Schematic of a wax precipitation curve.
process equipment to reveal the location where wax deposition can poten-
tially occur. This task can be completed using a wax deposition model that
usually contains a flow simulator to determine the hydraulic and heat trans-
fer characteristics in the pipe. The location where wax deposition starts to
occur is usually where the inner wall temperature of the pipe drops to below
the wax appearance temperature, as shown in Figure 1.7.
In addition to identifying the potential location of wax deposition in a
pipeline, a wax deposition model can evaluate the severity of wax deposition
by determining the growth rate of the wax deposits. To complete this task,
Introduction 7
Onset location of
wax deposition
FIGURE 1.7
Schematic of the inner wall temperature determined by a wax deposition model to indicate the
onset location for wax deposition.
2.1Introduction
Lets begin with the attempt to answer the first question of wax deposi-
tion, Do we have a problem for this oil field? While a definitive answer to
this question might require sophisticated fluid testing procedures, in many
industrial practices, one can often get a general idea by just looking at a few
parameters, for example, the wax content as well as the WAT. The values of
these parameters can often be found in the fluid testing conducted on the
stock tank oil sample of the production fluid. This testing is usually car-
ried out during the early stage of the flow assurance analysis for a field. An
example of the stock tank oil analysis is shown in Table 2.1.
9
10 Wax Deposition
TABLE 2.1
Example of Stock Tank Oil Analysis Results for Petroleum Fluid
Stock Tank Oil Analysis
Drilling fluid content <0.8 wt% STO
API gravity 32.1 at 60 F
Wax content 4.5 wt%
Asphaltene content 1.8 wt%
Sulfur content 0.9 wt%
Wax appearance temp. (CPM) 95 F
Pour point 30 F
Total acid number 0.07 mg KOH/g
Intuitively, we would think that an oil with more than 10 wt% wax con-
tent is more likely to cause wax deposition problem than another oil with
less than 2 wt% of wax, and that an oil with a WAT of 120F (49C) is likely
to be more difficult to manage than another oil with a WAT of 50F (10C).
Generally speaking, crude oils with more than 2% wax content and a WAT
higher than the ocean floor temperature (39F or 4C) could have wax depo-
sition risks. However, such thresholds only provide a very rough and often
overly conservative estimation. One can possibly find cases that an oil with a
WAT higher than 39F does not exhibit wax deposition during production. In
addition, two oils sharing similar wax content can have a great difference in
their deposition tendency. More importantly, one might find that two char-
acterizations on the same stock tank oil can show drastically different wax
contents or WAT values. In this case, how do we appropriately interpret the
test results and identify which one is more representative of the field?
In this chapter, we will go through a variety of testing methods in great
detail in order to introduce a comprehensive understanding on the labo-
ratory characterizations that would provide the data that are crucial to
reliably determine if a petroleum field is prone to the risk of wax depo-
sition. These characterizations focus on the two key thermodynamic
parameters: the WAT and the WPC. In Chapter 3, theoretical modeling of
wax deposition will be introduced as an alternative approach to study wax
precipitation.
Detect by:
First several Optical approaches;
precipitated wax Rheological approaches;
Controlled crystals Thermal approaches;
cooling Other novel approaches
Liquid oil
FIGURE 2.1
General strategy of WAT detection.
2.2.1Visual Inspection
Solid wax crystals scatter light, and thus the oilwax suspension appears
opaque. Using visual inspection, one can determine the WAT to be the tem-
perature at which cloudiness of the oil can first be observed due to wax
precipitation. Consequently, the WAT is also known as the cloud point.
Standard testing methods, such as ASTM-D2500, are developed based on
the visual inspection of wax solids by the naked eye (ASTM International,
2011). Figure 2.2 shows the experimental apparatus for the determination of
the WAT using ASTM-D2500.
Thermometer
Sample cell
Cooling bath
FIGURE 2.2
Apparatus for cloud point test using ASTM-D2500.
Experimental Characterization of Wax Thermodynamic Properties 13
A sample of crude oil is placed in a test jar and then submerged in a cool-
ing bath at a constant temperature. The decrease in the sample temperature is
monitored by a thermometer. After the temperature of the sample is decreased
by 1C, the test jar is removed from the cooling bath and inspected for signs of
solid precipitation. It should be noted that excessive temperature differences
between the cooling bath and the sample can lead to supercooling of the sample
and thus delay the formation of solid crystals (Bhat & Mehrotra, 2004; Coutinho
& Daridon, 2005). Therefore, it might appear to make sense to utilize computer-
programmed controllers to impose more carefully controlled slow cooling steps
to alleviate the effects of supercooling on wax crystal formation. However, even
with an infinitesimally slow cooling rate, the exact onset of wax precipitation is
still experimentally inaccessible due to the aforementioned detection limits of
the apparatus. Moreover, when the cooling rate approaches zero, measurement
of the WAT is unreasonably long and tedious. Practically speaking, when the
sample is cooled at a finite cooling rate ranging from 0.1C/min to 10C/min
(Hansen, Larsen, Pedersen, Nielsen, & Rnningsen, 1991; Rnningsen et al.,
1991). Within this practical range of cooling rates, supercooling is likely to occur,
and thus the measured WAT will vary with the cooling rate.
Regardless of the WAT detection methods used, wax precipitation must be
first induced by cooling the sample. Therefore, uncertainties in the WAT caused
by supercooling can be considered as a universal limitation in WAT detec-
tion techniques. In general, the measured WAT decreases as the cooling rate
increases because of the increasing extent of supercooling. Table 2.2 summa-
rizes the effects of the cooling rate on the WAT when different techniques are
applied.
It should be noted that the effect of the cooling rate on the WAT measured
by ASTM is not listed in Table 2.2 because the cooling rate in the ASTM is
not controlled.
The reliability of the WAT determined by ASTM-D2500 will vary for oil
with different waxy components. For example, for a model n-C12 (oil) and
n-C16 (wax) system that is transparent and colorless, the WAT measured by
ASTM-D2500 correlates well within ~1C difference of the values provided
by DSC and viscometry (Bhat & Mehrotra, 2004). Crude oils are usually dark
in color due to the light absorbance of the many other components such as
aromatics. Therefore, the determination of the first haziness during cooling
TABLE 2.2
Summary of the Effect of Cooling Rate on the Variation in Measured WAT
Technique Variation in Cooling Rates Decrease in WAT Source of Studies
CPM 0.1C/min to 0.5C/min ~3C/min Rnningsen et al. (1991)
DSC 1C/min to 3C/min <1C/min Alghanduri, Elgarni,
Daridon, and Coutinho
(2010)
Viscometry 0.03C/min to 0.2C/min <1.5C/min Rnningsen et al. (1991)
14 Wax Deposition
can be more difficult than in the case of model mixtures. The crude oil WAT
determined from ASTM-D2500 seems to be consistently ~5C lower than the
values determined by DSC and viscometry (Claudy, Letoffe, Neff, & Damin,
1986; Kruka, Cadena, & Long, 1995). Moreover, ASTM-D2500 is not applicable
when the color of the sample is darker than the ASTM-D1500 color of 3.5.
ASTM-D3117 was developed from the ASTM-D2500 based on the moti-
vation to facilitate visual observation of haziness in dark samples. In
ASTM-D3117, a smaller sample cell with an inner diameter (I.D.) of 0.8
in. is used in place of the ~1.2 in. I.D. sample cell used in ASTM-D2500.
Using a thinner sample cell, the thickness of the liquid sample film is
reduced, and the crude oil will appear less dark and thus allow bet-
ter observation of the first haziness. However, the improvement from
ASTM-D2500 to ASTM-D3117 might still not be sufficient to determine
the WAT of dark crude samples. Ijeomah et al. (2008) attempted to use
ASTM-D3117 to measure the WAT of a synthetic crude and an Alaska
North Slope (ANS) pipeline crude. In the study by Ijeomah et al., it was
reported that ASTM-D3117 is only applicable to the GTL sample because
it is transparent and light in color, while the determination of the WAT of
the ANS crude is impossible due to its dark color.
As is presented in this section, ASTM-D2500 and ASTM-D3117 provide
two rapid approaches to determine the WAT. However, the applicability of
these two methods in the determination of the crude oil WAT is greatly lim-
ited due to the dark color of the crude oil samples. In addition, the uncon-
trolled cooling in these two testing methods is likely to cause delay in wax
precipitation and underestimation of the WAT.
2.2.2CPM Techniques
CPM has the following advantages over ASTM-D2500 and ASTM-D3117:
The procedure to detect the WAT using CPM techniques is similar to that
of visual inspection: A drop of liquid oil sample is placed on a glass slide
and covered with a top glass to form a ~50 m liquid film. The liquid film is
then cooled at a cooling rate of usually 0.5 to 1C/min, and the first sign of
crystal formation is monitored by the microscope. Using a microscope, the
object can be magnified by ~100-fold, and one can observe solids as small
as ~0.5 to 1 m. In addition, the contrast between the wax solid and the liq-
uid is enhanced using cross-polarized light instead of unpolarized light.
Under a cross-polarized microscope, the liquid oil appears black because
Experimental Characterization of Wax Thermodynamic Properties 15
Vertical Horizontal
Plane of
polarizing polarizing
polarization
filter filter
rotated
FIGURE 2.3
Interactions between the cross-polarized light and (a) homogeneous liquid and (b) crystalline
material.
16 Wax Deposition
FIGURE 2.4
Photomicrograph of wax crystals under cross-polarized microscope; wax crystals appear as
white flocs. (From Rnningsen, H. P. et al., Energy & Fuels, 5, 895908, 1991.)
suspended in the liquid bulk in order for the change in viscosity to be detect-
able (Rnningsen et al., 1991).
It should be noted that in ASTM-D2500, ASTM-D3117, and the CPM
method, detection of the WAT is often dependent on the operators judg-
ment call of observing the first cloudiness or observing the first crys-
tals. Consequently, the WAT determined by different observers can vary.
Rnningsen et al. (1991) pointed out that the measurements made by one
operator normally are repeatable within ~1C, whereas the results obtained
by different operators agree within 1.5C2C. In order to reduce the subjec-
tivity in the measurements, Kok et al. (1996) utilized a photomonitor to track
the variation in the intensity of transmitted light through the sample, reduc-
ing the subjectivity in the WAT measurements due to different operators.
The CPM technique for the determination of the WAT overcomes the major
drawbacks in ASTM-D2500 and ASTM-D3117 and has become the most
applied technique in the determination of the WAT of crude oils. However, it
should be noted that both ASTM-D2500/D3117 and CPM methods are based
on the operators visual detection of the first cloudiness. Therefore, measure-
ments by different operators can vary, reducing the reliability of these two
TABLE 2.3
Summary of Experimental Studies on the Determination of WAT Using
CPMTechniques
Studies Sample Origin Number of Samples
Rnningsen et al. (1991) North Sea 17
Erickson, Niesen, and Brown (1993) Gulf of Mexico 15
Monger-McClure, Tackett, and Gulf of Mexico, Trinidad, 13
Merrill (1999) Michigan, and Oklahoma
Cazaux, Barre, and Brucy (1998) Undisclosed 2
Paso et al. (2009) Undisclosed 2
Experimental Characterization of Wax Thermodynamic Properties 17
FIGURE 2.5
CH2rocking vibrational modes in n-alkane chains.
A720
FIGURE 2.6
IR absorption spectrum of long-chain hydrocarbons with a characteristic absorbance at =
720 cm1.
18 Wax Deposition
A730
A722
A
FIGURE 2.7
IR spectrum of crystalline long-chain hydrocarbons with the original 720 cm1 absorption split
into two bands at 730 and 722 cm1.
A730
A722 Shaded area
represents I
A
FIGURE 2.8
Definition of the intensity of IR absorbance as a quantitative variable to correlate to the amount
of crystalline wax in oil.
Experimental Characterization of Wax Thermodynamic Properties 19
Intensity (I )
WAT
Temperature (T )
FIGURE 2.9
A drastic increase in the intensity as wax crystals form.
TABLE 2.4
Summary of Experimental Studies on the Determination of WAT Using
FT-IRTechniques
Literature Sample Origin Number of Samples
Roehner and Hanson (2001) Gulf of Mexico, Utah-Grand County, 3
and Alaska North Slope crude oil
Monger-McClure et al. Gulf of Mexico, Trinidad, Michigan, 13
(1999) and Oklahoma
Alcazar-Vara and Gulf of Mexico 3
Buenrostro-Gonzalez (2011)
mixture and crude oils that are dark in color. The sensitivity of the detec-
tion of the WAT by FT-IR can be comparable to the CPM method. However,
it is sometimes difficult to determine the inflection point in the intensity
temperature curve and thus difficult to determine the WAT. The uncertain-
ties in the determination of the inflection point prevail when the crude oil
of interest has a relatively wide n-alkane distribution or a low total n-alkane
content (Coutinho & Daridon, 2005).
2.2.4Viscometry
Above the WAT, the crude oil behaves as a Newtonian fluid, and the viscos-
ity of oil as a function of temperature can be described using the following
Arrhenius-type equation:
= AeEa/RT (2.2)
Yan & Luo, 1987). Determination of the WAT based on viscometry utilizes
the change in the slope in the viscosity vs. temperature curve. In a typical
experimental measurement of viscosity using a rheometer, the oil sample is
cooled at a rate of 0.03C2C/min (Rnningsen et al., 1991). The change of
viscosity as a function of temperature is measured under a constant shear
rate between 30 and 300 s1 (Rnningsen et al., 1991). In terms of reproduc-
ibility, the viscosity in the Newtonian region may have an error bar within
0.5%1%, while the viscosity in the non-Newtonian region can achieve an
error bar within 10% (Rnningsen et al., 1991).
In the viscositytemperature curve of a fluid, extrapolation from the
Newtonian region is performed to determine the WAT, as shown in Figure 2.10.
The fitting can be done as an exponential fitting of as a function T using a
computer, shown in Figure 2.10a or a semilog paper with ln() as a function of
1/T, shown in Figure 2.10b. Table 2.5 summarizes some of the representative
studies in which the WAT of the crude oils is determined based on viscosity
measurements.
It should be noted that when the temperature is further lowered below the
WAT, a substantial amount of wax can precipitate and remain suspended in
ln()
Actual viscosity
Extrapolated based on
Arrhenius equation
FIGURE 2.10
(a) Plot of a representative crude oil viscosity as a function of temperature with the dashed line
showing the extrapolation from the Newtonian region of the curve. (b) Plot of ln() as a func-
tion of 1/T with the dashed line showing extrapolation from the Newtonian region.
TABLE 2.5
Several Representative Studies Reporting the Measurement of WAT Based on
Viscosity Measurements
Literature Sample Origin Number of Samples
Rnningsen et al. (1991) North Sea 17
Kruka et al. (1995) Middle East 1
Kok et al. (1996) Undisclosed 15
Cazaux et al. (1998) Undisclosed 2
Kok et al. (1999) Undisclosed 8
Elsharkawy, Al-Sahhaf, and Fahim (2000) Middle East 8
Alcazar-Vara and Buenrostro-Gonzalez (2011) Gulf of Mexico 3
Erickson et al. (1993) Gulf of Mexico 15
Experimental Characterization of Wax Thermodynamic Properties 21
ln()
Non-Newtonian
.
. 1 Newtonian
2
.
3
T
Abnormal point WAT
FIGURE 2.11
Change of crude oil viscosity during cooling processes. (From Li, H., & Zhang, J. Fuel, 82,
13871397, 2003.)
the crude oil, causing the shear-thinning characteristics of the crude oil, as
shown in Figure 2.11.
In this section, the determination of the WAT based on viscometry was dis-
cussed. The sensitivity of the viscometry technique is lower than CPM and
FT-IR. In most cases, an amount of 0.30.4 wt% solid material needs to have
precipitated in order to cause a detectable deviation in the viscosity, whereas
the wax crystals only need to grow to ~1 m to be detectable by CPM, cor-
responding to ~0.1 wt% solid material in the liquid (Rnningsen et al., 1991).
Therefore, the WAT determined based on viscosity measurements can be
0C10C lower than the values determined by CPM (Erickson et al., 1993;
Rnningsen et al., 1991). Depending on the composition of the precipitating
n-alkane, the change in the slope of the viscositytemperature curve near the
WAT can be either sharp or gradual. If the precipitating n-alkane only covers a
narrow range of carbon numbers, the change in the slope is likely to be sharp
because n-alkanes precipitate in a narrow temperature range. Contrarily, if
the precipitating n-alkane covers a wide range of carbon numbers, the change
in the slope is likely to be gradual as the n-alkanes precipitate over a wider
temperature range. Sometimes, the change in the slope can become difficult
to identify when very little wax is precipitating, which will cause a significant
underestimation of the WAT (Kok et al., 1996; Kruka et al., 1995).
2.2.5Thermal TechniquesDSC
Wax precipitation is an exothermic process. Thermal techniques, e.g.,DSC,
capture the heat of crystallization released by the sample upon wax precipi-
tation and thus detect the onset of wax precipitation. It should be noted that
due to the thermal resistance of the sample and thus the finite rate of heat
transfer, the actual temperature of the sample does not respond instanta-
neously to changes in the temperature. When the sample is cooled at a rate
22 Wax Deposition
Baseline
Heat flow
Inflection
point
T WAT
FIGURE 2.12
The definition of baseline and the WAT based on a typical DSC thermogram.
Experimental Characterization of Wax Thermodynamic Properties 23
10
0
140 40 60
Temperature (C)
FIGURE 2.13
DSC thermogram of North Sea crude oil with narrow temperature range of wax precipitation.
(From Hansen, A. B. et al., Energy & Fuels, 5, 914923, 1991.)
WAT for this method (Standard Test Method for Measurement of Transition
Temperatures of Petroleum Waxes by DSC, 2013). It should be noted that the
shape of the thermogram greatly depends on the composition of the wax.
A sharp peak on the thermogram can be observed when the precipitation
of wax occurs in a narrow temperature range, as can be seen from Figure
2.13. The heat release peak can be gradual if the wax precipitation occurs in
a wide temperature range. Due to the complex composition of the wax in
the crude oil, its DSC thermogram can sometimes be irregular, as shown
inFigure 2.14 (Hansen et al., 1991).
10
Heat flow (mcal/sec)
5
WAT = 31.5C
0
140 40 60
Temperature (C)
FIGURE 2.14
DSC thermogram of North Sea crude oil with irregular shape. (From Hansen, A. B. et al.,
Energy & Fuels, 5, 914923, 1991.)
24 Wax Deposition
TABLE 2.6
Summary of Studies Reporting Measurement of WAT Using DSC
Literature Sample Origin Number of Samples
Hansen et al. (1991) North Sea 17
Kok et al. (1996) Undisclosed 15
Monger-McClure et al. (1999) Gulf of Mexico 13
Kok et al. (1999) Undisclosed 8
Elsharkawy et al. (2000) Middle East 8
Alghanduri et al. (2010) Libya 5
Alcazar-Vara and Buenrostro-Gonzalez (2011) Gulf of Mexico 3
de Oliveira et al. (2012) Brazil 2
Experimental Characterization of Wax Thermodynamic Properties 25
(c)
(d)
4 0 4 8 12 16 20 24 28 32 C
FIGURE 2.15
DSC thermogram obtained by cooling a crude oil sample at (a) 10C/min, (b) 7C/min, (c) 4C/
min, and (d) 1C/min. (From Jiang, Z. et al., Fuel, 80, 367371, 2001.)
The DSC method has two unique advantages over other detection meth-
ods. Firstly, DSC analysis can identify different temperature regions where
wax precipitates (Kruka et al., 1995). Hammami and Mehrotra (1995) have
demonstrated that using DSC, one can observe the precipitation of differ-
ent n-alkanes occurring at different temperature ranges. Figure 2.16 shows a
DSC thermogram of the C50 + C44 binary mixture.
Exothermal (arbitrary units)
C50 fraction
1.00
0.75
0.50
0.25
0.00
60 80 100 120
Temperature (C)
FIGURE 2.16
DSC thermogram of the C50 + C44 binary mixture at 1C/min cooling rate. The C50 fraction in
the binary mixture varies from 1.00 to 0.00. (From Hammami, A., & Mehrotra, A. K., Fluid Phase
Equilib., 111, 253272, 1995.)
26 Wax Deposition
It is seen from Figure 2.16 that two distinct peaks can be identified from
the thermogram. The peak in the higher temperature region is associated
with the crystallization of C50, and the other peak in the low temperature
region is due to the crystallization of C44.
Secondly, DSC can be adapted to determine the WAT of live crude oils that
often contain a significant amount of light ends. The light ends in live crude
often have a bubble point below atmospheric pressure. The effect of the light
end composition on the crude oil WAT can be profound (Juyal, Cao, Yen, &
Venkatesan, 2011). In this case, the WAT increases as the amount of dissolved
light end decreases. In the case of a North Sea gas condensate, an increase as
large as 10C in the WAT is observed due to depressurization and the loss of
light end (Daridon, Coutinho, & Montel, 2001). During a DSC test, one can
maintain the dissolved gas in live crude by regulating the sample pressure
over a range of 11000 bars using a specific high-pressure testing cell. The
WAT of live crudes can thus be measured (Juyal et al., 2011; Vieira, Buchuid,
& Lucas, 2010).
In this section, the determination of the WAT by DSC was discussed. The
sensitivity of DSC in determining the WAT is comparable to that of viscom-
etry but is less than that of CPM and FT-IR. The DSC method has two advan-
tages over other techniques: it can detect different temperature ranges of
wax precipitation and can be adapted to measure the WAT of live crude oils
at elevated pressure.
TABLE 2.7
Four Representative Reviews Comparing the WAT Values Obtained from
Different Detection Techniques
Literature CPM DSC Viscometry FT-IR
Erickson et al. (1993)
Kruka et al. (1995)
Monger-McClure et al. (1999)
Coutinho and Daridon (2005)
Experimental Characterization of Wax Thermodynamic Properties 27
80
WATDSC > WATCPM
(16 examples)
60
WATDSC (C)
40
20
WATDSC < WATCPM
(35 examples)
0
0 20 40 60 80
(a) WATCPM (C)
80
WATviscometry > WATCPM
(8 examples)
60
WATviscometry (C)
40
20
WATviscometry < WATCPM
(35 examples)
0
0 20 40 60 80
(b) WATCPM (C)
FIGURE 2.17
Comparison between the WAT measured by different methods: (a) DSC and CPM; (b) viscom-
etry and CPM.
It is seen that, statistically, the WAT determined by CPM is higher than the
WAT determined either by DSC or by viscometry, indicating that CPM provides
the most conservative value (the highest estimation) for the WAT in most cases.
Different from the DSC and viscometry measurements, which generally
provide lower WATs than CPM, FT-IR provides WAT values that are com-
parable with those measured by CPM, as shown in Figure 2.18 (Monger-
McClure et al., 1999).
Meanwhile, the viscometry and DSC method usually provide comparable
WAT values, as can be seen in Figure 2.19. The number of cases where the
WAT (DSC) is greater than the WAT (viscometry) is similar to the number of
cases where the WAT (DSC) is smaller than the WAT (viscometry), indicating
that both techniques often share similar sensitivity.
In summary, for the WATs determined by the four techniques, one might
wonder which is the correct WAT. Actually, none of these four WAT values
28 Wax Deposition
80
60
WATFT-IR (C)
40
20
0
0 20 40 60 80
WATCPM (C)
FIGURE 2.18
Comparison between the WAT obtained from FT-IR and the WAT obtained from CPM.
20
16
Number of examples
12
0
<10C 10C 5C 5C0C 0C5C 5C10C >10C
TDSCTviscometry
FIGURE 2.19
Comparisons between the WAT measured by the viscometry method and the DSC method.
is correct because none of them represents the exact onset of wax precipita-
tion. The WATs provided by CPM and FT-IR are the closest to the exact onset
at the thermodynamic WAT, whereas the WATs provided by DSC and viscos-
ity are usually further below the thermodynamic WAT. The WATs by DSC
and viscometry thus provide an estimation of temperature below which
severe wax precipitation could occur.
Experimental Characterization of Wax Thermodynamic Properties 29
T = T0 t + AT sin(t)
dT (2.3)
= + AT cos(t)
dt
The term represents the constant cooling rate, and the term ATcos(t)
represents the additional oscillation in the cooling rate imposed by TMDSC.
In theory, the additional oscillation in the cooling rate enhances the rate of
crystallization near the WAT. Unfortunately, the results reported by Jiang et
al. have not yet demonstrated a significant improvement over conventional
DSC (Coutinho & Daridon, 2005).
Other techniques to measure the WAT include NIR scattering, x-ray CT,
densitometry, and filter plugging. Details regarding the principles and detec-
tion limits of these techniques are summarized in Table 2.8. Sample studies
applying these techniques for WAT detection are also provided in Table 2.8.
TABLE 2.8
Summary of Principles, Detection Limits, and Sample Studies of NIR Scattering,
X-Ray CT, Densitometry, and Filter-Plugging for WAT Detection
Technique Principle Detection Limit Sample Study
NIR scattering Solid wax crystals can Theoretically more Paso et al. (2009)
scatter NIR light. sensitive than CPM and
can detect wax crystals
with size ~55nm
X-ray CT A solid wax crystal has a Comparable to ASTM Karacan, Demiral,
higher density than the but can detect the WAT and Kok (2000)
liquid bulk. for dark crude oils,
which cannot be
analyzed using ASTM
Densitometry Comparable to DSC and Alcazar-Vara and
viscometry Buenrostro-
Gonzalez (2011)
Filter-plugging Solid wax crystals larger Comparable to CPM Monger-McClure
than ~5 m in size clog et al. (1999)
a filter with pore size
~5m, resulting in a
detectable increase in
the differential pressure.
30 Wax Deposition
qWATT
wwax , WATT = (2.4)
H crystallization
WAT
Temperature (K)
FIGURE 2.20
A representative example of the WPC of a crude oil.
Experimental Characterization of Wax Thermodynamic Properties 31
TABLE 2.9
Summary of Different Experimental Techniques to Characterize the Wax
Precipitation Curve
Property Directly Relation to the Amount
Technique Measured of Precipitated Solid Wax Examples
DSC Heat flow due to The heat flow released Hansen et al. (1991);
precipitation during precipitation is Elsharkawy et al. (2000);
directly proportional to Martos et al. (2008, 2010);
the amount of wax Coto, Martos, Espada,
precipitated. Robustillo, Pea et al.
(2011); Coto, Martos,
Espada, Robustillo,
Merino-Garca et al. (2011)
Nuclear Resonance peak The peak intensity is Pedersen, Hansen, Larsen,
magnetic intensity directly proportional to Nielsen, and Rnningsen
resonance the number of hydrogen (1991)
(NMR) atoms in the solid phase.
FT-IR Intensity of FT-IR The intensity of the Alcazar-Vara and
absorbance characteristic FT-IR Buenrostro-Gonzalez
absorbance is directly (2011); Roehner and
proportional to the Hanson (2001)
amount of solid wax.
Baseline
T < WAT
WAT
T
FIGURE 2.21
Representative DSC thermogram of waxy crude oil. The cumulative heat released due to wax
precipitation can be determined by integrating the area enclosed by the baseline and the
thermogram.
32 Wax Deposition
Quick decay of
solid phase A long tail due
NMR intensity
Time
FIGURE 2.22
Representative NMR FID of a mixture with solid wax suspended in liquid oil.
34 Wax Deposition
t1: I ( s, t1 ) + I (l, t1 ) = I (T , t1 )
t2 t (2.5)
t2: I ( s, t1 )exp + I (l, t1 )exp 2 = I (T , t2 )
s l
As the intensities of the NMR signals, I(s, t1), I(l, t1), I(s, t2), and I(l, t2) are pro-
portional to the amount of solid and liquid, one can determine the relative
amount of solid and liquid by solving Equation 2.5 simultaneously.
Pedersen, Hansen et al. (1991) used the NMR technique to measure the
WPC of 17 crude oils. Detailed procedures to perform this NMR measure-
ment can be found in their published study.
In this section, the measurement of WPC by the NMR is discussed. Using
the NMR technique, measuring the solid content at a particular temperature
requires ~10 s. NMR measurements are expected to have a ~0.2% absolute
standard deviation in the measured solid wt% of replicates of samples.
14.5
Directly proportional to the
Intensity (cm1) 14.0 amount of solid
12.0
0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0
Temperature (C)
FIGURE 2.23
Intensity of FT-IR absorbance of a Gulf of Mexico crude oil as a function of temperature. (From
Roehner, R. M., & Hanson, F. V., Energy & Fuels, 15, 756763, 2001.)
extrapolation of the intensity from above the WAT to below the WAT, as
denoted by a dashed line, provides an estimation of the liquid-phase con-
tribution at temperatures below the WAT. The difference between the mea-
sured intensity and the extrapolation is directly proportional to the amount
of solid present in the mixture. This difference can thus be used to back-
calculate the amount of solid precipitated. Alcazar-Vara and Buenrostro-
Gonzalez (2011) and Roehner and Hanson (2001) provide two representative
studies regarding the application of FT-IR on WPC measurements.
Three techniques for the determination of the WPC, NMR, FT-IR, and DSC
have been discussed so far. At this point, the readers might be interested
in seeing a comparison between the WPC determined from different tech-
niques. Unfortunately, studies reporting the application of these techniques
are conducted individually with different crude oils, and there does not exist
a comprehensive study using a single crude oil in order to make comparison
of the WPC measured by these three techniques.
because of the following reason: Precipitation of wax molecules causes wax crys-
tals to grow into an interlocking network that can entrap liquid oil. Therefore, the
weight of the solid cake actually represents the total weight of the true
precipitated solid and the entrapped liquid.
Being aware of the issue that liquid oil can be trapped in the solid cake,
one can improve the aforementioned procedures by either performing fur-
ther removal of the entrapped oil or additional characterization of the pre-
cipitated solid phase to quantify its true solid content. In combination with
proper removal of the entrapped oil or characterization of the solid cake, the
proposed procedure based on separation of the solid and liquid phases can
be used to obtain WPCs with acceptable accuracy. Table 2.11 summarizes
representative studies where the oil precipitation curve is obtained based on
separation of the solid cake from the liquid bulk.
TABLE 2.11
Summary of Representative Studies on WPC Characterization Using Separation-
Based Techniques
Further Removal
Separation of Entrapped Characterization of the
Literature Techniques Liquid by Solid Phase
Burger, Perkins, and Centrifugation Washing by Not performed
Striegler (1981) toluene
Coto, Martos, Pea, Filtration Washing by Not performed
Espada, and acetone
Robustillo (2008)
Martos et al. (2008) Filtration Washing by Characterize the composition
acetone by HTGC, and calculate
C15+ content to be the solid
content
Roehner and Hanson Centrifugation Not performed Characterize the composition
(2001) of the solid by HTGC, and
define the precipitated wax
as the components with an
n-alkane-to-nonalkane ratio
in excess of the same
component in oil
Dauphin, Daridon, Filtration Not performed Characterize the composition
Coutinho, Baylre, and of the solid by HTGC, and
Potin-Gautier (1999) apply a mass balance to
Pauly, Dauphin, and calculate the true solid
Daridon (1998) content
Pauly, Daridon, and
Coutinho (2004)
Han, Huang, Senra, Centrifugation Not performed Characterize the composition
Hoffmann, and of the solid by HTGC, and
Fogler (2010) apply a mass balance to
calculate the true solid
content
Experimental Characterization of Wax Thermodynamic Properties 37
Attempts have been made to remove the entrapped oil by washing the
solid phase with toluene or acetone (Burger et al., 1981; Coto et al., 2008).
However, one cannot guarantee that washing by toluene or acetone com-
pletely removes the entrapped oil. Therefore, the WPC obtained after wash-
ing might still be an overestimation compared to the true WPC. Figure
2.24 shows a comparison of the WPC measured using fractional precipita-
tion followed by acetone wash and the WPC of the sample crude measured
by DSC.
It can be seen that the amount of solid precipitation measured by the frac-
tional precipitation method is significantly higher than the values from the
DSC method, indicating that acetone wash cannot completely remove the
entrapped oil.
In order to provide a more accurate WPC, Martos et al. (2008) further char-
acterized the composition of the solid cake using HTGC after washing with
acetone. n-Paraffins with carbon numbers less than 15 have melting points
lower than the room temperature. As a result, C15 paraffins usually do not
precipitate nor contribute to wax deposition. Martos et al. calculated the C15+
content of the solid cake based on the HTGC and used the calculated C15+
content to represent the true solid content of the solid cake. Martos et al.s
approach assumes that an n-paraffin with a carbon number less than 15 only
exists in the entrapped liquid phase, while C15+ paraffins exist in the solid
phase. Assuming that C15 is the lightest paraffin in the solid phase is purely
empirical. Roehner and Hanson (2001) proposed a more rigorous method
of determining the lightest n-paraffin component in the solid phase. They
determined whether or not an n-paraffin component is in the solid phase
by comparing its mass fraction in the solid cake with its mass fraction in the
original oil. The precipitating paraffin components are enriched in the solid
cake and thus have a higher mass fraction in the solid cake than in the crude
12.0
10.0
Filtration followed
8.0 by acetone wash
Weight % wax
6.0
DSC
4.0
2.0
0.0
250 260 270 280 290 300
Temperature (K)
FIGURE 2.24
Comparison between the WPC determined by DSC and filtration followed by acetone wash.
(From Martos, C. et al., Energy & Fuels, 22, 708714, 2008.)
38 Wax Deposition
oil. Therefore, Roehner and Hanson count the amount of Ci+ as the true
solid content, where i is the lowest carbon number component whose mass
fraction in the cake exceeds its mass fraction in the original oil. Both Martos
et al. and Roehner and Hanson assume that a particular carbon number
component exists only in the solid or liquid phase. Light components in the
entrapped oil tend to exist in the liquid phase, whereas heavy components
tend to exist in the solid phase. However, even heavy n-paraffins in the solid
cake can exist in both liquid and solid phase: Heavy n-paraffins in the solid
cake can remain dissolved in the entrapped oil. The heavy paraffins dis-
solved in the entrapped oil should not be counted as contributing to the solid
phase. Therefore, the methods of Martos et al. and Roehner and Hanson
might still overestimate the WPC. In order to account for the amount of dis-
solved heavy paraffin in the entrapped liquid phase, Dauphin et al. (1999),
Pauly et al. (1998, 2004), and Han et al. (2010) proposed a more mathemati-
cally involved method where they solved a set of mass balance equations for
the solid content of the solid cake.
In this section, separation-based techniques for the determination of the
WPC are discussed. Without proper characterization of the solid content of
the cake, these techniques tend to give overestimations of the WPC due to
entrapment of oil in the solid cake. With proper characterization of the cake
solid content, separation-based methods provide WPCs that are comparable
with the precipitation curves measured by DSC, FT-IR, and NMR.
TABLE 2.12
Summary of Experimental Techniques for the Characterization of the Composition
and Structure of Wax Deposit
Expected
Experimental
Outcome in the Case
of Deposit Being Rich
Analytical in Macrocrystalline Representative
Methods Measurable Wax Studies
DSC Heat of crystallization High heat of fusion Alghanduri et al.
and melting point and narrow melting (2010)
temperature range
XRD Diffraction peaks Can observe sharp Alghanduri et al.
diffraction peaks (2010)
1 HNMR Numbers of methyl A low methyl-to- Musser, Kilpatrick,
groups (CH3) and methylene ratio due and Carolina (1998)
methylene groups to fewer branches
(CH2)
13 CNMR Numbers of primary ~70% of the carbon Musser et al. (1998);
carbon (CH3), atoms exist in the Alghanduri et al.
secondary carbon form of primary (2010)
(CH2), and tertiary (CH3) and secondary
carbon (CH) carbon (CH2). Low
Number of carbon content of tertiary
atoms on aromatic carbon and carbon on
rings and cyclic rings rings
FT-IR Intensity of absorption >60% straight Musser et al. (1998)
at ~720cm1 methylene
associated with long
straight-chain
methylene and the
content of straight
methylene
Elemental analysis The ratio of the High H/C at ~2 Musser et al. (1998);
number of hydrogen Alghanduri et al.
atoms to the number (2010)
of carbon atoms
(H/C)
Mass spectroscopy Mass distribution of MS peaks cover a Musser et al. (1998)
wax range of 350600 a.u.,
and two neighboring
peaks are 14 a.u. (the
molar mass of one
methylene unit) apart
2.5Summary
In this chapter, experimental techniques to characterize the following key
wax precipitation information were discussed:
The WAT
The WPC
The WAT provides information regarding the onset location of wax depo-
sition in a transport pipeline. The WPC characterizes the solubility limit of
wax in crude oil and thus affects the concentration driving force for wax
deposition. The WAT and the WPC are both extremely important input
parameters required for wax deposition modeling. These experimental char-
acterizations are also commonly used as references for benchmarking wax
thermodynamic models. The characterization of the composition and of the
molecular structure of the precipitated wax was also briefly introduced, as
this characterization provides a piece of information regarding the mechani-
cal properties of the deposit.
It should be noted that the characterization of wax precipitation informa-
tion is not always possible and becomes increasingly difficult and sometimes
virtually impossible in actual field operations. When experimental charac-
terization of wax precipitation is difficult, modeling strategies are usually
applied to predict the wax precipitation characteristics. In Chapter 3, thermo-
dynamic modeling of the wax precipitation characteristics will be discussed.
3
Thermodynamic Modeling
of Wax Precipitation
3.1Introduction
In Chapter 2, we discussed different experimental approaches to character-
ize wax precipitation. The information of wax precipitation is not only indis-
pensable to identify if wax deposition would be a problem for a field but also
important to quantify how severe the problem could be.
However, experiments are known to have their own disadvantages. First,
cost or time can be an issue. For example, the characterization of wax pre-
cipitation of a live crude oil that is rich in light ends would need to be per-
formed in pressurized cells, which are often not readily accessible and the
experiments are expensive to perform. More importantly, reliability of the
results can be called into question. For example, for crude oils with low wax
content, the uncertainties in the measured WAT and WPC become large and
the results can be misleading. Although not quite frequently, sample mis-
handling or inappropriate design of experimental conditions can, in fact,
occur. Therefore, if we only rely on laboratory measurements, we might find
ourselves facing two results from similar testing methods showing quite dif-
ferent values. In this case, how would we gauge the validity of each result
without spending significant time and cost on a third test, which could also
be prone to uncertainties and errors? Sometimes experimental methods are
just not possible with todays technology: how can we conveniently take a
sample of wax deposit from an offshore subsea pipeline in the field to deter-
mine its composition, which could be an important piece of information for
one to determine the yield stress of the deposit and thus to optimize the pig-
ging frequency?
Consequently, we need theoretical modeling as an alternative to provide
us answers that are not prone to the errors and uncertainties associated with
experimental measurements. The results of experimental data, when bench-
marked with theoretical modeling, often provide the most reliable answers
for one to evaluate wax deposition concern for a field. In this chapter, we will
41
42 Wax Deposition
Liquid
Wax + liquid
Pressure
Liquid + vapor
Wax + liquid + vapor
Temperature
FIGURE 3.1
Typical phase diagram of a waxy crude oil. (From Leontaritis, K. J. Fuel Sci. Technol. Intl., 14,
1339, 1996.)
Among the 3N + 3 variables, the amount of the solid phase, nS, and the solid
phase composition, (s1, s2, , si, sn), are of the greatest interest in regard to
wax precipitation characterization. Compared to the live crude oils, the stock
tank oils (STOs) are usually deficient in light ends due to the gas removal
at a separator upstream of the stock tank. Therefore, wax thermodynamic
modeling of STO sometimes requires the simulation of the liquidsolid equi-
librium, as the gas phase is not present. Figure 3.2 shows the key steps in
performing wax thermodynamic modeling.
Construction of the
thermodynamic
equations
Simplification of the
thermodynamic
equations
FIGURE 3.2
Flowchart of the steps involved in wax thermodynamic modeling.
44 Wax Deposition
The first two steps as shown in Figure 3.2 will be introduced in detail in
the remainder of this chapter.
Vapor
f Vi
f Si f Li
Wax Liquid
FIGURE 3.3
Three-phase equilibrium among the vapor, liquid, and precipitated wax. The fugacity of each
component i in different phases should equal when phase equilibrium is achieved.
Thermodynamic Modeling of Wax Precipitation 45
The fugacity of species i, fiV, fiL, and fiS, depends on temperature, T, pres-
sure, P, and the molar composition of the phases. Equation 3.2 emphasizes
the dependency of fiV, fiL, and fiS on T, P, and the molar compositions:
Equating the fugacities of species in different phases, fiV = fiL = fiS, one
generates the phase equilibrium equations with respect to the molar compo-
sitions in different phases xi s, y i s, and si s:
nV y i + nL xi + nS si = niF (3.4)
iN y i = 1
iN xi = 1 (3.5)
iN si = 1
RT amixture
P= (3.6)
V bmixture V (V + bmixture ) + bmixture (V bmixture )
The mixing rules developed by Chueh and Prausnitz (1967) are used to cal-
culate amixture and bmixture for the mixture based on the as and bs of the pure
components, and the mixing rules typically involve empirical constants to
adjust the model to fit experimental measurements.
Since no EOS is available for the solid phase, the solidliquid phase equilib-
rium equation needs to be constructed in a different manner. We will go through
a step-by-step procedure for the development of solidliquid phase equilibrium
equations. Firstly, the fugacity of a pure solid is calculated based on the fugacity
of a pure liquid and the Gibbs free energy change from pure liquid to pure solid:
P
fiL, pure H i T Cpi Tf Cpi Ti f Vi
ln
fiS, pure
=
RT
1
Ti f
+
R
1 i
T
+
R
ln
T
+
RT d p (3.7)
P0
P
fL x L H i T Cpi Tf Cpi Ti f Vi
ln i S = ln i Si +
fi si i RT
1
Ti f
+
R
1 i
T
+
R
ln
T
+
RT d p
P0
(3.9)
Thermodynamic Modeling of Wax Precipitation 47
Cooling
140
20
Ttr,C23 Tm,C23
0
260 280 300 320 340
Temperature (K)
FIGURE 3.4
Change of enthalpy of n-C23 H48 subject to cooling. (From Dirand, M. et al., J. Chem. Thermodyn.,
34, 12551277, 2002.)
can coexist below the WAT, further complicating the thermodynamic mod-
eling of the liquidsolid phase equilibrium. Kravchenkos rule of thumb on
solid phase miscibility of n-alkanes suggests that total solid phase miscibil-
ity between two n-alkanes with different carbon numbers can only occur
when the difference in the carbon number is less than 6% the carbon number
of the alkanes (Kravchenko, 1946). An experimental study by Snyder, Conti,
Strauss, and Dorset (1993), Snyder, Goh, Srivatsavoy, Strauss, and Dorset
(1992), and Snyder et al. (1994) confirmed the existence of multiple immiscible
solid phases below the WAT. Figure 3.5 shows a phase diagram of a binary
mixture of C30 and C36.
As can be seen from the phase diagram, a hexagonal solid phase first
forms as the temperature decreases. The hexagonal phase starts to transit to
orthorhombic phase as the temperature is further decreased. The hexagonal
and orthorhombic phases can coexist during this transition depending on
the molar composition of the system. At temperatures below approximately
55C, two orthorhombic phases can coexist with each orthorhombic phase
rich in one of the two alkanes.
The existence of multiple solid phases poses further mathematical com-
plexity on wax thermodynamic modeling, as the fugacity of a particular
component i in each solid phase is equal, resulting in more phase equilib-
rium equations to be solved.
As discussed in this section, the following characteristics are important
for wax thermodynamic modeling:
80
L
L+H
70
60 H H+O2
O1 O2
50
Temperature (C)
40
30 Spontaneous
demixing at
room temperature
20
10
O1+O2
0
0.0 0.2 0.4 0.6 0.8 1.0
C30 X(C36) C38
FIGURE 3.5
Phase diagram of a binary mixture of C30 and C36 with varying concentration of C30: L: liquid;
H: hexagonal crystal; O1: orthorhombic crystal rich in C30; O2: orthorhombic crystal rich in C36.
(From Snyder, R. G. et al., J. Phys. Chem., 96, 1000810019, 1992.)
TABLE 3.1
Summary of Different Thermodynamic Models Treatments Regarding Phase Nonidealities and Existence of Multiple Solid Phases
Consider Liquid- Consider Solid- Consider Multiple Consider Secondary Incorporated in Commercial
Model Phase Nonideality? Phase Nonideality? Solid Phases? Phase Transition? Software Packages?
Conoco model No No Yes, empirically No Grand Unified
Thermodynamic Simulator
(GUTS)
Wons model Yes, by regular Yes, by regular No No No
solution theory solution theory
Pedersens model Yes, by regular Yes, by regular No No PVTsim
solution theory solution theory
Lira-Galeana et al.s Yes, by EOS No Yes, by phase No No
model stability analysis
Coutinhos model Yes, by Flory free Yes, by Wilson Yes, by UNIQUAC Yes Multiflash
volume theory and model or by model
by UNIQUAC UNIQUAC model
Functional-Group
Activity Coefficients
(UNIFAC) method
Wax Deposition
Thermodynamic Modeling of Wax Precipitation 51
model (Coutinho, 1998, 1999; Coutinho & Daridon, 2001; Coutinho, Edmonds,
Moorwood, Szczepanski, & Zhang, 2006; Coutinho & Stenby, 1996; Daridon,
Coutinho, & Montel, 2001; Pauly, Daridon, & Coutinho, 2004).
si Si H i T
ln = 1 f (3.12)
xi Li RT Ti
Vi ( i )2
ln i = (3.13)
RT
Example
The contribution to Gibbs free energy from the Poynting correction,
P
T
V d p, and the contribution from enthalpy change, H
P0
i i 1 ,
Ti f
will now be calculated. The material physical properties given in Table
3.2 were collected to carry out these calculations.
52 Wax Deposition
TABLE 3.2
Parameters Used for the Calculation of Phase Transition of n-C20H44
Material Property Value
Melting point (Tf) 310 K
Heat of fusion (H) 70 kJ/mol
Change of volume at phase transition (V) ~3.6 102 L/mol
Gas constant (R) 8.314 J/mol/K
P
L J
V d p V (p p ) = 3.6 10
P0
0
2
mol
101 kPa = 3.6
mol
(3.14)
T kJ 298 K J
H 1 = 70 1 = 2710 (3.15)
T f
mol 310 K mol
The contribution from the Poynting correction, 3.6 J/mol, is only 0.13%
of the contribution from the enthalpy of crystallization, 2710 J/mol.
However, the contribution from the Poynting correction increases with
pressure and can eventually become greater than the contribution from
the enthalpy of crystallization. Figure 3.6 shows the comparison of the
4.0
of fusion: H 1 T
3.0 Tf
Contribution from
2.0 Poynting correction:
p
~ ( p p0)V
p0Vdp
1.0
0.0
1 10 100 1000
Pressure (atm)
FIGURE 3.6
Comparison between the contribution to Gibbs free energy change in phase transition from
heat of fusion and the contribution from the Poynting correction.
Thermodynamic Modeling of Wax Precipitation 53
Example
Wons model will be used to calculate the amount of precipitated solid
at room temperature, i.e., T = 298 K, based on a (n-C10H44 + n-C24H50 +
n-C26H54) ternary mixture with the composition as shown in Table 3.3.
The goal of the thermodynamic simulation is to determine the number
of moles of the solid and liquid phase, nS and nL, respectively, as well as
the molar composition of each component in the two phases: s10, s24, s26,
x10, x24, and x26, where si and xi represent the mole fraction of component i
in the solid phase and the liquid phase, respectively. In order to calculate
these variables of interest, the physical properties of n-C10H44, n-C24H50,
and n-C26H54, including their heats of crystallization and melting point,
are needed as input parameters for thermodynamic modeling. These
parameters are collected from the National Institute of Standards and
Technology (NIST) database and are summarized in Table 3.4.
TABLE 3.3
Composition of the Ternary Mixture Investigated
in the Sample Thermodynamic Modeling
Component Content
n-C10H44 0.8 mol
n-C24H50 0.1 mol
n-C26H54 0.1 mol
54 Wax Deposition
TABLE 3.4
Material Properties of n-C10H44, n-C24H50, and n-C26H54
as Input Parameters to the Thermodynamic Model
s10 H 10 T
ln = 1 f (3.16)
x10 RT T10
s24 H 24 T
ln = 1 f (3.17)
x24 RT T24
s26 H 26 T
ln = 1 f (3.18)
x26 RT T26
The mole fractions of each species in the liquid and solid phases should
also satisfy the following constitutive equations:
nL = 0.822 mol
nS = 0.178 mol
s10 = 0.071
s24 = 0.427
(3.24)
s26 = 0.503
x10 = 0.958
x24 = 0.029
x26 = 0.013
As can be seen from these numerical results, the mole fraction of the
solid phase at room temperature is 0.178. The liquid phase at room tem-
perature is mainly composed of n-C10H44, whereas the solid phase is
mainly composed of the two heavy n-paraffins: n-C24H50 and n-C26H54.
Pedersen et al. improved Wons model by
100 100
Oil 1 Experiment (NMR) Oil 2
80 Model 80
Weight % wax
Weight % wax
60 60
40 40
20 20
0 0
200 250 300 350 200 250 300 350
Temperature (K) Temperature (K)
100 100
Oil 3 Oil 4
80 80
Weight % wax
Weight % wax
60 60
40 40
20 20
0 0
200 250 300 350 200 250 300 350
Temperature (K) Temperature (K)
FIGURE 3.7
Overprediction of the WPCs potentially due to the assumption of ideal solid phases. (From
Pedersen, K. S. et al., Energy & Fuels, 5, 924932, 1991.)
simulates the existence of multiple solid phases using a staged phase calculation
procedure: It calculates the partition of wax between the solid and liquid phase
at 1.8F temperature increments while assuming that the solid phase formed at
temperatures greater than T + 1.8F does not affect the phase equilibrium at T.
Different from the Conoco model, the thermodynamic model by Lira-
Galeana et al. treats the formation of multiple pure solid phases using phase
stability criteria. Thus, the Lira-Galena et al. and Conoco approaches model
precipitation as an expanding core of solids, somewhat analogous to an onion.
5
NMR
DSC
Wax precipitated (wt%)
4 WAT microscopy
WAT DSC
3 Wax 6 Wilson EOS
Wax 5 UNIQUAC EOS
2 Wax 4
Wax 3
1 Wax 2
Wax 1
0
0 10 20 30 40 50 60
Temperature (C)
FIGURE 3.8
Comparison between the WPCs by the Wilson model and the UNIQUAC model. (From
Coutinho, J. A. P. et al., Energy & Fuels, 20, 10811088, 2006.)
58 Wax Deposition
100
Experiment
80 Coutinho
Weight % wax
Won
60
Pedersen
40
20
0
250 270 290 310 330
Temperature (K)
FIGURE 3.9
Comparison of predicted WPC by different wax thermodynamic models with experimental
measurements. (From Pauly, J. et al., Fluid Phase Equilib., 149, 191207, 1998.)
study covering all mainstream models except the Conoco model and Lira-
Galeana et al.s model based on their predictions of the WPC of the same
model mixture. The comparison between different model predictions is
shown in Figure 3.9.
As can be seen from Figure 3.9, Coutinhos thermodynamic model pro-
vides the most accurate prediction for the WPC for this model mixture. All
other models slightly overpredict the WPC potentially due to the assump-
tion of ideal solid or liquid phases.
The single carbon number (SCN) distribution: A SCN fraction with car-
bon number i includes branched paraffin, cyclic paraffin, aromatics,
and n-paraffin with carbon number i. Therefore, the molecular for-
mula of the components in SCN with carbon number i can be writ-
ten as CiHj, where j = 2i + 2 represents saturated components and j
2i + 2 represents cyclic paraffins, aromatics, and other unsaturated
components.
The n-paraffin distribution: the mole/mass fraction of only n-paraffin,
n-CiH2i+2, as a function of carbon number i.
Postprocessing: tuning of
model against measured
WAT/WPC
FIGURE 3.10
Flowchart of performing wax thermodynamic modeling using software packages.
60 Wax Deposition
Integrated peak
Method A: Method B:
from valley to valley from baseline
FIGURE 3.11
Two different methods of integration to obtain SCN and n-paraffin distribution from oilgas
chromatograph. (From Coto, B., Coutinho, J. A. P., Martos, C., Robustillo, M. D., Espada, J. J., &
Pea, J. L., Energy & Fuels, 25, 11531160, 2011.)
5
By DSC
Precipitated paraffin (wt%)
Method B: from
4 baseline
By fractional
2 precipitation
1
Method A: from
valley to valley
0
50 40 30 20 10 0 10 20 30 40
Temperature (C)
FIGURE 3.12
Comparison of the WPC prediction based on n-paraffin distribution obtained from different
methods of integration: from valley to valley and from baseline with experimental measure-
ments by DSC and fractional precipitation. (From Coto, B., Coutinho, J. A. P. et al., Energy &
Fuels, 25, 11531160, 2011.)
Thermodynamic Modeling of Wax Precipitation 61
102
100
101
102
0 10 20 30 40
Carbon number
FIGURE 3.13
Composition of the crude oil including n-paraffin distribution of oil S measured by HTGC.
(From Zheng, S. et al., Energy & Fuels, 27, 73797388, 2013.)
the amount of the component with the highest carbon number in an HTGC
analysis is usually reported as one number including all C36+ components,
which is called a plus fraction. The wax precipitation characteristics in terms
of the WPC and the composition of the precipitated solid phase are sensitive
to the composition of the plus fraction because heavy components are the
first to precipitate as temperature is decreased below the WAT (Erickson et
al., 1993). Consequently, the mole fractions of the C36+, i.e., components with
carbon numbers greater than 36, and the corresponding n-paraffin fractions
need to be properly estimated and included in thermodynamic modeling
so that model predictions better represent the actual crude oil precipitation
characteristics.
Based on field experience, Pedersen, Thomassen, and Fredenslund (1985)
suggested extrapolating the SCN and the n-paraffin distribution of the heavy
end to a certain limiting carbon number by assuming that the amount of
hydrocarbon component CiHj and n-paraffin n-Ci-H2i+2 decreases logarithmi-
cally as the carbon number i increases. Pedersens method of extrapolation
has provided satisfactory modeling results for multiple crude oils (Hansen,
Fredenslund, Pedersen, & Rningsen, 1988). A sensitivity analysis showed
that for this particular crude oil, including C50+ alkane components in ther-
modynamic modeling affects the predicted precipitation curve by less than
0.1%. Figure 3.14 shows the composition of the oil with the mole fractions of
C36 C50 components extrapolated using Pedersens method.
In some other cases, the SCN and the n-paraffin distribution are extrapo-
lated to a carbon number as high as 75 (Erickson et al., 1993). In order for the
software packages to automatically extrapolate the plus fraction composition
to an appropriate limiting carbon number, additional information such as
62 Wax Deposition
102
100 Extrapolation
101
102
103
0 20 40 60
Carbon number
FIGURE 3.14
Extrapolated composition of the crude oil S including SCN distribution and n-paraffin distri-
bution. (From Zheng, S. et al., Energy & Fuels, 27, 73797388, 2013.)
the plus fraction density and the molecular weight should be supplied to the
software packages. If the plus fraction density and the molecular weight are
unknown, a trial-and-error method suggested by Calsep can be used to
adjust the plus fraction molecular weight and density until the flash calcula-
tion at stock tank condition predicts the reported STO molecular weight and
density.
Modeling the thermodynamic equilibrium of each n-paraffin and non-n-
paraffin components requires solving one vaporliquid and one liquidsolid
equilibrium equation for each component. With this crude oil, the compo-
sitional analysis resolved the weight distribution over approximately 100
components (both paraffin and nonparaffin). Therefore, approximately 200
equations need to be solved in order to determine the molar composition of
each component in each phase. The actual scenario is worse than it already
sounds. Each of the 200 equations contains the concentration of components
in the same phase. Therefore, each of the 200 equations contains the 100
unknown compositions. Solving this system of equation requires an enor-
mous amount of computer time. In order to reduce computational intensity,
components with neighboring carbon numbers are usually lumped to form
pseudo-components. Default PVTsim settings lumps real components into
12 pseudo-components (PVTsim 19 method documentation, 2009). After lump-
ing, only 24 equilibrium equations instead of some 200 equations need to be
solved.
After lumping, each pseudo-component represents a group of real compo-
nents, and thus, its physical properties should be an average of its constituent
real components. As important parameters in the thermodynamic modeling,
Thermodynamic Modeling of Wax Precipitation 63
the averaged critical pressure Pc, the averaged critical temperature Tc, and the
averaged acentric factor need to be calculated using empirical correlations.
Empirical correlations usually estimate Pc, Tc, and based on the molecular
weight of the pseudo-components.
The preceding three steps
TABLE 3.5
Summary of Tuning Options in Mainstream Thermodynamic Modeling Software
Packages and Their Shortcomings
Mainstream
Software Tuning Option Tune to Match Shortcoming
PVTsim Wt% of wax in oil WAT or wax Fluid might not represent oil
content or both composition after tuning
GUTS Wt% of wax at the WAT Inconsistent with lab testing
precipitation onset techniques
Multiflash n-Paraffin heat of fusion WAT Inconsistent with the physical
properties of the material
DSC thermogram
Integration
Pure component
properties
WPC
Correlation
n-Alkane-like
composition
SLE
New integration
WPC DSC thermogram
Correlation
SLE n-Alkane-like
composition
FIGURE 3.15
An iterative algorithm to improve the accuracy of DSC-measured WPC based on thermo
dynamic modeling. (From Coto, B. et al., Fuel, 89, 10871094, 2010.)
66 Wax Deposition
3.8Summary
In this chapter, different wax thermodynamic models were introduced.
Among these wax thermodynamic models, Coutinhos model is the most
rigorous model, while other models implement different assumptions for
model simplification. The assumption of ideal solidliquid phase behavior
usually leads to overestimation of the WAT and the WPC.
Wax thermodynamic models can predict the WAT, the WPC, and the com-
position of precipitated wax. These predictions provide important information
regarding wax precipitation characteristics to help assess the possibility and
severity of wax deposition issues. It should be noted that in order to obtain
reliable predictions from wax thermodynamic models, one should follow
certain established procedures. The industrial practice for using wax ther-
modynamic models was also introduced in this chapter. Reliable WAT and
WPC predictions by thermodynamic models can be used as input param-
eters for wax deposition modeling when the corresponding experimental
characterizations are not available. The application of wax thermodynamic
models in wax deposition modeling will be introduced in Chapters 4 to 6.
4
Wax Deposition Modeling
67
68 Wax Deposition
Pipe centerline
Flow direction
Dissolved wax
Bulk oil molecules
Wall
Deposit
Wall
FIGURE 4.1
Schematic of the two types of wax that could potentially cause wax deposition.
We can see that for the first mechanism, molecular diffusion, the subjects of
deposition are the dissolved wax molecules in the oil. The depositing mate-
rials in the other three mechanisms (shear dispersion, Brownian diffusion,
and gravity settling) are the suspended wax particles that have precipitated
from the oil. In order to identify which mechanism plays the primary role in
wax deposition, researchers have looked at the cause of each of these mecha-
nisms to see if the conditions where these mechanisms would occur are con-
sistent with what happens typically in a pipe flow that is subject to external
cooling (which would resemble the transportation of petroleum fluids in the
subsea pipelines). And several mechanisms were quickly identified to have
insignificant impact on wax deposition: The Brownian diffusion of wax par-
ticles is not likely to be valid because the wall temperature is lower than the
oil bulk temperature, resulting in more precipitated wax particles at the wall
than there are in the bulk oil. Consequently, the impact of Brownian diffu-
sion is to transport these precipitated particles from the wall toward the bulk
oil, instead of moving them toward the wall and deposit (Singh et al., 2000).
Gravity settling was also believed to be insignificant as there has not been a
study claiming that wax deposits are generally thicker at the bottom than at
the top of the pipe wall during single-phase oil flow conditions.
The remaining two mechanisms are shear dispersion and molecular dif-
fusion. At early times, the mechanism of shear dispersion was advocated
by several researchers (Burger et al., 1981; Hsu & Brubaker, 1995) and was
in fact incorporated in early deposition models (such as the wax deposition
model in the OLGA industrial flow simulator). However, another early study
by Bern et al. (1980) found that the wax deposition rate did not increase with
the increasing shear rate of the fluid, casting the first doubt on this mecha-
nism for wax deposition. In fact, many findings in the field of particle fluid
dynamics, both theoretical and experimental, seem to suggest that the wax
Wax Deposition Modeling 69
particles in the oil flow are not necessarily to be dispersed toward the wall:
The study by Saffman (1965) and Cleaver and Yates (1973) (based on particle
fluid mechanics theory) and those by Jimenez et al. (1988), Urushihara et al.
(1993), and Garcia et al. (1995) (based on particle velocimetry) all found that
particles located in the viscous layer near the wall tend to be re-entrained
into the bulk flow by a lifting force generated by the turbulent flow (known
as the Saffman lift force). These findings suggested that the mechanism
of shear dispersion is not significant or even possible. Based on the exten-
sive experimental observations on wax deposition over the past few decades
(Hunt Jr, 1962; Eaton & Weeter, 1976; Brown et al., 1993; Creek et al., 1999;
Singh et al., 2000; Hoffmann & Amundsen, 2010; Jemmett et al., 2012), it is
now generally accepted that the molecular diffusion is the main mechanism
for wax deposition. The precipitated wax crystals formed away from the wall
do not seem to contribute to the formation of the wax deposit but rather flow
with the oil to form a suspension in the fluid mixture. As of today, the mech-
anism of molecular diffusion has been used by many wax modeling studies
(Solaimany Nazar et al., 2001; Roehner & Fletcher, 2002; Hernandez et al.,
2003; Venkatesan, 2004; Edmonds et al., 2007; Akbarzadeh & Zougari, 2008;
Ismail et al., 2008; Lee, 2008; Merino-Garcia & Correra, 2008; Han et al., 2010;
Phillips et al., 2011; Huang et al., 2011a; Lu et al., 2012). Consequently, this
mechanism will be discussed in detail in the rest of this chapter.
Oil Oil
Deposit Deposit
Wall Wall
Step 1 Step 2
Oil Oil
Deposit Deposit
Wall Wall
Step 4 Step 3
FIGURE 4.2
Schematic of molecular diffusion as the wax deposition mechanism.
and form crystals. Precipitation of the waxy components can occur both at the
bulk oil and on the pipe wall as long as the temperature at that particular loca-
tion is below the WAT, as shown in Figure 4.2 (step 1). As discussed previously,
the precipitated wax crystals formed in the bulk are believed to flow with the
oil and not deposit on the pipe wall. Therefore, it is the precipitation of the
waxy components at the wall that forms the incipient layer of the wax deposit.
0.4
0.35
Wax content (wt. fraction)
0.3
0.25
0.2
0.15
0.1
0.05
0
0.0 1.0 2.0 3.0 4.0 5.0
Time (days)
FIGURE 4.3
Evolution of wax fraction in the deposit observed by Singh et al. (2000). The wax fraction in the
original oil is 0.0067. (Singh, P., Venkatesan, R., Fogler, H. S., & Nagarajan, N.: Formation and
aging of incipient thin film wax-oil gels. AIChE Journal. 2000. 46. 10591074. Copyright Wiley-
VCH Verlag GmbH & Co. KGaA. Reproduced with permission.)
Interlamellar disorders
nc
Growth direction Growth direction
FIGURE 4.4
Schematic of the crystal structure in the wax deposit proposed by Dirand et al. (1998) and
Singh et al. (2000). (Singh, P., Venkatesan, R., Fogler, H. S., & Nagarajan, N.: Formation and
aging of incipient thin film wax-oil gels. AIChE Journal. 2000. 46. 10591074. Copyright Wiley-
VCH Verlag GmbH & Co. KGaA. Reproduced with permission.)
0.27 mm 0.27 mm
FIGURE 4.5
Microscopic observation of the network structure of a waxy gel by Venkatesan et al. (2005).
(Reprinted from Chemical Engineering Science, 60, Venkatesan, R., Nagarajan, N. R., Paso, K., Yi,
Y.-B., Sastry, A. M., & Fogler, H. S., The strength of paraffin gels formed under static and flow
conditions, 35873598, Copyright 2005, with permission from Elsevier.)
Input
Temperature
increased insulation
profile
FIGURE 4.6
General algorithm for wax deposition modeling.
4.4.1Model Dimensions
The dimensions that a wax deposition model accounts for are often neglected
and are actually one of its most important aspects. Because wax deposition
occurs as a result of the radial heat loss of the petroleum fluids to the sur-
roundings, the radial profiles of temperature and the concentration of the
dissolved waxy components are essential for a wax deposition model to accu-
rately determine the growth rate of wax deposits. However, solving for the
transport characteristics in the radial direction with satisfactory resolution
often require significant computational intensity. Such intensity can become
a burden for industrial flow simulators, which are mostly discretized only
in the axial direction to look for characteristics of the fluids in the bulk flow.
An example of the temperature profile during single-phase oil flow is shown
in Figure 4.7.
With the current 1-D industrial flow simulators, a single value of tempera-
ture is used to represent the bulk oil with another value for the inner wall of
the pipe. A laminar sublayer is used to represent the region close to the wall,
where the temperature decreases linearly.
For the MWP (an academic wax deposition model), a 2-D (axial and radial)
discretization is used, as shown in Figure 4.8. Compared with the approach
of 1-D discretization, the 2-D discretization can account for radial variations
Pipe centerline
Temperature
profile in the oil
Bulk oil Toil
Boundary layer
Twall
Wall
FIGURE 4.7
Discretization of temperature profiles in the 1-D wax deposition models.
Wax Deposition Modeling 77
Pipe centerline
Temperature
profile in the oil
Bulk oil T = f (r,z)
Boundary layer
Wall
FIGURE 4.8
2-D discretization of temperature profiles in the MWP.
TABLE 4.1
Summary of Empirical Hydrodynamic Correlations for Different Wax Models
Huang, Lu,
Hoffmann,
Matzain Lindeloff and Amundsen,
Wax Deposition (1997); Rygg Krejbjerg Edmonds Hernandez and Fogler
Models et al. (1998) (2002) et al. (2007) et al. (2003) (2011)
Software OLGA wax PVTsim; FloWax TUWAX MWP
implementation DepoWax
Hydrodynamics Bendiksen, Mukherjee and Brill (1987) Xiao, Van Driest
correlation Maines, Brill (1985); Shonham, (1956)
Moe, and Bendiksen and Brill
Nuland etal. (1991) (1990);
(1991) Kaya,
Sarica, and
Brill (1999)
Multiphase Yes Yes Yes Yes Uses averaged
capability properties of
fluids as one
single-phase
fluid
2
Rpipe oil CpU dToil = 2 Rpipe hinternal (Toil Twall ) dz (4.1)
where
Rpipe = the radius of the pipe (m)
oil = density of the oil (kg/m3)
Cp = heat capacity of the oil (J/K/kg)
U = average velocity of the oil (m/s)
Wax Deposition Modeling 79
Tambient
U
Toil
Twall
Wall, insulation, etc. Heat loss
dz
Tambient
FIGURE 4.9
Schematic of heat transfer characteristics in a single-phase oil flow.
T = temperature (K)
z = axial coordinate (m)
hinternal = internal heat transfer coefficient based on the pipe inner diameter
(W/m2/K)
hexternal = external heat transfer coefficient (including the pipe wall, the pipe
insulation, and the ambient environment) based on the pipe inner
diameter (W/m2/K)
In this case, the temperatures of the bulk oil and at the inner wall, Toil and
Twall, are the two variables that are needed to be solved using Equations 4.1
and 4.2. The heat transfer coefficients can be estimated using the empirical
correlations summarized in Table 4.2. It should be noted that for multiphase
flows, most wax deposition models still use single-phase heat transfer cor-
relations with averaged properties of the phases to estimate the temperature
profiles.
The empirical correlation for the internal heat transfer coefficient does
not need to be used for the MWP, which has a 2-D discretization scheme.
Instead, using the energy balance in both the radial and axial directions
(Equation4.3) with appropriate boundary conditions (Equation 4.4), both the
radial and axial temperature profiles can be solved:
T (r , z) 1 k T (r , z)
Vz (r ) = r thermal + oil (4.3)
z r r oil Cp r
80 Wax Deposition
TABLE 4.2
Summary of Heat Transfer Correlations for Different Wax Deposition Models
Rygg et al. and
Matzain et al. Lindeloff and
(as Krejbjerg (as Edmonds
Implemented Implemented etal. (as
in OLGA in PVTsim Implemented
Wax) DepoWax) in FloWax) TUWAX MWP
Heat transfer Sieder and Tate (1936), Gnielinski (1976), or Sieder and Tate Deen (1998)
correlation Dittus and Blter (1985) (1936) oneddy
(single-phase thermal/
or oilwater mass
flow with diffusivities
averaged (single-phase)
Multiphase Uses averaged properties of fluids as one properties of Uses averaged
capability single-phase fluid fluids) properties of
Kim et al. (1999) fluids as one
(Gasliquid, single-phase
flow fluid
regime
dependent)
T
= 0, at r = 0
r
(4.4)
T
hexternal (Tambient Twall ) = koil , at r = Rpipe
r
Here, r(m) is the radial coordinate; thermal (m2/s) is the eddy thermal dif-
fusivity, which is determined by the von Krmn correlation (Wilkes, 2005);
and koil (W/m/K) is the thermal conductivity of the oil. The details of all the
transport equations in this model have been reported in the study of Huang,
Lee et al. (2011).
the precipitated wax particles in the bulk oil do not contribute to wax depo-
sition on the pipe wall. Instead, it is the radial diffusion of dissolved wax
molecules that results in the buildup of a wax deposit. The magnitude of
diffusion can be described in terms of radial mass fluxes (usually with the
unit of kg/m2/s), which are shown as JA and JB in Figure 4.10. These two
parameters can be further related to the deposition process based on a series
of mass balance equations, as shown in Equations 4.5 and 4.6, that are based
on the study of Singh et al.
Mass balance at the oildeposit interface is given by
d deposit
deposit Fwax = ( J A JB ) (4.5)
dt
( 2
deposit Rpipe (Rpipe deposit )2 dF )
wax
= JB (4.6)
2 (Rpipe deposit ) dt
where deposit (kg/m3) represents the density of the deposit; Fwax is the wax
mass fraction (the summation of all the waxy components) in the wax deposit
(both dissolved in the entrapped oil and precipitated as crystals); deposit (m)
represents deposit thickness; Rpipe (m) represents the radius of the pipe; and
t represents time (s).
Equation 4.5 represents a mass balance on dissolved waxy components at
the oildeposit interface, whereas Equation 4.6 is the mass balance on the
dissolved waxy components in the deposit, with a radially averaged value of
wax fraction in the deposit, Fwax. From these balance equations, it can be seen
Pipe centerline
Fluid flow
Bulk oil
Rpipe deposit
Rpipe
Wall
FIGURE 4.10
Schematic of the mass fluxes of the dissolved waxy components.
82 Wax Deposition
that the difference of the two mass fluxes (JA JB) corresponds to the growth
of the deposit layer, while the internal diffusion mass flux, JB, results in an
increase in the wax fraction within the deposit. The increase in the wax frac-
tion in the deposit with time, dFmax/dt, is also known as the aging rate of the
wax deposit, which is physically related to the hardening of the wax deposit
as time increases.
d deposit
deposit Fwax = J A (4.7)
dt
where deposit is the porosity of the wax deposit. Consequently, the mathe-
matical representation of molecular diffusion neglecting aging is shown in
thefollowing:
d deposit
deposit (1 deposit ) = J A (4.9)
dt
where
C = the concentration of dissolved waxy components (kg/m3)
Dwax = the diffusion coefficient of wax in oil (m2/s)
T = temperature (K)
mass transfer = the thickness of the mass transfer layer (m)
Pipe centerline
Concentration
profile of dissolved
wax molecules
Bulk oil Coil
Wall
FIGURE 4.11
Schematic of concentration profile for the 1-D wax deposition models.
84 Wax Deposition
T 1.47 B 10.2
Dwax = AHM
HaydukMinhas 0.71
, = 0.791 (4.11)
VA VA
( M )0.5 T
Dwax = BWC B B 0.6 (4.12)
WilkeChang
BVA
The values of the diffusion parameters AHM and BWC depend on the unit
ofthe diffusion coefficient. When the diffusion coefficient has a unit of m2/s,
the default values for AHM and BWC in these two correlations are 7.4 1012
and 13.3 1012, respectively. If the diffusion coefficient has a unit of cm2/s,
the default values of AHM and BWC becomes 7.4 10 8 and 13.3 10 8 respec-
tively. However, it should be noted that this correlation is originally designed
based on binary dilute solutions. In some wax deposition studies of petro-
leum fluids with a polydispersed composition, A and B become adjustable
parameters that are to be benchmarked by lab-scale deposition experiments
(Kleinhans, Niesen, & Brown, 2000).
Dissolved wax
Cwall
Twall Toil
Temperature Temperature
FIGURE 4.12
Conversion from the wax precipitation curve to the wax solubility curve.
(Coil Cwall) to the temperature difference (Toil Twall) between the bulk oil and
at the wall. The temperature difference is known to be the thermal driving
force for wax deposition. In many wax deposition studies, this thermal driv-
ing force is used as the indicator to cross-compare deposition experiments
at different temperatures (Creek et al., 1999; Singh et al., 2000; Paso & Fogler,
2004; Jennings& Weispfennig, 2005). However, it should be noted that two
assumptions are made in these analyses using the thermal driving force:
It can be seen that these two conditions are not generally valid, and thus
the use of thermal driving force often has its limitations. Such limitations
will be discussed in detail in Chapter 6. As to the second assumption, some
of the academic studies have tried to incorporate precipitation kinetics to
their models. For example, the TUWAX model includes an additional kinetic
precipitation step on the wall, so that the concentration at the wall is not nec-
essary at thermodynamic equilibrium, depending on the value of the surface
precipitation rate constant in the study of Hernandez et al. (2003). However,
in the MWP, instead of applying precipitation kinetics on the pipe wall, it
is applied on the boundary layer of the fluid, where the oil is expected to
experience drastic cooling while flowing in the pipeline. This implementa-
tion was achieved by using a 2-D mass transfer equation (axial advection and
radial diffusion), as shown in the following equation:
C 1 C
Vz = r( mass + Dwax ) k precipitation ( C C(eq )) (4.13)
z r r r
86 Wax Deposition
102
2.5
Dimensionless thickness, /ri kprecipitation = 0 (no precipitation)
2.0
1.0
Experiment
0.5
kprecipitation = 108 s1
(instantaneous precipitation)
0.0
0 5 10 15 20 25
Time (h)
FIGURE 4.13
The impact of wax precipitation kinetics in the oil flow on wax deposition predicted by the
MWP in the study of Huang, Lee et al. (2011).
It is seen that the second term on the right-hand side of Equation 4.13,
a kinetic precipitation term kprecipitation(C C(eq)), is used to account for the
kinetics of the precipitation of the dissolved wax molecules when their
concentration is above the solubility limit (Lee, 2008). Depending on the
value of the kinetic precipitation rate constant, kprecipitation, the concentra-
tion of dissolved wax molecules in the boundary layer can follow thermo-
dynamic equilibrium (kprecipitation ) or slow precipitation (kprecipitation
0s1). As theprecipitated wax crystals are found to most likely flow with
the oil instead of depositing on the pipe wall, a higher precipitation rate
would likely reduce the amount of dissolved waxy component in the
boundary layer potentially available for diffusion toward the pipe wall.
Consequently, a greater precipitation rate in the boundary layer actually
leads to less deposition on the pipe wall, as the precipitated wax crystals
often flow with the oil instead of depositing on the pipe wall. A compar-
ison with various experimental studies showed that a reasonable value
of kprecipitation can be on the order of 1.4 s1 (Huang, Lee et al., 2011) and is
shown in Figure 4.13.
Cbulk Cwall Sh
J A = Dwax = kmass transfer (Cbulk Cwall ) = Dwax (Cbulk Cwall ) (4.14)
mass transfer d
dpipe
mass transfer = diffusion = (4.15)
Sh
Equation 4.16. In this case, it is assumed that all the supersaturated wax mol-
ecules precipitated instantaneously in the oil:
dC Tbulk Twall
J A = Dwax (4.16)
dT T thermal
wall
In Equation 4.17, thermal represents the thickness of the layer where the
decrease in temperature can be assumed to be linear. Using a similar deri-
vation to the diffusion-equivalent layer shown previously, thermal can be
referred to as the conduction-equivalent layer. Analogous to Equation
4.15, thermal is equal to the pipe diameter divided by the Nusselt number, as
shown in the following equation:
d pipe
thermal = (4.17)
Nu
Similar to the clarification made previously, one would like to point out
the fact that this conduction-equivalent layer should not be mistaken as the
thermal boundary layer, as the thermal boundary layer often refers to the
layer where the temperature experiences a 99% decrease.
Comparing Equations 4.15 and 4.17, it can be seen that the conversion from
the diffusion-equivalent layer to the conduction-equivalent layer is merely
the ratio between the Sherwood number and the Nusselt number, as shown
in the following equation:
thermal Sh
= (4.18)
mass Nu
Using the correlation of Dittus and Boelter (1985) and the Chilton
Colburn analogy between heat and mass transfer in turbulent flows
(Chilton & Colburn, 1934), the Sherwood and Nusselt numbers can be
replaced by the Reynolds number, the Schmidt number (for mass trans-
fer), or the Prandtl number (for heat transfer), as shown in the following
equations:
Nu = 0.023Re0.8Pr0.3 (4.19)
Sh = 0.023Re0.8Sc0.3 (4.20)
Combining Equations 4.18 through 4.20, it is seen that the conversion from
the conduction-equivalent layer thickness to the diffusion-equivalent layer
Wax Deposition Modeling 89
thickness can be simply made via the Lewis number, which depends only on
the properties of the oil:
0.3
thermal Sh Sc koil
= = = Le0.3 = (4.21)
mass Nu Pr oilCp Dwax
4.5Summary
In this chapter, various aspects of wax deposition modeling are introduced,
and we start with the mechanism of wax deposition. Despite multiple mecha-
nisms proposed in the 1980s, it is now generally believed that molecular diffu-
sion is the primary mechanism for wax deposition. Based on this mechanism,
wax particles that precipitated in oil do not adhere to the wall and form a
deposit. Instead, wax deposition is mainly caused by the radial diffusion of
dissolved molecules of waxy components in the oil boundary layer toward the
wall (or toward the deposit surface once a layer of deposit has formed). After
the arrival of these dissolved waxy components at the deposit surface, they
can either precipitate on the surface (resulting in the further increase in the
deposit thickness) or continue to diffuse into the deposit and eventually pre-
cipitate to form a network structure (leading to the hardening of the deposit).
The mechanism of molecular diffusion has been widely used in many indus-
trial and academic wax deposition models. These models utilize the input of
operating conditions as well as the oil properties. Among these oil properties,
the most important parameter is the wax precipitation curve, which is often
obtained from the output of a wax thermodynamic model. With these two
pieces of information, a wax deposition model predicts the rate of the deposi-
tion on a pipe wall by performing the following transport calculations:
Most of the wax deposition models apply empirical correlations for trans-
port calculations. While different wax deposition models apply various
transport correlations, the main goal among these models is too similar,
which is to determine the radial mass flux of wax molecules, JA.
The key parameters in the mass flux, JA include the diffusivity of wax in
oil, the concentration difference, as well as the mass transfer thickness. Due
to the uncertainties associated with these parameters, wax deposition mod-
els are often benchmarked to identify the optimum model configurations for
field predictions.
5
Introduction to Wax Deposition Experiments
91
92 Wax Deposition
Heat
Oil storage tank Pump exchanger
Inlet
temperature
+ measurement
Test section
Differential
Coolant dp pressure
Heating storage tank
measurement
Removable
section
Flow
Pump (optional)
meter
Heat +
exchanger Outlet
temperature
measurement
FIGURE 5.1
Schematic of a typical flow-loop wax deposition apparatus.
Heating
FIGURE 5.2
Principles of estimating the depletion effect to determine the appropriate size of the oil storage
tank for wax deposition experiments.
Introduction to Wax Deposition Experiments 93
5.2.2Test Section
The test section is the most important part of a flow loop where wax deposition
is designed to take place. In the test section, oil is cooled by a cooling annulus.
Water is often selected as the coolant due to its high heat capacity (Hoffmann &
Amundsen, 2010). In some other studies, water is mixed with glycol to decrease
its freezing point in order to increase its cooling capacity (Hernandez, 2002). The
inlet temperature of the coolant is often controlled by a chiller.
A flow development section should be placed upstream of the test sec-
tion to allow the velocity profile to be fully developed. During the course
of wax deposition, the differential pressure across the test section increases
due to wax depositing on the wall and thereby decreasing the effective pipe
diameter. The effective diameter can be calculated from the pressure drop,
which can be measured with differential pressure transducers/transmitters
installed at the inlet and outlet of the test section.
Complications arise for the temperature characterization in the test sec-
tion. Temperature measurements within the test section are usually not as
satisfactory because the measuring device, such as a thermal couple, is likely
to interfere with the oil flow field as well as having wax depositing on its sur-
face, which greatly compromises its accuracy. Therefore, thermal couples are
usually installed upstream of the test section inlet and downstream of the
test section outlet to monitor the bulk temperature of the oil and the coolant.
It would be ideal to be able to measure the temperature at the inner pipe
wall or deposit surface because this temperature represents the local temper-
ature condition under which wax deposition occurs. However, the tempera-
ture close to a solid surface varies greatly within the boundary layer adjacent
to the surface (an example is highlighted in Figure 5.3). The thickness of the
boundary layer is typically less than 1 mm. Consequently, if the probe of a
thermal couple is not placed exactly at the inner wall surface, the boundary
layer of the flow can be interrupted, and the wall temperature might not be
correctly measured.
94 Wax Deposition
Flow centerline
Bulk temperature: 20C
Wall temperature: 12C
Re ~13,600
Pr ~32
Flow direction
Boundary layer
~0.3 mm
Wall
Temperature (C)
10. 12. 14. 16. 18. 20.
FIGURE 5.3
An example of the thickness of the boundary layer in normal turbulent pipe flow conditions.
5.3Deposit Characterization
At the end of a flow-loop wax deposition experiment, a layer of wax deposit
is formed on the inner pipe wall surface. The characterization of the deposit
should be carefully carried out, as the measurement results often provide
critical information for theoretical validation and model bench-marking.
Two parameters are usually measured: the thickness of the deposit and the
composition of the waxy components in the deposit. In this section we will
discuss different methods for these investigations.
which was used in the deposition study of Hoffmann and Amundsen (2010).
This correlation to calculate the effective pipe diameter is shown by
2
1.11
pipe 6.9
fDarcy = 1.8 log 10 + ( 4000 < Re < 108 ) (5.1)
3.7 dpipe Re
2
8 fDarcy oil LpipeQoil
Ppipe = (5.2)
5
2 dpipe
It can be seen from Equation 5.1 that the friction factor and the pressure
drop depend on the roughness of the pipe, pipe, and the Reynolds number of
the flow, which are often associated with a degree of uncertainties.
Several caveats should be noted when using this correlation. First, as the
wax deposit forms a layer of coverage on the pipe wall, the roughness of the
deposit surface instead of that of the pipe should be used in Equation 5.1.
However, the surface morphology of the deposit surface depends on mul-
tiple factors, including the composition of the waxy components as well as
the operating conditions. Therefore, the roughness of the deposit varies from
test to test. Based on the above Haaland correlation, the impact of the rough-
ness is shown in Figure 5.4. The y-axis represents the ratio between the two
10,000
1.11
100
pipe 6.9 Re = 2 106
/
3.7dpipe Re
10
Re = 2 105
1
Deposit
Re = 2 104 thickness as
roughness
0.1
1.E-03 1.E-02 1.E-01
Relative pipe roughness, pipe/dpipe
FIGURE 5.4
Contribution of friction pipe roughness on the friction factor at different values of Re.
96 Wax Deposition
terms inside the bracket on the right-hand side of Equation 5.1. This ratio
reflects the impact of the roughness on the friction factor. It can be seen
that such impact increases with increasing Reynolds number. For exam-
ple, based on a 2-in. steel pipe for laboratory wax deposition experiments,
the roughness of the bare pipe surface is approximately 50 m. This value
corresponds to a relative roughness (pipe/dpipe) of 2 10 3. As shown in
the dash line on the left of Figure 5.4, the ratio for a Reynolds number of
20,000 is less than 1, indicating that the friction factor (and thus the pres-
sure drop) does not significantly depend on the roughness at this range.
On the other extreme, the conclusion would be quite different if one
prefers to use the thickness of the deposit as the pipeline roughness in
this correlation, as was the approach in the study of Noville and Naveira
(2012). In this case, the value of the roughness parameter can be up to a
few millimeters. This value corresponds to a relative roughness of around
0.07 and is highlighted in the dash line on the right of Figure 5.4. The cor-
responding ratio for the same Reynolds number (20,000) is more than 50,
indicating that the parameter of roughness has a significant influence on
the pressure drop. In this case, as the thickness of the deposit is used as
the roughness of the pipe, the pressure drop can be rather sensitive to the
thickness of the deposit. Intuitively, the approach of using the deposit
thickness as the pipeline roughness for the pressure dropdiameter cor-
relation should be used only when the wax deposits are not expected to
have complete coverage of the test section (for example, in some of the
multiphase flow conditions where the oil does not have contact with part
of the pipe wall).
For wax deposition experiments in single-phase oil flow where the depos-
its form a complete coverage in the test section, Hoffmann and Amundsen
(2010) developed a technique that introduced short-term variations in the
flow rate to estimate the roughness of the wax deposit surface at different
operating conditions, and a range of 540 m was reported.
The second parameter that is often associated with uncertainty is the
viscosity of the oil, which is often used to determine the Reynolds num-
ber of the flow. Most of the available friction factor correlations are devel-
oped based on Newtonian fluid, where the viscosity used in determining
the Reynolds number is independent of shear rates. However, the precipi-
tated wax crystal in the oil flow can cause non-Newtonian behavior in the
flow. For non-Newtonian flow, the viscosity depends on the shear rate in the
flow field and thus is dependent on the oil flow rate. In addition, for a wax
deposition experiment with a certain oil flow rate, the variation in the oil
viscosity in the radial direction needs to be accounted for, as the shear rate
greatly changes radially in the flow field. Therefore, methods such as mul-
tidimensional computational fluid dynamics are needed to incorporate the
rheological characteristics of the non-Newtonian flow into the 2-D (axial and
radial) momentum equation to evaluate the pressure drop across the test sec-
tion in order to eventually determine the deposit thickness (Crochet, Davies,
Introduction to Wax Deposition Experiments 97
1 douter 1 d d
= + outer ln outer
U overall deffective inner hinternal 2 k pipe dinner
(5.5)
d dinner 1
+ outer ln +
2 kdeposit deffective inner hcoolant
It can be seen from Equation 5.5 that the overall heat transfer resistance,
1/Uoverall, represents the summation of each thermal resistance in the radial
direction (the oil phase, the deposit, the pipe wall, and the coolant). The ther-
douter 1
mal resistance of the oil phase, , depends on the heat transfer
deffective inner hoil
Oil
deffective inner
Heat loss dinner
douter
Deposit
Q
Pipe wall
Coolant
Tcoolant,inlet Tcoolant,outlet
Annulus centerline
FIGURE 5.5
Schematic of radial heat transfer during the course of wax deposition.
Introduction to Wax Deposition Experiments 99
coefficient in the oil phase, hoil, which can be obtained by empirical correla-
tions such as the Sieder and Tate correlation for turbulent flow:
0.14
koil koil
hoil = Nu oil = 0.027 Re0.8 Pr 1/3 centerline
deffective inner deffective inner Tdeposit interface
Lpipe
> 10
deffective inner
0.14
Tcenterline
It can be seen that a viscosity ratio, , is used to account for
Tdeposit interface
the effect of nonuniform radial viscosity due to radial heat transfer. The vis-
cosity of oil at the oildeposit interface can be difficult to determine. In some
cases, the viscosity ratio is assumed to be unity (Lund, 1998). The thermal
1
resistance of the coolant phase, , can be determined in a similar way
hcoolant
by calculating the heat transfer coefficient of the coolant side using empirical
correlations.
The other two thermal resistances are from the deposit layer and from the
pipe wall, and their values depend on the parameters related to the size and
the thermal properties of the deposit and the pipe wall. Among these param-
eters, the effective diameter of the oil flow is directly related to the thickness
of the deposit, as shown in Equation 5.7. Consequently, by solving Equation
5.6 for the effective diameter of the oil flow, defective inner, the deposit thickness
can be simply determined:
It should be noted that the heat transfer method suffers from a few disad-
vantages: First, for a laboratory-scale test section that is on the range from a
few meters to tens of meters, the difference in the temperatures between the
inlet and the outlet might not be significant. In this case, the uncertainty in
the thermal couple measurements can greatly compromise the accuracy of the
deposit thickness. In addition, the thermal conductivity of the deposit must be
known a priori for this measurement. However, the deposit represents a mix-
ture of liquid oil and wax crystals, and the wax crystal can have a thermal con-
ductivity that is 50%100% greater than that of the oil (Singh, 2000). Therefore,
the thermal conductivity of the deposit greatly depends on the oil fraction of
the deposit, which is not possible to determine during a wax deposition exper-
iment. Even after a wax deposition experiment, challenges still exist where the
100 Wax Deposition
oil fraction could change due to the change in temperature during the transfer
of the deposit from the pipe to the measuring device. Consequently, the heat
transfer method is not widely used to determine the deposit thickness.
Pipe wall
Mechanical support
FIGURE 5.6
Schematic of the laser technique for deposit thickness measurement. (Reprinted with per-
mission from Hoffmann, R., & Amundsen, L., Energy & Fuels, 24, 10691080. Copyright 2010
American Chemical Society.)
Introduction to Wax Deposition Experiments 101
using the pressure-drop technique and weight technique, and the laser tech-
nique showed satisfactory agreement and consistency with other methods.
However, due to interference with the flow field and device stability issues,
this method can only be performed after the wax deposition experiment is
completed and thus can only measure the final thickness of the wax deposit.
0.025 2h
12 h
1 day
0.015 2 days
3 days
C16C20
5 days
C15
0.005
C21
C22
C23
C24
C25
C26
C27
C28
Fi(t) Fi(0)
C29
C30
C31
C32
C33
C34
C35
C36
C37
0.005
Carbon number
0.015
0.025
0.035
FIGURE 5.7
Enrichment of heavy waxy components during the course of wax deposition discovered in the
study of Singh et al. (2000). (Singh, P., Venkatesan, R., Fogler, H. S., & Nagarajan, N.: Formation
and aging of incipient thin film wax-oil gels. AIChE Journal. 2000. 46. 10591074. Copyright
Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission.)
6
Fluid
Deposit with Twall = 7.5C
5 Deposit with Twall = 12.5C
Deposit with Twall = 17.5C
4
Increase in the critical
carbon number with
Weight %
0
10 15 20 25 30 35 40 45 50 55 60
Carbon number
FIGURE 5.8
Increase in the CCN in the wax deposit for deposition experiments with increasing wall tem-
perature. (Reprinted with permission from Hoffmann, R., & Amundsen, L., Energy & Fuels, 24,
10691080. Copyright 2010 American Chemical Society.)
Introduction to Wax Deposition Experiments 103
Chiller
TC
TC
Water bath
Heater
FIGURE 5.9
Schematic of the cold-finger wax deposition apparatus.
104 Wax Deposition
Shear cell
Mechanical
L mounting
stand
Variable speed
DC motor
FIGURE 5.10
Schematic of the organic solid deposition cell in the study of Akbarzadeh and Zougari (2008).
(Reprinted with permission from Akbarzadeh, K., & Zougari, M., Industrial and Engineering
Chemistry Research, 47, 953963. Copyright 2008 American Chemical Society.)
5.5.2Oil Characterization
The properties of oil are of great importance to comprehend the wax depo
sition testing results. Consequently, the oil should be well characterized.
The composition of the oil (especially the n-paraffin distribution) should be
measured, and wax precipitation tests should be carried out (detailed wax
thermodynamic characterizations and modeling are introduced in Chapters
2 and 3, respectively). The viscosity of the oil should be measured at dif
ferentshear rates and temperatures to assess the potential non-Newtonian
behavior that might be caused by the precipitated wax crystals. These pre
deposition characterizations provide essential information on waxoil thermo
dynamics and are an integral part to benchmarking wax deposition models
with the deposition experiments.
5.6Summary
In Chapter 2, we discussed different experimental techniques to study wax
precipitation. In this chapter, we introduced several typical wax deposition
apparatuses and a general procedure to carry out flow-loop wax deposition
experiments. The information as well as the referenced studies in these two
chapters provide a comprehensive description on the various tools for those
who look forward to developing a laboratory facility or extending an exist-
ing laboratory to analyze wax precipitation/deposition. In Chapter 6, we will
discuss the results of a series of wax deposition experiments and reveal some
of the most influential parameters governing the physics of wax deposition
through an in-depth theoretical analysis on the experimental results.
6
Applying Wax Deposition Models
to Flow Loop Experiments
6.1Introduction
The goal of wax deposition modeling is to achieve rigorous and reliable pre-
dictions on the severity of wax deposition for the pipelines in the field. To
achieve this goal, developing a comprehensive understanding of the per-
formance and limitations of wax deposition models is essential. While field
data would be the ideal candidate to directly test the wax deposition models,
it is usually not available nor as well controlled and monitored.
Contrarily, laboratory wax deposition experiments are able to provide
unique opportunities to benchmark wax deposition models due to ease of
experimental design and characterization. The results from laboratory exper-
iments can often help us develop educated guesses for the pipelines in the
field. The process of applying wax deposition models to well-characterized
flow-loop experiments helps one select the most appropriate deposition
model, to understand the governing physics for the deposition phenomenon,
and to identify the optimum configuration of the adjustable parameters in a
wax deposition model. In this chapter, the process of conducting wax deposi-
tion experiments and the analyses of these experiments with wax deposition
models will be applied.
107
108 Wax Deposition
A wax deposition experiment often provides two main results: the thick-
ness of the deposit as well as the composition of the wax deposit as a func-
tion of time. The thickness of the deposit is measured in nearly all of the wax
deposition experiments, and it is often selected as the variable to benchmark
wax deposition models. The composition of the deposit provides the most
valuable insights as to the diffusion of different waxy components during
the course of deposition. Unfortunately, it is not always measured and ana-
lyzed. The theoretical analysis on these two types of results will be discussed
extensively in the rest of this chapter.
6.3.1.2Test Oil
The oil used for the experiments is a waxy gas condensate from the North
Sea with a density of 809 kg/m3 and a WAT of 20C measured by DSC.
Measurements on the condensate using acetone precipitation (UOP 4664
method) revealed a wax content of approximately 4.5 wt%. The viscosity
of the oil was characterized at different shear rates using a Physica MCR
301rheometer, as shown in Figure 6.2.
The wax precipitation curve was measured by several methods, as
shown in Figure 6.3. The details of these methods were previously intro-
duced in Section 2.2. The method of characterization using centrifugation
and HTGC was previously reported by Han, Huang, Senra, Hoffmann, and
Fogler (2010). Another method utilizes the DSC based on an averaged value
of the heat of crystallization of 200 J/g. Considering that the actual heat
Applying Wax Deposition Models to Flow Loop Experiments 111
Removable part of
the test section Oil
reservoir
(4 m3)
dp cell
FIGURE 6.1
Schematic of the flow loop in Statoils Herya Research Center. (From Hoffmann, R., &
Amundsen, L. Energy & Fuels, 24, 10691080, 2010; Huang, Z. et al., Energy & Fuels, 25, 5180
5188, 2011.)
102
100 1/s
200 1/s
400 1/s
1000 1/s
2000 1/s
Viscosity (mPa s)
101
100
0 5 10 15 20 25 30 35 40
Temperature (C)
FIGURE 6.2
Viscosity curves of the North Sea condensate at different shear rates. (Reprinted with per-
mission from Hoffmann, R., & Amundsen, L. Energy & Fuels, 24, 10691080. Copyright 2010
American Chemical Society.)
of crystallization of the waxy components often lies between 300 and 100
J/g (Hansen et al., 1991), the methods of centrifugation + HTGC and DSC
appear to be consistent with one another.
6.3.1.3Experimental Conditions
Three series of wax deposition experiments were carried out: one series with
different oil temperatures, another with coolant temperatures, and the third
112 Wax Deposition
1.5
0.5
0
0 10 20 30 40
Temperature (C)
FIGURE 6.3
The precipitation curve of the North Sea condensate characterized with different methods.
(Reprinted with permission from Han, S., Huang, Z., Senra, M., Hoffmann, R., & Fogler, H. S.
Energy & Fuels, 24, 17531761. Copyright 2010 American Chemical Society.)
one with different oil flow rates. The detailed configurations of these series
of experiments are summarized in Table 6.1.
The thickness of the wax deposit was determined by three methods: the
pressure-drop technique, the weight technique, and the laser technique. The
pressure-drop measurement was nonintrusive and was carried out contin-
uously during the course of wax deposition, while intrusive methods like
TABLE 6.1
Summary of the Series of Wax Deposition Experiments Carried Out in Statoils
Herya Research Center
Oil (Inlet) Coolant (Inlet)
Temperature Temperature Oil Flow Coolant Flow
(C) (C) Rate (m3/h) Rate (m3/h)
Varying oil 15 5 20 5
temperature 20
25
35
Varying coolant 20 5 5
temperature 10
15
Varying oil flow 20 10 5
rate 10
15
21
25
Applying Wax Deposition Models to Flow Loop Experiments 113
weight and laser measurements were carried out at the end of the experiment.
A comparison of these three methods is highlighted in Figure6.4, indicating
the consistency among these measurement techniques. The reproducibility
of the experiments is also confirmed by repeating one of the experimental
conditions listed in Table 6.1 four times (oil temperature: 20C; coolant tem-
perature: 10C; oil flow rate: 5 m3/h; and coolant flow rate: 5 m3/h), as shown
in Figure 6.5.
2.5
Pressure drop
Weight
Laser
2
Wax thickness (mm)
1.5
0.5
0
0 50 100 150 200
Time (h)
FIGURE 6.4
Comparison of the deposit thickness measured by different methods. (Reprinted with per-
mission from Hoffmann, R., & Amundsen, L. Energy & Fuels, 24, 10691080. Copyright 2010
American Chemical Society.)
2
Thickness (mm)
0
0 50 100 150 200
Time (h)
FIGURE 6.5
Reproducibility of the wax deposition experiments.
114 Wax Deposition
3
Deposit thickness (mm)
2
Qoil = 5 m3/h
1 Qoil = 10 m3/h
Qoil = 15 m3/h
Qoil = 21 m3/h
Qoil = 25 m3/h
0
0 20 40 60 80
(a) Time (h)
3
Deposit thickness (mm)
Qoil = 5 m3/h
2 Qoil = 10 m3/h
Qoil = 15 m3/h
Qoil = 21 m3/h
Qoil = 25 m3/h
0
0 20 40 60 80
(b) Time (h)
FIGURE 6.6
Comparison of deposit thickness as a function of oil flow rate, Qoil, between (a) the experi-
mental results and (b) the model prediction by the MWP. The oil temperature, Toil, was main-
tained constant at 20C, and the coolant temperature Tcoolant, was maintained constant at 10C.
(Reprinted with permission from Lu, Y., Huang, Z., Hoffmann, R., Amundsen, L., Fogler, H. S.,
& Sheng, Z. Energy & Fuels, 26, 40914097. Copyright 2012 American Chemical Society.)
Applying Wax Deposition Models to Flow Loop Experiments 115
3
Deposit thickness (mm)
1 Toil = 15C
Toil = 20C
Toil = 25C
Toil = 35C
0
0 10 20 30 40 50 60
(a) Time (h)
3
Deposit thickness (mm)
Toil = 15C
2
Toil = 20C
Toil = 25C
Toil = 35C
0
0 10 20 30 40 50 60
(b) Time (h)
FIGURE 6.7
Comparison of deposit thickness as a function of oil temperature, Toil, between (a) the experi-
mental results and (b) the model prediction by the MWP. The oil flow rate, Qoil, was maintained
constant at 20 m3/h, and the coolant temperature Tcoolant, was maintained constant at 5C.
(Reprinted with permission from Huang, Z., Lu, Y., Hoffmann, R., Amundsen, L., & Fogler,
H.S. Energy & Fuels, 25, 51805188. Copyright 2011 American Chemical Society.)
116 Wax Deposition
3
Deposit thickness (mm)
Tcoolant = 5C
2
Tcoolant = 10C
1 Tcoolant = 15C
0
0 20 40 60 80 100 120
(a) Time (h)
3
Tcoolant = 5C
Tcoolant = 10C
Deposit thickness (mm)
2
Tcoolant = 15C
0
0 20 40 60 80 100 120
(b) Time (h)
FIGURE 6.8
Comparison of deposit thickness as a function of coolant temperature, Tcoolant, between (a) the
experimental results and (b) the model prediction by the MWP. The oil flow rate, Qoil, was
maintained constant at 5 m3/h, and the oil temperature Toil, was maintained constant at 20C
(Reprinted with permission from Huang, Z., Lu, Y., Hoffmann, R., Amundsen, L., & Fogler,
H.S. Energy & Fuels, 25, 51805188. Copyright 2011 American Chemical Society.)
TABLE 6.2
Summary of the Experimental Trends and That Predicted by the MWP
Operating Conditions Experimental Trend MWP
Increased flow rate Decreased thickness Decreased thickness
Increased oil temperature Decreased thickness Decreased thickness
Increased coolant temperature Decreased thickness Decreased thickness
Applying Wax Deposition Models to Flow Loop Experiments 117
It should be noted that the MWP was not used to exactly match the value
of the thickness of the wax deposit for each individual experiment. It is
believed that it is neither practical nor would it provide significant insight
to reveal the physics of wax deposition by tweaking one or two adjustable
parameters to match all the experimental data, given the degree of uncer-
tainties associated with wax deposition modeling that has been discussed
in Section 6.2.3. However, the MWP was able to predict all the experimen-
tal trends, and the importance of this achievement cannot be overstated:
First, the change in the oil flow rate and in the oil temperature is essen-
tial for the model scale-up for field pipelines, as the oil flow rate in these
laboratory experiments ranges from 5 to 25 m3/h; while the flow rate in the
field pipelines often ranges from 10,000 to 30,000 barrels per day (bpd) (i.e.,
66200 m3/h). While the oil temperature remains virtually constant in the
5.5 m long test section for the laboratory flow-loop experiments, the actual
oil temperature in the field pipelines can change from a value above 70C
down to close to the ocean floor temperature of about 4C. Consequently,
being able to capture the correct direction of change in the wax deposit
thickness during the change in these operating conditions is critical to the
development of reliable predictions to the field. To the authors knowledge,
we know of no other wax deposition model that has been reported to have
this success. Because of this success, more in-depth heat and mass transfer
analysis will be presented in Section 6.4 to help reveal the most influen-
tial parameters for wax deposition based on the studies of Huang, Lu et al.
(2011) and Lu et al. (2012).
force often refers to the temperature difference between the bulk oil and the
coolant. Wax deposition studies were often carried out with different ther-
mal driving forces. In the study by Jennings and Weispfennig (2005), a crude
oil from the Gulf of Mexico was used in a cold-finger wax deposition appa-
ratus. It was found that experiments with higher cold-finger temperatures
(smaller thermal driving force) yielded smaller amounts of deposited wax.
Creek et al. (1999) found similar results with a flow loop in which the deposit
thickness decreased when the coolant temperature increased (decreasing
thermal driving force). Consequently, it was often believed that the deposit
thickness would decrease when the thermal driving force decreases.
However, exceptions were, in fact, observed in other studies. In the cold-
finger study by Paso and Fogler (2004), it was observed that lower oil temper-
atures (a smaller thermal driving force) lead to thicker wax deposits. In the
study by Bidmus and Mehrotra (2009) using a mini flow loop, different trends
were observed: When the oil temperature is greater than the WAT (referred
to as hot flow), the amount of wax deposit decreases with increasing ther-
mal driving force. When the oil temperature is below the WAT (referred to as
cold flow), the amount of wax deposit decreases with decreasing thermal
driving force.
Given the many different experimental observations reported, it would
first appear that there is no established conclusion on the behavior of wax
deposition with varying temperatures. Indeed, with different oils and differ-
ent operating conditions, the observation in one set of experiments from one
oil sample might not necessarily be appropriate for another. However, the
fundamental heat and mass transfer equations for these experiments are the same or
rather similar, giving us some hope that the transport theories can be utilized
to identify a simple rule that can explain these different observations.
6.4.2Theoretical Analysis
In this section, an in-depth mass transfer analysis will be made to identify
the underlying physics that governs the behavior of wax deposition. To
achieve this goal, one of the classic approaches in transport phenomena
dimensionless analysiswill first be applied.
C 1 C
Vz = r( mass + Dwax ) (6.1)
z r r r
Applying Wax Deposition Models to Flow Loop Experiments 119
Centerline
Cinlet
Oil Concentration
profile Rinterface
Vz
Deposit Cinterface
deposit
Pipe wall
Coolant (water)
FIGURE 6.9
Schematic of the wax concentration profile at normal pipe flow conditions.
at z = 0, C = Cinlet
C
at r = 0, =0 (6.2)
r
at r = Rinterface , C = Cinterface
C Cinterface V z r dd2U
= , v= , = , = , Gz = (6.3)
Cinlet Cinterface U L rd (
L mass + Dwo )
1 4
v = (6.4)
Gz
at = 0, = 1
at = 0, = 0 (6.5)
at = 1, = 0
120 Wax Deposition
dC C
J A = Dwax,interface = Dwax,interface
dr interface r interface
In this case, the radial mass flux JA can be seen as a multiplication of two
parameters: The first one, noted as Jwax, appears to have the same unit of a
mass flux and is called the characteristic mass flux for wax deposition, while
the second parameter, the dimensionless concentration gradient, ,
interface
depends mainly on the oil flow rate.
Once we have Twall, we can accurately determine Cwall(eq) from the wax
solubility curve and Dwax,wall from the HaydukMinhas correlation (Hayduk
& Minhas,1982). It should be noted that in most of the lab-scale flow-loop
studies, no significant axial variation in the bulk concentration was observed
for the parameters due to the short length of the pipe in the flow-loop test
section. Therefore, Cinlet (eq) is referred to as Coil (eq) for subsequent quantita-
tive analyses.
TABLE 6.3
Comparison of the Parameters Involved in the Characteristic Mass Flux of Wax
Deposition, Jwax, among the Deposition Experiments with Different Toil while
Qoil and Tcoolant Remained Constant
Parameters Value
Toil (C) 15.3 20.3 25.3 35.4
Tcoolant (C) 5.0
Qoil (m3/h) 20.0
Twall (C) 9.5 12.0 14.7 20.5
Dwax,wall 1010 (m2/s) 2.11 2.49 2.89 3.64
Coil(eq) (wt%) 1.09 1.26 1.34 1.44
Cwall(eq) (wt%) 0.65 0.89 1.07 1.26
Coil(eq) Cwall(eq) (wt%) 0.45 0.36 0.27 0.17
Jwax 1010 (m/swt%) 19.25 16.20 13.18 10.12
Source: Huang, Z. et al., Energy & Fuels, 25, 51805188, 2011.
Applying Wax Deposition Models to Flow Loop Experiments 123
TABLE 6.4
Comparison of the Parameters Involved in the Characteristic Mass Flux of
Wax Deposition, Jwax, among the Deposition Experiments with Different
Tcoolant while Qoil and Toil Remained Constant
Parameters Value
Toil (C) 20.2
Tcoolant (C) 5.0 10.0 15.0
Qoil (m3/h) 5.0
Twall (C) 8.1 12.1 16.1
Dwax,wall 1010 (m2/s) 1.93 2.44 2.97
Coil(eq) (wt%) 1.26 1.26 1.26
Cwall(eq) (wt%) 0.48 0.89 1.13
Coil(eq) Cwall(eq) (wt%) 0.78 0.37 0.13
Jwax 1010 (m2/swt%) 31.64 16.54 7.92
Source: Huang, Z. et al., Energy & Fuels, 25, 51805188, 2011.
characteristic mass flux, Jwax, and its associated parameters for this change in
Tcoolant are shown in Table 6.4. It can be seen that the increase in Tcoolant leads
to the increase in Twall, which increased the values of only two parameters:
Dwo and Cwall(eq). Here, the concentration of wax in the bulk oil, Coil(eq), did
not change as the oil temperature was not varied, as in the previous case.
The increase in Dwax,wall tends to increase the Jwax. However, the increase in
Cwall(eq) from 0.48% to 1.13% caused a significant decrease in the mass trans-
fer driving force, [Coil(eq) Cwall(eq)], from 0.78% to 0.13%. This decrease in the
mass transfer driving force overcame the increase in the diffusivity, leading
to a decrease in the characteristic mass flux with increasing Tcoolant, as shown
in the last row of Table 6.4. This decrease in the characteristic mass flux
explains the reduced experimental deposit thickness, as shown in Figure 6.8.
In summary, the change in the characteristic mass flux, Jwax, in these exper-
iments with varying temperature conditions is shown to be consistent with
the change in the deposit thickness. It was further observed that the mass
transfer driving force, [Coil(eq) Cwall(eq)], is the most influential parameter
in determining the characteristic mass flux, and it should be ultimately the
indicator for the impact of temperature on wax deposition.
were found in this study: When Toil is above the WAT, the deposit thickness
decreases with increasing Toil. When Toil is below the WAT, the deposit thickness
increases with increasing Toil. While no mass transfer explanation was given in
the original study to help explain such observations, from our previous analy-
sis with the North Sea condensate, an investigation using the mass transfer
driving force, [Coil(eq) Cwall(eq)], as well as the characteristic mass flux for
wax deposition, Jwax, might generate some important insights. Consequently,
the same solvent, Norpar13 from Imperial Oil (Ontario, Canada), and the same
wax sample, Parowax, from Conros Corp. (Ontario, Canada), were purchased
(Huang, Lu et al., 2011). The measured wax precipitation curve is shown in
Figure 6.10. Comparing Figure 6.10 with Figure 6.3, one notes that the precipi-
tation curve for the model oilwax system is almost linear during the range of
experimental conditions, while the solubility curve for the North Sea conden-
sate exhibits more precipitation at lower temperatures.
The deposition experiments in the study of Bidmus and Mehrotra (2009)
are summarized in Figure 6.11, and the trend will now be analyzed with the
approach using the characteristic mass flux, Jwax. Using the same method of
analysis as for the North Sea condensate, the values of Jwax and its associated
parameters for the model waxoil experiments are summarized in Table 6.5.
It can be seen that the characteristic mass flux, Jwax, is again consistent with
the amount of deposit observed in Figure 6.11 in that they both increase with
increasing oil temperature until the oil temperature reaches the WAT. After
reaching the WAT, both Jwax and the amount of deposit decrease with increasing
temperature. The remarkable consistency in the relative magnitude of these two
parameters further confirms the effectiveness of the use of the characteristic
mass flux and the mass transfer driving force as predictors for the amount of
wax deposit on the wall.
3.0
2.5
Precipitated wax (wt%)
2.0
1.5
1.0
0.5
0.0
25 27 29 31 33 35
Temperature (C )
FIGURE 6.10
Solubility curve of wax for the model oil in the study of Huang, Lu et al. (2011). (Reprinted with
permission from Huang, Z., Lu, Y., Hoffmann, R., Amundsen, L., & Fogler, H.S., Energy & Fuels,
25, 51805188. Copyright 2011 American Chemical Society.)
Applying Wax Deposition Models to Flow Loop Experiments 125
2
1.8 Tcoolant = 25C WAT
1.6
wall surface area (kg/m2)
Deposit mass per inner
1.4
1.2
1
0.8
0.6
0.4
0.2
0
26.5 30 32.5 35 38.5
Toil (C)
FIGURE 6.11
Comparison of the amount of deposit among experiments with different oil temperatures for the
study of Bidmus and Mehrotra (2009). (From Bidmus, H.O., & Mehrotra, A.K., Energy & Fuels, 23,
31843194, 2009. Reprinted with permission from Huang, Z., Lu, Y., Hoffmann, R., Amundsen,
L., & Fogler, H. S., Energy & Fuels, 25, 51805188. Copyright 2011 American Chemical Society.)
TABLE 6.5
Comparison of the Parameters for Jwax among Deposition Experiments with
Different Toil while Qoil and Tcoolant Are Maintained Constant in the Study of
Bidmus and Mehrotra (2009) Using a Model OilWax System
Parameters Value
Toil (C) 26.5 29.0 33.0 35.0 38.5
Tcoolant (C) 25.0
Qoil (m3/h) 0.4
Twall (C) 25.4 25.9 27.1 28.0 29.6
Dwax,wall 1010 (m2/s) 2.81 2.94 3.27 3.53 4.07
Coil(eq) (wt%) 4.19 4.80 5.74 6.00 6.00
Cwall(eq) (wt%) 3.94 4.05 4.34 4.55 4.93
Coil(eq) Cwall(eq) (wt%) 0.25 0.75 1.40 1.45 1.07
Jwax 1010 (m2/swt%) 27.66 86.81 180.24 201.52 171.28
Source: Huang, Z. et al., Energy & Fuels, 25, 51805188, 2011.
TABLE 6.6
Comparison of the Changes in Toil, Twall, Coil(eq), and Cwall(eq) for the
Experiments with the North Sea Oil A and the Model OilWax
Mixturewith Varying Toil
North Sea Condensate Model OilWax Mixture
Change in Toil (C) 15.335.4 26.535.0
Toil (C) 20.1 8.5
Twall (C) 11.0 2.6
Toil > Twall Yes Yes
Coil(eq) (wt%) 0.35 1.81
Cwall(eq) (wt%) 0.61 0.61
Coil(eq) > Cwall(eq) No Yes
the increases in Toil for both studies. The major difference is that for the model
oilwax mixture, the gradient of its solubility curve (Figure 6.10) is virtually
constant (the solubility curve is close to a straight line), so that the changes
in Coil(eq), and Cwall(eq) simply reflected the changes in Toil and Twall. These
changes eventually lead to the increase in the mass transfer driving force,
[Coil(eq) Cwall(eq)], and in the mass flux, Jwax, for the model mixture when Toil
increased. However, the gradient of the solubility curve for the North Sea con-
densate decreased with increasing temperature because the solubility curve is
concave, as shown in Figure 6.3. This concave shape of the solubility curve for
the North Sea condensate results in the change in [Coil(eq) Cwall(eq)] to deviate
from the change in [Toil Twall]. This finding further demonstrates the advan-
tage of using the mass transfer driving force in comparison to the thermal driv-
ing force as it includes the impact of the solubility curve on wax deposition.
12.0
10.0
Weight percent (wt%)
8.0
6.0
4.0
2.0
0.0
20 30 40 50 60 70 80
(a) Carbon number
12.0
10.0
Weight percent (wt%)
8.0
6.0
4.0
2.0
0.0
20 30 40 50 60 70 80
(b) Carbon number
FIGURE 6.12
CND of the heavy components in (a) the model oilwax mixture and (b) the North Sea con-
densate. (Reprinted with permission from Huang, Z., Lu, Y., Hoffmann, R., Amundsen, L., &
Fogler, H. S. Energy & Fuels, 25, 51805188. Copyright 2011 American Chemical Society.)
concave shape of the solubility curve of the North Sea condensate at low
temperatures and eventually explains the contradictory observed trends of
growth of deposit thickness between Figures 6.7 and 6.11.
TABLE 6.7
Summary of Two Sets of Operating Conditions:
Aand B Resulting in Similar Deposit CND
Condition A Condition B
Qoil (m3/h) 5.00 21.00
Toil (C) 20.17 15.24
Tcoolant (C) 5.00 10.00
Twall (C) 8.39 12.49
Applying Wax Deposition Models to Flow Loop Experiments 129
50%
40% Experiment A
Ci/Ci (%)
30% Experiment B
20%
10%
0%
(a) 20 25 30 35 40 45
50%
Ci,deposit/Ci,deposit (%)
40%
30%
20%
10%
0%
20 25 30 35 40 45
(b) Carbon number
FIGURE 6.13
(a) Comparison of the concentration driving force distribution between conditions A and B and
(b) comparison of the deposit CND generated by two sets of operating conditions. (Reprinted
with permission from Zheng, S., Zhang, F., Huang, Z., & Fogler H. S., Energy & Fuels, 27, 7379
7388. Copyright 2013 American Chemical Society.)
The deposit CND under the particular set of the operating conditions
should resemble the concentration driving force distribution.
6.6Summary
In this chapter, the performance of a model developed at the University of
Michigan, known as the Michigan Wax Predictor (MWP), is highlighted
130 Wax Deposition
Based on this consistency, the transport equations in the MWP are carefully
examined to identify the reasons for such agreements. With the approach of a
dedimensionalized analysis, a representative parameter to highlight the degree
of wax deposition is proposed: the characteristic mass flux for wax deposition:
Within this parameter, the effect of temperature can be separated into two
categories: The first is on the transport properties of the oil (to affect the wax
diffusion coefficient at the wall temperature, Dwax,wall), and the second is on the
thermodynamic properties of the oil (to affect the wax equilibrium concentra-
tions at the bulk oil and at the wall). It was found that in most of the wax depo-
sition experiments, the effect of temperature on the thermodynamic properties
often outweighs the effect on the transport properties, so that the change in the
deposit thickness is often found to be consistent with the change in the concen-
tration difference, [Coil(eq) Cwall(eq)]. Such difference is considered to be the
mass transfer driving force for deposition, and it is an important indicator on
how the deposition behavior varies when the operating conditions are altered.
This methodology was further validated with a waxoil model mixture in a
different deposition flow loop. The consistency of this methodology for differ-
ent oils and different experimental apparatuses has strategic importance for
model scale-up for field predictions. In fact, in Chapter 7, it will be seen that
the MWP is found to be able to achieve satisfactory agreement with field data.
In addition to the comparison of the deposit thickness, an analysis is made
based on the comparison of the composition of the wax deposit between
experimental characterization and model predictions. In this case, wax
deposition modeling by the MWP is combined with a thermodynamic model
(Coutinho & Stenby, 1996; Coutinho & Ruffier-Mray, 1997; Coutinho, 1998;
Coutinho & Daridon, 2001; Coutinho, Edmonds, Moorwood, Szczepanski,
& Zhang, 2006) to resolve the mass transfer of alkanes with different chain
length. It was discovered that the evolution of deposit composition is also
dominated by the distribution of concentration driving forces [Coil(eq)
Cwall(eq)] among alkanes with different chain length.
7
Applying Wax Deposition
Models for Field Predictions
7.1Introduction
7.1.1Wax Control Strategies for the Field
In this chapter, thermodynamic modeling and transport modeling are used
to help develop remediate strategy characteristics for wax deposition in the
field. A general methodology of wax deposition evaluation and control is
highlighted in Figure 7.1.
The first step to assess the possibility of wax deposition is to determine the
wax appearance temperature (WAT). The WAT is usually compared with the
minimum wall temperature of the pipeline during operations. The wall tem-
perature of the pipeline can be calculated by most flow simulators that have
been used extensively for pipeline designs for the past decades. Depending
on the values of the WAT and the minimum wall temperature, the following
possible scenarios can occur:
Scenario 1: Wax deposition would probably not pose a risk for nor-
mal operations if the WAT is significantly lower than the mini-
mum wall temperature. However, one should be on the lookout for
situations where the flow rate of the fluids decreases (production
turndown), causing the pipeline wall temperature to fall below
the WAT.
Scenario 2: If the minimum wall temperature of the pipeline is close
to or only slightly below the WAT, there can be a risk of wax depo-
sition in the pipeline. In this case, the problem of wax deposition
might be addressed by increasing the insulation or installing heat-
ing facilities (for example, direct electrical heating or pipe-in-pipe
heating) over the section of the pipe where the wall temperature
might drop to below the WAT. For short pipelines, insulation and
heating methods might still be economically viable even if the mini-
mum wall temperature is significantly below the WAT.
131
132 Wax Deposition
Thermodynamic Deposition
characterizations experiments and
(Chapters 2 and 3) model benchmarking
(Chapters 4 through 6)
FIGURE 7.1
Summary of general wax deposition screening and control methodologies.
It should be noted that wax concerns in the subject of flow assurance are
not limited to deposition during production. During pipeline shutdowns,
the flow is stopped, and the fluid is being cooled to the ambient temperature.
In this case, if the cool down time is too long, the waxy components can also
precipitate from the oil and cause gelation to occur in the pipeline. Wax gela-
tion normally occurs during static conditions (pipeline shutdown) in which
the entire fluid (from the wall to the pipe centerline) can be gelled, while wax
deposition mainly occurs on the pipe wall during flow conditions. Although
Applying Wax Deposition Models for Field Predictions 133
Restart experiments
Rheological
and model
characterizations
benchmarking
FIGURE 7.2
Summary of general wax gelation screening and control methodologies.
not discussed in detail in this book, the control strategies of wax gelation are
analogous to those for wax deposition, as summarized in Figure7.2.
Subsea
PLEM
Gas16 in. 23 km
Oil12 in. 23 km
Central
processing
platform
FPSO
Export oil rate = 55,000 bpd
Export temp = 170F
FIGURE 7.3
Field layout in the study of Singh et al. (2011). (From Singh, A. et al., Flow assurance: Validation
of wax deposition models using field data from a subsea pipeline. In Offshore Technology
Conference [pp. 119]. Houston, TX: Offshore Technology Conference, 2011.)
TABLE 7.1
Summary of Operating Conditions in the Study of Singh et al. (2011)
Pipeline diameter 12 in.
Pipeline length 23 km
External heat transfer coefficient 22 W/m2/K
(insulation included)
Oil flow rate 300,000 bpd
Inlet temperature 75C
Outlet pressure 350 psig
Predeposition pressure drop 200 psi
Pigging frequency Weekly
Source: Singh A. et al., Flow assurance: Validation of wax deposition models
using field data from a subsea pipeline. In Offshore Technology Conference
(pp. 119). Houston, TX: Offshore Technology Conference, 2011.
Applying Wax Deposition Models for Field Predictions 135
10
All hydrocarbons
n-Alkanes
1
Wt%
0.1
0.01
0 10 20 30 40 50 60
Carbon number
FIGURE 7.4
Composition of the oil and its n-paraffin distribution. (From Singh, A. et al., Flow assurance:
Validation of wax deposition models using field data from a subsea pipeline. In Offshore
Technology Conference [pp. 119]. Houston, TX: Offshore Technology Conference, 2011.)
the WAT of 58C. The composition and the n-paraffin distribution of the oil
are shown in Figure 7.4. Based on this analysis, the oil has 17% n-C19+, which
is generally known as the waxy components. Other detailed properties of
the oil can be found in the original study of Singh et al. (2011).
20
10
Coutinho model
(untuned prediction)
Coutinho model
(tuned to WAT)
0
0 10 20 30 40 50 60 70
Temperature (C)
FIGURE 7.5
Comparison of the wax precipitation curves from different thermodynamic models.
80
70 Bulk oil
temperature
Temperature (C)
60
WAT
50
40
Inner wall Actual outlet
temperature temperature
30
~3 km
20
0 5000 10,000 15,000 20,000 25,000
Length (m)
FIGURE 7.6
Prediction of the temperature profiles in the pipeline.
any wax has deposited. In fact, it is very important to compare the predicted
bulk temperature of the fluid with field measurements if possible to confirm
the validity of the heat transfer calculation by the wax model. It is seen that
due to the insulation (external heat transfer coefficient of 22 W/m/K) of the
pipeline, the bulk and the wall temperature profiles are quite similar. The
temperature profiles are consistent with field measurements, which have
recorded an outlet temperature of 27C29.5C. The agreement between the
calculated temperature and the measured value is critical to develop reliable
deposition predictions. The location where the wall temperature drops to
below the WAT is around 3 km downstream of the inlet, which indicates that
even with insulation, approximately 20 km of the pipeline is subject to the
risk of wax deposition.
After gaining insights from the temperature profile along the pipe, we are
ready to investigate the growth of the wax deposit thickness. The visualiza-
tion of deposit thickness from field predictions is slightly different from that
during the modeling of laboratory deposition experiments: In previous dis-
cussions in Chapter 6 about the laboratory experiments, the evolution of axial
averaged deposition thickness (x-axis) as a function of time (y-axis) was shown.
Because for the laboratory experiments, the length of the pipe is much shorter,
the deposit profile throughout the pipe is expected to be uniform. This axial
averaged deposit thickness was compared with the experimental measure-
ments using the pressure drop method, i.e., the pressure-drop of the pipe was
used to determine the effective radius of the pipe. The difference of the effec-
tive radius and the pipe inner radius is the deposit thickness. This method
of plotting was also used in many other laboratory wax deposition studies
(Hernandez, 2002; Hoffmann & Amundsen, 2010; Huang, Lu, Hoffmann,
Amundsen, & Fogler, 2011; Lund, 1998; Singh et al., 2000; Venkatesan, 2004).
138 Wax Deposition
30
7 days
Deposit thickness (mm) 25
20
15
2 days
10
24 h
5
14 h
0
0 5000 10,000 15,000 20,000 25,000
Length (m)
FIGURE 7.7
Predictions of the deposit thickness along the pipeline at different times by the MWP.
Lpipe, i oilU 2
Ppipe, i = fDarcy (7.2)
deffective 2
Applying Wax Deposition Models for Field Predictions 139
Ppipe = P
i
pipe, i (7.3)
450
400
Pipeline pressure drop (psi)
350
300
250
200
150
100
50
0
0 2 4 6 8
Time (days)
FIGURE 7.8
Comparison of pipeline pressure drop during the course of wax deposition between MWP
predictions and field measurements.
140 Wax Deposition
30
7 days
Deposit thickness (mm) 25
20
15
10
5
14 h
0
0 5000 10,000 15,000 20,000 25,000
Length (m)
FIGURE 7.9
Prediction of deposit thickness along the pipe after 7 days.
For this study, it was eventually found that pigging could be performed
every 7 days without causing any operational concerns (Singh et al., 2011).
After 7 days, the axial distribution of the deposit thickness predicted by the
MWP is 27 mm, as previously shown in Figure 7.7. In other words, in this
study, if the pigging frequency is determined by using the deposition model
of MWP, the threshold for the maximum thickness allowed could be as much
as 27 mm. In this case, applying the general value of 4 mm seems conserva-
tive for this case, and the pigging frequency would change to every 14 h,
which in hindsight seems conservative (Figure 7.9).
Top side
FIGURE 7.10
Schematic of the production pipeline for the study of Noville and Naveira (2012).
TABLE 7.2
Summary of the Operating Conditions in the Study of
Noville and Naveira (2012)
Pipeline diameter 6 in.
Pipeline length 5 km
Oil flow rate 2617 bpd
Gas oil ratio (GOR) 618
Lift gas flow rate 5.8 mmscfd
Reservoir temperature 78C
Reservoir pressure 2217 psig
Topside pressure 155 psig
Pigging frequency 14 days
Source: Noville, I. & Naveira, L., Comparison between real
field data and the results of wax deposition simula-
tion. In SPE Latin America and Caribbean Petroleum
Engineering Conference (pp. 112). Mexico City, Mexico:
Society of Petroleum Engineers, 2012.
for wax deposition modeling. The detailed descriptions for the wax ther-
modynamic model in PVTsim can be found in Chapter 3, Section 3.3. The
detailed transport equations for the wax deposition model in OLGA can be
found in Chapter 4, Sections 4.2 and 4.3.
The next step is to import this information of the wax precipitation curve
to a wax deposition model. However, here emerges a problem about indus-
trial practicality: The wax deposition model in the industrial flow simula-
tor OLGA only accepts the input of a wax precipitation curve in a format
that is developed by selected industrial thermodynamic models that it sup-
ports. In this study, the program of PVTsim is used for this purpose. The
program of PVTsim generates a wax table that consists of the amount of
wax precipitated at different temperatures and pressures. This wax table
could be directly used by OLGA for wax deposition simulations. In this case,
the n-paraffin distribution is adjusted in a thermodynamic model order to
tune the predicted wax precipitation curve to match the DSC results. The
thermodynamic model used is based on the study by Pedersen, Thomassen,
and Fredenslund (1985). However, we can see that this information is not as
important here as the thermodynamic model was only used as a dummy.
The real source of the precipitation curve is the DSC experiment. The DSC
measurements as well as the thermodynamic model fitting are shown in
Figure 7.11. It is seen that the predicted precipitation curve matches the
experimental measurements from 0C to 17C. However, at the range of
17C35C, there seems to be a very small amount of wax precipitation that
was underpredicted by the thermodynamic model.
This method of thermodynamic characterization indicates that the method
of DSC was preferred in this study to study the wax precipitation curve,
which is different from the previous example (Example 1) where thermo-
dynamic prediction was carried out based on the n-paraffin distribution
measured by the HTGC. A comparison of these two strategies is high-
lighted in Table 7.3.
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
0 5 10 15 20 25 30 35 40 45 50 55 60
Temperature (C)
FIGURE 7.11
Experimental measurement of the wax precipitation curve and the thermodynamic fitting in
thestudy of Noville and Naveira (2012). (From Noville, I. & Naveira, L., Comparison between real
field data and the results of wax deposition simulation. In SPE Latin America and Caribbean Petroleum
Engineering Conference [pp. 112]. Mexico City, Mexico: Society of Petroleum Engineers, 2012.)
Applying Wax Deposition Models for Field Predictions 143
TABLE 7.3
Comparison of the Two Thermodynamic Characterization Strategies Used in
Examples 1 and 2
Example 1 Example 2
Experiments n-Paraffin distribution using HTGC Entire wax precipitation curve
carried out WAT using CPM using DSC
Purpose of using To predict a wax precipitation curve to To generate a wax precipitation
thermodynamic be used by a deposition model curve same as that measured
modeling by DSC for a deposition model
Parameters to A different thermodynamic model offers The n-paraffin distribution
adjust during different tuning options, for example,
tuning the the heat of crystallization and the
thermodynamic detection limit of the WAT measuring
model device. A complete list for various
models was discussed in Table 3.2
Parameters to match WAT measured using CPM Entire wax precipitation curve
during modeling measured by DSC
1600
1550
Field measurements OLGA wax prediction
pessimistic (Dwax 100)
Downhole pressure (psig)
1500
1450
1400
1350
OLGA wax prediction
optimistic (Dwax 22)
1300
FIGURE 7.12
Comparison of the increase in the downhole pressure between OLGA wax predictions and
field measurements. (From Noville, I. & Naveira, L., Comparison between real field data and
theresults of wax deposition simulation. In SPE Latin America and Caribbean Petroleum Engineering
Conference [pp. 112]. Mexico City, Mexico: Society of Petroleum Engineers, 2012.)
25
20 Well B
15
Well A
10
5 Well C
0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
Flowline length (m)
FIGURE 7.13
Predictions of deposit thickness using the OLGA wax based on the optimistic scenario (Dwax
multiplied by 22). (From Noville, I. & Naveira, L., Comparison between real field data and
theresults of wax deposition simulation. In SPE Latin America and Caribbean Petroleum Engineering
Conference [pp. 112]. Mexico City, Mexico: Society of Petroleum Engineers, 2012.)
Applying Wax Deposition Models for Field Predictions 145
The time for the predicted maximum deposit thickness along the
pipe to reach a certain threshold
The time for the predicted maximum mass of the wax deposit in the
entire pipeline to reach a certain threshold
These two thresholds are designed to ensure that the pig does not become
trapped in the pipeline during pigging operations. The final pigging fre-
quency is based on the maximum number of days that satisfies both of these
criteria. According to the original study, the threshold for the maximum
predicted deposit thickness is reported to be 11 mm, and that for the maxi-
mum mass of wax deposit acceptable is 750 kg. It was found that the time for
the first criterion is 7 days, and the time for the second criterion is 12 days.
Consequently, the proposed pigging frequency was 710 days. The field was
operating with this pigging frequency initially, and the pigging frequency
eventually increased to once every 14 days.
7.4Summary
7.4.1Wax Thermodynamic Characterizations
One of the differences in the two examples in this chapter is the method
used to characterize the wax precipitation curve. The first study by Singh
et al. (2011) is based on the n-paraffin distribution measured by the HTGC.
The wax precipitation curve is obtained from the prediction from a thermo-
dynamic model, which is tuned to match the WAT, which was measured
experimentally. In the second study by Noville and Naveira (2012), the wax
precipitation curve is directly measured, and the use of a thermodynamic
model is simply to generate the output of the precipitation curve for the input
of a deposition model. In this case, the n-paraffin distribution is adjusted to
fulfill the purpose. These two different methodologies reflect different lev-
els of confidence on the wax thermodynamic characterization techniques, as
summarized and shown in Table 7.4.
It should be noted that uncertainties are associated with the technolo-
gies used in these two approaches. For example, the n-paraffin distribu-
tion from HTGC is obtained by integration of the area under the peaks in
a chromatogram, which can depend on the integration method (either the
baseline method or the valley-to-valley method). The WAT detected by the
CPM can be sensitive to its detection limit and the DSC measurements have
been reported to have difficulties to differentiate between the noise and the
signal for low wax content oils (such as the one for Example 2). A detail of the
146 Wax Deposition
TABLE 7.4
Comparison of the Two Wax Thermodynamic Characterization Methods Used in
the Two Different Examples in This Chapter
n-Paraffin Parameter to Match Measuring Technologies That
Examples Distribution during Tuning the Method Relies On
1 Measured WAT HTGC, CPM
experimentally Wax thermodynamic models
2 Adjustable Wax precipitation DSC
parameters curve
7.5Future Outlook
Based on the above field applications, it can be seen that despite several
decades of wax deposition research, today, the application of wax modeling
for field design still faces significant uncertainties. Consequently, industrial
Applying Wax Deposition Models for Field Predictions 147
increase in the pressure drop of the pipeline during the course of wax depo-
sition, for Example 2, the agreement could not be reached without drastic
tuning. The main reason for this discrepancy is on the occurrence of mul-
tiphase oilgas flow. In addition to oilgas flow, oilwater two-phase flows
and oilgaswater three-phase flows can also occur in field pipelines. Very
few experimental studies attempt to understand the effects of the interplay
between phases on wax deposition from fundamental perspectives of view.
Moreover, the oilgaswater three-phase wax deposition study has not yet
been conducted. Lack of experimental characterization poses difficulties
on the development of fundamental wax deposition models for multiphase
flows. A variety of research directions are being undertaken at the University
of Michigan in order to understand and model wax deposition in multiphase
flows from first principles. As an example, in the case of wateroil flows,
emulsified water can greatly change the rheological behavior of waxy crude
oil, while it is still unclear how the rheology of waxy oil emulsion should be
considered in any wax deposition models. The existence of a second phase
in addition to the oil phase may also change the interaction between the flow
and the pipe wall. In the case of wateroil two-phase flows, the wax deposi-
tion characteristics might also depend on the hydrophilicity and hydropho-
bicity of the pipe wall surface. Research and developments to understand the
multiple perspectives of wax deposition characteristics in multiphase flows
are essential to improving the efficiency in the numerical modeling.
Bibliography
Agrawal, K. M., Khan, H. U., Surianarayanan, M., & Joshi, G. C. (1990). Wax deposi-
tion of Bombay high crude oil under flowing conditions. Fuel, 69, 794796.
Akbarzadeh, K., & Zougari, M. (2008). Introduction to a novel approach for mod-
eling wax deposition in fluid flows. 1. TaylorCouette system. Industrial and
Engineering Chemistry Research, 47, 953963.
Alana, J. D. (2003). Investigation of heavy oil single-phase paraffin deposition characteristics
(M.S. thesis). University of Tulsa, Tulsa, OK.
Alcazar-Vara, L. A., & Buenrostro-Gonzalez, E. (2011). Characterization of the wax
precipitation in Mexican crude oils. Fuel Processing Technology, 92, 23662374.
Alghanduri, L. M., Elgarni, M. M., Daridon, J.-L., & Coutinho, J. A. P. (2010).
Characterization of Libyan waxy crude oils. Energy & Fuels, 24, 31013107.
Apte, M. S., Matzain, A., Zhang, H.-Q., Volk, M., Brill, J., & Creek, J. L. (2001).
Investigation of paraffin deposition during multiphase flow in pipelines and
wellboresPart 2Modeling. Journal of Energy Resources Technology, 123,
150157.
Ashford, J. D., Blount, C. G., Marcou, J. A., & Ralph, J. M. (1990). Annular packer
fluids for paraffin control: Model study and successful field application. SPE
Production Engineering, 5, 351355.
ASTM International. (2008). Standard test method for temperature calibration of differential
scanning calorimeters and differential thermal analyzers (ASTM Standard E967-08).
doi:10.1520/E0967-08.
ASTM International. (2010). Standard test method for measurement of transition tempera-
tures of petroleum waxes by differential scanning calorimetry (DSC) (ASTM Standard
D4419-90). doi:10.1520/D4419-90R10.
ASTM International. (2011). Standard test method for cloud point of petroleum products
(ASTM Standard D2500-11). doi:10.1520/D2500-11.
Bai, C., & Zhang, J. (2013a). Effect of carbon number distribution of wax on the
yield stress of waxy oil gels. Industrial and Engineering Chemistry Research, 52,
27322739.
Bai, C., & Zhang, J. (2013b). Thermal, macroscopic, and microscopic characteristics of
wax deposits in field pipelines. Energy & Fuels, 27, 752759.
Bai, Y., & Bai, Q. (2012). Subsea engineering handbook (1st ed.). Houston, TX: Gulf
Professional Publishing.
Barry, E. G. (1971). Pumping non-Newtonian waxy crude oils. Journal of the Institute
of Petroleum, 57, 7485.
Bendiksen, K. H., Maines, D., Moe, R., & Nuland, S. (1991). The dynamic two-fluid
model OLGA: Theory and application. SPE Production Engineering, 6, 171180.
Bern, P. A., Withers, V. R., & Cairns, R. J. R. (1980). Wax deposition in crude oil pipe-
lines. In European Offshore Petroleum Conference & Exhibition (p. 571). London:
Earls Court.
Berne-Allen, A., Jr., & Work, L. T. (1938). Solubility of refined paraffin waxes in petro-
leum fractions. Industrial and Engineering Chemistry, 30, 806812.
149
150 Bibliography
Bhat, N. V., & Mehrotra, A. K. (2004). Measurement and prediction of the phase
behavior of waxsolvent mixtures: Significance of the wax disappearance tem-
perature. Industrial and Engineering Chemistry Research, 43, 34513461.
Bidmus, H. O., & Mehrotra, A. K. (2009). Solids deposition during cold flow of
waxsolvent mixtures in a flow-loop apparatus with heat transfer. Energy &
Fuels, 23, 31843194.
Bilderback, C. A., & McDougall, L. A. (1969). Complete paraffin control in petroleum
production. Journal of Petroleum Technology, 21, 11511156.
Bokin, E., Febrianti, F., Khabibullin, E., & Perez, C. E. S. (2010). Flow assurance and sour
gas in natural gas production.
Brill, J. P. (1987). Multiphase flow in wells. Journal of Petroleum Technology, 39, 1521.
Brown, T. S., Niesen, V. G., & Erickson, D. D. (1993). Measurement and prediction of
the kinetics of paraffin deposition. ProceedingsSPE Annual Technical Conference
and Exhibition, 353368.
Bruno, A., Sarica, C., Chen, H., & Volk, M. (2008). Paraffin deposition during the flow
of water-in-oil and oil-in-water dispersions in pipes. In ACTE 2008 Proceedings:
Proceedings of SPE Annual Technical Conference and Exhibition. Denver, CO: Society
of Petroleum Engineers, SPE 114747, September 2124, 2008.
Burger, E. D., Perkins, T. K., & Striegler, J. H. (1981). Studies of wax deposition in the
Trans Alaska Pipeline. Journal of Petroleum Technology, 33, 10751086.
Cazaux, G., Barre, L., & Brucy, F. (1998). Waxy crude cold start: Assessment through
gel structural properties. In SPE Annual Technical Conference and Exhibition. New
Orleans, LA: Society of Petroleum Engineers, SPE 49213.
Chilton, T. H., & Colburn, A. P. (1934). Mass transfer (absorption) coefficients pre-
diction from data on heat transfer and fluid friction. Industrial and Engineering
Chemistry, 26, 11831187.
Chueh, P. L., & Prausnitz, J. M. (1967). Vaporliquid equilibria at high pressures:
Calculation of partial molar volumes in nonpolar liquid mixtures. AIChE
Journal, 13, 10991107.
Claudy, P., Letoffe, J. M., Neff, B., & Damin, B. (1986). Diesel fuels: Determination of
onset crystallization temperature, pour point and filter plugging point by dif-
ferential scanning calorimetry. Correlation with standard test methods. Fuel, 65,
861864.
Cleaver, J., & Yates, B. (1973). Mechanism of detachment of colloidal particles from
a flat substrate in a turbulent flow. Journal of Colloid and Interface Science, 44,
464474.
Coto, B., Coutinho, J. A. P., Martos, C., Robustillo, M. D., Espada, J. J., & Pea, J. L. (2011).
Assessment and improvement of n-paraffin distribution obtained by HTGC to
predict accurately crude oil cold properties. Energy & Fuels, 25, 11531160.
Coto, B., Martos, C., Espada, J. J., Robustillo, M. D., Merino-Garca, D., & Pe, L. (2011).
A new DSC-based method to determine the wax porosity of mixtures precipi-
tated from crude oils. Energy & Fuels, 25, 17071713.
Coto, B., Martos, C., Espada, J. J., Robustillo, M. D., & Pea, J. L. (2010). Analysis
of paraffin precipitation from petroleum mixtures by means of DSC: Iterative
procedure considering solidliquid equilibrium equations. Fuel, 89, 10871094.
Coto, B., Martos, C., Espada, J. J., Robustillo, M. D., Pea, J. L., Coutinho, J. A. P.,&
Pe, L. (2011). Study of new methods to obtain the n-paraffin distribution of
crude oils and its application to flow assurance. Energy & Fuels, 25, 487492.
Bibliography 151
Coto, B., Martos, C., Pea, J. L., Espada, J. J., & Robustillo, M. D. (2008). A new method
for the determination of wax precipitation from non-diluted crude oils by frac-
tional precipitation. Fuel, 87, 20902094.
Coutinho, J. A. P. (1998). Predictive UNIQUAC: A new model for the description
of multiphase solidliquid equilibria in complex hydrocarbon. Industrial and
Engineering Chemistry Research, 37, 48704875.
Coutinho, J. A. P., & Daridon, J. L. (2001). Low-pressure modeling of wax formation
in crude oils. Energy & Fuels, 15, 14541460.
Coutinho, J. A. P., & Daridon, J.-L. (2005). The limitations of the cloud point mea-
surement techniques and the influence of the oil composition on its detection.
Petroleum Science and Technology, 23, 11131128.
Coutinho, J. A. P., Edmonds, B., Moorwood, T., Szczepanski, R., & Zhang, X. (2006).
Reliable wax predictions for flow assurance. Energy & Fuels, 20, 10811088.
Coutinho, J. A. P., & Ruffier-Mray, V. (1997). Experimental measurements and ther-
modynamic modeling of paraffinic wax formation in undercooled solutions.
Industrial and Engineering Chemistry Research, 36, 49774983.
Coutinho, J. A. P., & Stenby, E. H. (1996). Predictive local composition models for
solid/liquid equilibrium in n-alkane systems: Wilson equation for multicompo-
nent systems Industrial and Engineering Chemistry Research, 35, 918925.
Coutinho, P. (1999). Predictive local composition models: NRTL and UNIQUAC and
their application to model solidliquid equilibrium of n-alkanes. Fluid Phase
Equilibria, 158, 447457.
Couto, G. (2004). Investigation of two-phase oil-water paraffin deposition (M.S. thesis).
University of Tulsa, Tulsa, OK.
Creek, J., Lund, H. J., Brill, J. P., & Volk, M. (1999). Wax deposition in single phase
flow. Fluid Phase Equilibria, 158160, 801811.
Crochet, M. J., Davies, A. R., & Walters, K. (1991). Numerical simulation of non-Newtonian
flow (3rd ed.). New York: Elsevier B.V.
Daridon, J., Coutinho, J. A. P., & Montel, F. (2001). Solid-liquid-vapor phase bound-
ary of a North Sea waxy crude: Measurement and modeling. Energy & Fuels, 15,
730735.
Dauphin, C., Daridon, J., Coutinho, J., Baylre, P., & Potin-Gautier, M. (1999). Wax con-
tent measurements in partially frozen paraffinic systems. Fluid Phase Equilibria,
161, 135151.
Deen, W. M. (1998). Analysis of transport phenomena (2nd ed.). Oxford University Press.
DeepStar. (2011). DeepStar 10203 wax prediction confidence and pigging risk.
Retrieved from http://www.deepstar.org/PhaseXProjects/.
de Oliveira, M. C. K., Teixeira, A., Vieira, L. C., de Carvalho, R. M., de Carvalho,
A. B.M., & do Couto, B. C. (2012). Flow assurance study for waxy crude oils.
Energy & Fuels, 26, 26882695.
Ding, J., Zhang, J., Li, H., Zhang, F., & Yang, X. (2006). Flow behavior of Daqing waxy
crude oil under simulated pipelining conditions. Energy & Fuels, 20, 25312536.
Dirand, M., Bouroukba, M., Briard, A.-J., Chevallier, V., Petitjean, D., & Corriou, J.-P.
(2002). Temperatures and enthalpies of (solid + solid) and (solid + liquid) transi-
tions of n-alkanes. The Journal of Chemical Thermodynamics, 34, 12551277.
Dirand, M., Chevallier, V., & Provost, E. (1998). Multicomponent paraffin waxes
and petroleum solid deposits: Structural and thermodynamic state. Fuel, 77,
12531260.
152 Bibliography
Dittus, F., & Boelter, L. (1985). Heat transfer in automobile radiators of the tubular
type. International Communications in Heat and Mass Transfer, 12, 322.
Eaton, P. E., & Weeter, G. Y. (1976). Paraffin deposition in flow lines. In 16th National
Heat Transfer Conference. St. Louis, MO.
Economic analysis methodology for the 5-year OCS Oil and Gas Leasing Program for 2012
2017. (2011). Washington, DC.
Edmonds, B., Moorwood, T., Szczepanski, R., & Zhang, X. (2007). Simulating wax
deposition in pipelines for flow assurance. Energy & Fuels, 22, 729741.
Elliott, J. R., & Lira, C. T. (2012). Introductory chemical engineering thermodynamics.
Upper Saddle River, NJ: Prentice Hall.
Elsharkawy, A. M., Al-Sahhaf, T. A., & Fahim, M. A. (2000). Wax deposition from
Middle East crudes. Fuel, 79, 10471055.
Erickson, D. D., Niesen, V. G., & Brown, T. S. (1993). Thermodynamic measurement and
prediction of paraffin precipitation in crude oil. In SPE Annual Technical Conference
and Exhibition (pp. 353368). Houston, TX: Society of Petroleum Engineers.
Esbensen, K. H., Halstensen, M., Tnnesen Lied, T., Svalestuen, J., de Silva, S., & Hope,
B. (1998). Acoustic chemometricsFrom noise to information. Chemometrics and
Intelligent Laboratory Systems, 44, 6176.
Garcia, M. C. (2001). Paraffin deposition in oil production. In SPE International
Symposium on Oilfield Chemistry (pp. 17). Houston, TX: Society of Petroleum
Engineers.
Garcia, M., Lopez, F., & Nino, Y. (1995). Characterization of near-bed coherent struc-
tures in turbulent open channel flow using synchronized high-speed video and
hot-film measurements. Experiments in Fluids, 19, 1628.
Gluyas, J. G., & Underhill, J. R. (2003). The Staffa Field, Block 3/8b, UK North Sea.
Geological Society, London, Memoirs, 20, 327333.
Gnielinski, V. (1976). New equations for heat and mass transfer in turbulent pipe and
channel flow. International Journal of Chemical Engineering, 16, 359368.
Golczynski, T. S., & Kempton, E. (2006). Understanding wax problems leads to deep-
water flow assurance solutions. World Oil, D7D10.
Green, D. W. (2008). Perrys chemical engineers handbook (8th ed.). New York:
McGraw-Hill.
Haaland, S. E. (1983). Simple and explicit formulas for the friction factor in turbulent
pipe flow. Journal of Fluids Engineering, 105, 8990.
Halstensen, M., Arvoh, B. K., Amundsen, L., & Hoffmann, R. (2013). Online estima-
tion of wax deposition thickness in single-phase sub-sea pipelines based on
acoustic chemometrics: A feasibility study. Fuel, 105, 718727.
Hammami, A., & Mehrotra, A. K. (1995). Liquid-solid-solid thermal behaviour of
n-C44+n-C50 and n-C25+n-C28 paraffinic binary mixtures. Fluid Phase Equilibria,
111, 253272.
Han, S., Huang, Z., Senra, M., Hoffmann, R., & Fogler, H. S. (2010). Method to deter-
mine the wax solubility curve in crude oil from centrifugation and high tem-
perature gas chromatography measurements. Energy & Fuels, 24, 17531761.
Hansen, A. B., Larsen, E., Pedersen, W. B., Nielsen, A. B., & Rnningsen, H. P. (1991).
Wax precipitation from North Sea crude oils. 3. Precipitation and dissolution
of wax studied by differential scanning calorimetry. Energy & Fuels, 5, 914923.
Hansen, J. H., Fredenslund, A., Pedersen, K. S., & Rningsen, H. P. (1988). A thermo-
dynamic model for predicting wax formation in crude oils. AIChE Journal, 34,
19371942.
Bibliography 153
Hayduk, W., & Minhas, B. (1982). Correlations for prediction of molecular diffusivi-
ties in liquids. Canadian Journal of Chemical Engineering, 60, 295299.
Hernandez, O. C. (2002). Investigation of single-phase paraffin deposition characteristics
(M.S. thesis). University of Tulsa, Tulsa, OK.
Hernandez, O. C., Hensley, H., Sarica, C., Brill, J. P., Volk, M., & Delle-Case, E. (2003).
Improvements in single-phase paraffin deposition modeling. In SPE Annual
Technical Conference and Exhibition (pp. 19). Denver, CO: Society of Petroleum
Engineers.
Hoffmann, R., & Amundsen, L. (2010). Single-phase wax deposition experiments.
Energy & Fuels, 24, 10691080.
Hsu, J. J. C., & Brubaker, J. P. (1995). Wax deposition measurement and scale-up mod-
eling for waxy live crudes under turbulent flow conditions. In International
Meeting on Petroleum Engineering (pp. 110). Beijing, China: Society of Petroleum
Engineers.
Huang, Z., Lee, H. S., Senra, M., & Fogler, H. S. (2011). A fundamental model of wax
deposition in subsea oil pipelines. AIChE Journal, 57, 29552964.
Huang, Z., Lu, Y., Hoffmann, R., Amundsen, L., & Fogler, H. S. (2011). The effect of
operating temperatures on wax deposition. Energy & Fuels, 25, 51805188.
Huang, Z., Senra, M., Kapoor, R., & Fogler, H. S. (2011). Wax deposition modeling of
oil/water stratified channel flow. AIChE Journal, 57, 841851.
Hunt, E. B., Jr. (1962). Laboratory study of paraffin deposition. Journal of Petroleum
Technology, 14, 12591269.
Ijeomah, C. E., Dandekar, A. Y., Chukwu, G. A., Khataniar, S., Patil, S. L., & Baldwin,
A. L. (2008). Measurement of wax appearance temperature under simulated
pipeline (dynamic) conditions. Energy & Fuels, 22, 24372442.
Ismail, L., Westacott, R. E., & Ni, X. (2008). On the effect of wax content on paraf-
fin wax deposition in a batch oscillatory baffled tube apparatus. Chemical
Engineering Journal, 137, 205213.
Jemmett, M. R., Deo, M., Earl, J., & Mogenhan, P. (2012). Applicability of cloud point
depression to cold flow. Energy & Fuels, 26, 26412647.
Jennings, D. W., & Weispfennig, K. (2005). Effects of shear and temperature on wax
deposition: Coldfinger investigation with a Gulf of Mexico crude oil. Energy &
Fuels, 19, 13761386.
Jiang, Z., Hutchinson, J., & Imrie, C. (2001). Measurement of the wax appearance tem-
peratures of crude oils by temperature modulated differential scanning calorim-
etry. Fuel, 80, 367371.
Jimenez, J., Moin, P., Moser, R., & Keefe, L. (1988). Ejection mechanisms in the sub-
layer of a turbulent channel. Physics of Fluids, 31, 13111313.
Juyal, P., Cao, T., Yen, A., & Venkatesan, R. (2011). Study of live oil wax precipitation with
high-pressure micro-differential scanning calorimetry. Energy & Fuels, 25, 568572.
Karacan, C. O., Demiral, M. R. B., & Kok, M. V. (2000). Application of x-ray CT imag-
ing as an alternative tool for cloud point determination. Petroleum Science and
Technology, 18, 835849.
Kaya, A. S., Sarica, C., & Brill, J. P. (1999). Comprehensive mechanistic modeling
of two-phase flow in deviated wells. In SPE Annual Technical Conference and
Exhibition. Houston, TX: Society of Petroleum Engineers.
Kim, D., Ghajar, A. J., Dougherty, R. L., & Ryali, V. K. (1999). Comparison of twenty two-
phase heat transfer correlations with seven sets of experimental data, including
flow pattern and tube inclination effects. Heat Transfer Engineering, 20(1), 1540.
154 Bibliography
Kleinhans, J., Niesen, V., & Brown, T. (2000). Pompano paraffin calibration field trials.
In SPE Annual Technical Conference and Exhibition (pp. 115). Dallas, TX: Society
of Petroleum Engineers.
Kok, M. V., Jean-Marie, L., Claudy, P., Martin, D., Garcin, M., Vollet, J., Ura, C.
(1996). Comparison of wax appearance temperatures of crude oils by DSC, ther-
momicroscopy and viscometry. Fuel, 75, 787790.
Kok, M. V., Letoffe, J. M., & Claudy, P. (1999). DSC and rheometry investigations of
crude oils. Journal of Thermal Analysis and Calorimetry, 56, 959965.
Kravchenko, V. (1946). No title. Acta Physicochimica, 21, 335.
Kruka, V. R., Cadena, E. R., & Long, T. E. (1995). Cloud-point determination for crude
oils. Journal of Petroleum Technology, 47, 681687.
Labes-Carrier, C., Rnningsen, H. P., Kolnes, J., & Leporcher, E. (2002). Wax deposi-
tion in North Sea gas condensate and oil systems: Comparison between opera-
tional experience and model prediction. In SPE Annual Technical Conference and
Exhibition. San Antonio, TX.
Lee, H. S. (2008). Computational and rheological study of wax deposition and gelation in
subsea pipelines (Ph.D. thesis). University of Michigan.
Leontaritis, K. J. (1996). The asphaltene and wax deposition envelopes. Fuel Science
and Technology International, 14, 1339.
Li, H., & Zhang, J. (2003). A generalized model for predicting non-Newtonian viscos-
ity of waxy crudes as a function of temperature and precipitated wax. Fuel, 82,
13871397.
Lindeloff, N., & Krejbjerg, K. (2002). A compositional model simulating wax deposi-
tion in pipeline systems. Energy & Fuels, 16, 887891.
Lira-Galeana, C., Firoozabadi, A., & Prausnitz, J. M. (1996). Thermodynamics of wax
precipitation in petroleum mixtures. AIChE Journal, 42, 239248.
Lu, Y., Huang, Z., Hoffmann, R., Amundsen, L., Fogler, H. S., & Sheng, Z. (2012).
Counterintuitive effects of the oil flow rate on wax deposition. Energy & Fuels,
26, 40914097.
Lund, H.-J. (1998). Investigation of paraffin deposition during single-phase liquid flow in
pipelines (M.S. thesis). University of Tulsa.
Martos, C., Coto, B., Espada, J. J., Robustillo, M. D., Gmez, S., & Pea, J. L. (2008).
Experimental determination and characterization of wax fractions precipitated
as a function of temperature. Energy & Fuels, 22, 708714.
Martos, C., Coto, B., Espada, J. J., Robustillo, M. D., Pea, J. L., & Merino-Garcia, D.
(2010). Characterization of Brazilian crude oil samples to improve the prediction
of wax precipitation in flow assurance problems. Energy & Fuels, 24, 22212226.
Matzain, A. (1997). Single phase liquid paraffin deposition modeling (M.S. thesis).
University of Tulsa.
Matzain, A., Apte, M. S., Zhang, H.-Q., Volk, M., Redus, C. L., Brill, J. P., & Creek,
J.L. (2001). Multiphase flow wax deposition modeling. In ETCE 2001: Petroleum
Production Technology Symposium. Houston, TX.
Merino-Garcia, D., & Correra, S. (2008). Cold flow: A review of a technology to avoid
wax deposition. Petroleum Science and Technology, 26, 446459.
Monger-McClure, T. G., Tackett, J. E., & Merrill, L. S. (1999). Comparisons of cloud
point measurement and paraffin prediction methods. SPE Production & Facilities,
14, 416.
Monrad, C. C., & Pelton, J. G. (1942). Heat transfer by convection in annular spaces.
Transations of AIChE, 38, 593611.
Bibliography 155
Mukherjee, H., & Brill, J. P. (1985). Pressure drop correlations for inclined two-phase
flow. Journal of Energy Resources Technology, 107, 549554.
Musser, B. J., Kilpatrick, P. K., & Carolina, N. (1998). Molecular characterization of
wax isolated from a variety of crude oils. Energy & Fuels, 59, 715725.
Nazar, A. R. S., Dabir, B., Vaziri, H., & Islam, M. R. (2001). Experimental and math-
ematical modeling of wax deposition and propagation in pipes transporting
crude oil. In SPE Production and Operations Symposium (pp. 111). Oklahoma
City, OK: Society of Petroleum Engineers.
Niesen, V. (2002). The real cost of subsea pigging. E&P Magazine, 9798.
Noville, I., & Naveira, L. (2012). Comparison between real field data and the results of
wax deposition simulation. In SPE Latin America and Caribbean Petroleum Engineer
ing Conference (pp. 112). Mexico City, Mexico: Society of Petroleum Engineers.
Paso, K. G., & Fogler, H. S. (2004). Bulk stabilization in wax deposition systems.
Energy & Fuels, 18, 10051013.
Paso, K., Kallevik, H., & Sjblom, J. (2009). Measurement of wax appearance tempera-
ture using near-infrared (NIR) scattering. Energy & Fuels, 23, 49884994.
Pauly, J., Daridon, J.-L., & Coutinho, J. A. P. (2004). Solid deposition as a function
of temperature in the nC10 + (nC24nC25nC26) system. Fluid Phase Equilibria,
224, 237244.
Pauly, J., Dauphin, C., & Daridon, J. L. (1998). Liquidsolid equilibria in a decane+
multi-paraffins system. Fluid Phase Equilibria, 149, 191207.
Pedersen, K. S., Skovborg, P., & Rnningsen, H. P. (1991). Wax precipitation
from North Sea crude oils. 4. Thermodynamic modeling. Energy & Fuels, 5,
924932.
Pedersen, K. S., Thomassen, P., & Fredenslund, A. (1985). Thermodynamics of petro-
leum mixtures containing heavy hydrocarbons. 3. Efficient flash calculation
procedures using the SRK equation of state. Industrial & Engineering Chemistry
Process Design and Development, 24, 948954.
Pedersen, W. B., Hansen, A. B., Larsen, E., Nielsen, A. B., & Rnningsen, H. P. (1991).
Wax precipitation from North Sea crude oils. 2. Solid-phase content as function
of temperature determined by pulsed NMR. Energy & Fuels, 5, 908913.
Petitjean, D., Schmitt, J. F., Laine, V., Bouroukba, M., Cunat, C., & Dirand, M. (2008).
Presence of isoalkanes in waxes and their influence on their physical properties.
Energy & Fuels, 22, 697701.
Phillips, D. A., Forsdyke, I. N., McCracken, I. R., & Ravenscroft, P. D. (2011). Novel
approaches to waxy crude restart: Part 2: An investigation of flow events fol-
lowing shut down. J. Pet. Sci. Eng., 77, 286304.
PVTsim 19 method documentation. (2009). Calsep.
Reistle, C. E. (1928). Methods of dealing with paraffin troubles encountered in producing
crude oil.
Reistle, C. E. (1932). Paraffin and gongealing-oil problems.
Roehner, R., & Fletcher, J. (2002). Comparative compositional study of crude oil solids
from the Trans Alaska Pipeline System using high-temperature gas chromatog-
raphy. Energy & Fuels, 16, 211217.
Roehner, R. M., & Hanson, F. V. (2001). Determination of wax precipitation tempera-
ture and amount of precipitated solid wax versus temperature for crude oils
using FT-IR spectroscopy. Energy & Fuels, 15, 756763.
Rnningsen, H. P. (2012). Production of waxy oils on the Norwegian continental shelf:
Experiences, challenges, and practices. Energy & Fuels, 26, 41244136.
156 Bibliography
Rnningsen, H. P., Bjamdal, B., Hansen, A. B., & Pedersen, W. B. (1991). Wax
precipitationfrom North Sea crude oils. 1. Crystallization and dissolution tem-
peratures, and Newtonian and non-Newtonian flow properties. Energy & Fuels,
5, 895908.
RPSEA. (2007). RPSEA 1201 controlling wax deposition in the presence of hydrates
technology evaluation. Retrieved from http://www.netl.doe.gov/kmd/RPSEA
_Project_Outreach/07121-1201-University_of_Utah_Phase1.pdf.
Rygg, O. B., Rydahl, A. K., & Rnningsen, H. P. (1998). Wax deposition in offshore
pipeline systems. In BHRGroup 1st North American Conference Multiphase
Technology. Baniff, Canada.
Saffman, P. G. (1965). The lift on a small sphere in a slow shear flow. Journal of Fluid
Mechanics, 22, 385400.
Sandler, S. I. (2006). Chemical, biochemical, and engineering thermodynamics. Hoboken,
NJ: Wiley.
Senra, M. J. (2009). Assessing the role of polydispersity and cocrystallization on crystallizing
n-alkanes in n-alkane solutions (Ph.D. thesis). University of Michigan, Hoboken,
NJ.
Sieder, E. N., & Tate, G. E. (1936). Heat transfer and pressure drop of liquids in tubes.
Industrial and Engineering Chemistry, 28, 14291435.
Singh, A., Lee, H., Singh, P., & Sarica, C. (2011). Flow assurance: Validation of wax
deposition models using field data from a subsea pipeline. In Offshore Technology
Conference (pp. 119). Houston, TX: Offshore Technology Conference.
Singh, P. (2000). Gel deposition on cold surfaces (Ph.D. thesis). University of Michigan.
Singh, P., & Venkatesan, R. (2001). Morphological evolution of thick wax deposits
during aging. AIChE Journal, 47, 618.
Singh, P., Venkatesan, R., Fogler, H. S., & Nagarajan, N. (2000). Formation and aging
of incipient thin film wax-oil gels. AIChE Journal, 46, 10591074.
Smith, B. (1999). Infrared spectral interpretation, a systematic approach. New York: CRC
Press.
Snyder, R. G., Conti, G., Strauss, H. L., & Dorset, D. L. (1993). Thermally induced mix-
ing in partially microphase segregated binary n-alkane crystals. The Journal of
Physical Chemistry, 97, 73427350.
Snyder, R. G., Goh, M. C., Srivatsavoy, V. J. P., Strauss, H. L., & Dorset, D. L. (1992).
Measurement of the growth kinetics of microdomains in binary n-alkane
solid solutions by infrared spectroscopy. The Journal of Physical Chemistry, 96,
1000810019.
Snyder, R. G., Hallmark, V. M., Strauss, H. L., & Maroncelli, M. (1986). Temperature
and phase behavior of infrared intensities: The poly(methy1ene) chain. The
Journal of Physical Chemistry, 94720, 56235630.
Snyder, R. G., Srivatsavoy, V. J. P., Cates, D. A., Strauss, H. L., White, J. W., & Dorset,
D. L. (1994). Hydrogen/deuterium isotope effects on microphase separation
in unstable crystalline mixtures of binary n-alkanes. The Journal of Physical
Chemistry, 98, 674684.
Streitwieser, A., & Heathcock, C. (1976). Introduction to organic chemistry. New York:
Macmillan Publishing.
Suppiah, S., Ahmad, A., Alderson, C., Akbarzadeh, K., Gao, J., Shorthouse, J.,
Jamaluddin, A. (2010). Waxy crude production management in a deepwater sub-
sea environment. In SPE Annual Technical Conference and Exhibition (pp.118).
Florence, Italy: Society of Petroleum Engineers.
Bibliography 157
Urushihara, T., Meinhart, C. D., & Adrain, R. J. (1993). Investigation of the logarithmic
layer in pipe flow using particle image velocimetry. In Near Wall Turbulent Flows,
(pp. 433446). Amsterdam, Netherlands: Elsevier.
Van Driest, E. R. (1956). On turbulent flow near a wall. J. Aeronaut. Sci., 23, 10071011.
Venkatesan, R. (2004). The deposition and rheology of organic gels (Ph.D. thesis).
University of Michigan.
Venkatesan, R., & Fogler, H. S. (2004). Comments on analogies for correlated heat and
mass transfer in turbulent flow. AIChE Journal, 50(7), 16231626.
Venkatesan, R., Nagarajan, N. R., Paso, K., Yi, Y.-B., Sastry, A. M., & Fogler, H. S.
(2005). The strength of paraffin gels formed under static and flow conditions.
Chemical Engineering Science, 60, 35873598.
Vieira, L. C., Buchuid, M. B., & Lucas, E. F. (2010). Effect of pressure on the crystal-
lization of crude oil waxes. I. Selection of test conditions by microcalorimetry.
Energy & Fuels, 24, 22082212.
Wilke, C. R., & Chang, P. (1955). Correlation of diffusion coefficients in dilute solu-
tions. AIChE Journal, 1, 264270.
Wilkes, J. O. (2005). Fluid mechanics for chemical engineers (2nd ed.). Upper Saddle
River, NJ: Prentice Hall.
Won, K. W. (1986). Thermodynamics for solid solution-liquid-vapor equilibria: Wax
phase formation from heavy hydrocarbon mixtures. Fluid Phase Equilibria, 30,
265279.
Xiao, J. J., Shonham, O., & Brill, J. P. (1990). Comprehensive mechanistic model for
two-phase flow in pipelines. In SPE Annual Technical Conference and Exhibition,
New Orleans, LA: Society of Petroleum Engineers.
Yan, D., & Luo, Z. (1987). Rheological properties of Daqing crude oil and their appli-
cation in pipeline transportation. SPE Production Engineering, 2, 267276.
Zheng, S., Zhang, F., Huang, Z., & Fogler, H. S. (2013). Effects of operating conditions
on wax deposit carbon number distribution: Theory and experiment. Energy &
Fuels, 27, 73797388.
Zhu, T., Walker, J. A., & Liang, J. (2008). Evaluation of wax deposition and its control dur-
ing production of Alaska North Slope OilsFinal Report.
Appendix A: Nomenclature
159
160 Appendix A
163
164 Index
Diffusion coefficient, 143. See also Wax in Statoils Herya Research Center,
diffusion coefficient 110113, 111f113f, 112t
of waxy components in the oil, 70, 77, test section, 93, 94f
8384 Fourier transform infrared
Diffusion-equivalent mass transfer spectroscopy (FT-IR)
layer (diffusion), 87 for WAT determination, 1718,
Dissolved waxy components 17f19f, 19t, 26t, 28f
concentration of, 8486, 85f, 86f for WPC determination, 3435, 35f
on deposit surface, 71 Free induction decay (FID), 33, 33f
internal diffusion of, 7172, 72f73f FT-IR. See Fourier transform infrared
precipitation of, 6970 spectroscopy (FT-IR)
radial concentration gradient of, 70 Fugacity, 45
calculation of, 46
vaporliquid, 46
E
Enthalpy of crystallization, 32, 47, 52 G
Equation of state (EOS)
PengRobinson, 46 Gibbs free energy, 46, 47, 51, 52f
vaporliquid fugacity, 46 Gravity settling, and wax deposition, 67
Experiments, wax deposition, 91105
deposit characterization, 94102
H
effect of temperature on, 117118
flow loop. See Flow loop wax Haaland correlation, 9495
deposition experiment HaydukMinhas correlation, 120121,
heat transfer/mass transfer analysis 143
of, 117127 Heat transfer
reproducibility of, 113f analysis, of wax deposition
selection of, 110113, 111f113f, 112t experiments, 117127
equations/correlations, 7880, 79f,
80t
F and wax deposit thickness, 97100,
FID. See Free induction decay (FID) 98f
Field, wax deposition in, 131133, 132f, Herya Research Center, Statoil ASA
133f flow loop wax deposition
examples of experimentin, 110113,
gas/oil flow, 140145 111f113f, 112t
single-phase pipe flow, 134140 High-temperature gas chromatography
predictions, Michigan Wax Predictor (HTGC), 59, 60f, 61, 101
for, 136140, 137f140f n-paraffin distributions using, 133
Filter-plugging Hot flow, 118
for WAT detection, 29t Hydrodynamics, in wax deposition
FloWax software package, 75, 78t, 80t model, 7778, 78t
Flow loop wax deposition experiment,
9193, 92f
I
apparatus setup, 104
conditioning system, 9193, 92f Industrial practice, wax thermodynamic
oil characterization, 105 model, 5864, 59f
pumping system, 9193, 92f fluid characterization, 5863, 59f62f
Index 165
K N
Kinetic precipitation rate constant, 86 n-alkane
solidliquid phase transition, 47
Near-infrared (NIR) scattering, 11
L
for WAT detection, 29t
Lab-scale wax deposition experiments, 91 Newtonian fluid, 96
Laser technique, for deposit thickness NMR. See Nuclear magnetic resonance
measurement, 100101, 100f (NMR)
Lewis number, 89 n-paraffins, 1, 5, 37
Liquid displacement-level detection carbon number of, 101, 102f
technique, 100 distribution of, 5962, 59f, 60f, 62f,
133, 135, 135f
thermodynamic modeling for,
M
135136, 136f
Macrocrystalline wax, 38 Nuclear magnetic resonance (NMR)
Mass balance equations, 45, 81 for WPC determination, 3334, 33f
Mass diffusivity for wax, 108, 119, 130, Nusselt number, 88
143, 144, 146
Mass flux
O
characteristic of wax deposition,
120121, 122t123t Offshore petroleum production
determination of, 83 shifted from onshore toward, 12, 2f
Mass transfer analysis, of wax Oil flow rate, 93, 96, 115, 117, 120, 130
deposition experiments, Oil temperature, 2, 117, 118
117127 impact on wax deposition, 123125,
Mass transfer driving force, 121 124f125f, 125t
Mass transfer layer (mass transfer) OLGA wax commercial/wax deposition
conduction-equivalent, 8789 model, 74, 140141, 142, 144,
diffusion-equivalent, 87 144f
thickness of, 8689 1-D discretization, in wax deposition
Michigan Wax Predictor (MWP), 75, model, 83, 83f
114f116f, 114117, 116t heat transfer equations, 7879
defined, 129 temperature profiles in, 76, 76f
for field predictions, 136140, 137f140f Onshore petroleum production
2-D discretization, 7677, 77f, 79 shifted toward offshore production,
Microcrystalline wax, 38 12, 2f
Molecular diffusion, and wax and wax deposition problem, 1
deposition, 67, 68 Operating temperature
dissolved waxy components impact on wax deposition, 122127
on deposit surface, 71
internal diffusion of, 7172, 72f73f
P
precipitation of, 6970
radial concentration gradient of, 70 Pedersens wax precipitation model, 55
mathematical representation of, PengRobinson EOS
8082, 81f for fugacity calculation, 46
steps in, 6973, 70f Phase equilibrium equations, 4445, 44f
MWP. See Michigan Wax Predictor construction of, 46
(MWP) solidliquid, 46
166 Index