Minimizing Power Consumption Related to BOG Reliquefaction in an

Jun 14, 2016 - In this work, we study BOG reliquefaction at a typical regasification ...... (4) Park, C.; Song, K.; Lee, S.; Lim, Y.; Han, C. Retrofit...
0 downloads 0 Views 6MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article pubs.acs.org/IECR

Minimizing Power Consumption Related to BOG Reliquefaction in an LNG Regasification Terminal Harsha N. Rao, Khai H. Wong, and Iftekhar A. Karimi* Department of Chemical & Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, Singapore 117585 ABSTRACT: The generation of boil-off gas (BOG) is inevitable in liquefied-natural-gas (LNG) regasification terminals. It can be a safety concern and, in most cases, demands BOG reliquefaction before regasification. This study presents a fundamental thermodynamic analysis of the BOG reliquefaction−LNG regasification process, which consumes significant power. We propose an evolving series of schemes that exploit the cryogenic energy of the send-out LNG to reduce the total power consumption. We synergistically integrate options such as BOG precooling, recompression, interstage cooling, and recondensation in our proposed schemes. We develop a nonlinear programming (NLP) formulation to minimize the power required for each scheme for a given amount of BOG. Under some idealistic assumptions, the optimal solutions to these NLPs give lower bounds on the total power consumption, which can be used to evaluate the cost-effectiveness of terminal operations. A novel scheme consisting of BOG precooling before compression and two-stage recondensation by direct mixing has the least power consumption. For a case-study terminal, it offers 11.9% lower power consumption than the best-reported design in the literature for 15% BOG.

1. INTRODUCTION Natural gas (NG), the cleanest fossil fuel, has been the fastestgrowing energy resource for more than two decades and will be so for the foreseeable future. It offers the highest efficiency in power generation with the lowest carbon emissions. Although pipelines are the best for transporting NG over short distances, today’s diverse and long supply chains necessitate global NG transport over long distances. For this purpose, liquefied natural gas (LNG) is the most common and economical option that has overcome geopolitical constraints on the global gas supply.1 NG from a gas reservoir is pretreated prior to its energyintensive liquefaction at an onshore plant or offshore vessel. The resulting LNG is stored in atmospheric storage tanks at a loading terminal before being shipped in specially insulated LNG tankers to a receiving terminal. The latter also stores the incoming LNG in atmospheric tanks and regasifies it continuously to supply NG to gas distribution grids or customers. From the time it is liquefied to the time it is regasified, LNG is continuously maintained at atmospheric pressure and cryogenic temperatures (about −160 °C) through autorefrigeration. Despite extensive insulation, heat leakage into the LNG during this entire duration is unavoidable. This causes part of the LNG to vaporize throughout its supply chain to form what is known as boil-off gas (BOG). Managing BOG at various points (storage, loading, transport, unloading, and regasification) in the LNG supply chain2,3 is critical, because transporting or storing it as a gas is uneconomical. Unless reliquefied and/or used in some form, BOG represents a loss of valuable product, and its flaring is becoming increasingly unacceptable. The loss of BOG also © XXXX American Chemical Society

changes the LNG composition and quality over time. This affects the LNG trade, as LNG is priced based on its energy content.3 Most regasification terminals deliver NG at high pressures (30−80 bar) to their gas grids or customers. Compressing atmospheric BOG to such high pressures is also uneconomical. Therefore, unless an opportunity for direct, gainful, complete, and inexpensive use of the BOG exists without much energy penalty, most regasification terminals employ a system to reliquefy the BOG. Efficient BOG management systems are thus very important. In fact, the power/energy required to manage BOG is a major cost for regasification terminals. Thus, minimizing the energy consumption at LNG regasification terminals is of great interest and a worthwhile undertaking that has received some attention in the literature.4,5 BOG arises from several sources in a regasification terminal. The heat leaking continuously into the LNG storage tanks is an obvious source, but loading and unloading activities at the terminal are even larger sources. The tanks are typically far from the jetties, and the transfer lines connecting the two must be continuously kept cold to avoid substantial and periodic expansion/contraction of pipes and other infrastructure, which could cause damage. This line cooling is typically achieved by recirculating the LNG from the tanks continuously through the transfer lines. The recirculating LNG sweeps heat from the Received: April 7, 2016 Revised: June 8, 2016 Accepted: June 14, 2016

A

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

Figure 1. Schematic of a typical regasification terminal.

scenarios. For each scheme, we present a nonlinear programming (NLP) formulation to obtain its minimum power usage. Finally, we demonstrate significant power reductions from our schemes for the terminal considered by Park et al.4

surroundings, pumps, turbulence, and friction. In some terminals, this LNG is dumped back into the storage tanks, which generates a major portion of the total BOG. Most literature on LNG regasification has focused on exploiting the cryogenic energy of LNG to produce power by integrating various cycles such as Organic Rankine, Brayton, etc.6,7 Most of these have focused on exergetic and/or energetic efficiency, and very few8−10 have considered cost. In contrast, the literature on BOG management topics such as estimation,2,11−13 minimization,2,14,15 reliquefaction,4,5,16−18 direct compression,13 and utilization19 is limited. Wordu and Peterside12 estimated the BOG arising from the heat leak into an LNG storage tank. Hasan et al.2 estimated BOG generation during loading, laden voyage, unloading, and ballast voyage of an LNG tanker and performed simulations to find the optimal heel for maritime transport. Adom et al.11 studied the effect of pressure and heat leakage on BOG generation. Park et al.14 simulated and optimized the line-cooling recirculation in an LNG receiving terminal. Lee et al.15 suggested an unloading procedure for minimizing BOG generation in an operation with a combination of above-ground and in-ground LNG tanks. Jang et al.17 presented an algorithm for an optimal operating schedule of BOG compression depending on the dynamic behavior of an LNG tank. Shin et al.18 proposed a mixed integer linear problem (MILP) formulation to optimize BOG compressor operation based on the dynamics of tank pressure. Liu et al.5 studied BOG reliquefaction by direct, multistage mixing with pressurized LNG. However, they did not use highpressure LNG (HP-LNG) for cooling, which can reduce the power usage significantly. Later, Park et al.4 and Li et al.16 studied the retrofit design of a BOG reliquefaction system. They used HP-LNG for some cooling, but did not maximize its use and did not consider all options for power reduction. Querol et al.13 reviewed the various recondensation and direct compression options and proposed a combined heat and power system. Recently, Park et al.19 integrated LNG regasification with a gas-to-liquid (GTL) plant that uses the BOG as a chemical feedstock. In this work, we study BOG reliquefaction at a typical regasification terminal (Figure 1). Using selected design/ operation guidelines, we evolve several novel reliquefaction schemes that maximize the use of LNG cold energy in different ways to substantially reduce power consumption. We begin with a precise problem statement, develop our reliquefaction schemes, and then evaluate them under some idealistic

