CO2 Pipeline Transportation Analysis
CO2 Pipeline Transportation Analysis
of thermodynamics
Vol. 35(2014), No. 1, 117–140
DOI: 10.2478/aoter-2014-0008
ANDRZEJ WITKOWSKI∗
MIROSŁAW MAJKUT
SEBASTIAN RULIK
∗
Corresponding Author. E-mail: [Link]@[Link]
118 A. Witkowski, M. Majkut and S. Rulik
1 Introduction
are based on the physical properties of CO2 containing impurities and were
obtained using real gas equations of state (EOS) with the Lee and Kesler
equation modified by Plocker, Knapp, and Prausnitz (the LKP equation of
state) [2,8] and the Peng-Robinson equation of state with Boston-Mathias
modifications (PRBM) [2,14]. All these equations are included in the AS-
PEN PLUS V 0.7 design process simulator [2].
Moreover, Fig. 2 shows that CO2 compressibility is nonlinear in the
range of pressures common for pipeline transport and is highly sensitive to
any impurities as predicted by the Peng-Robinson equation of state [10].
Thus, it is necessary to use accurate representations of the phase behavior,
density, and viscosity of CO2 while designing a pipeline. To reduce difficul-
ties in design and operation, it is generally recommended that the pressure
in a CO2 pipeline should be higher than 8.6–10 MPa because the abrupt
changes in CO2 compressibility can then be avoided across the range of
temperatures that may be encountered in the pipeline system [5,10].
Figure 2. Nonlinear compressibility of CO2 in the range of pressures common for pipeline
transport as predicted by the Peng-Robinson equation of state [10].
Analysis of pipeline transportation systems. . . 121
3 Pipeline configuration
A complete CO2 capture and sequestration (CCS) system requires safe, reli-
able and cost-efficient solutions for transmission of CO2 from the capturing
facility to the permanent storage site. The CO2 pipeline facilities from the
power plant flue stack, through the separation, compression, transportation
and pumping systems to the injection well head are shown schematically
in Fig. 3. For sequestration purposes, CO2 is generally injected to depths
exceeding 1000 m [11]. At greater depths, CO2 increases in density and
below 800 m it becomes a supercritical fluid. A large compressing system
is required to compress the source CO2 under nearly atmospheric pressure
conditions. Figure 3 provides an example of an eight-stage integrally geared
compressor system, which was selected in [16] as the most available, reliable
and efficient compression technology and injection configuration. Compres-
sor stations in a CO2 pipeline system can be subdivided into two classes:
stations positioned at the pipeline inlet, and booster stations located along
the pipeline to compensate for the pressure drop due to friction and eleva-
tion losses.
certain distance. Otherwise, it will transform into the liquid state. Mean-
while, some authors propose that transporting liquid CO2 at a relatively
low temperature is preferred in terms of a reduction in the pressure drop
along the pipeline [19]. In the present paper the CO2 working area is as-
sumed to be either in the liquid or in the supercritical state and the results
of these two states are compared.
For subcooled liquid CO2 transportation, a facility cooling CO2 to 15 o C
or less is needed so that it can be kept below its critical temperature down
the line. Conceptually, refrigeration could be added along the pipeline using
CO2 as the working medium, but this will obviously increase the capital and
operating costs as well as reduce the overall energy efficiency. Generally,
CO2 transportation in the subcooled liquid state has some advantages over
the supercritical state transport, most importantly because of the liquid
lower compressibility and higher density within the pressure range consid-
ered here, which permits smaller pipe sizes and generates lower pressure
losses. To reduce difficulties in design and operation, it is generally recom-
mended that the pipeline should operate at pressures higher than 8.6 MPa.
Then, the abrupt changes in CO2 compressibility and specific heat can be
avoided across the range of temperatures (Fig. 2) that may be encountered
in the pipeline system [5,10].
Most new pipelines are laid underground, despite the higher initial costs,
for environmental, security and safety reasons. Underground temperatures
are much more stable than surface temperatures. Hence the operating tem-
peratures of CO2 pipelines are generally dictated by the temperature of the
surrounding soil.
The pressure keeps dropping along the pipeline until CO2 evaporates and
the pipeline may eventually be blocked. This means that there is a max-
imum safe transport distance. If there is a need to transport CO2 farther
than over this maximum distance, boosting pump station are required along
the pipeline.
A pipeline segment is defined as a length of pipeline for which the inlet
pressure and the minimum outlet pressure values are specified, e.g., a length
of pipeline between two compressor stations.
pmop d2
t= , (1)
2S E F
in which pmop is the maximum operating pressure of the pipeline (Pa), d2
is the outside pipe diameter (m), S is the specified yield stress for the pipe
material (Pa), E is the longitudinal joint factor, and F is the design factor.
Analysis of pipeline transportation systems. . . 125
For the purpose of estimating the pipe wall in this case study, the maximum
operating pressure is assumed at 153 MPa, the longitudinal joint factor is
1.0, and the design factor is 0.72 (US Code of Federal Regulations – CFR,
[13]). The minimum yield stress is specified as 483 MPa, which corresponds
to the API (American Petroleum Institute) for 5L X-70 steel line pipe [1].
Alternatively, a thermal insulation layer was assumed on the pipeline exter-
nal surface. The flow process of CO2 along the pipeline is mainly influenced
by three factors: friction forces, heat exchange through the pipe wall be-
tween the soil and the thermal insulation layer, and the change in elevation.
Based on the data given above, thickness of the pipe wall, t, was calculated
to be 10 mm. Since the critical temperature of CO2 is predominately higher
than the normal temperature of the soil or the ambient temperature, after
CO2 is compressed and boosted to above 10 MPa, the thermal insulation
will slow down the CO2 temperature decrease process, increasing the pres-
sure drop in the pipeline. Therefore, the thermal insulation layer should
not be laid on the pipeline external surface in Poland. In the present work,
the CO2 working area was assumed to be either in the liquid or in the
supercritical state.
a)
b)
Figure 4. Comparison of the pressure drop (a) and density changes (b) along the pipeline
for adiabatic and isothermal conditions at different inlet CO2 temperatures.
CO2 at the inlet is in the liquid state because the liquid is incompressible. It
also shows that a negligible difference exists between adiabatic and isother-
mal transmission when the CO2 inlet state is liquid at 15 o C. The difference
becomes significant when CO2 reaches the supercritical fluid state (35 o C).
Under isothermal conditions and with an initial subcooled liquid, there is
a sufficient heat transfer to the ground, and, therefore, the transmission dis-
tance is shortened. In adiabatic conditions, there is no heat transfer across
the pipe wall and the transmission distance is longer. As it was accepted
earlier, the lowest transport pressure of the CO2 pipeline is 9 MPa. In this
Analysis of pipeline transportation systems. . . 127
case, the maximum safe transport distance depends on the thermal condi-
tions, changing from 469 km at the CO2 subcooled temperature of 0 o C to
339 km at the supercritical temperature of 35 o C, and under isothermal con-
ditions. For the adiabatic case, there is a discontinuity in the density profile.
This discontinuity corresponds with the saturation state when CO2 changes
density quickly, e.g., from the subcooled liquid state to a two-phase state.
Even for supercritical fluids where there is no phase change, the density
variation has a very nonlinear region. With initial temperatures above the
supercritical point, the CO2 density changes abruptly within the pipeline
once the temperature reaches the saturation point, and reaches two-phase
flow conditions. Moreover, Fig. 4 also indicates that if the inlet temperature
is higher than critical, CO2 repressurization is needed after a much shorter
distance, whether adiabatic or isothermal conditions prevail.
a)
b)
c)
Figure 5. Maximum safe transport distances (a), density (b) and velocity (c) with differ-
ent internal pipe diameters in adiabatic and isothermal conditions.
Analysis of pipeline transportation systems. . . 129
the increase in capital and maintenance costs. The pipeline needs not be
insulated when advantage can be taken of cold ground conditions, which
maintain liquid conditions. In the present work, the CO2 working area was
assumed to be either in the liquid or in the supercritical state. In the latter
case, thermal insulation was assumed on the external surface of the pipe
(Fig. 6). Engineering experience indicates that long distance pipelines are
usually placed underground at a depth of 1.2–1.5 m. It is also found that the
annual lowest and highest soil temperatures at a depth of 1.5 m in Poland
are between 5 and 16 o C (Fig. 7).
Figure 6. Cross section of the underground and insulated pipeline. r1 – internal radius
of pipe, r2 – external radius of pipe, re – external wall radius of the thermal
insulation layer, t – pipewall thickness, tsurf ace – ground surface temperature,
tsoil – soil temperature.
Two cases – with the lowest and the highest soil temperatures – are
considered in the pipeline design to ensure that the pipeline can operate
well over the whole year.
A two-dimensional heat conduction formula [7] can be used to calculate
the overall heat transfer coefficient between the ground and the CO2 in the
pipeline (Fig. 6):
1
k= r1 r2 r1 re r1 r1 1 , (2)
λpw ln r1 + λti ln r2 + λsoil ln 2z
re + z αag
where heat conductivities of the pipe wall and of the insulation layer mate-
rials are λpw = 25 W/mK and λti = 0.058 W/mK, respectively.
The thermal conductivity of the soil is assumed at λsoil = 1.21 W/(mK),
and the distance between the ground surface and the pipe center is 1.225 m.
The air convection heat transfer coefficient is αg = 5 W/m2 K. According
130 A. Witkowski, M. Majkut and S. Rulik
Figure 7. Annual soil temperatures at different depths and ambient temperatures, tsoil
– soil temperature.
to [18] and [19], the thermal resistance of the convective thermal transfer
between CO2 and the inner pipe wall is much smaller than that between the
pipe wall and the heat insulation layer, so it is assumed that the tempera-
ture of the inner pipe wall is equal to the temperature of CO2 on the same
cross section. The value of the overall heat transfer coefficient k between
the ground and the CO2 calculated from Eq. (2) for a pipeline with a 0.05 m
and 0.03 m heat insulation layer is 0.7387 and 0.912 W/m2 K, respectively,
and for a pipeline without insulation – 2.11 W/m2 K. Resistance to the heat
transfer from the tubing and casing is ignored, as the conductivity of the
steel used in the tubing and casing is at least an order of magnitude larger
than any other conductivity in the system.
An increase in the pressure drop means higher operating costs and possibly
the need to introduce compressor stations. The pressure drop along the
pipeline is dependent on the flow velocity, ambient temperature, as well as
on geometric characteristics of the pipeline such as length, elevation changes
etc. The pressure drop and temperature changes along the pipeline reduce
the CO2 density and increase velocity, which, in turn, increases the pressure
drop and ultimately leads to a choking condition at a certain distance. The
maximum safe pipeline length to prevent choking is preliminarily selected
at a value about 10% smaller than the choking point [18]. The dependence
Analysis of pipeline transportation systems. . . 131
Figure 8. Dependence of the safe pipeline length and of the choking point on CO2 pipe
inlet temperatures (20 o C, 30 o C, 35 o C, 45 o C) at energy balance with sur-
roundings.
Figure 9. Dependence of the safe pipeline length and of the choking point on ambient
temperatures (0 o C, 10 o C, 20 o C, 30 o C).
132 A. Witkowski, M. Majkut and S. Rulik
of the safe distance on the pipeline inlet and ambient temperatures for the
reference case study is presented in Figs. 8 and 9, respectively. It can be
seen that both the CO2 temperature (Fig. 8) and the ambient temperature
(Fig. 9) have a significant impact on the pipeline maximum safe length.
In order to avoid the choking condition, recompression of CO2 becomes
necessary.
a)
Figure 10. Comparison of the pressure drop (a), the change in temperature (b), in density
(c), and in velocity (d) along a pipeline with and without thermal insulation
for the conditions of energy balance with surroundings at different ambient
temperatures tamb : 0o C, 15o C, 20o C, 30o , tCO2 – inlet CO2 temperature, k –
overall heat transfer coefficient; LKP equation of state.
b)
c)
d)
Figure 10b-d (continued).
Analysis of pipeline transportation systems. . . 135
a)
b)
c)
Figure 11a-c. For caption see next page.
136 A. Witkowski, M. Majkut and S. Rulik
d)
Figure 11. Comparison of the pressure drop (a), the change in temperature (b), in density
(c), and in velocity (d) along a pipeline with and without thermal insulation
for the conditions of energy balance with surroundings at different ambient
temperatures: 0o C, 15o C, 20o C, 30o ; PRBM equation of state.
Figure 12. Maximum safe length of a pipeline with (k = 0.738) and without (k = 2.11) an
insulating layer; distance to subsequent booster station at different ambient
temperatures; LKP equation of state.
8 Conclusions
The paper presents the analysis of the influence of multiple factors, including
pipe diameter, pipeline inlet temperature, ambient temperature or thermal
Analysis of pipeline transportation systems. . . 137
Figure 13. Maximum safe length of a pipeline with and without an insulating layer; dis-
tance to subsequent booster station at different ambient temperatures; PRBM
equation of state.
References
[1] American Petroleum Institute: Spec 5L-Specification for Line Pipe, 43rd Edn. Wash-
ington 2004.
[2] Aspen Plus, Version 7.0. Computer program, 2008.
[3] Beggs H.D., Brill J.P.: A Study of two-phase flow in inclined pipes. J. Petrol
Technol., 25(1973), 607–617.
[4] Chmielniak T., Łukowicz H.: Condensing power plant cycle – assessing possibil-
ities of improving its efficiency. Arch. Thermodyn. 31(2010), 3, 105–113.
[5] Farris C.B.: Unusual design factors for supercritical CO2 pipelines. Energy
Progress 3(1983), 3, 150–158.
[6] Gottlicher Gerold: The Energetics of Carbon Dioxide Capture in Power Plants.
VDI Verlag, Dusseldorf 1999.
[7] Incropera F.P., DeWitt D.P.: Introduction to heat transfer, 3rd Edn. John
Wiley & Sons, Inc., New York 1996.
[8] Ludtke K.H.: Process Centrifugal Compressors. Springer, 2004.
[9] Łukowicz H., Dykas S., Rulik S., Stępczyńska K.: Thermodynamic and
economic analysis of a 900 MW ultra-supercritical power unit. Arch. Thermodyn.
32(2011), 3, 231–244.
[10] McCoy S.T., Rubin E.S.: An engineering-economic model of pipeline transport of
CO2 with application to carbon capture and storage. Int. J. Greenhouse Gas Contr.
2(2008), 2, 219–229.
[11] Mohitpour M., Seevam P., Botros K.K., Rothwell B., Ennis C.: Pipeline
Transportation of Carbon Dioxide Containing Impurities. ASME, New York 2012.
[12] Polish Committee for Standardization: Thermal Insulation of Ducts, Fittings and
Installations, PN-B-02421, Jul. 2000.
[13] Transportation. Title 49 Code of Federal Regulations Pt. 195, 2005 Edn. 170–171.
[14] Rao B.: Multiphase Flow Models Range of Applicability. CTES, L.C, May 18, 1998.
[15] Wayne C., Edmister Byung Ik Lee: Applied Hydrocarbon Thermodynamics Vol.
1. Gulf Publishing Company, London Houston, London, Paris, Tokyo 1984.
[16] Witkowski A., Majkut M.: The impact of CO2 compression systems on the com-
pressor power required for a pulverized coal power plant in post-combustion carbon
sequestration. Arch. Mechanical Eng., LIX(2012), 3, 343–360.
[17] Witkowski A., Rusin A., Majkut M., Rulik S., Stolecka K.: Comprehensive
Analysis of the Pipeline Transportation Systems for CO2 Sequestration. Thermody-
namics and Safety Problems. In: Proc. 3rd Int. Conf. Contemporary Probl. Thermal
Eng. CPOTE 2012, 18-20 September, 2012, Gliwice.
140 A. Witkowski, M. Majkut and S. Rulik
[18] Zhang Z.X., Wang G.X. Massarotto P., Rudolph V.: Optimization of
pipeline transport for CO2 sequestration. Energ. Conver. Manage. 47(2006), 702–
715.
[19] Zhang D., Wang Z., Sun J., Zhang L., Zheng L.: Economic Evaluation of
CO2 pipeline transport in China. Energ. Conver. Manage. 55(2012), 127–135.