2. PROBLEM STATEMENT Consider the seawater-based LNG regasification terminal in Figure 1. On a long-term basis, it receives and gasifies F t/h (throughput) of LNG, on average, and supplies F t/h (tph) of NG (send-out rate) to a local gas distribution grid at specified conditions. Let βF tph (0 < β < 1) of BOG exit the tank tops at (P1, T1) and (1 − β)F tph of LNG be available for send-out at (p1, t1). Normally, BOG is superheated (compared to tank conditions) as a result of the heat ingress into the tanks and/or downstream pipelines. We allow the user the flexibility to give BOG and LNG specifications independently based on real terminal data. With this setup, we address the following problem in this work. Given (1) the LNG throughput or gas send-out rate F t/h (tph) (with no loss of generality, we assume F = 1 tph); (2) the average composition and conditions (p1, t1) of the LNG available for send-out; (3) the average BOG composition, generation fraction β, and conditions (P1, T1); (4) the send-out temperature (TS) and pressure (PS); and (5) the ratio (psw) of the power required to pump seawater to the rate of heating required for regasification Obtain the minimum power consumption (kW) of the terminal per tonne of send-out NG for any given β, and the scheme of process steps to achieve the same Assuming that (1) the flows, compositions, temperatures, pressures, and pressure drops of streams are long-term averages (essentially, the dynamics of the terminal operation arising due to loading, unloading, delivery, recirculation cooling, vessel cooling, etc., are neglected by averaging over a long time); (2) the LNG and BOG compositions are in their normal ranges found in the industry (i.e., the nitrogen content is sufficiently low to make BOG reliquefaction possible); B

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

represent the LNG exiting the tanks, and let (PS, HS, TS) represent the NG sent out to the gas grid or customers. Note that the LNG and BOG need not be in equilibrium, which is the reality in regasification terminals, contrary to the popular assumption of equilibrium.5 We require a series of processing steps that transform the BOG at (P1, H1, T1) and the LNG at (p1, h1, t1) to the send-out NG at (PS, HS, TS) using minimal power. We exploit the following basic thermodynamic facts: (1) Work for isentropic compression/pumping is given by the enthalpy change on a P−H diagram. (2) Work for isentropic compression/pumping increases with the feed-fluid temperature. (3) Pumping needs much less energy than compression for a given pressure ratio. (4) The latent heat of vaporization decreases with pressure. Obviously, we must heat and pressurize β tph of the BOG and (1 − β) tph of the LNG to obtain 1 tph of the send-out gas. As gas compression is much more energy-demanding than liquid pumping, we must maintain the LNG’s liquid state and avoid compressing the BOG needlessly. Therefore, it is better to recondense the BOG first and then pump it to PS rather than compressing it all the way from P1 to PS. Now, recondensing the BOG at P1 is impossible in the absence of an external coolant. Hence, we must compress it to some P2 > P1 and then recondense it using any available LNG streams. To this end, we have two LNG streams. One is the LNG at (p1, h1, t1). We can use it as such, if it can recondense all of the BOG while remaining liquid. If it cannot, then we must pump it to some pL > p1 to make it a viable coolant (cooling potential increases with pressure). We call this pressurized LNG LP-LNG at (pL, hL, tL) on the P−H diagram. This is our primary coolant for recondensing the BOG at P2. The second cold stream, called HP-LNG, is the cold ungasified version of the send-out NG. HP-LNG at (PS, hH, tH) is a mixture of two LNG streams. One is the recondensed BOG (or RBOG), and the other is LP-LNG after it has recondensed the BOG. To lower the energy

(3) the pumps and compressors are isentropic with 100% efficiency; (4) the heat exchangers are countercurrent with a minimum temperature approach of zero; (5) the pressure drops through the equipment and pipelines are zero; and (6) the net positive suction head (NPSH) for all pumps is zero.

3. BOG RECONDENSATION SCHEMES Let Figure 2 represent the pressure−enthalpy (P−H) diagram for the LNG feed. The thick solid line in green represents the

Figure 2. P−H diagrams for one-stage compression schemes.

vapor−liquid (V−L) phase boundary, and the thin solid lines in yellow represent the isentropic curves. In this article, we use uppercase labels for the gaseous streams and lowercase labels for the liquid streams. Thus, let (P1, H1, T1) on this P−H diagram represent the BOG exiting the tanks, let (p1, h1, t1)

Figure 3. Process flow diagram for scheme 1.1. C

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research consumption in reliquefaction, we must reduce the BOG compression energy (thus P2) as much as possible. To recondense the superheated BOG from (P2, H2, T2), we must first cool it to its dew-point temperature (DPT2) and then recondense it to its bubble-point temperature (BPT2). Therefore, we heat a fraction y1 of HP-LNG from (pS, hH, tH) to (pS, hH1, tH1) to cool the BOG from (P2, H2, T2) to (P2, H21, tH). This heat exchange must respect the relations y1(hH1 − hH) = β(H2 − H21)

(1)

DPT2 ≤ t H ≤ t H1 ≤ T2

(2)

(1 − β)[h(pL , t L2) − h(pL , t L1)] = β[H(P2 , t H) − H(P2 , DPT2)]

To recondense the BOG over the interval [DPT2, BPT2], we use LP-LNG at (pL, tL). Therefore (1 − β)[h(pL , t L1) − hL] (6)

(8)

BPT2 ≥ t L

(9)

(14)

(15)

This completes our scheme 1.1 (first of the one-stage compression schemes) for recondensing the BOG indirectly with one-stage compression. To determine the minimum PC, we can solve the following nonlinear programming (NLP) problem (M1.1) with P2 and PL as independent variables. M1.1: Minimize PC (eq 15) subject to eqs 1−12 and 14. We consider the Korean regasification terminal studied by Park et al.4 for illustration. Although we study the terminal in full detail later for all values of β, for now, we use the specific case of β = 0.15 just to illustrate the effectiveness of our regasification schemes and the impact of process modifications. For β = 0.15 and psw = 0.0049, scheme 1.1 requires a power of 12.33 kW per tonne of send-out gas. 4.2. Direct Recondensation. In scheme 1.1, LP-LNG and RBOG are eventually mixed to form HP-LNG. This is in contrast to the prevalent practice of recondensing BOG by directly mixing it with LP-LNG.4 Since BOG is a result of vaporization, it must be slightly lighter in composition than the LP-LNG. From the BOG perspective, mixing it with the slightly heavier LP-LNG would increase its bubble-point temperature and decrease its latent heat of vaporization. As a result, its condensation load would reduce, which would essentially mean a lower P2 and lower BOG compression energy. From a practical viewpoint, direct recondensation eliminates the resistances inherent in indirect heat transfer. However, it requires

(5)

t L1 ≤ bpt L

(13)

− h(pL , t L2)]

(4)

(7)

DPT2 ≥ t L1

− H(P2 , BPT2)] + (1 − β)[hL − h1 + h(PS , t L3)

(3)

t L ≤ t L1 ≤ DPT2

(12)

PC = psw [hS − h(PS , t H1)] + β[H2 − H1 + h(PS , T22)

where h(p,t) denotes the LNG enthalpy at (p, t) and bptL is the bubble-point temperature of LP-LNG at pL. Clearly, the following relation must hold for complete recondensation of BOG

= β[H(P2 , DPT2) − h(P2 , BPT2)]

t L2 ≤ bpt L

Next, pumping the seawater to regasify the HP-LNG from (PS, tH1) to (PS, TS) requires psw[hS − h(PS, tH1)] kW. Adding the power required for compressing the BOG from P1 to P2 in one isentropic stage and also pumping the LNG streams isentropically, we obtain the total power consumption for the terminal as

where H(P,T) denotes the BOG enthalpy at (P, T). Since LPLNG at (pL, tL) is our primary coolant and must not be allowed to vaporize, its maximum cooling capacity is

Q LP ≥ Q BOG

(11)

hH = (1 − β)h(PS , t L3) + βh(PS , T22)

4. ONE-STAGE COMPRESSION Two possibilities exist for recondensing the BOG at P2. We can mix it directly with LP-LNG (direct recondensation) as in existing terminals, or we can condense it indirectly using LPLNG without mixing (indirect recondensation). Figure 2 shows the P−H diagram for all of the one-stage compression schemes discussed below. 4.1. Indirect Recondensation. We first exhaust the cooling duty of HP-LNG (y1 = 1 in eq 1) to cool the BOG from (P2, T2) to (P2, tH). Then, we cool the BOG to DPT2 and recondense it over the interval [DPT2, BPT2] in a series of countercurrent heat exchangers with LP-LNG (Figure 3). This requires the following cooling and condensation duty in total

Q LP = (1 − β)[h(pL , bpt L) − hL]

t L1 ≤ t L2 ≤ t H

where tL2 is the exit temperature of the LP-LNG. Equations 8 and 12 prevent the LP-LNG from vaporizing. Equations 7 and 13 provide the same upper bound for tL1. We now pump the LP-LNG at (pL, tL2) from the above cooler to (PS, tL3), pump the RBOG from (P2, BPT2) to (PS, T22), and then mix the two to get HP-LNG at (PS, hH, tH). Hence

Equation 2 assumes DPT2 ≤ tH, which holds for low P2, because HP-LNG is a mixture of RBOG and heated LP-LNG. Throughout this work, we exclude the possibility that HP-LNG can recondense the BOG and limit its use for sensible cooling only. We can compress the BOG from (P1, H1, T1) to (P2, H2, T2) in either one or more stages. For simplicity, we exclude three or more stages, but the methodology extends easily.

Q BOG = β[H(P2 , t H) − H(P2 , BPT2)]

(10)

where tL1 is the exit temperature of LP-LNG from the recondenser. To cool the BOG to DPT2 before recondensation, we use the LP-LNG at (pL, tL1) exiting the recondenser

P2 = pL

(16)

(1 − β)hL + β[H(P2 , t H)] = h(P2 , tR )

(17)

where tR is the recondenser temperature. Now, we can pressurize the recondensed liquid to PS to get HP-LNG and D

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

Figure 4. Process flow diagram for scheme 1.2.

Figure 5. Process flow diagram for scheme 1.3.

4.3. Direct Recondensation with HP Precooling before Compression. In schemes 1.1 and 1.2, we compressed the BOG from the tanks directly. We can reduce the compression work, if we cool the BOG before compressing it (Figure 5). Hence, we now use a second portion (y2) of the HP-LNG to precool the BOG to tH, in addition to that (y1) already used to postcool the BOG after compression. For this precooling, we have

then follow the remaining operations in scheme 1.1. Then, the total power consumption is PC = psw [hS − h(PS , t H1)] + β[H2 − H1] + (1 − β)[hL − h1] + [h(PS , t H) − h(P2 , t R )]

(18)

For this scheme (Figure 4), we can solve NLP problem M1.2 with P2 as the only independent variable to obtain the minimum power consumption. M1.2: Minimize PC (eq 18) subject to eqs 1, 2, 16, and 17. For the case-study scenario defined in scheme 1.1, scheme 1.2 requires 2.11% lower power. This explains why all existing regasification terminals use direct instead of indirect recondensation. E

β[H1 − H(P1 , t H)] = y2 (hH2 − hH)

(19)

DPT1 ≤ t H ≤ t H2 ≤ T1

(20)

y1 + y2 = 1

(21)

0 ≤ y1 , y2 ≤ 1

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

Figure 6. Process flow diagram for scheme 1.4.

The two HP-LNG portions exiting the precooler and the postcooler can now be recombined to form the feed to regasification. Therefore y1hH1 + y2 hH2 = h(PS , t H12)

PC = psw [hS − h(PS , t H12)] + β[H2 − H(P1 , t H)]

T11 ≥ DPT1

(27)

(1 − x)(1 − β)hL + βH(P2 , t H) = [β + (1 − x)(1 − β)]h(P2 , t R )

(28)

(23)

The liquid at (P2, tR) from the recondenser should now be pressurized to (pM, tMR) to mix with the MP-LNG at (pM, tH). Thus

In this case, we can solve NLP problem M1.3 to minimize the power usage. M1.3: Minimize PC (eq 23) subject to eqs 1, 2, 16, 17, and 19−22). M1.3 has P2 and y1 as independent variables. For β = 0.15, this scheme uses 3.98% lower power than scheme 1.2. To the best of our knowledge, no work in the literature has considered the option of using HP-LNG to precool the BOG. 4.4. Direct Recondensation with HP and MP Precooling before Compression. In scheme 1.3, HP-LNG could precool the BOG only to tH. Since tH ≥ DPT1, we can reduce the compression work further by cooling the BOG to T11 such that DPT1 ≤ T11 < tH. LP-LNG is the only viable option for this cooling, and we can use a fraction x (0 ≤ x ≤ 1) of the LP-LNG to achieve this. However, this fraction would reduce the muchneeded flow of LP-LNG to the recondenser, and hence increase P2. Thus, we must minimize x to keep P2 low. For the minimum value of x, the LP-LNG must exit the precooler at tH (zero-minimum approach). However, tH might exceed bptL, and the LP-LNG might vaporize. Hence, to maintain the LPLNG exiting the precooler as a liquid, we pump it from (pL, tL) to some (pM > pL, tM) so that bptM ≥ tH. We call this new pressurized fraction of LP-LNG medium-pressure LNG (MPLNG). We can now write

h(pM , tMC) = [β + (1 − x)(1 − β)]h(pM , tMR ) + x(1 − β)h(pM , t H)

(29)

where tMC is the temperature of the resultant mixture. This mixture is now pumped to (PS, tH) to form HP-LNG. With the other terminal operations as in the last scheme, we obtain the total power consumption as PC = psw [hS − h(PS , t H12)] + β[H2 − H(P1 , T11)] + (1 − β)(hL − h1) + x(1 − β)[h(pM , tM) − hL] + [β + (1 − x)(1 − β)][h(pM , tMR ) − h(P2 , t R )] + [h(PS , t H) − h(pM , tMC)]

(30)

This gives us NLP problem M1.4 with four independent variables (P2, pM, T11, and y1) to obtain the minimum power consumption. M1.4: Minimize PC (eq 30) subject to eqs 1, 2, 16, 19−22, and 24−29. Again, for the case-study scenario, scheme 1.4 (Figure 6) uses 6.13% lower power than scheme 1.3, which is the effect of MP precooling alone. The use of MP-LNG and MP precooling is absent in the regasification terminals known to us and is a novel option proposed in this study. Overall, precooling the BOG with both HP-LNG and MP-LNG reduced the power

x(1 − β)[h(pM , t H) − h(pM , tM)] = β[H(P1, t H) − H(P1, T11)] (24)

T11 ≥ tM

(26)

where bptM is the bubble-point temperature of MP-LNG at pM. Equation 27 is necessary to avoid liquid in the BOG entering the compression stage. Furthermore, we must rewrite eq 17 to use the decreased LP-LNG flow for recondensation

(22)

Then, the total power requirement for this scheme (Figure 5) is

+ (1 − β)[hL − h1] + [h(PS , t H) − h(P2 , t R )]

tM ≤ t H ≤ bpt M

(25) F

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research consumption by 9.86% compared to scheme 1.2, which has no precooling. All of the schemes considered so far assumed single-stage compression. If P2/P1 exceeds 4−5, however, then two or more stages of compression become necessary. Therefore, we now consider two stages.

β[H(Pi , Ti ) − H(Pi , t H)] = y3 [h(PS , t H3) − hH]

(31)

DPTi ≤ t H ≤ t H3 ≤ Ti

(32)

y1 + y2 + y3 = 1

5. TWO-STAGE COMPRESSION With two compression stages, interstage cooling becomes possible. Therefore, we now modify scheme 1.4 by introducing a second stage of compression with intercooling. Figure 7 shows the P−H diagram for the two-stage compression schemes discussed below.

0 ≤ y1 , y2 , y3 ≤ 1

y1hH1 + y2 hH2 + y3 hH3 = hH13

(33) (34)

where Ti is the BOG temperature from stage 1, y3 is the portion of the HP-LNG entering the intercooler, tH3 is its exit temperature, and hH13 is the enthalpy of the regasifier feed. The other terminal operations are the same as in scheme 1.4. The total power consumption for this scheme (Figure 8) is given by PC = psw (hS − hH13) + β[Hi(Pi , Ti ) − H(P1 , T11)] + β[H2 − H(Pi , t H)] + (1 − β)(hL − h1) + x(1 − β)[h(pM , tM) − hL] + [β + (1 − x)(1 − β)] [h(pM , tMR ) − h(P2 , t R )] + [h(PS , t H) − h(pM , tMC)] (35)

In this case, NLP problem M2.1 with six independent variables (P2, Pi, pM, T11, y1, and y2) gives the minimum power consumption. M2.1: Minimize PC (eq 35) subject to eqs 1, 2, 16, 19, 20, 24−29, and 31−34. Two-stage compression reduced the power consumption by 1.47% compared to scheme 1.4 for the case-study scenario. 5.2. Direct Recondensation, HP and MP Precooling, and HP and MP Intercooling. A logical extension of the last scheme is to add MP-LNG intercooling. As discussed earlier, using the MP-LNG for this would decrease the amount of LPLNG available for recondensation, and P2 would increase accordingly. We thus use a fraction z (0 ≤ z ≤ 1) of MP-LNG for this additional intercooling and heat it to tH. Recall that the MP-LNG is a fraction x of the LP-LNG pressurized from pL to pM. For this case, we have

Figure 7. P−H diagrams for two-stage compression schemes.

5.1. Direct Recondensation, HP and MP Precooling, and HP Intercooling. Let us compress the BOG from P1 to Pi in stage 1 and from Pi to P2 in the second stage. Also, we use HP-LNG to intercool the BOG from stage 1 to tH before it enters stage 2 (Figure 8). Then, we have

Figure 8. Process flow diagram for scheme 2.1. G

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

Figure 9. Process flow diagram for scheme 2.2.

Figure 10. Process flow diagram for scheme 2.3.

The resultant mixture is then pressurized to (PS, tH) and regasified as described in the previous cases. The total power requirement for this case is

zx(1 − β)[h(pM , t H) − h(pM , tM)] = β[H(Pi , t H) − H(Pi , Ti1)] (36)

tM ≤ Ti1 ≤ t H ≤ bpt M

(37)

Ti1 ≥ DPTi

(38)

PC = psw (hS − hH13) + β[H(Pi , Ti ) − H(P1 , T11)] + β[H2 − H(Pi , Ti1)] + (1 − β)(hL − h1) + x(1 − β)[h(pM , tM) − hL] + [β + (1 − x)(1 − β)]

(1 − z)x(1 − β)[h(pM , t H) − h(pM , tM)] = β[H(P1 , t H) − H(P1 , T11)]

[h(pM , tMR ) − h(P2 , t R )] + [h(PS , t H) − h(pM , tMC)]

(39)

(41)

where the BOG enters the second compression at Ti1 and DPTi is its dew-point temperature at Pi. The two MP-LNG fractions from MP precooling and MP intercooling are mixed with the recondensed liquid pumped to pM. The energy balance for this mixing is

For this scheme (Figure 9), we have NLP problem M2.2 with seven independent variables (P2, Pi, pM, T11, Ti1, y1, and y2) to obtain the minimum power consumption. M2.2: Minimize PC (eq 41) subject to eqs 1, 2, 16, 19, 20, 25−28, 31−34, and 36−40. For the case-study scenario, MP intercooling was not beneficial, as it increased power consumption compared to that of scheme 2.1. Thus, the increased process complexity from MP intercooling is not justified.

h(pM , tMC) = [β + (1 − x)(1 − β)]h(pM , tMR ) + zx(1 − β)h(pM , t H) + (1 − z)x(1 − β)h(pM , t H) (40) H

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research 5.3. Two-Stage Direct Recondensation, HP and MP Precooling. In schemes 2.1 and 2.2, we condensed the BOG after two stages of compression. Having two compression stages allows the freedom to recondense a part of the BOG after the first stage and then recondense the rest after the second stage using the liquid from the first recondenser. This gives a scheme (Figure 10) with two-stage recondensation. In this scheme, all of the BOG is compressed to a lower pressure than before, and only the uncondensed gas is compressed to a higher pressure. In addition, the saturated vapor feed to the second compression stage reduces the compression work. Figure 11 shows the P−H diagram for this scheme. In this

(Li + Vi )h(pM , tMR ) + x(1 − β)h(pM , t H) = h(pM , tMC) (48)

The remaining steps are as in scheme 2.1, and the total power use is PC = psw (hS − hH13) + β[H(Pi , Ti ) − H(P1 , T11)] + Vi [H(P2 , TR3) − H(Pi , t R1)] + (1 − β)(hL − h1) + x(1 − β)[h(pM , tM) − hL] + Li[h(P2 , t R2) − h(Pi , t R1)] + (Li + Vi )[h(pM , tMR ) − h(P2 , t R )] + [h(PS , t H) − h(pM , tMC)]

(49)

Finally, we solve NLP problem M2.3 with six independent variables (P2, Pi, pM, T11, y1, and y2) to obtain the minimum power consumption. M2.3: Minimize PC (eq 49) subject to eqs 1, 2, 19, 20, 24−27, 31−34, and 42−48. This scheme (Figure 10) uses 18.77% lower power than scheme 1.2 for the case-study scenario, and it is the best of all schemes. This novel and thermodynamically superior regasification scheme is the most important contribution of this work. More than two condensation stages would reduce power consumption even further, but at the expense of process complexity and capital cost.

6. SOLUTION STRATEGY Solving NLP problems M1.1−M2.3 in an algebraic environment such as GAMS is inherently approximate, as obtaining accurate analytical expressions for the various thermophysical properties and transformation steps is a challenge. Therefore, we preferred to use an enumerative but ingenious approach based on rigorous process simulations using Aspen HYSYS with the Peng−Robinson model from the HYSYS databank as the fluid package. We carefully prioritized the optimization variables for each scheme and used physical insights to fix all but one or two of these variables. We then enumerated the combinations of these variables using Case Study in Aspen HYSYS to find the optimal solution. We now describe our procedure separately for each scheme and present the numerical results for the illustrative case study with β = 0.15, as was used in sections 4 and 5. 6.1. Scheme 1.1. M1.1 has P2 and pL as its two independent variables. Consider varying P2 for a given pL while ensuring full BOG recondensation (the BOG leaves the recondenser at BPT2). Because power consumption increases with P2, we decrease P2. The minimum acceptable P2 is such that the minimum of the three end approaches in the LP-LNG cooler and recondenser becomes zero. We can compute this P2 value easily using Adjust and Spreadsheet in Aspen HYSYS. Thus, we obtained the minimum P2 for the given pL. Next, we used Case Study in Aspen HYSYS to simulate the above procedure for various pL values. For very low pL values, we observed that the LP-LNG exiting the LP cooler vaporized, which is not acceptable. Therefore, we considered only the acceptable pL values. Figure 12a shows plots for the minimum P2 and corresponding total power as functions of pL for the illustrative case study. Values of pL ≤ 5.5 bar are clearly infeasible, and the total power increases slightly for pL ≥ 5.5 bar. Therefore, pL = 5.50 bar, P2 = 5.90 bar, and power =12.33 kW represent the optimal solution for M1.1.

Figure 11. P−H diagram for scheme 2.3.

scheme, all constraints except the following for the first recondenser remain the same as in scheme 2.1 (1 − x)(1 − β)hL + βH(Pi , t H) = Lih(Pi , t R1) + VH i (Pi , t R1) (42)

Li + Vi = (1 − x)(1 − β) + β

(43)

PL = Pi

(44)

The liquid from the first recondenser at (Pi, tR1) is pumped to (P2, tR2), whereas the uncondensed BOG is compressed from (Pi, tR1) to (P2, TR3). A straightforward extension from the single-stage recondensation schemes is to use an HP cooler after the second compression stage and before the second recondensation stage to cool the compressed BOG from (P2, TR3) to (P2, tH). For this purpose, we use a fraction y3 of HPLNG, which is heated from (pS, hH, tH) to (pS, hH3, tH3). The corresponding energy balance and temperature approach constraints are β[H(P2 , TR3) − H(P2 , t H)] = y3 (hH3 − hH)

(45)

DPT2 ≤ t H ≤ t H3 ≤ TR3

(46)

The liquid from the first recondenser pumped to (P2, tR2) and the BOG from the HP cooler are mixed in the second recondenser to condense all of the BOG. In other words Lih(P2 , tR 2) + VH i (P2 , t H) = (Li + Vi )h(P2 , t R )

(47)

Subsequently, the fully condensed liquid from the second recondenser is pressurized to (pM, tMR) before being mixed with the MP-LNG from the MP cooler I

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

Figure 12. Plots for the case study. J

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research 6.2. Scheme 1.2. P2 is the only variable in M1.2. Therefore, we simply used Case Study to vary pL = P2. Figure 12b shows plots of the vapor flow from the recondenser and the total power consumption versus P2 for the illustrative case study. P2 ≤ 5.58 bar is unacceptable because of incomplete recondensation, and total power increases with P2. Therefore, pL = P2 = 5.58 bar and power =12.07 kW is the optimal solution for M1.2. Alternatively, we could use Adjust to find the lowest feasible P2 value directly. 6.3. Scheme 1.3. P2 and y1 are the two independent variables in M1.3. We fixed y1 by using just enough HP-LNG to ensure zero approach in the HP cooler. We did this automatically in Aspen HYSYS to reduce M1.3 to a singlevariable NLP problem. Then, using the procedure for M1.2, we obtained the plots of vapor flow from the recondenser and total power versus P2 shown in Figure 12c. Therefore, P2 = 5.58 bar, y1 = 0.094, and power =11.59 kW is the optimal solution for M1.3. 6.4. Scheme 1.4. M1.4 has two more variables in addition to P2 and y1, namely, pM and T11. We fixed a value of T11 such that T11 ∈ [DPT1, tH]. As in scheme 1.3, we fixed y1 to force the minimum approach to zero in the HP cooler. Then, we fully utilized the cooling potential of the MP-LNG by fixing pM such that the MP-LNG exits the MP cooler at tH = bptM. This is the optimal pM, because the MP-LNG in this case is similar to the LP-LNG in scheme 1.1. As pM increases beyond its minimum feasible value, the power increases slightly. Now, we can find the lowest P2 that recondenses the BOG fully by using Adjust as in scheme 1.2. This fixes all variables for any T11 value. By repeating this procedure for many values of T11 (Figure 12d), we found that the power consumption decreases with T11. Hence, T11 = DPT1 gives the minimum total power for M1.4. 6.5. Scheme 2.1. M2.1 has Pi, P2, pM, T11, y1, and y2 as the independent variables. We fixed a Pi > P1. As for scheme 1.4, we fixed y1 and y2 to make the minimum approaches zero in the HP coolers, pM to get the MP-LNG to exit the MP cooler at tH = bptM, and T11 = DPT1. We then found the lowest value P2 as in scheme 1.4. This fixes all variables for any Pi. We repeated this procedure to obtain a plot of total power versus Pi, as shown in Figure 12e. Table 2 (below) gives the optimal Pi values and other variables. 6.6. Scheme 2.2. M2.2 has seven independent variables, namely, Pi, P2, pM, Ti1, T11, y1, and y2. We fixed a Pi > P1 and a Ti1 ∈ [DPTi, tH]. As in scheme 2.1, we fixed y1, y2, pM, T11, and P2. Then, we used Case Study to enumerate all combinations of Pi and Ti1. Figure 12f shows that the total power and P2 increase as Ti1 decreases and that Ti1 = tH gives the minimum power. This result holds for all Pi. For Ti1 = tH, scheme 2.2 reduces to scheme 2.1; hence, M2.2 has the same solution as M2.1. 6.7. Scheme 2.3. M2.3 has Pi, P2, pM, T11, y1, and y2 as independent variables. We fixed a Pi > P1. We also fixed y1, y2, pM, T11, and P2 as in scheme 2.2. This fixed all variables for any Pi. We repeated this procedure for various Pi to obtain Figure 12g and the optimal Pi. 6.8. Scheme of Park et al.4 This scheme has Pi, P2, and y1 as independent variables. For a fixed Pi, we fixed y1 and P2 as before. Then, we enumerated Pi to get the minimum power. Figure 12h shows the resulting plot.

1 lists the relevant composition data and process conditions for LNG, BOG, and send-out NG. We solved NLP problems Table 1. Stream Data for the Regasification Terminal from Park et al.4 pressure (bar)

temperature (°C)

LNG

1.12

−160

BOG sendout NG

1.12 74.53

−120 0

stream

composition 89.26% CH4, 8.64% C2H6, 1.44% C3H8 0.35% n-C4H10, 0.04% N2, 0.27% iC4H10 99.37% CH4, 0.04% C2H6, 0.59% N2 90.63% CH4, 7.50% C2H6, 1.25% C3H8 0.30% n-C4H10, 0.08% N2, 0.23% iC4H10

M1.1−M2.3 for this terminal for 0.00 ≤ β ≤ 0.15. We varied β just to illustrate its effects for a given BOG composition. The trend should hold for any other BOG, as our model is perfectly valid for any BOG composition. Table 2 reports the minimum power and optimal operating conditions for each scheme. As expected, the BOG compression pressure and total power increase with β for all schemes. Figure 13 compares the power as a function of β for all of our proposed schemes as well as that proposed by Park et al. (Figure 14).4 To ensure a fair comparison, we simulated the scheme proposed by Park et al.4 using our assumptions and optimized the results for each β value. Consider scheme 1.1 with indirect recondensation. The optimal solution occurs at PL and P2 values that satisfy tL1 = DPT2, tL2 = bptL, and QBOG = QLP for any β. The first condition ensures a zero minimum temperature approach for both cooling and condensation, and the second forces the LP-LNG exiting the cooler to be a saturated liquid. Table 2 suggests pL < P2 for all optimal solutions. For scheme 1.2 with direct recondensation, the minimum power corresponds to the lowest P2 = pL at which the mixing of the LP-LNG and the BOG gives a saturated liquid. As β increases, eq 16 is satisfied at higher and higher P2 values. Thus, the power increases with β. In comparison to the indirect recondensation, the direct mixing of the BOG with the heavier LP-LNG reduces the latent heat of the BOG. Therefore, P2 and the power demand are slightly lower for direct versus indirect recondensation. As an illustration, consider β = 0.15 and P2 = 5.58 bar. At this pressure, the BPTs of LP-LNG, BOG, and the BOG−LNG mixture are −133.77, −136.96, and −134.24 °C, respectively. In other words, direct mixing increases the BPT of the BOG by roughly 2.72 °C and decreases that of the LP-LNG by 0.47 °C. However, the positive effect of the former (in reducing the BOG recondensation load) outweighs the negative effect of the latter (in reducing the cooling potential of the LP-LNG). The net effect is to reduce P2 by 5.4% and, thus, the power by 2.1%. Given the other practical advantages of direct recondensation mentioned before, the prevalent practice of direct recondensation at LNG terminals is well justified. The addition of HP precooling to direct recondensation as in scheme 1.3 yields a significant reduction in power, as the lower BOG temperatures reduce the compression work without any change in the optimal recondensation pressure from scheme 1.2. This is particularly useful for higher β values. For example, for β = 0.15, the power is reduced by 3.98% . Therefore, it is clearly advisable to precool the BOG with the abundant HPLNG before the compression stage. To our knowledge, this

7. CASE STUDY To evaluate our schemes, we consider the Pyeongtaek LNG receiving terminal in South Korea studied by Park et al.4 Table K

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research Table 2. Detailed Optimization Results for Various BOG Generation Scenarios in the Terminal Case Study pressure (bar) β

a

Pi

P2

0.00 0.03 0.06 0.09 0.12 0.15

− − − − − −

1.12 1.59 2.32 3.26 4.45 5.90

0.00 0.03 0.06 0.09 0.12 0.15

− − − − − −

1.12 1.55 2.22 3.09 4.20 5.58

0.00 0.03 0.06 0.09 0.12 0.15

− − − − − −

1.12 1.55 2.22 3.09 4.20 5.58

0.00 0.03 0.06 0.09 0.12 0.15

− − − − − −

1.12 1.56 2.28 3.28 4.65 6.44

0.00 0.03 0.06 0.09 0.12 0.15

1.12 1.40 2.00 2.50 3.30 4.30

1.12 1.56 2.28 3.28 4.65 6.44

0.00 0.03 0.06 0.09 0.12 0.15

1.12 1.40 2.00 2.50 3.30 4.30

1.12 1.56 2.28 3.28 4.65 6.44

0.00 0.03 0.06 0.09 0.12 0.15

1.12 1.30 1.60 2.10 2.60 3.30

1.12 1.57 2.31 3.36 4.86 6.89

0.00 0.03 0.06 0.09 0.12 0.15

1.12 1.12 1.14 1.41 1.80 2.20

1.12 1.55 2.22 3.09 4.20 5.58

pL

temperature (°C) pM

T11

power (kW)

Ti1

pump

compressor

SW pump

Scheme 1.1: Indirect Recondensation − 4.45 0.00 0.95 − 4.54 0.24 0.93 − 4.62 1.03 0.90 − 4.70 2.38 0.87 − 4.79 4.27 0.84 − 4.86 6.67 0.80 Scheme 1.2: Direct Recondensation 1.12 − − − 4.45 0.00 0.95 1.55 − − − 4.53 0.22 0.93 2.22 − − − 4.62 0.97 0.90 3.09 − − − 4.70 2.25 0.87 4.20 − − − 4.79 4.06 0.84 5.58 − − − 4.87 6.40 0.80 Scheme 1.3: Direct Recondensation with HP Precooling 1.12 − − − 4.45 0.00 0.95 1.55 − − − 4.53 0.17 0.93 2.22 − − − 4.62 0.79 0.90 3.09 − − − 4.70 1.91 0.87 4.20 − − − 4.79 3.61 0.84 5.58 − − − 4.87 5.92 0.81 Scheme 1.4: Direct Recondensation with HP and MP Precooling 1.12 − − − 4.45 0.00 0.95 1.56 1.84 −159.20 − 4.54 0.16 0.93 2.28 2.68 −159.20 − 4.63 0.74 0.90 3.28 3.86 −159.20 − 4.72 1.75 0.87 4.65 5.45 −159.20 − 4.82 3.22 0.84 6.44 7.55 −159.20 − 4.91 5.16 0.81 Scheme 2.1: Two-Stage Recondensation with HP and MP Precooling and HP Intercooling 1.12 − − − 4.45 0.00 0.95 1.56 1.84 −159.20 − 4.54 0.16 0.93 2.28 2.68 −159.20 − 4.63 0.73 0.90 3.28 3.85 −159.20 − 4.72 1.72 0.87 4.65 5.45 −159.20 − 4.82 3.14 0.84 6.44 7.55 −159.20 − 4.91 5.00 0.81 Scheme 2.2: Two-Stage Recondensation with HP and MP Precooling, and Intercooling 1.12 − − − 4.45 0.00 0.95 1.56 1.84 −159.20 −152.40 4.54 0.16 0.93 2.28 2.68 −159.20 −146.70 4.63 0.73 0.90 3.28 3.85 −159.20 −140.60 4.72 1.72 0.87 4.65 5.45 −159.20 −134.20 4.82 3.14 0.84 6.44 7.55 −159.20 −127.60 4.91 5.00 0.81 Scheme 2.3: Two-Stage Recondensation with HP and MP Precooling 1.12 − − − 4.45 0.00 0.95 1.57 1.85 −159.20 − 4.54 0.12 0.93 2.31 2.72 −159.20 − 4.63 0.55 0.90 3.36 3.95 −159.20 − 4.73 1.33 0.87 4.86 5.70 −159.20 − 4.84 2.48 0.85 6.89 8.08 −159.20 − 4.96 4.04 0.81 Scheme from Park et al.4 1.12 − − − 4.45 0.00 0.95 1.55 − − − 4.53 0.17 0.93 2.22 − − − 4.62 0.79 0.90 3.09 − − − 4.70 1.89 0.87 4.20 − − − 4.79 3.45 0.84 5.58 − − − 4.87 5.46 0.81 1.12 1.55 2.30 3.10 4.30 5.50

− − − − − −

− − − − − −

total

reductiona (%)

5.40 5.70 6.55 7.96 9.89 12.33

0.00 −0.33 −1.05 −1.70 −2.10 −2.15

5.40 5.68 6.48 7.82 9.69 12.07

0.00 0.00 0.00 0.00 0.00 0.00

5.40 5.63 6.30 7.49 9.24 11.59

0.00 0.85 2.78 4.30 4.69 3.98

5.40 5.63 6.26 7.34 8.88 10.88

0.00 0.90 3.41 6.17 8.38 9.86

5.40 5.63 6.26 7.31 8.80 10.72

0.00 0.90 3.53 6.56 9.20 11.18

5.40 5.63 6.26 7.31 8.80 10.72

0.00 0.90 3.53 6.56 9.20 11.18

5.40 5.58 6.08 6.93 8.17 9.81

0.00 1.73 6.23 11.38 15.72 18.77

5.40 5.63 6.30 7.46 9.08 11.13

0.00 0.85 2.78 4.62 6.31 7.79

Reduction (%) in power compared to scheme 1.2. L

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research

is clearer as T11 = DPT1 gives the minimum power (Figure 12d). Note that the MP-LNG exits its cooler at bptM. The reduction in power reaffirms the impact of BOG cooling prior to compression. Pumping a portion of the LP-LNG to MPLNG made it possible to limit the increase in P2. We now analyze the two-stage compression schemes (schemes 2.1−2.3). As expected, having two stages with HP intercooling (scheme 2.1) reduces the power demand, particularly for higher β values. The effect of the second stage is negligible for low β values. For instance, for β = 0.15, two-stage compression reduces the power consumption by only about 1.47% compared to that of scheme 1.4. Again, P2 remains unchanged irrespective of the compression stages, and the distribution of work between the two stages is crucial in reducing the power use (Figure 12e). Since the first compression is at a lower temperature as a result of the MP precooling, the optimal pressure ratio is higher for the first stage than for the second for all β (Table 2). Thus, two-stage compression with HP intercooling is beneficial for higher β values. The loss of LP-LNG to BOG precooling and intercooling in scheme 2.2 increases P2 (Figure 12f) and the total power; hence, Ti1 = tH gives the optimal solution. In other words, scheme 2.2 simplifies to scheme 2.1. However, MP intercooling might still benefit real designs with nonideal compression stages. The final scheme demonstrates the significant impact of twostep recondensation. Scheme 2.3 has the lowest power for any β value. For β = 0.15, the total power is 18.77% lower than for scheme 1.2 and 8.49% lower than for scheme 2.1 (Table 2). As explained before, this reduction occurs because only a portion of the BOG is compressed to P2 and the feed BOG to the second compressor is the coldest possible. To judge the effectiveness of our proposed schemes, we simulated and optimized the best regasification scheme (Figure 14) reported by Park et al.4 using the idealistic assumptions in our work. Their scheme consists of direct recondensation, twostage compression, HP intercooling, and HP cooling prior to recondensation. Table 2 and Figure 13 also include the optimized results for this scheme. Clearly, schemes 2.1 and 2.3 need much lower power than their scheme for any β value. As can be seen in Figure 13, the power reduction increases with β.

Figure 13. Power consumption vs β for all schemes.

very useful feature seems to be absent from existing LNG terminal designs. Scheme 1.4 further lowers the power use by achieving additional precooling from the MP-LNG prior to compression. Because of its much lower temperature, the MP-LNG is able to cool the BOG much further toward its DPT, which reduces the compression work significantly. For instance, for β = 0.15, the power consumption is decreased by 9.86% compared to that of scheme 1.2 and by 6.13% compared to that of scheme 1.3. However, as discussed before, the decreased flow of the LPLNG available for recondensation increased P2 by 15.53% for β = 0.15. The reduction in compression work outweighed the power increase due to the higher P2 value. Although this result is consistent across all β values, the effect of precooling is more profound for higher β values. The powerful effect of precooling

Figure 14. BOG recondensation scheme proposed by Park et al.4 M

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research For β = 0.15, the power consumed by scheme 2.3 is 11.86% lower than that for Park et al.4 and that consumed by scheme 2.1 is 3.68% lower. Although the reduction in scheme 2.1 is only due to BOG precooling before compression, that in scheme 2.3 is enhanced by two-stage recondensation.

P = gas-phase pressure, bar p = liquid-phase pressure, bar psw = ratio of seawater pump power to regasification heat duty T = gas-phase temperature, °C t = liquid-phase temperature, °C V = mass flow rate of vapor stream exiting recondenser 1 in scheme 2.3, t/h y = fraction of HP-LNG x = fraction of LP-LNG z = fraction of MP-LNG β = ratio of BOG to send-out natural gas

8. CONCLUSIONS A thermodynamic analysis of BOG reliquefaction at a typical regasification terminal to minimize its power use is presented. Our analysis yielded a series of progressively evolving schemes, where the solution of an NLP problem gave the minimum power consumption. Using a case study from Park et al.,4 we established that direct recondensation is better than indirect recondensation, because the latent heat of vaporization of the BOG decreases as a result of direct mixing with the heavier LNG. Precooling the BOG before compression is extremely useful in reducing power consumption. However, introducing a second stage of recondensation has an even greater impact. Therefore, the scheme with the lowest power demand involved two stages of compression with BOG precooling by HP-LNG and MP-LNG and two stages of recondensation with HP cooling prior to each stage. This scheme reduced the power demand by as much as 18.77% compared to that of the simplest scheme (scheme 1.2) with single-stage compression, singlestage recondensation, and HP cooling before recondensation. It amounts to nearly 11.9% lower power than the best scheme reported by Park et al.4 The idealistic scenarios assumed in this work provide lower bounds for the power used in regasification terminals. Clearly, all schemes need not be cost-effective, and some might be too complex. However, our work provides useful benchmarks for evaluating existing regasification terminals.



Subscripts

1 = streams at inlet pressures P1 or p1 2 = streams at pressure P2 L = streams at pressure pL i = streams at pressure Pi M = streams at pressure pM H = streams at pressure pS (high pressure) S = streams at send-out pressure PS Abbreviations



AUTHOR INFORMATION

NG = natural gas BOG = boil-off gas LNG = liquefied natural gas LP-LNG = low-pressure liquefied natural gas MP-LNG = medium-pressure liquefied natural gas HP-LNG = high-pressure liquefied natural gas NLP = nonlinear programming NPSH = net positive suction head RBOG = recondensed boil-off gas

REFERENCES

(1) Mokhatab, S.; Mak, J. Y.; Valappil, J. V.; Wood, D. A. Dynamic Simulation and Optimization of LNG Plants and Import Terminals. In Handbook of Liquefied Natural Gas; Elsevier: Amsterdam, 2014; Chapter 8, pp 321−357. (2) Hasan, M. M. F.; Zheng, A. M.; Karimi, I. A. Minimizing Boil-Off Losses in Liquefied Natural Gas Transportation. Ind. Eng. Chem. Res. 2009, 48 (21), 9571−9580. (3) Dobrota, Đ.; Lalić, B.; Komar, I. Problem of Boil-off in LNG Supply Chain. Transactions on Maritime Science 2013, 02, 91−100. (4) Park, C.; Song, K.; Lee, S.; Lim, Y.; Han, C. Retrofit design of a boil-off gas handling process in liquefied natural gas receiving terminals. Energy 2012, 44 (1), 69−78. (5) Liu, C.; Zhang, J.; Xu, Q.; Gossage, J. L. ThermodynamicAnalysis-Based Design and Operation for Boil-Off Gas Flare Minimization at LNG Receiving Terminals. Ind. Eng. Chem. Res. 2010, 49 (16), 7412−7420. (6) Gómez, M. R.; Garcia, R. F.; Gómez, J. R.; Carril, J. C. Thermodynamic analysis of a Brayton cycle and Rankine cycle arranged in series exploiting the cold exergy of LNG (liquefied natural gas). Energy 2014, 66, 927−937. (7) Ferreiro García, R.; Carbia Carril, J.; Romero Gomez, J.; Romero Gomez, M. Power plant based on three series Rankine cycles combined with a direct expander using LNG cold as heat sink. Energy Convers. Manage. 2015, 101, 285−294. (8) Szargut, J.; Szczygiel, I. Utilization of the cryogenic exergy of liquid natural gas (LNG) for the production of electricity. Energy 2009, 34 (7), 827−837. (9) Choi, I.-H.; Lee, S.; Seo, Y.; Chang, D. Analysis and optimization of cascade Rankine cycle for liquefied natural gas cold energy recovery. Energy 2013, 61, 179−195. (10) Xue, X.; Guo, C.; Du, X.; Yang, L.; Yang, Y. Thermodynamic analysis and optimization of a two-stage organic Rankine cycle for

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was funded in part under Energy Innovation Research Programme (EIRP) Award NRF2014EWTEIRP003008, administrated by the Energy Market Authority (EMA). The EIRP is a competitive grant call initiative driven by the Energy Innovation Programme Office and funded by the National Research Foundation (NRF) of Singapore. This work was also funded in part by the National University of Singapore through a seed grant (R261-508-001-646/733) for CENGas (Center of Excellence for Natural Gas). We acknowledge AspenTech Inc. for allowing the use of HYSYS under an academic license provided to the National University of Singapore.



NOMENCLATURE BPT = bubble-point temperature of BOG, °C bpt = bubble-point temperature of LNG, °C DPT = dew-point temperature of BOG, °C F = mass flow rate of send-out natural gas, t/h H = gas-phase mass enthalpy, kJ/t h = liquid-phase mass enthalpy, kJ/t L = mass flow rate of liquid stream exiting recondenser 1 in scheme 2.3, t/h N

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Article

Industrial & Engineering Chemistry Research liquefied natural gas cryogenic exergy recovery. Energy 2015, 83, 778− 787. (11) Adom, E.; Islam, S. Z.; Ji, X. Modelling of Boil-Off Gas in LNG Tanks: A Case Study. Int. J. Eng. Technol. 2010, 2 (4), 292−296. (12) Wordu, A. A.; Peterside, B. Estimation of Boil-off-Gas BOG from Refrigerated Vessels in Liquefied Natural Gas Plant. Int. J. Eng. Technol. 2013, 3 (1), 44−49. (13) Querol, E.; Gonzalez-Regueral, B.; García-Torrent, J.; GarcíaMartínez, M. J. Boil off gas (BOG) management in Spanish liquid natural gas (LNG) terminals. Appl. Energy 2010, 87 (11), 3384−3392. (14) Park, C.; Lee, C.-J.; Lim, Y.; Lee, S.; Han, C. Optimization of recirculation operating in liquefied natural gas receiving terminal. J. Taiwan Inst. Chem. Eng. 2010, 41, 482−491. (15) Lee, C.-J.; Lim, Y.; Park, C.; Lee, S.; Han, C. Synthesis of Unloading Operation Procedure for a Mixed Operation of AboveGround and In-Ground Liquefied Natural Gas Storage Tanks Using Dynamic Simulation. Ind. Eng. Chem. Res. 2010, 49 (17), 8219−8226. (16) Li, Y.; Chen, X.; Chein, M.-H. Flexible and cost-effective optimization of BOG (boil-off gas) recondensation process at LNG receiving terminals. Chem. Eng. Res. Des. 2012, 90 (10), 1500−1505. (17) Jang, N.; Shin, M. W.; Choi, S. H.; Yoon, E. S. Dynamic simulation and optimization of the operation of boil-off gas compressors in a liquefied natural gas gasification plant. Korean J. Chem. Eng. 2011, 28 (5), 1166−1171. (18) Shin, M. W.; Shin, D.; Choi, S. H.; Yoon, E. S.; Han, C. Optimization of the Operation of Boil-Off Gas Compressors at a Liquified Natural Gas Gasification Plant. Ind. Eng. Chem. Res. 2007, 46 (20), 6540−6545. (19) Park, S.; Park, C.; Lee, U.; Jung, I.; Na, J.; Kshetrimayum, K. S.; Han, C. Comparative Study of Process Integration and Retrofit Design of a Liquefied Natural Gas (LNG) Regasification Process Based on Exergy Analyses: A Case Study of an LNG Regasification Process in South Korea. Ind. Eng. Chem. Res. 2014, 53 (37), 14366−14376.

O

DOI: 10.1021/acs.iecr.6b01341 Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX