Nationwide, Regional, and Statewide Energy Supply Chain

Sep 5, 2013 - The Energy Information Administration (EIA) projects that the increase ... (57) studied a real-world natural gas supply chain from the g...
1 downloads 0 Views 8MB Size
Subscriber access provided by NEW MEXICO STATE UNIV

Article

Nationwide, Regional, and Statewide Energy Supply Chain Optimization for Natural Gas to Liquid Transportation Fuel (GTL) Systems Josephine Anastasia Elia, Richard C Baliban, and Christodoulos A. Floudas Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/ie401378r • Publication Date (Web): 05 Sep 2013 Downloaded from http://pubs.acs.org on September 14, 2013

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Industrial & Engineering Chemistry Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Nationwide, Regional, and Statewide Energy Supply Chain Optimization for Natural Gas to Liquid Transportation Fuel (GTL) Systems Josephine A. Elia, Richard C. Baliban, and Christodoulos A. Floudas∗ Department of Chemical and Biological Engineering Princeton University, NJ 08544, USA August 14, 2013

Abstract An optimization-based supply chain framework is proposed for the nationwide, regional, and statewide analyses of natural gas to liquids (GTL) systems for the United States. Using optimized GTL refineries of differing capacities (i.e., 1, 5, 10, 50, and 200 thousand barrels per day) and fuel product ratios (i.e., unrestricted, maximization of diesel, maximization of kerosene, and commensurate with the United States demand), the optimal nationwide, regional, and statewide supply chains are obtained by solving a large-scale mixed-integer linear optimization (MILP) model that minimizes the total cost of fuel production. The mathematical formulation includes the locations of natural gas in the United States discretized by county, the delivery locations of fuel products, the transportation costs of every input and output of the refinery, the material balances of each GTL refinery, water resources, electricity requirement/production of the supply chain, and the CO2 sequestration capacities in the United States. Solutions of the proposed MILP optimization model provide useful information for the strategic locations of GTL refineries, the allocations of feedstocks and products in the supply chain, as well as a quantitative basis in evaluating each cost contributing factor. Comprehensive ∗ Corresponding

author; Tel: (609) 258-4595; Fax: (609) 258-0211; E-mail: [email protected].

1

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

analyses on the effects of modifying the geographical scope of the supply chain problem are completed and the economic performances of GTL supply chains in different regions are compared. Results suggest that GTL supply chains can produce highly competitive liquid fuels in United States and the Southwest and Central regions of the United States are the most profitable areas for these supply chains.

1

Introduction

The Energy Information Administration projects that the increase in the transportation sector’s energy demand will be fulfilled by “non-petroleum” derived sources 1 . The total transportation fuel demand in 2010 was 13,466 thousand barrels per day (kBPD), 72.1% of the nationwide petroleum consumption 2 , and this demand is projected to reach 14,540 kBPD by 2035 1 . If these fuels can be produced using domestic feedstocks, the United States’ dependence on crude oil imports will decrease, the nation can enhance its energy independence, and the effects of fluctuating crude oil prices will be reduced in the national energy and fuel market. Recently, a review has highlighted the process design alternatives that can produce liquid hydrocarbon fuels and other fuel additives using coal, biomass, and natural gas feedstocks, or any combination of the three carbon based feedstocks, which are available domestically 3 . Coal is an attractive feedstock due to its abundant supply in the United States and its low delivered cost ($2.0-$2.5/MM Btu) 1 , but its high carbon content may result in high production of CO2 that needs to be vented, sequestered, or recycled 4–8 . Biomass is an important feedstock that can reduce the life cycle greenhouse gas emissions of liquid fuels from the CO2 uptake during cultivation 9–13 . An increase in the national crop production, however, needs to take place to generate sufficient amount of sustainable feedstock for fuel production 14–16 . In addition, processes that convert coal and biomass to liquid hydrocarbon fuels using a synthesis gas (syngas) intermediate require high capital costs and extensive syngas cleaning steps to remove the sulfur- and nitrogen-containing gas components. The gasification and cleaning processes for conventional coal or biomass plants approximately account for 50%-65% of the entire plant cost 17,18 . Technologies that convert natural gas feedstock to liquid transportation fuels have existed and studied for decades, but were not economically favorable in the United States due to the high

2

ACS Paragon Plus Environment

Page 2 of 79

Page 3 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

natural gas prices. Commercial natural gas to liquids (GTL) plants exist around the world such as Shell’s Bintulu GTL plant in Malaysia and Shell’s Pearl GTL plant in Qatar. However, recent discoveries of shale gas in the United States, coupled with the horizontal drilling and fracking technologies that enable shale gas production, have driven the cost of natural gas down and made it a very attractive feedstock for liquid fuels and energy production 1 . Sasol has announced plans to construct the first GTL plant in the United States in Louisiana 19 . In hydrocarbon production, the high hydrogen to carbon ratio in the composition of natural gas is highly advantegous in increasing the yield of carbon in the fuel products. This, in turn, will reduce the capital investment required and mitigate the CO2 production from the process. GTL processes also do not require steps to remove sulfur from the synthesis gas since the pipeline requirement already enforces a very low amount of sulfur in the transported gas. Consequently the resulting liquid fuels have close to zero sulfur content. Various studies have proposed hybrid systems where the advantages of coal, biomass, and natural gas can be combined in a synergistic way 3 . Examples of hybrid systems include coal and natural gas to liquids 20–28 , biomass and natural gas to liquids (BGTL) 29–33 , and coal, biomass, and natural gas to liquids (CBGTL) 5–8,34,35 . However, the development of a single feedstock natural gas to liquids (GTL) refinery can take advantage of certain locations that are abundant in natural gas, but far away from coal and biomass resources. For example, GTL technologies can tap into unutilized stranded gas resources, which make up a significant portion of the United States gas supply 36–39 . Stranded gas is low value natural gas that is co-produced in small amount with oil production (associated gas) or gas in wells that are inaccessible from the current pipeline infrastructure. These resources are currently untapped, or in the case of associated gas, flared to the atmosphere since the its amount is too small to justify building an infrastructure to utilize it 36–39 . GTL technologies that can take advantage of stranded gas resources can bring economic value from adding fuels to the market and reduce the environmental damage from flaring. Efficient management of the upstream and downstream operations of GTL plants will determine their profitability. In the literature, extensive research has been done on biomass supply chains for fuel production to determine biorefinery locations with respect to the distributed nature of biomass resources 3,40–57 . Natural gas supply chain studies generally focus on the operational aspect of natural gas production complexes and the pipeline distribution systems that include re3

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

fineries, compression stations, and mixing of products in the pipelines 58,59 . Hamedi et al. 58 studied a real-world natural gas supply chain from the gas and oil wells, refineries, compression stations, local gas companies, city gate stations, and consumers. They also included a literature review of studies on the operations of natural gas supply chains. Contesse et al. 60 studied the natural gas supply chain system that involves the producers, transportation companies, and local distribution companies to optimize the daily purchasing and transportation contracts between the companies. ¨ Ozelkan et al. 61 proposed a framework that analyzes the design of liquefied natural gas (LNG) terminals for the LNG supply chain. dos Santos et al. 62 used a simulation-based approach to develop a management system for the transportation logistics of natural gas. In the decisions to built new GTL refineries in the United States, it is important to analyze what strategic placements of GTL refineries will maximize their profitability in addition to the operational constraints. The strategic placements of GTL refineries have to ensure that the feedstock can be delivered sustainably to the plant over a long time horizon and that the fuel products can be distributed efficiently to the market. Environmental limitations such as freshwater availabilities and CO2 sequestration capacities in the local region must also be fulfilled. These logistical challenges must be considered simultaneously, and an optimization framework is an important tool that can maximize the profitability of the entire GTL supply chain. Previous studies have been done on a nationwide hybrid feedstock supply chain optimization for the coal, biomass, and natural gas to liquids system where the locations of the hybrid plants are optimized with respect to the feedstocks, fuel demand locations, plant investment costs, and environmental constraints such as water, electricity, and sequestration requirements 40,41 . In this paper, an optimization framework is considered for the nationwide, regional, and statewide GTL supply chains for the United States, using novel GTL processes that have been globally optimized in 63 . The supply chain begins at the natural gas feedstock locations, ends at the demand locations, and is optimized to yield the strategic placements of GTL facilities with minimum overall cost of fuel production. The supply chain also addresses the electricity supply, water supply, and the allocation of sequestered CO2 to injection sites. The rest of the paper is organized as follows.

4

ACS Paragon Plus Environment

Page 4 of 79

Page 5 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

2

Natural Gas Resources

The first set of parameter inputs into the energy supply chain optimization model is the natural gas feedstock parameters. This category of parameters include (i) the locations of natural gas sources, (ii) the availabilities at each location, and (iii) the purchase price of natural gas at each location. The supply chain takes into consideration the geographical variabilities of these parameters, which influence the profitability of the GTL refineries across the country. The same type of GTL refinery in two different locations will have different operational costs due to the logistics and configuration of the upstream and downstream operations. In this study, one generic composition for natural gas is used based on an average composition of natural gas wellhead productions in the United States 63,64 . Note that the uniformity of natural gas feedstock is a reasonable assumption due to the quality requirement of natural gas for pipeline transportation that requires above 90% of methane content and very low sulfur and CO2 contents 65 . The data on natural gas resources across the United States are accessed from various statewide reports of gas wells production and these data are based on the 2010 natural gas production level. Natural gas production from oil and gas wells, as well as shale gas and coal-bed methane are included in the total production level on a county basis. When data on the well-specific production are available, the figures are grouped based on the county of the wells. Table 1 lists the sources of natural gas resources data that are used as inputs in this paper and Figure 1 shows the geographical distribution of natural gas production in the United States. The states that have major natural gas productions are Alabama, Arkansas, California, Colorado, Kansas, Kentucky, Louisiana, Michigan, Mississippi, Montana, New Mexico, Ohio, Oklahoma, Pennsylvania, Texas, Utah, Virginia, West Virginia, and Wyoming. A total of 813 counties in the United States have natural gas production and these counties serve as the beginning points of the GTL supply chain as well as the candidate locations for potential GTL refineries. The breakdown of the number of counties and the 2010 natural gas production level for each state can also be seen in Table 1. The county-specific, statewide natural gas productions are compared with the EIA’s statewide production figures to make sure that the data corroborate each other. The wellhead costs of natural gas are taken from the EIA website, which are discretized by state and include all costs prior to shipment from the lease, including the costs for gathering, compression, and processes to prepare the feeds

5

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

for pipeline quality. Since the county-specific natural gas prices are not available, it is assumed that all natural gas within a state has a uniform price as given by the EIA website. The amount of natural gas that is available for the GTL supply chain is assumed to come from further expansions from the 2010 natural gas production level. In other words, the GTL supply chain does not interrupt the current usage of natural gas in various energy, industrial, or residential sectors. In the supply chain case studies, 25%, 50%, 75%, and 100% expansion levels from the 2010 natural gas production are considered to analyze the amount of liquid fuels that can be produced using the optimized GTL refineries 63 .

3

Gas to Liquids (GTL) Refineries

The design and topology of the GTL refineries that may exist in the supply chain are optimized using the process synthesis methods developed in Baliban et al. 63 . For GTL systems, typical studies focus on fixed process designs and topologies where simulations are used to determine the heat and mass balances for the process 20,21,24,25,32,66–74,74–82 . However, recent developments in process synthesis strategies for hybrid feedstock systems 7,8,12,13,33–35,43,44,63,83–89 have shown that rigorous and computationally efficient methodologies can be developed and used to simultaneously analyze thousands of process designs to extract the optimal process design(s) and establish trade-offs among them. The superstructure optimization method can directly compare the technoeconomic and environmental benefits of GTL processes in a singular mathematical model, simultaneously analyzing several existing or novel processes via a process superstructure to determine the optimal topology which will have either the lowest cost or highest net present value. Combining an environmental objective function with the optimization of refinery designs in a multiobjective optimization framework has also been proposed to measure the tradeoff between economic and environmental performances of energy systems 43,44,86–89 . Yue et al. 86 proposed a life cycle optimization framework based on the functional unit of energy systems, which can be the mass or volumetric amount of the product or in a multi product system, the sum of the properties such as heating value or market value of the product streams. They proposed a four-stage steps of the life cycle assessment, begining with the definition of the scope and boundaries of the study (e.g., wellto-wheel, well-to-gate, etc.), inventory analysis, impact assessment, and interpretation, the last step 6

ACS Paragon Plus Environment

Page 6 of 79

Page 7 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 1: Summary of the states that produce natural gas, the number of counties, and the total production level in 2010. The statewide sources of the natural gas data are included for reference.

State

Region

Number of counties

Total natural gas production (billion SCF/yr)

Alabama

Southeast

21

240.71

Arkansas

Southwest

32

1,085.71

California

Western

27

262.88

Colorado

Central

37

1,917.81

Kansas

Central

56

333.18

Kentucky

Southeast

38

135.15

Louisiana

Southwest

60

2,142.41

Michigan

Midwest

43

143.83

Mississippi

Southeast

37

59.53

Montana

Central

30

105.29

New Mexico

Southwest

9

1,397.35

Ohio

Midwest

45

79.88

Oklahoma

Southwest

63

1,706.70

Pennsylvania

Northeast

33

671.11

Texas

Southwest

197

6,499.12

Utah

Central

10

493.26

Virginia

Northeast

7

147.26

West Virginia

Northeast

48

287.24

Wyoming

Central

20

2,522.94

Total

Nationwide

813

20,231.36

7

Data source

Geological Survey of Alabama State Oil and Gas Board www.gsa.state.al.us Arkansas Oil and Gas Commission www.aogc2.state.ar.us State of California Department of Conservation Division of Oil, Gas & Geothermal Resources www.conservation.ca.gov Colorado Oil and Gas Conservation Commission cogcc.state.co.us Kansas Geological Survey www.kgs.ku.edu Kentucky Geological Survey www.kgs.uky.edu State of Louisiana Department of Natural Resources: Sonris Database sonris.com State of Michigan Department of Environmental Quality www.michigan.gov/deq Mississippi State Oil & Gas Board www.ogb.state.ms.us State of Montana Department of Natural Resources & Conservation dnrc.mt.gov New Mexico Oil & Gas Association www.nmoga.org Ohio Department of Natural Resources www2.ohiodnr.gov Oklahoma Geological Survey www.ogs.ou.edu Pennsylvania Department of Environmental Protection www.depweb.state.pa.us Texas Oil & Gas Association www.txoga.org Utah Department of Natural Resources Division of Oil, Gas & Mining oilgas.ogm.utah.gov Virginia Department of Mines, Minerals and Energy www.dmme.virginia.gov West Virginia Oil and Natural Gas Association www.wvonga.com Wyoming Oil and Gas Conservation Commission wogcc.state.wy.us

ACS Paragon Plus Environment

660

880

LA

1,100

1,320

1,540 Miles

FL

GA MS AL

AR

184004757.000001 - 422199351.000000

82904348.470001 - 184004757.000000

25070843.000001 - 82904348.470000

1.000000 - 25070843.000000

Amount

Natural Gas Sources

422199351.000001 - 1197584292.000000

0 110 220

NM AZ

440

OK

TX

KY TN

KS CO

CA

OR

NV

ID

WY UT

NE

IA

MO

IL

IN

MI

OH

NC SC

PA NJ WV MDDE VA

VTNH NY CTMA MI WI SD

MN ND MT WA

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 79

ME

Industrial & Engineering Chemistry Research

Figure 1: The geographical layout of natural gas resources for the GTL supply chain network in the United States. The figures are based on the 2010 natural gas production on a county basis. There are a total of 19 states that produce natural gas in the United States. The counties highlighted in the figure are also candidate locations for GTL refineries.

8

ACS Paragon Plus Environment

Page 9 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

being coupled with the multiobjective optimization framework. The method has been applied to superstructure optimization of biorefineries that convert cellulosic biomass to gasoline and diesel via gasification 87 , hybrid poplar feedstock to gasoline and diesel range biofuels via fast pyrolysis and hydroprocessing steps 89 , and algae to biodiesel via lipid extraction and algal oil processing with CO2 from power plant flue gas 88 . In Baliban et al. 63 , a rigorous global optimization strategy 34,90–101 is used to mathematically guarantee that the GTL process design selected by the framework will have an overall cost (or profit) that is within a small percentage of the best value possible. The key process alternatives include (i) the inclusion and mathematical modeling of steam reforming of natural gas, auto-thermal reforming of natural gas, direct conversion of natural gas to methanol via partial oxidation, and direct conversion of natural gas to olefins via oxidative coupling as conversion technologies, (ii) the direct usage of natural gas in the fuel combustor unit to provide process heat and in the gas turbine for electricity production, (iii) synthesis gas conversion via Fischer-Tropsch (FT) or methanol synthesis, (iv) methanol conversion via methanol-to-gasoline (MTG) or methanol-to-olefins (MTO), and (v) hydrocarbon upgrading via ZSM-5 zeolite catalysis, olefin oligomerization, or boiling point fractionation and subsequent treatment, (vi) different product compositions (i.e., gasoline, diesel, and kerosene) considered, namely the unrestricted composition (U), maximization of diesel (D), maximization of kerosene (K), and compositions commensurate with the United States demand ratio (R), and (vii) calculations of the life cycle emissions of GTL systems compared to petroleumbased processes and natural gas-based electricity production 63 . The framework includes a simultaneous heat and power integration 6,7 using an optimization-based heat-integration approach 102 and a series of heat engines that can convert waste heat into electricity 5–7 . A comprehensive wastewater treatment network 8 that utilizes a superstructure approach 103–106 to determine the appropriate topology and operating conditions of process units is utilized to minimize wastewater contaminants and freshwater intake. Finally, recycling of gases and CO2 from and to various sections of the refinery (i.e., reformers, Fischer-Tropsch, reverse water gas shirt reactor) is allowed to maximize fuel production and limit CO2 emissions 7,8,34,35 . The life cycle analysis of the refineries is done on a well-to-wheel basis and an environmental constraint is included in the optimization model to limit the GTL refinery emissions to be at most equal to petroleum-based processes or natural gas-based electricity production (see Section 9) 63 . 9

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For this study, five different refinery capacities are considered, namely plants that produce 1 kBPD, 5 kBPD, 10 kBPD, 50 kBPD, and 200 kBPD of total liquid fuels to accommodate geographical locations that vary in natural gas resources. In addition to the varying capacities, four different modes of fuel production apply to the GTL refineries: (i) the unrestricted fuel composition (U) where only gasoline in produced with byproduct LPG, (ii) maximization of diesel (D) where 75 vol% diesel and 25 vol% gasoline are produced, (iii) maximization of kerosene (K) where 75 vol% kerosene and 25 vol% gasoline are produced, and (iv) fuel ratios commensurate with the 2010 United States demand (R) that produce 67 vol% gasoline, 22 vol% diesel, and 11 vol% kerosene 107 . The 20 combinations of capacity and fuel ratios are optimized via the process synthesis method and the outputs include (i) natural gas requirement, (ii) electricity production, (iii) butane requirement, (iv) water requirement, (v) investment cost, (vi) liquid fuel flow rates, (vii) CO2 vented, (viii) byproduct propane/liquefied petroleum gas (LPG) flow rates, and (ix) CO2 sequestered. The results for the 20 GTL refineries in the supply chain can be seen in Table 2 From Table 2 it can be seen that no butane requirement is needed except for the (R) refineries, and (K) refineries do not produce byproduct LPG. Economies of scale are observed as the capacities increase from 1 kBPD to 200 kBPD in the normalized investment cost figures. Electricity is produced in all cases, which excludes the need of consuming additional fossil fuels to produce electricity and add to the overall greenhouse gas emissions.

4

Liquid Fuel Delivery Locations

The 2010 demand data for petroleum based products are extracted from the Energy Information Administration (EIA) website, which are available on a statewide basis. Among the list of products, motor gasoline, distillate fuel oil, and jet fuel represent a large majority of the total product volumes 2 . In this study, the end destinations for the fuel products are the existing oil refineries in the United States to be distributed across the nation or to be blended with petroleum-based fuels before distribution. In other words, the GTL supply chain will take advantage of the hydrocarbon pipeline infrastructure connected to the refineries. Current oil refineries are viable options to send the GTL products to because they can help improve the quality of petroleum-based fuels via blending. The GTL products are high-quality liquid fuels that have close to zero sulfur content 10

ACS Paragon Plus Environment

Page 10 of 79

Page 11 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 2: Overall material balance, electricity output, and total investment cost for the 20 types of GTL refineries that may exist in the supply chain. The inputs to the GTL refineries are natural gas, butane, and water, while the outputs include gasoline, diesel, kerosene, LPG, sequestered CO2 , vented CO2 , and electricity. The four modes of fuel production are (i) unrestricted (U), where only gasoline is produced, (ii) maximization of diesel (D), where gasoline and diesel are produced, (iii) maximization of kerosene (K), where gasoline and kerosene are produced, and (iv) fuel ratios commensurate with the United States demand (R), where all three fuels are produced.

Natural gas (mscf/hr) Butane (kBD) Water (kg/s) Gasoline (kBD) Diesel (kBD) Kerosene (kBD) LPG (kBD) Seq. CO2 (tonne/hr) Vented CO2 (tonne/hr) Electricity Investment cost (MM $) Normalized inv. cost ($/bpd)

Natural gas (mscf/hr) Butane (kBD) Water (kg/s) Gasoline (kBD) Diesel (kBD) Kerosene (kBD) LPG (kBD) Seq. CO2 (tonne/hr) Vented CO2 (tonne/hr) Electricity Investment cost (MM $) Normalized inv. cost ($/bpd)

U-1

U-5

U-10

U-50

0.36 1.92 1.00 0.09 1.36 4.03 -2 138 138127

1.82 9.59 5.00 0.45 6.80 20.14 -12 478 95656

3.61 14.01 10.00 0.90 10.98 40.99 -25 798 79808

18.53 74.21 50.00 4.50 84.64 204.30 -141 3354 67072

K-1

K-5

K-10

K-50

0.36 1.42 0.25 0.75 1.35 4.15 -4 143 142978

1.81 7.09 1.25 3.75 6.76 20.75 -19 501 100208

3.54 19.85 2.50 7.50 10.00 40.14 -32 812 81192

18.31 85.89 12.50 37.50 73.66 214.71 -218 3393 67866

11

GTL Refineries U-200 D-1 73.57 400.25 200.00 19.78 347.18 796.88 -460 11685 58427

0.39 1.80 0.25 0.75 0.04 2.34 3.40 -2 139 138799

GTL Refineries K-200 R-1 77.28 300.19 50.00 150.00 587.99 808.07 -878 11732 58658

0.36 0.03 1.68 0.67 0.22 0.11 0.02 2.02 3.86 -4 145 145227

ACS Paragon Plus Environment

D-5

D-10

D-50

D-200

1.93 9.01 1.25 3.75 0.20 11.70 17.02 -11 464 92756

3.76 17.34 2.50 7.50 0.37 16.53 35.64 -25 800 79997

18.89 90.89 12.50 37.50 1.52 91.88 176.63 -134 3299 65977

73.85 406.92 50.00 150.00 6.06 265.71 706.08 -475 11384 56922

R-5

R-10

R-50

R-200

1.81 0.13 8.42 3.36 1.08 0.56 0.08 10.12 19.29 -19 510 102027

3.62 0.26 16.84 6.72 2.15 1.13 0.16 20.25 38.59 -38 834 83393

18.02 1.14 63.37 33.60 10.77 5.63 0.74 114.88 174.36 -141 3527 70538

72.20 4.29 370.23 134.39 43.10 22.51 3.01 463.60 696.58 -563 12347 61736

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and mixing them with petroleum-derived fuels can reduce the overall sulfur content and meet the imposed fuel standards. To specify the demand amount at each location, the statewide demand data are distributed over these oil refineries. The EIA database provides the maximum operating capacities of United States operating petrochemical refineries 108 . We assumed that the demand amount or the amount of GTL fuels delivered to each oil refinery does not exceed their operating capacities. The mixtures between gasoline, diesel, and kerosene in the delivered fuels are unspecified.

5

Transportation Costs

The modes of transportation for each commodity in the supply chain for both the upstream and downstream operations are defined. Natural gas is transported by pipeline, water is transported by pipeline, and fuel products are transported by a combination of truck, pipeline, and barge 40,41 . Natural gas and liquid fuels transportation costs by pipeline are calculated based on average values on a per mile per volumetric flow rate basis derived from the informational postings of 10 largest natural gas pipeline and 10 largest hydrocarbon pipeline companies in the United States 40 . These postings consist of information on pipeline origin and end points, and usage fees, which can be used to estimate the average transportation costs for natural gas and liquid fuel products. Connections between GTL refinery locations, ports (origin and destination), and demand locations are allowed for transportation by pipeline and barge, and the transportation costs for barge transportation follow the assumptions in Elia et al. 40 where a list of 50 United States ports was extracted from the World Port Source 109 consisting of the 49 highest ranked ports according to the total quantity of imports received in a year (excluding Virgin Islands and Puerto Rico) 110 and 1 port in Anchorage, AK. Port-to-port distances are determined from an online World Ports Distances calculator 111 and the ports that do not have an entry in the World Ports Distances database are listed in Table 7 of Elia et al. 40 , along with their respective reference ports (closest port location with available distance information) and the adjustment distance (i.e., the known distance in nautical miles between the unaccounted ports and their corresponding reference ports). Transportation costs by truck, barge, and water pipeline are calculated using Eq. 1, where DFC is the distance fixed cost, DVC is the distance variable cost, Distance is the distance traveled, and DM is the distance multiplier. Distance is the straight line distance calculated from the latitude and longitude coordinates of two 12

ACS Paragon Plus Environment

Page 12 of 79

Page 13 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

points, and to account for path curvatures, DM is equal to 1.1 for truck transportation and water pipeline 112 . For fuel products transportation by truck, the DFC value is $3.318/bbl and the DVC value is $0.124/bbl-mi; for fuel products transportation by barge, the DFC value is $1.218/bbl and the DVC value is $0.0063/bbl-mi; and for water transportation by pipeline, the DFC value is $0.003/kg and the DVC value is $5E-6/kg-mi 112 . The CO2 pipeline transportation costs are discussed in further detail in Section 8.

Transportation Cost = DFC + DVC · Distance · DM

6

(1)

Electricity

When a GTL refinery produces electricity, it can be placed anywhere in the country and a profit will be taken from selling electricity to the grid. However, when a plant requires electricity, this requirement will add to the electrical power demand for the country and three options are available, namely drawing electricity from the grid, producing electricity from solar resources or wind resources. The 2009 net electricity generation capacity in the United States was reported to be 1025.4 GW in the summer and 1063.8 GW in the winter 113 , and most of this electrical power relies on fossil fuels. It is important to cap the total electricity usage in the supply chain network to prevent excessive expansion of the current power supply and limit the greenhouse gas emissions from power generation. On the single plant approach, heat and power integration is completed for each GTL refinery topology such that the combined heat and power requirement is minimized and electricity is recovered from waste heat 6,7,35,63 . On the supply chain network approach, a total cap on the grid expansion for the whole country is imposed, which is set to be at most 10% of current generation, amounting to 102.5 GW 41 . Electricity generated from renewable sources (i.e., solar and wind) is incorporated into the supply chain in the following way. The investment costs of each electricity option are accounted in the model, and the optimization model will select one or a combination of electricity sources for a given location. Electricity from solar and wind is assumed to be produced from solar and/or wind plants built on the site of the GTL refinery. The maximum capacity is set to 40 MW for a wind plant and 20 MW for a solar plant. Only counties that have adequate solar irradiation and wind 13

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

speed are allowed to have solar and wind power generation. Solar irradiation data are obtained from the National Solar Radiation Data Base (NSRDB) 1991-2005 Update 114 , which contain the longitude and latitude specifications of the station sites with solar irradiation. The longitude and latitude data are used to calculate the distance between solar sites and candidate locations of the GTL refineries. If a county is within 10 miles of a solar station site, and the data indicate that there is adequate irradiation surrounding the county, then a solar plant may be built with the GTL refinery 41 . Wind profiles across the United States are obtained from the National Renewable Energy Laboratory (NREL) Wind Integration Dataset 115 , namely the Eastern Wind Dataset and the Western Wind Datasets. The Western Wind Datasets contain 32,043 locations with the highest wind energy density in the western part of the United States. From the Eastern Wind Datasets, sites from the land based and Midwest ISO (MISO) areas are used in this analysis. The longitude and latitude data of the station sites are used to calculate the distance between wind sites and the candidate locations of the GTL refineries. If a county is within 10 miles of a wind site, and the data indicate that there is adequate wind speed surrounding the county, then a wind farm may be built with the GTL refinery. Note that due to the onsite power generation, the costs for transmission and electrical grid expansion are not included 41 .

7

Water Resources

To prevent over-stressing the regional water resources, the water supply chain is included in the GTL supply chain formulation. On the single plant level, the water network superstructure is optimized simultaneously with the process synthesis framework to minimize the freshwater input and wastewater discharge from the GTL refineries 8,35,63 . On the supply chain level, a regional constraint is included such that the water requirement for any GTL refinery must be satisfied locally within a 5-mile radius. Water usage data in the United States are obtained from the United States Geological Survey (USGS) database 116 , where estimates of freshwater consumption for domestic, agricultural, and industrial sectors are reported for each county. The available water for the GTL supply chain is calculated by taking the minimum of 1.5 times of current industrial use or 15% of total freshwater use in the county to prevent excessive expansion in the infrastructure. The water 14

ACS Paragon Plus Environment

Page 14 of 79

Page 15 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

transportation costs by pipeline is included in the total cost of the GTL supply chain 41 .

8

CO2 Sequestration

The optimized GTL refineries include the selection of CO2 capture and sequestration in the optimal topologies (see Section 3). In the supply chain formulation, it is important to include the constraint that ensures viable locations of sequestration sites with enough capacities to accommodate the amount produced by the GTL refineries. Currently, potential geological storage sites are being explored in the United States and the National Energy Technology Laboratory (NETL) published estimates on the capacities of underground CO2 storage 117 . In this study, these estimates from the NETL 2010 Carbon Sequestration Atlas of the United States and Canada (Third Edition, Atlas III) 117 project are used to determine the end points and capacities of sequestered CO2 . Data from Appendix C of the document 117 are used to calculate the total estimated storage resources (low capacity scenario) by state and the sequestration sites are assumed to be at the center of the state from the purpose of calculating the distance of CO2 transportation. Note that the approach taken in the optimization model can be readily used where more discretized information on storage sites become available. The CO2 to be sequestered exits the GTL refineries at 150 bar and is transported via pipeline to the storage sites. The levelized costs of investment for CO2 pipelines are calculated using a linearization of the pipeline equations in Ogden 118 as functions of the CO2 flow rates and transportation distances, which are calculated explicitly in the model. The cost of injection and levelized cost for new wells also follow the approach in Ogden 118 and are minimized in the energy supply chain optimization model.

9

Life Cycle Analysis

Baliban et al. 63 described the methods of calculating the life cycle greenhouse gas emissions from the optimized GTL refineries. An environmental constraint is included in the process synthesis model to prevent excessive emissions from the process. For each GTL refinery, the total GHG emission target was set to be at most equal to that for petroleum-based production of liquid fuels 15

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

or natural gas-based production of electricity. The GHG emission rates (in kg CO2 eq/s) for the eight major point sources in the refinery are given in Baliban et al. 63 and include (a) acquisition and transportation of the natural gas and butane feeds, (b) transportation and use of the gasoline, diesel, kerosene, and LPG, (c) transportation and sequestration of any CO2 , and (d) venting of any process emissions. In other words, the system boundaries start at the feedstock acquisition and end at the end use of the liquid fuels 35,40,41,63 . The GHG emissions for feedstock acquisition and transportation in (a), product transportation in (b), and CO2 transportation in (c) are calculated from the GREET model for well-to-wheel emissions 119 , with the given average distances traveled for each mode of transportation. The GHG emissions from product use in (b) are calculated assuming that each product will be completely combustion to generate CO2 that is simply vented to the atmosphere. The lifecycle GHG emissions (LGHG) was set to be the sum of the total emissions from each stage of the process. The GHG emissions avoided from liquid fuels (GHGAF) is equivalent to the total energy of fuels produced multiplied by a typical petroleum-based emissions level (i.e., 91.6 kg CO2eq /GJLHV ) while the GHG emissions avoided from electricity (GHGAE) is equivalent to the energy produced by electricity multiplied by a typical natural gas-based emissions level (i.e., 101.3 kg CO2eq /GJ) 120 . Note that since the GTL refineries are multi-product systems, the total fuel production is levelized to the total heating value of gasoline, diesel, and kerosene produced. The GHG emissions index (GHGI) represents the division of LGHG by the sum of GHGAF and GHGAE, and values less than unity are indicative processes with superior lifecycle GHG emissions than current processes. In the refinery topology, the CO2 management includes the option to vent CO2 to the atmosphere, capture and sequester CO2 , and recycling of CO2 into various units in the process (i.e., reformers, reverse water gas shift reactor, Fischer-Tropsch reactors) via the reverse water gas shift reaction to produce CO and liquid fuels subsequently 35,63 . Since the GTL refineries do not consume renewable carbon feedstock, the only negative contribution in the emissions is the CO2 sequestration. Baliban et al. 63 showed that for each of the optimized GTL refinery, the overall life cycle emissions is equal to petroleum-based processes (i.e., GHGI exactly equal to 1, according to the environmental constraint). The majority of the emissions are attributed to the liquid fuels consumption (i.e., 70% of total emissions). Natural gas is a particularly GHG intensive feedstock due to the small amount of methane that is leaked to the atmosphere during extraction from the 16

ACS Paragon Plus Environment

Page 16 of 79

Page 17 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

ground. Nevertheless, it is still economical to develop GTL processes that can have appropriate GHG emissions targets. On the supply chain level, the transportation distances for the feedstocks and products will not be exactly the same as the assumed GREET parameters in Baliban et al. 63 . Some feedstocks and products will travel a greater distance and some will travel less distance, and the GTL refineries in the supply chain will have slight variations on the total life cycle emissions. Note that an environmental constraint on the life cycle emissions of the GTL supply chain may be included in the model formulation. However, since the GTL supply chain does not involve any renewable source and the emissions are comparable to petroleum-derived fuels, this constraint is not included in this study.

10

Energy Supply Chain Optimization Model

The previous subsections detail the parameter inputs required in the energy supply chain optimization model. These parameter inputs include: (i) the location, availabilities, and purchase prices of natural gas feedstock, (ii) the transportation costs of natural gas, water, and fuel products, (iii) investment costs, feedstock requirement, electricity requirement/output, fuel product amounts, and other input-output material balances for the 20 GTL refineries (see Table 2), (iv) locations of fuel product destination (i.e., oil refineries), and (v) price of utilities (i.e., water and electricity). The GTL supply chain optimization model is formulated as a mixed-integer linear optimization (MILP) model and is solved to give the (a) the strategic locations of the GTL refineries on a nationwide, regional, and statewide level, (b) the capacity of the selected GTL refinery at each location, (c) the supply chain topology from the feedstocks to product destination, and (d) the costs breakdown associated with each segment of the supply chain. Continuous variables are used to represent the levelized investment cost at a location (CostlI ), the total electricity requirement (EllT ), electricity produced (EllP ), solar electricity required (EllS ), wind electricity required (EllW ), grid electricity required (EllG ), the water requirement (W Fl ), the CO2 sequestration amount (SFl ), the flow of feedstock from the producing locations to the plant facility (x f ,c,l,m ), the required flow rate of natural gas feedstock to a plant facility (FR f ,l ), the flow of products transported from the GTL refineries to the demand locations (z p,l,c,m ), the flow of 17

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 79

freshwater to the facilities (wc,l ), and the sequestration amount from a facility (seqc,l ). Note that the model allows for electricity input in the supply chain, but since all of the GTL refineries have net electricity output, the EllT , EllS , EllW , and EllG values are zero in all cases. Binary variables (y f ,l,t,q ) represent the selection of a specific GTL refinery at location l with natural gas feedstock ( f ) ∈ FC, capacity t, and fuel product ratio q. The complete variable list can be found in the Notation Section. At most one GTL refinery may be selected at a given candidate loccation, which is specified at the county centroid where natural gas feedstock exists (Eq. 2). Equations 3-5 impose lower and upper limits in the number of selected GTL refineries for the entire network and per plant capacity. The typical bound for the overall number of refineires used in this paper is N = 1000, with Ntmax and Ntmin customized for each case study. For the nationwide supply chain case studies, the Ntmin values for the 1 kBPD and 5 kBPD refineries are set to be 5 and 3, respectively, to represent demonstration-scale GTL refineries in the early stage of the supply chain development. All other Ntmin values are set to zero. In the regional and statewide case studies, the lower bound for the 5 kBPD refineries is removed. The value of Ntmax is set to be 3 for the 200 kBPD refineries to prevent the model to excessively favor the large capacity due to economies of scale, since large risks are associated with installing 200 kBPD refineries. All other Ntmax values are set to be 500. y f ,t,l,q ≤ 1



∀l ∈ LF

(2)

( f ,t,l,q)∈FL

y f ,l,t,q ≤ N



(3)

( f ,t,l,q)∈FL

∑ ( f ,t,l,q)∈FL

∑ ( f ,t,l,q)∈FL

y f ,l,t,q ≤ Ntmax ∀t ∈ T

(4)

y f ,l,t,q ≥ Ntmin ∀t ∈ T

(5)

The correct investment cost parameters and the natural gas feedstock requirements associated with the selected GTL facilities are determined by Equations 6-7.

∑ ( f ,t,l,q)∈FL

FR f ,l =

y f ,l,t,q LC f ,t,q = CostlI ∑

( f ,t,l,q)∈FL

∀l ∈ LF

FRGf,t,q y f ,l,t,q ∀ f ∈ F G , l ∈ LF

18

ACS Paragon Plus Environment

(6) (7)

Page 19 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Equation 8 constrains the flow rates from each natural gas feedstock source locations not to exceed the amount available at that location. The availability of natural gas is determined by the level of natural gas expansion as described in Section 2. The total natural gas flow rates arriving at a particular GTL refinery has to match the feedstock requirement of the selected GTL refinery (Eq. 9). Equation 10 constrains the product ratios exiting each GTL refinery to match either the unrestricted fuel ratio, maximum diesel production, maximum kerosene production, or 2010 United States demands for gasoline, diesel, and kerosene (see Section 3). Finally, the total product flow rates arriving at each demand location must be less than the known maximum capacities of the oil refineries and less than the total demand amount within the state. All product flows in the supply chain must amount to TotalFuel, which is determined by the maximum amount that can be produced by the available natural gas feedstock via a parametric analysis (Eq. 11-12).



x f ,c,l,m ≤ FA f ,c

∀( f , c) ∈ CF

(8)

x f ,c,l,m = FR f ,l

∀ f ∈ F, l ∈ LF

(9)

y f ,l,t,q PR p,t,q ∀l ∈ LF , p ∈ P

(10)

( f ,c,l,m)∈FT

∑ ( f ,c,l,m)∈FT

z p,l,c,m =

∑ (p,l,c,m)∈PT

∑ ( f ,t,l,q)∈FL

z p,l,c,m ≤ DM p,c

∀(p, c) ∈ CP

(11)

z p,l,c,m = TotalFuel

∀(p, c) ∈ CP

(12)

∑ (p,c,l,m)∈PT

∑ (p,c,l,m)∈PT

The constraints associated with electricity requirement and production by the selected facilities are defined in Equations 13-18. Produced electricity is assumed to be sold to the grid at $0.07/kWh. Equations 13 and 14 define the electricity requirement and produced for each selected facility, respectively. The electricity requirement is fulfilled in a combination of three ways, namely through solar, wind, or grid electricity (Eq. 15). The amount of solar and wind electricity is capped at 20 MW and 40 MW, respectively, and only apply for locations that receive enough solar irradiation and wind speed for electricity generation (Eq. 16 and 17). Finally, a maximum for the total grid electricity expansion for the whole country is imposed (Eq. 18). In this paper, CapG is equal to

19

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 79

102.5 GW (i.e., about 10% increase of the grid capacity in 2009 113 ). y f ,l,t,q EL f ,t,q = EllT ∀l ∈ LF

(13)

y f ,l,t,q EPf ,t,q = EllP ∀l ∈ LF

(14)

EllT = EllS + EllW + EllG

∀l ∈ LF

(15)

EllS ≤ CapSl

∀l ∈ LS

(16)

EllW ≤ CapW l

∀l ∈ LW

(17)

∑ EllG ≤ CapG

∀l ∈ LF

(18)

∑ ( f ,t,l,q)∈FL

∑ ( f ,t,l,q)∈FL

l∈LF

Note that since all GTL refineries described in Section 3 are net producers of all electricity, all variables associated with the electricity requirement for the supply chain will take the value of zero. Equation 19 defines the freshwater input requirement for each selected facility. This value has to be matched with all water flowrates from various sources (Eq. 20) and they cannot exceed the water availabilities in each source location (Eq. 21).



y f ,l,t,q FW f ,t,q = W Fl ∀l ∈ LF

(19)

∀l ∈ LF

(20)

∀c ∈ CW

(21)

( f ,t,l,q)∈FL

W Fl = ∑ wc,l c∈CW

∑ wc,l ≤ WAc l∈LF

Similarly, Equation 22 determines the amount of CO2 sequestered from each facility. They are delivered to the sequestration sites (Eq. 23), which have maximum capacities that cannot be exceeded (Eq. 24).



y f ,l,t,q SQ f ,t,q = SFl ∀l ∈ LF

(22)

∀l ∈ LF

(23)

∀c ∈ CSQ

(24)

( f ,t,l,q)∈FL

SFl = ∑ seqc,l c∈CSQ

∑ seqc,l ≤ SQCAPc l∈LF

The objective function includes the total overall cost of the energy supply chain network, which 20

ACS Paragon Plus Environment

Page 21 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

covers (i) the investment costs associated with the new GTL plants, (ii) the electricity costs/gains, (iii) feedstock purchase and transportation costs, (iv) product transportation costs, (v) freshwater purchase and transportation costs, and (vi) sequestered CO2 transportation and injection costs. Note that the electricity sales and purchases are treated separately in the objective function to allow a degree of flexibility for the system to sell or purchase electricity at different costs and at different locations.

∑F CostlI − EllPCost El,P

l∈L

+

∑F EllGCost El,G + EllSCost El,S + EllW Cost El,W

l∈L

+∑

∑F wc,l (CostcW P +Costc,lW T )

+∑

∑F seqc,l (Costc,lCO2,T +CostcCO2,In j )

+



c∈C l∈L

c∈C l∈L

x f ,c,l,m (Cost Ff,c +Cost FT f ,c,l,m )

( f ,c,l,m)∈FT

+



PT z p,l,c,mCost p,l,c,m

(25)

(p,c,l,m)∈PT

Equations 2-25 represent a large-scale mixed-integer linear optimization (MILP) model that can be solved using CPLEX. Note that the model is flexible and can be adapted to solve for nationwide, regional, and statewide GTL supply chains by modifying the input parameters. For a nationwide supply chain case study with 50% natural gas expansion, the model consists of 16,044 binary variables, 1,516,410 continuous variables, and 19,757 constraints. The size of the model varies with the geographical scope of the supply chain problem. The solution of the MILP gives the active binary variables that represent the existence and selection of GTL refineries, the flow rates of feedstocks and products, water, CO2 , electricity amounts in the supply chain. The computation is completed on a single computer containing 8 Intel Xeon 2.83 GHz processors and shared memory parallelization, and the optimization model is solved using CPLEX with eight parallel threads. The best incumbent solutions are reported and the optimality gap for all performed computational studies are below 1%, and the computational times are within 24 CPU hours for regional and statewide case studies, and within 48 hours for the nationwide case studies.

21

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

11

Computational Studies

The GTL supply chain optimization model is solved for (i) nationwide, (ii) regional, and (iii) statewide supply chain networks. In the nationwide supply chain studies, different levels of natural gas expansion are considered to show the maximum possible amount of liquid fuels that can be produced for the United States transportation sector via the optimized GTL technologies. For the regional supply chain studies, the optimization problem is solved for the six regions of the United States, namely the Northeast (NE), Southeast (SE), Midwest (MW), Southwest (SW), Central (CE), and Western (WE) regions. The regional divisions of the lower 48 states follow the delineations defined by EIA for natural gas pipelines 121 , and the states that belong in each region can be seen in Table 3. The specifications of the regional supply chain problem include (i) only natural gas resources within the specified region may be used, (ii) the facilities may only be selected within the region, and (iii) the liquid fuels can be delivered to other regions. Finally, the statewide supply chain is defined similarly as the regional case studies, with an additional specification that no inter-county natural gas transportation is allowed. In other words, the county where a GTL refinery is located must be able to satisfy the feedstock requirement of that refinery. The statewide supply chain case study is presented for Texas, which has an abundant resource of natural gas that will greatly enhance the profitability of the GTL refineries in the state. Note, however, that such an analysis can be applied for every state with natural gas resources. Table 3: State groupings and regional definitions of the United States. The grouping defined here is based on the 10 Federal Regions of the U.S. Bureau of Labor Statistics, as published on the EIA website. 121

Region

States

Northeast Southeast Midwest Southwest Central Western

CT, ME, MA, NH, RI, VT, NJ, NY, DE, DC, MD, PA, VA, WV AL, FL, GA, KY, MS, NC, SC, TN IL, IN, MI, MN, OH, WI AR, LA, NM, OK, TX IA, KS, MO, NE, CO, MT, ND, SD, UT, WY AZ, CA, NV, ID, OR, WA

The labeling convention for the supply chain case studies is given in Table 4. The first indicator denotes the geographical scope of the supply chain problem, either nationwide (N), regional (regional acronyms), or statewide (state acronyms, e.g., TX). The second indicator denotes the

22

ACS Paragon Plus Environment

Page 22 of 79

Page 23 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

percentage of the expansion level of natural gas production and takes the value of 25, 50, 75, or 100. Indicators (U) and (R) signify the fuel ratios produced in the supply chain, whether they are unrestricted (U) or constrained to produce minimum portions of fuel ratio (R). When the fuel ratio is unrestricted, any composition of gasoline, diesel, and kerosene may be produced by the supply chain, and when the fuel ratio is restricted, at least 55 vol% of the total fuels must be gasoline, at least 20 vol% is diesel, and at least 10 vol% is kerosene. These minimum ratios are imposed to represent the composition of the United States transportation fuel demands. Lastly, the case studies with “10” are those that only allow 10 kBPD GTL refineries to exist in the supply chain. These represent the scenario in which modular 10 kBPD plants are built across the country. The 10 kBPD GTL facilities are especially attractive economically since they (i) harness the economies of scale for using one train of process units, (ii) yield good economic results in terms of investment cost and levelized cost of fuels produced, and (iii) are at a manageable scale in terms of feedstock requirements 63 .

11.1

Nationwide GTL Supply Chains

The nationwide supply chain case studies were completed for different levels of natural gas expansion (see Section 2). Based on the 2010 natural gas production levels, the availability of natural gas feedstock for the GTL supply chain is assumed to come from expansion and not from the established production levels. Four (4) levels of natural gas expansion are considered, namely 25%, 50%, 75%, and 100% of the 2010 natural gas production amounts and the usage of these amounts for the GTL supply chain is maximized. In other words, the maximum amount of liquid fuels produced is targeted, rounded to the nearest 100 kBPD. Using this natural gas figures, the resulting total amounts of liquid transportation fuels that can be produced are 1,500 kBPD, 3,100 kBPD, 4,500 kBPD, and 6,200 kBPD, respectively. The lower bounds on the number of the 1 kBPD and 5 kBPD refineries are set to be 5 and 3, respectively, to account for the initial developments of GTL technologies and demonstration of financibility using plants with lower total investment costs. The upper bound is imposed for only the 200 kBPD refineries and is set to be 3, except in the case studies with 100% natural gas expansion level, where the number is increased to 6. The summary of results for all nationwide supply chain case studies can be seen in Table 5.

23

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 79

Table 4: Case study labeling conventions.

Case study label

Geographical scope

Expansion level

Fuel ratio

Facilities

N-100-U

Nationwide

100%

Unrestricted

All

N-75-U

Nationwide

75%

Unrestricted

All

N-50-U

Nationwide

50%

Unrestricted

All

N-25-U

Nationwide

25%

Unrestricted

All

N-100-R

Nationwide

100%

Ratio

All

N-50-R

Nationwide

50%

Ratio

All

N-50-U-10

Nationwide

50%

Unrestricted

10 kBPD

N-50-R-10

Nationwide

50%

Ratio

10 kBPD

NE-50-U

Northeast

50%

Unrestricted

All

NE-50-R

Northeast

50%

Ratio

All

SE-50-U

Southeast

50%

Unrestricted

All

SE-50-R

Southeast

50%

Ratio

All

MW-50-U

Midwest

50%

Unrestricted

All

MW-50-R

Midwest

50%

Ratio

All

SW-50-U

Southwest

50%

Unrestricted

All

SW-50-R

Southwest

50%

Ratio

All

CE-50-U

Central

50%

Unrestricted

All

CE-50-R

Central

50%

Ratio

All

WE-50-U

Western

50%

Unrestricted

All

WE-50-R

Western

50%

Ratio

All

TX-50-U

Texas

50%

Unrestricted

All

TX-50-R

Texas

50%

Ratio

All

TX-50-U-10

Texas

50%

Unrestricted

10 kBPD

TX-50-R-10

Texas

50%

Ratio

10 kBPD

24

ACS Paragon Plus Environment

Page 25 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 5: Summary of results for the nationwide GTL supply chain case studies. Case Study

N-100-U

N-75-U

N-50-U

N-25-U

N-100-R

N-50-R

N-50-U-10

N-50-R-10

BEOP ($/bbl)

$70.31

$67.34

$65.05

$61.87

$71.40

$66.25

$77.96

$80.25

Levelized cost ($/GJ)

$14.77

$14.38

$14.03

$13.67

$14.82

$13.92

$16.38

$16.35

128

103

71

35

129

71

310

311

1 kBPD

5

5

5

5

5

5

0

0

U

5

4

5

5

4

5

0

0

D

0

0

0

0

0

0

0

0

K

0

1

0

0

1

0

0

0

R

0

0

0

0

0

0

0

0

5 kBPD

3

3

3

3

5

3

0

0

U

3

3

3

3

2

3

0

0

D

0

0

0

0

3

0

0

0

K

0

0

0

0

0

0

0

0

R

0

0

0

0

0

0

0

0

10 kBPD

18

18

13

8

17

13

310

311

U

17

15

10

8

13

11

252

186

D

0

0

0

0

4

0

0

83

K

1

3

3

0

0

2

57

42

R

0

0

0

0

0

0

1

0

50 kBPD

96

74

47

16

96

47

0

0

U

59

53

35

14

62

30

0

0

D

0

0

0

0

16

9

0

0

K

37

21

12

2

18

8

0

0

R

0

0

0

0

0

0

0

0

200 kBPD

6

3

3

3

6

3

0

0

U

5

3

3

3

2

1

0

0

D

1

0

0

0

4

2

0

0

K

0

0

0

0

0

0

0

0

R

0

0

0

0

0

0

0

0

19.67

13.96

9.33

4.04

18.13

8.98

8.22

8.11

Total number of facilities

Electricity produced (GW) Total fuels (kBPD) Gasoline (kBPD) Diesel (kBPD)

6,200

4,500

3,100

1,500

6,200

3,100

3,100

3,110

4,655.00

3,689.25

2,627.50

1,425.00

4,283.00

2,147.50

2,669.21

2,172.50

150.00

0

0

0

1,241.25

637.50

2.15

622.50

Kerosene (kBPD)

1,395.00

810.75

472.50

75.00

675.75

315.00

428.63

315.00

Total LHV output

34.89

25.16

17.29

8.29

35.25

17.62

17.27

17.67

18.09

13.14

9.05

4.38

18.19

9.09

8.83

8.97

Natural gas purchase

5.76

5.81

5.80

5.85

5.73

5.72

5.68

5.63

Natural gas trans.

2.83

2.79

2.75

2.81

2.86

2.56

3.07

2.98

Butane purchase

0

0

0

0

0

0

0

0

Investment cost

6.77

6.89

6.86

6.76

6.66

6.69

8.29

8.12

Product trans.

1.19

0.77

0.61

0.56

1.26

0.64

0.69

0.70

Water purchase and trans.

0.01

0.01

0.01

0.01

0.01

0.01

0.01

0.01

CO2 trans. and injection

0.45

0.51

0.45

0.36

0.44

0.44

0.96

0.99

Propane sale

1.30

1.46

1.55

1.87

1.27

1.29

1.54

1.33

Electricity sale

0.95

0.93

0.91

0.82

0.86

0.86

0.80

0.77

(million GJ/day) Natural gas usage (trillion SCF/yr) Average cost breakdown ($/GJ)

25

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

11.1.1

Costs of Fuel Production

In this section, the base case for the nationwide GTL supply chain is case study N-50-U where 50% natural gas expansion produces a total of 3,100 kBPD of liquid fuels, approximately 23.5% of the total liquid fuel demand in the United States. The break even oil price (BEOP), which is obtained from deducting the refiner’s margin from the total cost for all gasoline, diesel, and kerosene produced and dividing the remainder with the total volume produced in the supply chain, is $65.05/bbl. This BEOP is the minimum price of crude oil at which the supply chain becomes competitive with petroleum based processes and is also equivalent to a levelized cost of $14.03/GJ LHV (lower heating value) of total fuels produced. The contributing factors of the fuel production cost for each case study are comprised of: (1) natural gas purchase cost, (2) natural gas transportation cost, (3) GTL refinery investment cost, (4) water purchase cost, (5) water transportation cost, (6) butane purchase cost, (7) electricity cost, (8) fuel product transportation cost, (9) CO2 transportation for sequestration, and (10) CO2 injection cost. Natural gas, water, butane, and electricity are inputs to the GTL refineries, and gasoline, diesel, and kerosene are the main outputs of the refineries. None of the GTL refineries requires a net input of electricity, meaning that no electricity purchase cost applies to all case studies. Instead, the GTL refineries are net producers of electricity, and the sale profits from selling electricity to the grid are used to offset the overall total cost of fuel production. The second byproduct of these GTL refineries is propane/liquefied petroleum gas (LPG), and the sale profits are accounted for in the total cost figure. For case study (N-50-U), the cost breakdown for the entire supply chain is shown in Figure 2. The highest contributing factor is the investment cost to build the GTL refineries ($6.86/GJ), amounting to 41.60% of the total cost, followed by the natural gas purchase cost ($5.80/GJ) at 35.20%, and natural gas transportation cost ($2.75/GJ) at 16.67%. The relative contributions of the costs are consistent with the cost breakdown of the individual GTL refineries in Baliban et al. 63 . The product transportation at $0.61/GJ (3.70%) contributes to the overall cost in a lesser degree than the natural gas transportation due to the locations of the GTL refineries that are close to the oil refineries (see Figure 4), the demand end points for the GTL fuels. Water purchase and transportation, and CO2 transportation and injection costs contribute minimally to the overall cost

26

ACS Paragon Plus Environment

Page 26 of 79

Page 27 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

figure at 0.09% and 2.74%, respectively. A total of 9.33 GW of electricity is produced, resulting in an amount of $0.91/GJ that offsets the overall cost, and an additional $1.55/GJ is deducted from the cost by the sales of byproduct LPG. The co-production of electricity is significant because no futher fossil fuels resources are added to satisfy the power requirement of the GTL refineries, implying that no additional greenhouse gas from electricity production is emitted. Since the optimization model relies on a large amount of data, the input parameters in this model are subject to uncertainties. The major sources of uncertainties are the prices of natural gas, the availabilities, and transportation costs. These uncertainties can be incorporated in a robust optimization framework for suppply chain design under uncertainty, which is the subject of future study. To illustrate the effect of changes in the price parameters, however, a parametric analysis is completed on case study (N-50-U) in which the natural gas prices are varied to 10% below and above the current nominal values. Note that each location has a different natural gas price parameter and the drop or increase in the price is done uniformly across the supply chain. For these case studies, the BEOP values are $61.78/bbl and $68.62/bbl, respectively, resulting in 5% change in the BEOP value from the nominal case study. The selection of the locations is generally uniform, with exchanges between the selection of one additional 50 kBPD refinery and five less 10 kBPD refineries, and vice versa. When the liquid fuels ratio is unrestricted, the supply chain produces a total of 2,627.50 kBPD of gasoline (84.76 vol%) and 472.50 kBPD of kerosene (15.24 vol%), with no production of diesel. The reason why no (D) GTL refineries that maximize the production of diesel is selected is because the natural gas feedstock requirements of the (D) refineries are higher than their unrestricted (U), maximization of kerosene (K), and United States fuel ratio (R) equivalents (see Table 2). In other words, to produce the same volume of liquid fuels, the (D) refineries require a greater amount of natural gas. In the supply chain formulation where it is desired to maximize the total volume of liquid fuels from the existing natural gas, it is intuitive that the (D) refineries are not selected even though their investment costs are comparatively lower (see Table 2). For the United States, however, some amount of diesel must be produced to fulfill the transportation fuel demands, which accounts for about 22 vol% of the total demand 107 . To impose the production of diesel, the case studies with minimum portions of liquid fuels are presented. Case study (N-50-R) produces a total of 2,147.50 kBPD of gasoline (69.27 vol%), 637.50 kBPD of 27

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Product   Trans   3.70%  

Water   Purchase   0.08%  

Water  Trans   CO2  Seq   2.67%   0.01%  

CO2  Inj   0.07%  

NG  Purchase   35.20%   Investment   cost   41.60%   NG  Trans   16.67%  

Figure 2: The cost contributions for the nationwide supply chain case study (N-50-U).

28

ACS Paragon Plus Environment

Page 28 of 79

Page 29 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

diesel (20.56 vol%), and 315.00 kBPD of kerosene (10.16 vol%), with a BEOP value of $66.25/bbl, slightly higher than case study (N-50-U). Note in Table 5 that even though the BEOP is higher compared to the unrestricted case study, the levelized cost is lower ($13.92/GJ) due to the higher total amount of LHV output when diesel is produced (17.62 million GJ/day as opposed to 17.29 million GJ/day). The levelized investment cost is $6.69/GJ, the natural gas purchase cost is $5.72/GJ, and the natural gas transportation cost is $2.56/GJ, accounting for 41.63%, 35.62%, and 15.97% of the total costs, respectively (see Figure 3). The total electricity output is reduced from 9.33 GW to 8.98 GW due to the fact that the (D) refineries produce less amounts of electricity output compared to the (K) refineries (see Table 2). Product   Water   Water  Trans   CO2  Seq   2.67%   Trans   Purchase   0.01%   CO2  Inj   0.08%   3.96%   0.07%  

NG  Purchase   35.62%   Investment   cost   41.63%   NG  Trans   15.97%  

Figure 3: The cost contributions for the nationwide supply chain case study (N-50-R).

In the case where 100% of natural gas expansion is assumed (N-100-U), the total fuels produced is 6,200 kBPD with a BEOP value of $70.31/bbl and a levelized cost value of $14.77/GJ. The contributions of capital investment, natural gas purchase, and natural gas transportation costs to the 29

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

overall cost are 39.79%, 33.84%, and 16.65%, respectively. Product transportation cost accounts for 6.98% to the overall cost, CO2 transportation and injection for sequestration account for 2.65%, and water purchase and transportation account for less for 0.08%. In this case, a small amount of diesel is produced (150 kBPD) due to the selection of one 200 kBPD refinery that maximizes the production of diesel (refinery D-200). The product compositions are 75.08 vol% gasoline, 2.42 vol% diesel, and 22.50 vol% kerosene. When the fuel ratio is enforced in case study (N-100-R), the BEOP increases to $71.40/bbl and the levelized cost to $14.82/GJ. The amount of electricity produced is reduced from 19.67 GW to 18.13 GW, and the supply chain produces 4,283.00 kBPD of gasoline (69.08 vol%), 1,241.25 kBPD of diesel (20.02 vol%), and 675.75 kBPD of kerosene (10.90 vol%). The investment cost accounts for 39.28% of the total cost with similar composition for the rest of the cost contributors. The effect of natural gas availabilities to the amount of fuels that can be produced can be seen in case studies (N-100-U), (N-75-U), (N-50-U), and (N-25-U). The BEOP decreases from $70.31/bbl, $67.34/bbl, $65.05/bbl, to $61.87/bbl, respectively, and the contributions of the investment cost are 39.79%, 41.05%, 41.60%, and 41.36%, respectively. The electricity outputs are 19.67, 13.96, 9.33, and 4.04 GW, and the propane byproduct sales grows more significant as the fuel production decreases. While the cost figures are similar, the topological network of the supply chains differ between the four case studies and the differences will be discussed in Section 11.1.2. Two additional case studies are presented with only 10 kBPD GTL refineries allowed in the supply chain, representing scenarios in which modular 10 kBPD plants are built across the country. The 10 kBPD GTL refineries are especially attractive economically since they maximize the economies of scale for using one train of process units, yield good economic results in terms of investment and levelized cost of fuels produced, and are at a scale where natural gas feedstock can be supplied locally. The unrestricted case study (N-50-U-10) has a BEOP of $77.96/bbl, 20% higher than case study (N-50-U), producing 2.669.21 kBPD gasoline, 2.15 kBPD diesel, and 428.63 kBPD kerosene. This increase is largely due to the investment cost contribution (44.27%) since the 10 kBPD refineries have higher levelized investment costs compared to the 50 kBPD and 200 kBPD refineries. The natural gas purchase contribution to the cost is decreased due to the lower natural gas usage overall. However, the transportation cost for natural gas is increased from $2.75/GJ to $3.07/GJ, due to the more dispersed layout of the GTL refineries locations. Lastly, when the fuel 30

ACS Paragon Plus Environment

Page 30 of 79

Page 31 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

ratio is enforced (N-50-R-10), the BEOP further increases to $80.25/bbl due to lower byproduct LPG and electricity production. 11.1.2

Supply Chain Network Topology

The base case study (N-50-U) selected a total of 71 GTL refineries, 5 of 1 kBPD capacity, 3 of 5 kBPD capacity, 13 of 10 kBPD capacity, 47 of 50 kBPD capacity, and 3 of 200 kBPD capacity. Note that the solution consistently hits the lower and upper bounds for the 1, 5, and 200 kBPD refineries, respectively, for all of the nationwide case studies. The selection of larger facilities are more favorable, with most of fuel productions taking place in 10, 50, and 200 kBPD refineries. Each GTL refinery of a given capacity can have four modes of fuel production, namely (i) unrestricted (U), where no fuel ratio restriction is imposed, (ii) maximization of diesel (D), with a ratio of 75 vol% diesel and 25 vol% gasoline, (iii) maximization of kerosene (K), with a ratio of 75 vol% kerosene and 25 vol% gasoline, and (iv) fuel ratios commensurate with the United States demand (R), with a ratio of 67 vol% gasoline, 22 vol% diesel, and 11 vol% kerosene. The breakdown of the selected refinery types for each case study can be seen in Table 5. In case study (N-50-U), only the (U) and (K) GTL refineries are selected. The 1, 5, and 200 kBPD refineries are (U) GTL refineries, and the 10 and 50 kBPD refineries are a mixture of (U) and (K) refineries. The selection of these plants are not only motivated by the costs of the plants, since on an individual plant basis, the (D) refineries have the lowest investment costs. Instead, it is also driven by the amount of natural gas feedstock required by the refineries, which in this case study amounts to 9.05 trillion SCF annually. Even though the (D) refineries have lower overall costs, they require a higher amount of natural gas per unit volume of total fuels produced (see Table 2). For example, the (D-10) refinery requires 3.76 MSCF/hr, higher than the 3.61 MSCF/hr required for (U-10), 3.54 MSCF/hr for (K-10), and 3.62 MSCF/hr for (R-10). In the supply chain case studies where we maximized the production of liquid fuels to use as much natural gas resources as possible, the selections of the (U) and (K) refineries dominate. Most of the selected 71 refineries are located in the Southwest and Central regions of the United States, with several facilities around the intersections of the Northeast, Midwest, and Southeast regions. Figure 4 shows the graphical representation of the supply chain network layout, with the blue boxes indicating where the natural gas resources are drawn from, the brown circles indicating 31

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

"""" " " " " "" "^ " _ "MT ND "" " " """ _ ^ " " " MN MI "" """^ _ "" " " ^ _ ME "" " " """ " " "" " " "" " " SD WI ID """ " " """ " VTNH " " " " "^ _ OR MI " " " " " " " "" NY "WY" _" ^ _" ^ """""" "^ _ """ """ ^ IA "" " _ "" ^ " " _ " " " " CTMA " " NE " " " " " " "" """ _" " " " "" " "" " "" _ ^ "" ""PA " _ " ^"^ " " " RI _ ^ " " " " _ " " ^ " " " " " " " " " " " " " " "" "" " "" " " NJ " "" " " "" """"" " _ "" IL IN OH " " "" " " " "" " "" "" " "^ " " " "" MD "" " " " UT " NV ^ " _ " " " "" " " " " " " " " " " " " " " CO " " DE "^ """" WV "" MO " " " " "" " "" " " _ "" ^ """" " "" " " "" "KS " " "" " " "" "" " " """" "_ " " " "" "" " " " " " KY " " """ "" " "" " " " " " " " "" "" "" " """ "" "" " "" VA " " "" " " " "" """ " " " " " " " " "" _ " ^ _ ^ " " " " " " " " "" " " " " " " " _ _ ^ ^ " " " " " " " "^ CA """ "" """ "" _ "" """ """" _^ ^ " "" " "" TN "" _"" ^ """ """" _ " " " OK NC _" " " ^ " "" " " " """" " " " " "" " " " "" _ ^ " " " " " " " """"" " " " AR "" " " _ ^ " " " " " " NM " " " " " AZ " " _ ^ " " " " " " " """" " " """"" " SC " " _ ^ " " " " " " " " " " " """AL _ " ^ " " """ "" " "" "" " " " "" "" " " _ ^ " """ GA " " " " " " " " MS " _ ^ " " " "" " _ ^ " " "" " "" "" "^ " """ " _ _ "" " " " "" " _ """ ^ " " " " "" "" """"^ " """ """ "" "" " " " " " " " " " _ "" ^ " " " " " " " " " " " " " " " " TX " " " " " " " " " " """ " " " " " " " " " " " " " " " LA " " " " " " "" "" " """ " " """ "" "" _ " "" " " "" "" _ ^ " " ""^ _" "" """"" ^ " " "" _"" ^ " " "" _ " ^ "" """ Total capacity 3,100,000 BPD Natural Gas Feed Source " "" _ ^ " " " "" " " _"" " """"" " ^ " " " " " " " " " " " Plant size (BPD) Flow " """"""" FL " "" 1000 _ ^ "" " 0.000000 - 22.131700 "" "" " " " 5000 22.131701 73.574100 " """ " 10000 " 73.574101 - 180.423600 180.423601 367.672100 " 50000 0 110 220 440 660 880 1,100 1,320 1,540 Miles " 367.672101 - 866.572800

WA

200000

_ ^

Fuel Products Delivery Locations

Figure 4: Graphical representation of the locations of selected facilities for unrestricted nationwide supply chain case study with 50% natural gas expansion (N-50-U). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

32

ACS Paragon Plus Environment

Page 32 of 79

Page 33 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

the locations of the selected refineries, and the red stars indicating the end locations of fuel delivery. The produced liquid fuels are delivered to oil refineries to be distributed to the final consumer points or to be blended with petroleum-based fuels before distribution. The GTL fuels are high-quality liquid fuels that have close to zero sulfur content, and can be used to improve the quality of the petroleum fuels via blending. The locations of the facilities align well with the locations of the fuel delivery locations, which correspond to the fact that the product transportation cost contributes in a lesser degree than the natural gas transportation cost. This is intuitive since there are more locations of natural gas resources than fuel delivery locations and multiple sources of natural gas are required to fulfill the requirement of the facility. Locating the GTL refineries farther from the oil refineries will not guarantee that the natural gas transportation cost will be reduced. When the nationwide supply chain produces a ratio of gasoline, diesel, and kerosene (case study N-50-R), the network topology consists of the same total amount of GTL refineries and the same breakdown for each capacity level. However, the types of the selected GTL refineries differ due to the minimum amount of diesel imposed. Diesel is produced by the 9 refineries of 50 kBPD capacity and 2 refineries of 200 kBPD capacity. The resulting total amount of natural gas usage increases from 9.05 to 9.09 trillion SCF/yr, and the modified topological layout of case study (N50-R) is shown in Figure 5. Note that no (R) refineries are selected in the supply chain due to their high investment costs compared to the (U), (D), and (K) refineries. When 100% of natural gas expansion is allowed (N-100-U), (U) and (K) refineries are selected except for 1 refinery of 200 kBPD capacity. A total of 18.09 trillion SCF/yr natural gas is required and more populated Southwest and Central regions are observed (see Figure 6). In case study (N-100-R), the selected (D) refineries selected are the 3 refineries of 5 kBPD capacity, 4 refineries of 10 kBPD capacity, 16 refineries of 50 kBPD capacity, and 4 refineries of 200 kBPD capacity (see Figure 7 and Table 5). The selection of more (D) facilities in the larger refineries is caused by the switch in the ranking of natural gas requirement between the modes of fuel production. At 200 kBPD, the (D) refinery requires 73.85 MSCF/hr of natural gas (D-200), while the (K-200) refinery requires 77.28 MSCF/yr, the (U-200) refinery requires 73.57 MSCF/yr, and the (R-200) refinery requires 72.20 MSCF/yr of natural gas. Thus, no (K-200) refinery is selected in any of the nationwide case studies and only (U-200) and (D-200) refineries are selected. The (R-200) refinery requires the lowest amount of natural gas feestock, but the portion of diesel produced is 33

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 79

""""""" " " "" "^ " _ MT """ " " """" " ND _ ^ MN MI "" " ME "" _ ^ _ ^ " " " "" " " " " " " " " " " " " SD WI ID " " " " " " VT " " " " " " " " " OR " _ MI ^ """ """" NY NH "" WY" _" ^ _" ^ "" " "" "" _ """ ^ "" """ ^ MA IA " " " " " " _ " " ^ _ " NE " " " " " " "" " _"" " """"^ " "" ^ " """ CT " " " " " _ " " PA " " _ " " ^ " " " _ ^ " " " _ ^ " " " " "" " " """"" "" NJ "" " " " "" " " " " " IL IN "OH _ ^ " " " "" """ " "" " " " " " " ""MDDE "" " " NV ^ " " """ " " _ UT " " CO " " " " " MO " " " " "" " " " "" " " WV " " " " " " KS " " " " " " " " " " " " _ " " ^ _ " " ^ " " " " """ " " "" _ """" "" ^ "" " " " "" """ "" " "" " " " """ " KY "" "" " " " " " " "" " " "" VA "" " "" "" "" " " " " " "" " " " " " "" " " " _ ^ " " " " " """ "" "" " "" "" " " "" " " " _ _ ^ ^ " " " CA " " " " " " " " " " " " " _ ^ " " """ """""""""" _ ^ " "" TN " "" "" " _" " ^ OK NC ""^ _ " " " _" " ^ "" "" "^ """" " "" "" " " " "" "" "" " """" AR " "_"""" " "" " """^ " _ ^ " " NM " " " " " " AZ " _ " ^ " " " SC _ " " " " " " " " " "" """ """ " " "" " """" "" GA """ " "" _ ^ " "" "" "" "" " _ " ^ AL " "" " " " _ " " ^ " _ MS ^ "" " " " " " "" "" " "" " _ "" "" ^ " " " " "" " "" _ _ " ^ ^ " " _ ^ " " " " " " " " " " " " " " " " "" " " " " " " " " " _ """ " " " ^ " " " " " " " " " " " TX " " " " " " """ " " "" " "" " "" "" " " """ " " " "LA "" " "" """ " "" "" " " "" " " _""" """ ^ " " "" " "" "" " " " " " """ " _ " ^ " " "" " " " _ ^ "" " " "" " " " " " " " " _ ^ _ ^ " " " _ ^ " " """ " " " " " " _ ^ "" """ """ "" " "" "" " "" " " " " FL " " "" " _ ^ " " " """" " Legend " WA

United States - Total Capacity: 3,100,000 BPD Natural Gas Feed Source Plant size (BPD) 1000 5000

10000 50000

200000

_ ^

Flow (mscf/day)

" " " " "

0 105210

420

630

840

1,050 1,260 1,470 Miles

0.000000 - 22.826400

22.826401 - 76.020200

76.020201 - 185.241200

185.241201 - 360.608300 360.608301 - 964.403200

Fuel Products Delivery Locations

Figure 5: Graphical representation of the locations of selected facilities for the nationwide supply chain case study with 50% natural gas expansion and minimum portions of liquid fuels produced (N-50-R). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

34

ACS Paragon Plus Environment

Page 35 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

not sufficient to fulfill the portion of diesel required for the supply chain. Consequently, along with the fact that the (R) refineries incur higher overall costs compared to the (U), (D), and (K) refineries, they are not selected in the supply chain.

"""" " " " " "" "^ " _ "MT ND "" " " """ _ ^ " " " MN " MI "" """^ _" " " ^ _ ME "" " " """ _ ^ " " "" " " " SD WI ID """ " " """ " " " VTNH " " " " " " _ OR ^ MI " " " " " " " "" NY "WY" _" ^ _" ^ """""" "^ _ """ """ ^ IA "" " _" "" ^ "" " CTMA _ " " "^ " " " NE " " " " " " " _ " ^ " " "" " " "" "" ""PA " " "" "_ "" "" "" _" " ^ "" " " RI _ ^ " " " " " _ " " ^ "" " " "" " " """ " "" " " " " "" " "" " " " " "" " _ ^ _NJ ^ " "" " IL IN "OH " " " " " " " " " " " " " " " " " " " " UT " " NV " _ " ^ " " " "" " " " "" " " "" " MDDE " " " """ " " "" " " "" " CO "" WV "" " " " " " "" " "" " " _ " ^ """^ _ "" MO " " " " " " "" " "KS " " "" " " "" "" " _ " ^ """" "^ """ " "" "" " " " _""" " "" " KY " " "" " VA "" "" "" " " " " " " "" " "" " " """ "" "" " " """" " " " " " " " " " " " " " " " " " " " "" _ " ^ _ ^ " """ "" """" "" "" "" _ " " """ "" _" ^ """ "" " "^ "^ " CA " "" " _ "" " """ """" _^ ^ " " " "" "" TN "" _"" ^ """ """" _ " " OK NC _" " " " ^ " "" " _" " " ^ " " " " " "" " " " " " " _ ^ " " " " " " " " " " " " " AR"" "" "" """"" " "" " """ _ ^ NM "" " " AZ " _ ^ " " " " _ " ^ " " " " " " SC " "" _" ^ " " "" " "" " "" " " """ " " """" " _ " ^ "" " " "" "" "" " " _ ^ "" "" AL " "" _ ^ " _ """ GA " " "" " " "" " MS "^ " " " " "" " " _ ^ " " "" " """ "" " "^ " """ " _ _ "" " ^ " " "" " _ """ ^ " " " " " " " " " " " " " " " " " " " " " " _ "" " ^ " _ ^ " " " " " " " " " " " " " " " " " " " " " TX " " " " " " " " """""" " " " " "" " " "" " " " "LA "^ _" """ """ "" "" " _"^ "" " """ " "" "" "" " _ " "" " " _ ^ " " ^ _ ^ " " " _ """"" ^ " " " " "" _ ^ " " " " "" _ " ^ _ ^ """ "" """ " Total capacity 6,200,000 BPD Natural Gas Feed Source "^ _ ^ _ "" "" " " " "" " " " _ _ ^ " " ^ " " " " """ _ ""^ "" Plant size (BPD) Flow "" FL _ "" "" "^ "" " " 1000 " _ ^ "" " 0.000000 - 48.366100 "" " " "" " " 5000 " 48.366101 - 150.694800 "" 10000 " 150.694801 - 318.361200 " 318.361201 - 721.216700 50000 0 110 220 440 660 880 1,100 1,320 1,540 Miles " 721.216701 - 1526.546200

WA

200000

_ ^

Fuel Products Delivery Locations

Figure 6: Graphical representation of the locations of selected facilities for unrestricted nationwide supply chain case study with 100% natural gas expansion (N-100-U). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

The described pattern also applies in case studies (N-50-U-10) and (N-50-R-10). In the former case, mostly (U) and (K) 10 kBPD refineries are selected, with 1 (R-10) refinery selected. In the latter case, (U), (D), and (K) refineries are selected to fulfill the minimum fuel ratio. Since the (D-10) refineries require more natural gas, the total natural gas usage increases from 8.83 to 8.97 MSCF/yr. A larger number of total plants are selected in these case studies, 310 and 311, respectively, and Figures 8 and 9 show that the Southwest and Central regions are the most

35

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 79

_ ^ _ WA ^

""""""" " " "" "^ " _ MT """ " " """" " ND _ ^ MN MI "" " ME "" _ ^ _ ^ " " " "" " " " " _ ^ " " " " " " " " SD WI ID " " " " " " VT " " " " " " " " " OR " _ MI ^ """ """" NY NH "" WY" _" ^ _" ^ "" " "" "" _ """ ^ "" """ ^ MA IA " " _ ^ " " " " _ " " ^ _ " NE " " " " " " "" " _"" " """"^ " "" " """ CT " " " " " _ " " PA " " _"^ " ^ " " " " " _ ^ " " " _ ^ " "" " " "" _NJ " " """"" "" "" " " "" " " " " " IL IN "OH _ ^ " " " " "" """ " "" " " " " " " ""MD^ "" " " NV ^ " " """ " " " "" _ UT " " " CO " " " " " MO " " " " "" " " " "" DE " " WV " " " " " KS " " " " " " " " " " " " " _ " " ^ _ " " ^ " " " " " " "" "" " " " "" """" "" _ ^ """ " " _ " " ^ " """ "" " "" " _ ^ " " " " " " " " KY " VA " " " " " " " " " " " " " " "" " " " " " " " " " " " " " " " " " " " " " " " _ ^ " """ "" """" "" " """ """ "" "" """" _"" " _" " " ^ CA ^ " "" "" " "" _ ^ " " " _ ^ " "" TN " "" "" "" " _" " ^ " OK NC ""^ _ " " ^ "" " " " " " " " " _ ^ " "" " " " _ " " ^ " " " " " " " " " " " _ " " " " " " " " " " " AR"" " """" "" "" """ " _ ^ " NM " """"""" " " " " AZ " _ ^ " SC " " " _ ^ _" ^ " " " " " "" """ """ " " "" " """" "" GA """ " "" _ ^ " "" "" "" "" " _ " ^ AL " "" " " " _ " " ^ " _ MS ^ "" " " " " " "" "" " "" " _ "" "" ^ " " " " "" " " "" _ _ " ^ ^ " " _ ^ " " " " " " " " _ ^ " " " " " " " " " " " " " " " " " " _ """ " " " ^ " " " " " " " " " " " TX " " " " " " " "" " "" " " "" "" " " """ " " " "" "" "LA "" " " "" """ " "" """ " " "" " " " _""" """ ^ " " " " " " " " " _ " ^ " " _ " ^ " " " " " " _ ^ " " " " " " " " " " " " " " _ """ "^ " _ "" "^ "^ _ _ "^ "" " " " _ _ ^ "" "^ " " "" _ ^ " "" " " " "" "" " "" _ " ^ " " " FL " " " " _ ^ " """ "" " Legend "" " United States - Total Capacity: 6,200,000 BPD Natural Gas Feed Source Plant size (BPD) 1000 5000

10000 50000 200000

_ ^

Flow (mscf/day)

" " " " "

0 105 210

420

630

840

1,050 1,260 1,470 Miles

0.000000 - 47.099100

47.099101 - 147.148300

147.148301 - 306.977300 306.977301 - 721.216700

721.216701 - 1526.546200

Fuel Products Delivery Locations

Figure 7: Graphical representation of the locations of selected facilities for the nationwide supply chain case study with 100% natural gas expansion and minimum portions of liquid fuels produced (N-100-R). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

36

ACS Paragon Plus Environment

Page 37 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

populated areas due to their abundant natural gas resources. A series of oil refineries along the southern border of Texas and Louisiana receive the fuels produced from the GTL facilities.

"" """" " "" ""^ ND _"MT " " " """" " ^ _ " MN "^ " MI " ME " "" _ _"" "" SD ^ "" " " " " "" " " " " " WI " ID " " " VT " "" " " OR " " " " " _ ^ ""MI """" NY NH WY" "" " _" ^ _ ^ " "" " "" " " """ ^ _ ^ IA " " " " " _ ^ " " CTMA " " " _ ^ ""^ " " _"" ^ " " " " " " " " "" NE " """ ^ " _""" " " _ PA " " "_ " " " " " _"^ ^ "" " " " " _ ^ " " " " " _ ^ " " "" OH " _ " " " " " NJ " " " " " " "" "" " "" ^ "" " " IL IN "" " " """" " " _ " " " "" " " " " " " NV ^ " " " " "" _ UT" "" " " " " " MO " " """ "" " " " ""CO " " " WV " " "" MDDE KS " " " " " " " " " " " " " " " " _ ^ " " _ ^ " " " " " " " " " " " " " " " " " " _ ^ " " " " " " " " " " """" " VA " " " " " "" " "" " " "" """ " """ " "" " """ " " " " ""KY "" " " " " "" " " " _ ^ _ ^ " " " " " "" "" " " "" " " "" " """ CA " " _ ^ " _ ^ " " """ " " " " " " " "" " " " " " " _ ^ " " TN " _ ^ " " " " " " " " OK NC " " " _ ^ " " " " _" " "" " "" " _ "NM ""^ ^ "" " " " _" " " " " " " "" " " " " " _"""" ^ "" " """ " "" " AR " " " "" " "^ " " _ ^ " " " AZ " " " " " "" " " _ ^ SC " " " " _ ^ " " " " " " " "" """ "" " " "" GA " " """ " """ " " _ ^ " "" "" "" " " " AL " _ ^ " " MS " " " " " " " " "" " " _ ^ "" " " " " " " _ ^ " "" " " "" " " " " _ ^ " " " " " " _ ^ _ ^ " " _ ^ " " " " " " " " " " " " " " " " " "" " " " " " " " " " " " " " " _ """ ^ _ ^ " "" "" " TX "" """ ""LA " " "" " "" " " "" """ " " "" "" " " "" " "_ " " " " " "" " " " " " " " " """ " " " " " "" "" " _ ^ " ""^ _ ^ " " " " " """" " "" _ ^ " " " " " " " _ ^ " "" " " " " " " " " " _ ^ " _ ^ _ ^ " " " _ ^ " _ ^ " " " " " "" "^ _ " " " _"" ^ "" " " "" "" " "" " " _ ^ " FL " " " " _ ^ " " """ "" Legend "" " WA

Plant size (BPD) 10000

_ ^

Fuel Products Delivery Locations

Natural Gas Feed Source Flow (mscf/day)

" " " " "

0.000000 - 8.220400

0 110 220

440

660

880

1,100

1,320

1,540 Miles

8.220401 - 23.111000

23.111001 - 41.819000

41.819001 - 68.173400 68.173401 - 86.614100

Figure 8: Graphical representation of the locations of selected facilities for unrestricted nationwide supply chain case study with 50% natural gas expansion and only 10 kBPD facilities allowed (N-50-U-10). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

11.2

Regional GTL Supply Chains

While the nationwide supply chain studies can address the feasibility question of whether the natural gas supply in the United States can replace petroleum derived fuels and by how much, the implementation of the GTL supply chains will more applicable in a regional or statewide scale since each region will primarily optimize in favor of its own interest and energy needs. Thus, the optimization problem is solved for the six regions of the United States and the natural gas 37

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 79

"" """" " "" ""^ ND _"MT " " " """" " ^ _ " MN "^ " MI " ME " "" _"" "" _ ^ "" " "" " " "" " " " SD " " WI " ID " " " VT " "" " " OR " " " " MI _ ^ "" """ "" NY NH WY" "" " "" _^ ^ _ " "" "" """""""" CTMA "" " " """ ^ _ ^ IA " " _ ^ ""^ " " " " " " " " " " " "" NE " """ " _""" " " _ PA " " "_ " "" " " " " _"^ ^ "" " " _ ^ " " _ ^ """ " "" OH " " " " " " NJ " " " " " " IN " " " " " " " IL " " " " " " _ ^ _ " " " " " " " ""MD^ NV ^ " " " """ " " " " " " _ UT" """""CO "" " " " " " " " DE " " " " " WV " " " " KS " " " " " " " " MO " " """"""" " " " " " "" _ ^ "" VA _ ^ """ "" " " " """ """ "" " " " _"" ^ " " " " " " "" " " _ ^ " " KY " "" " " " " "" " "" "" " "" " " " " "" " " " " "" " "" " " " " " " " " " _ ^ " " " " " " " " " " " " " " " " " " _ ^ " _ ^ " " CA " " " " " " " " " " " " " "" " " _ ^ " "" TN " _ "" " """ "" "" OK NC " _" ^ ""^ " " _ " " " "" " _" "NM " ^ "^ " " " " _ "^ " " " "" " " " " " _"""" "" " "" " "" "" " AR "" " " " " " " "^ " " " _ ^ " " AZ " " "^ " " " " " " _ ^ SC " " " " _" " " " " " " " " GA "" """" " " " " " " "" " """ " " _ ^ " "" "" "" " " " AL " _ ^ " " " MS " """ " " " " "" " " _ ^ "" " " " " " """" " " " _ ^ " " " " " " " " " _ ^ " " " " " " _ ^ _ ^ " " _ ^ " " " " " " " " " " " " " " " " " " " " " " " " " " " """ _ """"""" ^ " " "" " TX "" "" "" " " "" " " """" " " " "" " " " " LA " "" "^ " " " " " "" " " " " " " " " " """ " _ " " " " " "" " " _ ^ " ""^ _ ^ " " " " " """" " " _ _ ^ " " " " " " " _ ^ " "" " " " " " " " " " _ ^ " _ ^ _ ^ " " " _ ^ " _ ^ " " " " " "" "^ _ " " " _"" ^ "" " " "" "" " "" " " _ ^ " FL " " " " _ ^ " " """ "" Legend "" " WA

Plant size (BPD) 10000

_ ^

Fuel Products Delivery Locations

Natural Gas Feed Source Flow (mscf/day)

" " " " "

0.000000 - 7.707800

0 110 220

440

660

880

1,100

1,320

1,540 Miles

7.707801 - 21.682400

21.682401 - 41.350500 41.350501 - 66.817700 66.817701 - 90.338800

Figure 9: Graphical representation of the locations of selected facilities for the nationwide supply chain case study with 50% natural gas expansion, minimum portions of liquid fuels produced, and only 10 kBPD facilities allowed (N-50-R-10). The facilities are represented by dark brown circles with sizes that correspond to their respective capacities. The blue boxes indicate where the natural gas resources are drawn from and the red stars indicate the end locations of fuel delivery.

38

ACS Paragon Plus Environment

Page 39 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

expansion level is set to 50% for all case studies. The specifications of the problem are modified such that (i) only natural gas resources within the specified region may be used, (ii) the refineries may only be located within the region, and (iii) the liquid fuels can be delivered to other regions. The lower bound of the number of 5 kBPD refineries is removed to accommodate regions with less natural gas resources and the summary of results for the unrestricted and minimum fuel ratio regional case studies are shown in Tables 6 and 7. 11.2.1

Northeast Region

The Northeast region includes Connecticut, Maine, Massachusetts, New Hampshire, Rhode Island, Vermont, New Jersey, New York, Delaware, DC, Maryland, Pennsylvania, Virginia, and West Virginia. The unrestricted fuel scenario (NE-50-U) yields the selection of 14 refineries, 5 of 1 kBPD capacity, 1 of 5 kBPD capacity, 6 of 10 kBPD capacity, and 2 of 50 kBPD capacity, producing a total of 170 kBPD liquid fuels. The selected facilities are located in the western side of Pennsylvania, western side of Virginia, and the northern and southern sides of West Virginia. Figure 10 shows the topological layout of the Northeast GTL supply chain. The natural gas resources, depicted as blue boxes, are drawn from the northern and western parts of Pennsylvania and across the state of West Virginia, and the GTL refineries deliver fuels to oil refineries in New Jersey, Ohio, and West Virginia. The 170 kBPD of fuels are composed of 95 kBPD gasoline and 75 kBPD kerosene, with an overall BEOP value of $71.86/bbl. When the minimum portions of fuel ratio are imposed (NE-50-R), the supply chain network selects the same total number of refineries. However, the 10 kBPD refineries changed from being all (U-10) refineries to 2 (U-10) refineries, 1 (D-10) refineries, and 3 (K-10) refineries. The two 50 kBPD refineries that maximize kerosene production in case study (NE-50-U) are not selected and instead, 1 (U-50) and 1 (D-50) refineries are selected. The decrease in the total electricity production from 0.61 GW to 0.47 GW and an increase in the natural gas usage from 0.49 trillion SCF/yr to 0.50 trillion SCF/yr are observed due to the selection of (D) refineries, and the resulting fuel composition consists of 102.50 kBPD gasoline, 45 kBPD diesel, and 22.50 kBPD kerosene, increasing the BEOP value to $73.45/bbl. In this topology, the delivery of fuels goes to more destinations, including locations in Ohio, Pennsylvania, New Jersey, and West Virginia (see Figure 11). The differences in the cost breakdown are most marked in the natural gas transportation, where the cost 39

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 79

Table 6: Summary of results for the regional GTL supply chain case studies with unrestricted fuel ratios. Case Study

NE-50-U

SE-50-U

MW-50-U

SW-50-U

CE-50-U

WE-50-U

BEOP ($/bbl)

$71.86

$81.33

$81.02

$65.60

$72.46

$106.27

Levelized cost ($/GJ)

$14.64

$17.31

$17.25

$14.18

$14.84

$21.84

Total number of facilities

14

11

8

80

30

9

1 kBPD

5

5

5

5

5

5

U

5

5

5

5

4

5

D

0

0

0

0

1

0

K

0

0

0

0

0

0

R

0

0

0

0

0

0

5 kBPD

1

0

0

1

1

1

U

1

0

0

1

1

1

D

0

0

0

0

0

0

K

0

0

0

0

0

0

R

0

0

0

0

0

0

10 kBPD

6

6

3

54

13

3

U

6

6

3

40

13

3

D

0

0

0

0

0

0

K

0

0

0

14

0

0

R

0

0

0

0

0

0

50 kBPD

2

0

0

17

10

0

U

0

0

0

13

3

0

D

0

0

0

0

0

0

K

2

0

0

4

7

0

R

0

0

0

0

0

0

200 kBPD

0

0

0

3

1

0

U

0

0

0

3

0

0

D

0

0

0

0

0

0

K

0

0

0

0

0

0

R

0

0

0

0

1

0

0.61

0.16

0.09

5.57

2.86

0.10

Electricity produced (GW) Total fuels (kBPD) Gasoline (kBPD) Diesel (kBPD)

170

65

35

2,000

840

40

95.00

65.00

35.00

1,745.00

511.14

40.00 0

0

0

0

0

43.85

Kerosene (kBPD)

75.00

0

0

255.00

285.01

0

Total LHV output

0.97

0.36

0.19

11.13

4.79

0.22

0.49

0.19

0.10

5.80

2.43

0.11

(million GJ/day) Natural gas usage (trillion SCF/yr) Average cost breakdown ($/GJ) Natural gas purchase

5.61

6.39

6.42

5.97

5.09

6.24

Natural gas trans.

1.97

3.04

3.02

2.62

2.84

7.41

Butane purchase

0

0

0

0

0.11

0

Investment cost

7.56

8.83

9.23

7.08

6.95

9.33

Product trans.

0.82

0.41

0.20

0.52

0.99

0.11

Water purchase and trans.

0.01

0.01

0.01

0.01

0.01

0.01

CO2 trans. and injection

0.63

1.31

0.88

0.47

0.57

1.47

Propane sale

0.91

1.92

1.74

1.65

0.71

1.97

Electricity sale

1.06

0.77

0.77

0.84

1.00

0.76

40

ACS Paragon Plus Environment

Page 41 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 7: Summary of results for the regional GTL supply chain case studies with minimum portions of fuel ratios. Case Study

NE-50-R

SE-50-R

MW-50-R

SW-50-R

CE-50-R

WE-50-R

BEOP ($/bbl)

$73.45

$90.97

$87.83

$68.36

$73.34

$123.09

Levelized cost ($/GJ)

$14.99

$17.92

$17.38

$14.23

$14.99

$23.48

Total number of facilities

14

12

8

88

28

10

1 kBPD

5

5

5

5

5

5

U

5

1

4

5

1

1

D

0

1

0

0

4

1

K

0

3

1

0

0

3

R

0

0

0

0

0

0

5 kBPD

1

2

0

1

3

3

U

1

1

0

1

3

1

D

0

1

0

0

0

1

K

0

0

0

0

0

1

R

0

0

0

0

0

0

10 kBPD

6

5

3

64

12

2

U

2

2

1

49

12

1

D

1

2

1

0

0

1

K

3

1

1

15

0

0

R

0

0

0

0

0

0

50 kBPD

2

0

0

15

6

0

U

1

0

0

12

2

0

D

1

0

0

0

0

0

K

0

0

0

3

2

0

R

0

0

0

0

2

0

200 kBPD

0

0

0

3

2

0

U

0

0

0

0

0

0

D

0

0

0

3

1

0

K

0

0

0

0

0

0

R

0

0

0

0

1

0

0.47

0.17

0.10

5.51

2.39

0.11

Electricity produced (GW) Total fuels (kBPD)

170

65

35

2,000

840

40

Gasoline (kBPD)

102.50

35.75

19.25

1,325.00

513.54

22.00

Diesel (kBPD)

45.00

19.50

7.50

450.00

217.64

12.00

Kerosene (kBPD)

22.50

9.75

8.25

225.00

108.82

6.00

Total LHV output

0.97

0.37

0.20

11.41

4.82

0.23

0.50

0.19

0.10

5.80

2.43

0.12

(million GJ/day) Natural gas usage (trillion SCF/yr) Average cost breakdown ($/GJ) Natural gas purchase

5.66

6.19

6.18

5.83

5.06

6.08

Natural gas trans.

2.56

2.93

2.92

2.47

3.24

7.22

Butane purchase

0

0

0

0

0.16

0

Investment cost

7.48

8.70

8.89

6.93

6.70

9.36

Product trans.

0.68

0.36

0.19

0.58

0.85

0.12

Water purchase and trans.

0.01

0.01

0.01

0.01

0.01

0.01

CO2 trans. and injection

0.72

1.56

0.85

0.45

0.55

2.56

Propane sale

1.31

1.06

0.86

1.22

0.76

1.08

Electricity sale

0.81

0.77

0.80

0.81

0.83

0.80

41

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 79

MI

ME

WIWI WI

MI

IL

NY

IN

KY

" " " " " " " " " " " "" "" " " " " " " " PA _ ^ ^ " _ "" " " " " " OH " "" " " " _ ^ "" "" " " " MD " "" " "" " " " " " " "" DE DC " " ""WV " " _"" " " "" " " ^ "" " "" " " " VA "" " " " " "

TN AL

VT

GA

MA

CT

RI

NJ Northeast Region (170,000 BPD) Natural Gas Feed Source Plant size (BPD)

Flow

1000 5000

10000 50000

_ ^ 0 35 70

NC

NH

Fuel Products Delivery Locations

140

210

280

" " " " "

350

0.000200 - 5.380100

5.380101 - 14.197200

14.197201 - 29.013600 29.013601 - 52.496600

52.496601 - 156.820500

420

490 Miles

SC

Figure 10: Graphical representation of the supply chain case study for the Northeast region with unrestricted fuel ratio (NE-50-U). The Northeast region includes CT, ME, MA, NH, RI, VT, NJ, NY, DE, DC, MD, PA, VA, and WV. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

42

ACS Paragon Plus Environment

Page 43 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

in (NE-50-R) is higher than (NE-50-U) due to the increased natural gas requirement. The propane sale increases from the unrestricted to the minimum fuel ratio case study, since (K) refineries do not produce byproduct LPG, but the electricity sale decreases due to less electricity production in (D) refineries compared to (K) refineries. The levelized investment cost decreases from (NE-50-U) to (NE-50-R) due to the lower costs of the (D) refineries, and this trend is consistent for all regional case studies. Legend

Northeast Region (170,000 BPD) Plant size (BPD) 1000 5000

10000 50000

_ ^

Fuel Products Delivery Locations

Natural Gas Feed Source Flow (mscf/day)

" " " " "

0.000000 - 4.062100

4.062101 - 14.197200

14.197201 - 32.217900 32.217901 - 52.936900

52.936901 - 130.865500

ME VT

MI

" _ " " " " " _^ ^ " " " " " "" "" " " " """ " PA " _ ^ ^ _ "" " " " " " OH " " " "" " _ ^ " "" " " " " MD " "" " " " " " " "" " " DE DC " "" "" " " WV _"""" " "" " " ^ " " "" " " " KY VA """ 0 """ " " "

TN

NH

NY MA CT

RI

NJ

35 70

140

210

280

350

420

490 Miles

NC

Figure 11: Graphical representation of the supply chain case study for the Northeast region with minimum portions of liquid fuels produced (NE50-R). The Northeast region includes CT, ME, MA, NH, RI, VT, NJ, NY, DE, DC, MD, PA, VA, and WV. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

11.2.2

Southeast Region

The Southeast region includes Alabama, Florida, Georgia, Kentucky, Mississippi, North Carolina, South Carolina, and Tennessee. When the fuel ratio is unrestricted (case study SE-50-U), the solution yields the selection of 11 refineries, 5 of 1 kBPD capacity and 6 of 10 kBPD capacity, producing a total of 65 kBPD liquid fuels. The selected facilities are located in the eastern side of Kentucky, easterm side of Mississippi, and western part of Alabama, and deliver fuels to the 43

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 79

same three states. From Figure 12, it can be seen that the locations of the GTL refineries align with the oil refineries and two clusters of GTL supply chain networks are formed along the border of Mississippi and Alabama. The 65 kBPD of fuels is composed entirely of gasoline and the overall BEOP value of $81.33/bbl. MI

IA IN

IL KS

OK

TX

MO

""" "" "" " " "" "

PA

OH

LA

"" """ " " " "" " " " " "" " "" " _ "AL ^

" MS " "" " " "" " _" " " " ^ """ """" " " " "" """" " "" _ "" "" ^ "

NJ MD DC

WV

"" _ ^ " "" "" " " KY "" """ " " """ "" _ ^ "" " "

MA CT RI NY

DE

VA

TN AR

NY

NC SC GA Southeast Region (65,000 BPD) Plant size (BPD)

_ ^

1000

10000

FL

_ ^ 0 40 80

Fuel Products Delivery Locations

160

240

320

Natural Gas Feed Source Flow

" " " " " 400

0.000000 - 1.238300 1.238301 - 4.427800 4.427801 - 8.738000

8.738001 - 36.002100

36.002101 - 84.542700

480

560 Miles

Figure 12: Graphical representation of the supply chain case study for the Southeast region with unrestricted fuel ratio (SE-50-U). The Southeast region includes AL, FL, GA, KY, MS, NC, SC, and TN. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

When the minimum portions of fuel ratio are imposed (case study SE-50-R; Figure 13), the supply network selects 12 refineries, 5 of 1 kBPD capacity, 2 of 5 kBPD capacity, and 5 of 10 kBPD capacity, located in the same three states as in case study (SE-50-U). The 5 refineries of 1 kBPD capacity consist of 1 facility with unrestricted fuel production (U), 1 facility that maximizes diesel (D), and 3 facilities that maximize kerosene (K). The 2 facilities of 5 kBPD capacity consist of 1 facility with unrestricted fuel production and 1 facility that maximizes diesel. Finally, the 5 facilities of 10 kBPD capacity consist of 2 facilities with unrestricted fuel production, 2 that maximize diesel, and 1 that maximizes kerosene. The supply chain produces 35.75 kBPD gasoline, 19.50 kBPD diesel, and 9.75 kBPD kerosene and an increase in the total electricity production from 44

ACS Paragon Plus Environment

Page 45 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

0.16 GW to 0.17 GW is observed since the (D) and (K) refineries produce more electricity than the (U) refineries. The BEOP value increases significantly from $81.33/bbl to $90.97/bbl due to primarily the decrease in byproduct propane sale from $1.92/GJ to $1.06/GJ, as the (D) refineries produce less LPG byproduct than the (U) refineries, and the (K) refineries do not produce LPD at all. IA

WI

MO

OK

TX

""" " "" "" " " "" TN

AR

LA

Plant size (BPD) 1000 5000

10000

_ ^

Fuel Products Delivery Locations

CT RI

NY

DE

VA NC

" MS " "" " "" " " " " _" " " ^ """ """"" " " """ "" " "" _ " "" ^ " ^ _ "

Southeast Region (65,000 BPD)

NJ MD DC

WV

"^ " _ " " "" " " " " KY """"" " " """"" _ ^ " "" "

"" """ " " " "" " "" " " "" " "" _"AL ^

Legend

PA

OH

IN

IL

VTMA

NY

MI

SC GA

FL

Natural Gas Feed Source Flow (mscf/day)

" " " " "

0.000000 - 1.588600 1.588601 - 5.325900

5.325901 - 13.384700

13.384701 - 28.800200 28.800201 - 84.542700

0

45 90

180

270

360

450

540

630 Miles

Figure 13: Graphical representation of the supply chain case study for the Southeast region with minimum portions of liquid fuels produced (SE50-R). The Southeast region includes AL, FL, GA, KY, MS, NC, SC, and TN. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

11.2.3

Midwest Region

The Midwest region includes Illinois, Indiana, Michigan, Minnesota, Ohio, and Wisconsin, whose natural gas resources allow the production of 35 kBPD of total fuels. The unrestricted fuel scenario (MW-50-U) yields the selection of 8 refineries, 5 of 1 kBPD capacity and 3 of 10 kBPD capacity, all of which have unrestricted fuel production ratios (U) or only gasoline production. The selected facilities are located in the southeastern part of Michigan but they draw natural gas resources mostly from the northern part of Michigan and the western part of Ohio. Note that the natural 45

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 79

gas resources in this region are located in two clusters that are far away from each other, and the selected refinery locations are placed in between the two clusters (see Figure 14). For case study (MW-50-U), the overall BEOP value of $81.02/bbl with a high contribution from the investment cost of the refineries. In case study (MW-50-R), the supply network also selects 8 facilities with the same locations as in case study (MW-50-U), as shown in Figure 15. One refinery of 1 kBPD capacity maximizes kerosene, 1 refinery of 10 kBPD capacity maximizes diesel, and 1 refinery of 10 kBPD capacity maximizes kerosene. This case study produces a total of 19.25 kBPD gasoline, 7.50 kBPD diesel, and 8.25 kBPD kerosene. Similar to the Southeast region, the increase in the BEOP value from $81.02/bbl to $87.83/bbl is primarily due to the decrease is byproduct propane sale from $1.74/GJ to $0.86/GJ as a lesser number of (U) refineries are selected.

_ ^

Fuel Products Delivery Locations Natural Gas Feed Source

Midwest Region (35,000 BPD) Plant size (BPD) 10000

0 40 80

ND MN

IA

IL

TX

OK

320

400

1.180801 - 3.680700 3.680701 - 7.288000

7.288001 - 12.736500

12.736501 - 47.881300

480

" " " " "" " """" " " " "" " " " "" " "" " MI" " """ " " " " " "" " "^ _ " " " "_ " " ^ " " " "" " "" " " " "" " " " "" " "" " " " " OH " "" " IN "" " """ " " " ""

WI

KS

240

0.000100 - 1.180800

560 Miles

MI

SD

NE

160

Flow

" " " " "

1000

PA MD DC

WV

MO

AR

NY

KY TN MS

AL

GA

NJ DE

VA

SC

NC

Figure 14: Graphical representation of the supply chain case study for the Midwest region with unrestricted fuel ratio (MW-50-U). The Midwest region includes IL, IN, MI, MN, OH, and WI. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

46

ACS Paragon Plus Environment

Page 47 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

0

35 70

140

210

280

350

420

490 Miles

Legend

Midwest Region (35,000 BPD) Plant size (BPD) 1000

10000

_ ^

Fuel Products Delivery Locations

ND MN

MI WI

SD

IA

NE

IL KS

MO

Flow (mscf/day)

" " " " "

0.000100 - 1.025500 1.025501 - 3.139300 3.139301 - 6.160000

6.160001 - 12.736500

12.736501 - 48.278600

" """"" " """ " " " " "" " " " " " "MI" " " " "" " " " " " " " "" " "^ _ " " " "^ " _ " " "" " " """ " " " """ """ " " " " " OH " "" "" "" IN "" """""" " " " " WV KY

OK

Natural Gas Feed Source

NY PA MD

VA

Figure 15: Graphical representation of the supply chain case study for the Midwest region with minimum portions of liquid fuels produced (MW50-R). The Midwest region includes IL, IN, MI, MN, OH, and WI. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

47

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

11.2.4

Southwest Region

The Southwest region is the region with the most abundant natural gas resources, allowing the regional supply chain to produce 2,000 kBPD total fuels. The region includes Arkansas, Louisiana, New Mexico, Oklahoma, and Texas. In case study (SW-50-U), the optimization model yields the selection of 80 refineries, 5 of 1 kBPD capacity, 1 of 5 kBPD capacity, 54 of 10 kBPD capacity, 17 of 50 kBPD capacity, and 3 of 200 kBPD capacity. Most of the fuels produced are gasoline (1,745 kBPD) and the remaining 255 kBPD is kerosene produced by the 14 refineries of 10 kBPD capacity and 4 refineries of 50 kBPD capacity that maximize kerosene production. The overall BEOP value for the case study is $65.60/bbl, the lowest out of all six regions, with output electricity of 5.57 GW. Figure 16 shows that the locations of the 80 refineries span across the entire region. Clusters of GTL refineries form along the borders of the states, with the highest density along the border between Texas and Louisiana. These locations correspond with the oil refineries locations and the GTL liquid fuels are delivered to oil refineries across the five Southwestern states plus one location in Mississippi. Natural gas is mostly drawn from across Arkansas, Oklahoma, Louisiana, and Texas. When the minimum portions of fuel ratio are imposed (SW-50-R), the supply network selects 88 facilities, 5 of 1 kBPD capacity, 1 of 5 kBPD capacity, 64 of 10 kBPD capacity, 15 of 50 kBPD capacity, and 3 of 200 kBPD capacity. All 3 of the 200 kBPD refineries maximize the production of diesel, and 15 refineries of 10 kBPD capacity and 3 refineries of 50 kBPD capacity maximize kerosene. The fuel composition is composed of 1,325 kBPD gasoline, 450 kBPD diesel, and 225 kBPD kerosene, with a BEOP value of $68.36/bbl. The topological layout of the supply chain is largely similar to case study (SW-50-U) with notable shifts in certain clusters, namely the reduction of refineries around the Texas and New Mexico border and the increase of refineries along the southeast border of Texas (see Figure 17). 11.2.5

Central Region

The Central region, the second largest fuel-producing region, is made up of Iowa, Kansas, Missouri, Nebraska, Colorado, Montana, North Dakota, South Dakota, Utah, and Wyoming. The total fuels produced in this region amounts to 840 kBPD with a BEOP value of $72.46/bbl. Case study

48

ACS Paragon Plus Environment

Page 48 of 79

Page 49 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

CO

UT

_ "^ _" ^ AZ

KS

" NM

_ ^

KY

_ ^ "" " " " "" " " _ ^ " " " " " " " """"^ " " _ ^ " _ " _ " ^ " " " " "" " " " "^ TN _" " _" " OK" """" " "" " "" "" " "^ " " " " "" " " """ " " " " " " " " " " " "" "" " " " """ " """ " AR "" "" ^ _ "" " " " " " _" ^ " """" " " " " " "" " " " " " "" "_ " " _ ^ " """ " " """^ """"" AL _" " "" ^ " " _ ^ "" "" "" "" " " MS """ " _ "" "" "^ _"^ ^ " "" " "" " "^ " "" _ _ " " """" " ^ " _ " " " " " " " "" " " " " " " " " " "" " " " """ " " """" "" " " " " "" " " " TX " "" " "" " LA "" " "" " " """ """ """ " " ""_ FL " "" " """ " "" " " "^ " """ ""^ " " " _ ^ " " _ "" " _ ^ "" """" " " " " " " " _ ^ " _ ^ " " _ ^ " " " _ " " " "^ _ " "" " " " " " """ " " ^ "" " " " " " " "" " " "^ _ Southwest Region (2,000,000 BPD) Natural Gas Feed Source "" "" " Plant size (BPD) Flow _ ^ " " """ 1000 " 0.000200 - 21.616200 5000 " " "" " 21.616201 - 63.652600 10000 " "" " 63.652601 - 146.133700 50000 " " 146.133701 - 444.036800 200000 " 444.036801 - 763.273100

"

"

IN

IL

MO

"

_ ^

Fuel Products Delivery Locations

0 40 80

160

240

320

400

480

560 Miles

Figure 16: Graphical representation of the supply chain case study for the Southwest region with unrestricted fuel ratio (SW-50-U). The Southwest region includes AR, LA, NM, OK, and TX. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

49

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

UT

CO

_ "^

"

KS

"

Page 50 of 79

IL IN

MO

^ _ "" " " " " "" " " _ ^ " " " " " " """"^ " " " " _ " _" " " "^ ^ " " " "" _ " " "_" " ^ ""OK" "" " " """

KY

" "" " " " " " " " " "" " " " " """ " " "" " " " " " """ "" " " " " " AR """ " "" ^ _ "" " NM " " " AZ " _"" ^ " """" "" " "" " "" " " " """" " _ "^ " _ ^ " " " " " """" " " _" " " ^ _ ^ " " " " " " " " " """" MS """"""" _ _""^ ^ " " " "^ " "^ " "" _ _ """"" " ^ " " " " """" _ " " "" " " " " " " " " " " " " " " " "" " " " " " " """" _ ^ "" " " " "" " " " "" " " " TX " "" """" """"" LA"" "" " " " " " " " " " " _""""" " ^ " """ " " "^ _ " " _ """" ^ "" "" "" " _ ^ """ "" " " "^ " _ " _ _" ^ ^ " " " ^ "^ _ """ _ " "" "" " _" """ "" "" ^ " _ "^ _ " " " """ "^ "" " " "^ _ "" """ " _ ^ " """" " "" " Legend " "" Southwest Region (2,000,000 BPD) Natural Gas Feed Source "

^" _

"

Plant size (BPD) 1000 5000

10000 50000

200000

_ ^

Fuel Products Delivery Locations

TN

AL

Flow (mscf/day)

" " " " "

0.000200 - 21.712800

21.712801 - 60.439100

60.439101 - 141.612700

141.612701 - 444.753100

444.753101 - 1037.128200

0 35 70

140

210

280

350

420

490 Miles

Figure 17: Graphical representation of the supply chain case study for the Southwest region with minimum portions of liquid fuels produced (SW50-R). The Southwest region includes AR, LA, NM, OK, and TX. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

50

ACS Paragon Plus Environment

Page 51 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(CE-50-U) selects 30 refineries, 5 of 1 kBPD capacity, 1 of 5 kBPD capacity, 13 of 10 kBPD capacity, 10 of 50 kBPD capacity, and 1 of 200 kBPD capacity. The fuel composition consists of 511.14 kBPD gasoline, 43.85 kBPD diesel, and 285.01 kBPD kerosene, being the only region that produces diesel in the unrestricted case. The diesel portion is produced by the 1 refinery of 1 kBPD capacity that maximizes diesel (D-1) and the 1 refinery of 200 kBPD capacity that produces fuels in United States demand ratio (R-200). The kerosene is produced by the same (R-200) refinery in addition to 7 refineries of 10 kBPD capacity that maximize kerosene (K-10). The locations of the refineries are spread across Montana, Wyoming, Utah, Colorado, and Kansas, and the fuel delivery locations are not only located within the region, but also Nevada, New Mexico, and Texas in the Western and Southwestern regions (see Figure 18). When the minimum portions of fuel ratio are imposed in case study (CE-50-R), the supply network selects 28 refineries, 5 of 1 kBPD capacity, 3 of 5 kBPD capacity, 12 of 10 kBPD capacity, 6 of 50 kBPD capacity, and 2 of 200 kBPD capacity. Four refineries of 1 kBPD capacity and 1 refinery of 200 kBPD maximize diesel, 2 refineries of 50 kBPD capacity maximize kerosene, 2 refineries of 50 kBPD capacity and 1 refinery of 200 kBPD capacity produce fuels in United States demand ratio. The BEOP is $73.34/bbl due to the decrease in electricity production and increase in natural gas transportation cost. The fuel cpmposition is 513.54 kBPD gasoline, 217.64 kBPD diesel, and 108.82 kBPD kerosene. In this case study, less fuel production takes place in the eastern part of Utah and more production takes place in the eastern part of Kansas compared to the unrestricted case (see Figure 19). 11.2.6

Western Region

The Western region has fuel production with the highest BEOP value at $106.27/bbl. The states within the region are Arizona, California, Nevada, Idaho, Oregon, and Washington. A total of 9 refineries are selected, 5 of 1 kBPD capacity, 1 of 5 kBPD capacity, and 3 of 10 kBPD capacity, where all of them produce only gasoline (U refineries). In case study (WE-50-R), 10 refineries are selected, 5 of 1 kBPD capacity, 3 of 5 kBPD capacity, and 2 of 10 kBPD capacity, producing 22 kBPD gasoline, 12 kBPD diesel, and 6 kBPD kerosene. A significant increase in the BEOP value to $123.09/bbl is observed due to the fact that the natural gas resources in the Western region are all located in California and are dispersed from each other (see Figures 20 and 21). The increase 51

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

" "" " " " " " "" "^ " ND _"" MT " " "" " _ ^ "" " MN " ""^ _ " " " ^ "_ " " " " " SD _ "" " " "^ " WY " _" ^ " _" ^ " IA _" " ^ " NE _" " " _ ^ " _ "^ ^ " " " "" " " " "" " "" " "" " _ ^ "" " " " " "" "" UT " "" " " "" " CO " " KS "" " " " " " MO "" _" " " """ """"^ " " """""" " " " " " " " " "" "" " " " " "" " "" " " " " " " "" " " " "" _ ^

WA ID

OR

NV

_ ^

CA

_ ^

_ ^

AZ

_ ^

Plant size (BPD) 1000 5000

10000 50000 200000

0

Flow

" " " " "

65 130

OK

_ ^ NM TX

_ ^

Fuel Products Delivery Locations Natural Gas Feed Source

Central Region (840,000 BPD)

Page 52 of 79

0.000900 - 37.424700

MI MI WI

IL

KY

TN

AR MS LA

IN

AL FL

37.424701 - 126.180600 126.180601 - 238.007900 238.007901 - 418.203000 418.203001 - 1223.397200

260

390

520

650

780

910 Miles

Figure 18: Graphical representation of the supply chain case study for the Central region with unrestricted fuel ratio (CE-50-U). The Central region includes IA, KS, MO, NE, CO, MT, ND, SD, UT, and WY. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

52

ACS Paragon Plus Environment

Page 53 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

OR

""" " " " " " "" "^ " ND _" " MT " " """ _ ^ MN "" " " " ""^ _" " " "_^ "" " " SD _ "" " " "^ _" ^ " " WY " _" ^ " _" " ^ NE _" "" " _ ^ " _"^ ^ " " " " " " "" " "" " "^ " _ " "" "" " "" " CO" UT "" " "" " " " "" " "" " """ ""KS " " " _""" " """ """"^ """""""" " " _ ^ " " "" " "" " """ "" " " "" " """ """" """" " _"" ^

ID

_ NV^

_ ^

CA

_ ^

AZ

_ ^

OK

NM

_ ^

Legend

TX

Central Region (840,000 BPD) Plant size (BPD)

0 50100

200

300

400

500

600

MI WI IA

MO

Fuel Products Delivery Locations

Flow (mscf/day)

" " " " "

GA

AL

LA Natural Gas Feed Source

50000

_ ^

TN

MS

10000

OH

KY

AR

5000

200000

IN

IL

1000

700 Miles

MI

FL

0.000900 - 15.434000

15.434001 - 56.685700

56.685701 - 107.373900

107.373901 - 333.945500 333.945501 - 902.240800

Figure 19: Graphical representation of the supply chain case study for the Central region with minimum portions of liquid fuels produced (CE-50R). The Central region includes IA, KS, MO, NE, CO, MT, ND, SD, UT, and WY. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

53

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 54 of 79

in the overall cost is due to the increase in CO2 transportation and injection costs for sequestration and the decrease in byproduct propane sale. Note that the Western region is the only region that has an average natural gas transportation cost that is higher than the natural gas purchase cost.

WA

SD

ID

OR

"

WY

" "" """ " "" _ ^ " " " "" " CA "" " "

_" ^

UT

NV

_" ^ " _" ^ ^ " _ _ ^

Products

1000 5000

10000

_ ^

Fuel Products Delivery Locations

NE

CO

NM

AZ

"

Western Region (40,000 BPD)

ND

MT

Natural Gas Feed Source Flow

" " " " "

TX

0.000100 - 1.353900 1.353901 - 5.692900 5.692901 - 13.557500 13.557501 - 26.208600 26.208601 - 71.392000

0 55 110

220

330

440

550

660

770 Miles

Figure 20: Graphical representation of the supply chain case study for the Western region with unrestricted fuel ratio (WE-50-U). The Western region includes AZ, CA, NV, ID, OR, and WA. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

11.2.7

Regional Comparisons

The BEOP values for the six regional case studies for both the unrestricted fuel ratio case and the minimum portions of fuels are compared in Figure 22. The Southwest region has the most profitable regional supply chain, followed by the Central, Northeast, Midwest, Southeast, and Western regions. The Western region has a significantly higher BEOP value compared to the other regions due to the limited supply of natural gas, higher price of natural gas, and dispersed nature of the natural gas counties. The BEOP values of the (R) case studies are always higher than the (U) case studies but in varying degree across the six regions. Areas with abundant sources of natural gas, such as the Northeast, Southwest and Central regions have BEOP increases that 54

ACS Paragon Plus Environment

Page 55 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

WA

MT

ID

OR

WY

"

Legend

Western Region (40,000 BPD) Plant size (BPD) 1000 5000

10000

_ ^

Fuel Products Delivery Locations

" "" NV "" " "" _ " _^ ^ " "" " CA "" " " _" _" ^ ^ " "^ _ ^ Natural Gas Feed Source " _ _" ^ Flow (mscf/day)

" " " " "

UT

AZ

CO

NM

0.000100 - 1.353900 1.353901 - 3.799300 3.799301 - 9.260500

9.260501 - 22.833500

22.833501 - 86.614100

0 40 80

160 240 320 400 480 560 Miles

Figure 21: Graphical representation of the supply chain case study for the Western region with minimum portions of liquid fuels produced (WE50-R). The Western region includes AZ, CA, NV, ID, OR, and WA. The facilities and natural gas sources must be located within the region, and the liquid fuels can be delivered to other regions.

55

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

are less than $3/bbl, while the other regions experience a BEOP increase ranging from $6.81/bbl to $16.82/bbl. Note also that the regions with small increases in BEOP experience drops in the electricity production, but the regions with large increases in BEOP experience an increase in the electricity production. The reason for this phenomenon is because the Southwest, Central, and Northeast regional supply chains selected (K) refineries in the unrestricted case studies, which have higher electricity production than the (D) refineries. Thus, when the minimum portions of fuels are imposed and the (D) refineries are selected, the total electricity production is reduced. On the contrary, the Southeast, Midwest, and Western regional supply chains only selected (U) refineries in the unrestricted case studies. When the minimum portions of fuels are imposed, the combined selections of (U), (D), and (K) refineries produce a higher total amount of electricity for the same total volumes of liquid fuels. $140.00   $120.00   $100.00   $80.00   $60.00  

2,000,000.00  

$0.00  

Northeast   Southeast   Midwest   Southwest   Central   Region  and  Total  Volume  of  Fuels  (BPD)  

40,000.00  

Fuel  RaRo   840,000.00  

$20.00   35,000.00  

Unrestricted  

65,000.00  

$40.00  

170,000.00  

BEOP  ($/bbl)  

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 56 of 79

Western  

Figure 22: BEOP comparisons for the six regional supply chains for the unrestricted and minimum fuel ratio case studies.

The breakdown of the cost contributions for each regional supply chain is shown in Figures 23 and 24. It can be readily seen that the investment cost, natural gas purchase and transportation account for 90% of the total cost or more. Generally, the contribution from the natural gas transportation cost is less than the natural gas purchase cost except in the Western region due to the separated clusters of natural gas sources. The regional supply chains that locate the refineries very 56

ACS Paragon Plus Environment

Page 57 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

close to the oil refineries have less cost contributions from the product transportation segment. 100%   90%   80%  

CO2  injecBon  

70%  

CO2  trans  

60%  

Water  purchase  and  trans.  

50%  

Butane  purchase   Product  trans.  

40%  

Investment    

30%  

NG  trans.  

20%  

NG  purchase  

10%   0%   Northeast  Southeast   Midwest   Southwest   Central   Western  

Figure 23: Comparison of the cost contributions in the six regional supply chains for the unrestricted case studies.

11.3

Statewide GTL Supply Chains: Texas

The last set of supply chain case studies limits the geographical scope of the problem further from United States regions to individual states. In these statewide supply chain problems, similar restrictions as in the regional case studies are imposed, namely only natural gas resources within the state can be used, the facilities must be located within the state, and the fuels may be delivered to other states. However, one additional restriction is that only counties that can supply entirely the natural gas requirement of the selected refinery can be selected. In other words, no inter-county natural gas transportation is allowed. The statewide supply chain case studies are applied to the state of Texas, the state with the most abundant natural gas sources and the largest number of counties that produce natural gas (see Table 1). Table 8 summarizes the results from the four case studies for Texas, where 50% natural gas expansion level is applied in all cases. Very low BEOP values are observed for Texas, ranging from $55.58/bbl to $65.05/bbl. Using the restrictions described above, Texas can produce a total of 200 kBPD liquid fuels when no inter-county natural gas transportation is allowed. Note that the 57

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

100%   90%   80%  

CO2  injecBon  

70%  

CO2  trans  

60%  

Water  purchase  and  trans.  

50%  

Butane  purchase   Product  trans.  

40%  

Investment    

30%  

NG  trans.  

20%  

NG  purchase  

10%   0%   Northeast  Southeast   Midwest   Southwest   Central   Western  

Figure 24: Comparison of the cost contributions in the six regional supply chains for the minimum portions of liquid fuels case studies.

state can produce much more than 200 kBPD if all of its natural gas resources are used. Case study (TX-50-U), with a BEOP value of $55.58/bbl produces 200 kBPD of gasoline in 5 refineries of 1 kBPD capacity, 1 refinery of 5 kBPD capacity, 9 refineries of 10 kBPD capacity, and 2 refineries of 50 kBPD capacity. The number of facilities is the same in case study (TX-50-R) when 136.25 kBPD gasoline, 41.25 kBPD diesel, and 22.50 kBPD kerosene are produced. The BEOP increases to $60.21/bbl due to the higher investment costs of (D) and (K) refineries and the reduction in byproduct LPG sales. The geographical layout of the statewide supply chain can be seen in Figure 25 and the same locations are selected for case study (TX-50-R). The counties that can provide the natural gas resources without inter-county natural gas transportation are located in the northern and eastern part of Texas, and the fuels are delivered to oil refineries along the northern, eastern, and southeastern border of Texas. When only 10 kBPD refineries are allowed to exist in the supply chain, the BEOP values for the unrestricted and minimum portions of liquid fuels are $59.01/bbl and $65.05/bbl, respectively. In both cases, 20 refineries are selected, with all 20 producing gasoline in case study (TX-50-U10) and 12 refineries with unrestricted fuel production, 5 refineries that maximize diesel, and 3 refineries that maximize kerosene in case study (TX-50-R-10). Compared to the case studies that

58

ACS Paragon Plus Environment

Page 58 of 79

Page 59 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 8: Statewide GTL supply chain profiles for Texas. Case Study

TX-50-U

TX-50-R

TX-50-U-10

TX-50-R-10

BEOP ($/bbl)

$55.58

$60.21

$59.01

$65.05

Levelized cost ($/GJ)

$12.63

$12.83

$13.25

$13.72

Total number of facilities

17

17

20

20

1 kBPD

5

5

0

0

U

5

5

0

0

D

0

0

0

0

K

0

0

0

0

R

0

0

0

0

5 kBPD

1

1

0

0

U

1

0

0

0

D

0

1

0

0

K

0

0

0

0

R

0

0

0

0

10 kBPD

9

9

20

20

U

9

6

20

12

D

0

0

0

5

K

0

3

0

3

R

0

0

0

0

50 kBPD

2

2

0

0

U

2

1

0

0

D

0

1

0

0

K

0

0

0

0

R

0

0

0

0

200 kBPD

0

0

0

0

U

0

0

0

0

D

0

0

0

0

K

0

0

0

0

R Electricity produced (GW) Total fuels (kBPD)

0

0

0

0

0.53

0.55

0.51

5.52

200

200

200

200

200.00

136.25

200.00

140.00

Diesel (kBPD)

0

41.25

0

37.50

Kerosene (kBPD)

0

22.50

0

22.50

Total LHV output

1.10

1.14

1.10

1.14

0.58

0.58

0.57

0.58

Gasoline (kBPD)

(million GJ/day) Natural gas usage (trillion SCF/yr) Average cost breakdown ($/GJ) Natural gas purchase

6.08

5.90

5.99

0

0

0

0

Investment cost

7.89

7.61

8.36

8.13

Product trans.

0.76

0.73

0.74

0.72

Water purchase and trans.

0.01

0.01

0.01

0.01

CO2 trans. and injection

0.69

0.65

0.89

1.12

Propane sale

1.97

1.28

1.97

1.34

Electricity sale

0.82

0.81

0.77

0.77

Natural gas trans.

59

ACS Paragon Plus Environment

5.86

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

UT

CO

KS

NM

AZ

IL

MO

AR

_ ^ " " TX

"

0

MS

_ " "" _ ^ "^ " ""

" ^ _ " _ ^ " _ ^

KY TN

OK

" " _ ^ "

Page 60 of 79

40 80

AL

LA

FL

^ _ _ ^ "

160

240

320

400

480

560 Miles

Legend

Texas (200,000 BPD) - Unrestricted Fuel Ratio Natural Gas Feed Source Plant size (BPD) 1000 5000

10000 50000

_ ^

Fuel Products Delivery Locations

Flow (mscf/day)

" " " "

8.738000

8.738001 - 43.690000

43.690001 - 86.614100

86.614101 - 444.753100

Figure 25: Graphical representation of the supply chain case study for the state of Texas with unrestricted fuel ratio (TX-50-U). The facilities and natural gas sources must be located within the state, and the liquid fuels can be delivered to other states. The facilities must be located at a location where enough natural gas is available and no inter-county natural gas transportation is allowed.

60

ACS Paragon Plus Environment

Page 61 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

include all refinery capacities, these case studies have higher levelized investment costs and lower electricity production, but similar natural gas usage. Figure 26 displays the geographical layout of the statewide supply chain, showing that more facilities are selected in the northern and eastern parts of Texas as well as four additional refineries along the southwestern border of Texas. Note that case study (TX-50-R-10) selects the same locations for its 20 refineries. UT

CO

KS

NM

AR

_ ^ _ ^ "

"" "" TX

"" ""

"

MS

_ _ "" ^ "^ " ""

_ ^ Legend

KY TN

OK

" _ " ^ AZ

IL

MO

AL

LA

FL

_ ^

_ ^

Plant size (BPD)

_ ^

"

10000

_ ^

Fuel Products Delivery Locations

"

86.614100

Flow (mscf/day)

0

40 80

160

240

320

400

480

560 Miles

"

Figure 26: Graphical representation of the supply chain case study for the state of Texas with unrestricted fuel ratio and only 10 kBPD facilities allowed (TX-50-U-10). The facilities and natural gas sources must be located within the state, and the liquid fuels can be delivered to other states. The facilities must be located at a location where enough natural gas is available and no inter-county natural gas transportation is allowed.

11.4

Overall Comparison

The importance of the geographical scope of the supply chain studies is highlighted in the average cost or BEOP values obtained from the nationwide, regional, and statewide case studies. Comparing the nationwide case study (N-50-U) with the regional equivalent of the formulation (i.e., case studies NE-50-U, SE-50-U, MW-50-U, SW-50-U, CE-50-U, and WE-50-U), the BEOP value is lower for the nationwide supply chain compared to the individual regional case studies

61

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(i.e., $65.05/bbl compared to $65.60/bbl-$106.27/bbl). In the nationwide formulation, there is no restriction in the flow rates across regional and state boundaries and thus a more optimal network can be achieved. However, the benefit of analyzing the regional case studies lies in the practicality of implementing a regional management of potential GTL supply chains. Extending the pattern further to the statewide supply chain case studies, it can be observed that the cost of fuel production can vary from state to state or from region to region depending on the amount of local resources and the geographical layout of the natural gas counties and the existing oil refineries. The state of Texas is an ideal state for GTL refineries due to the synergistic combination between all of the supply chain components, yielding an overall BEOP ($55.58/bbl) that is lower than the nationwide average.

12

Conclusions

This paper proposed and reported results on the nationwide, regional, and statewide GTL supply chains for the United States transportation sector. Multiple segments of the supply chain from the acquisition of the natural gas feedstock, transportation costs of feedstock and fuel products, investment costs of GTL refineries, and environmental constraints such as water, electricity, and CO2 sequestration resources are taken into account simultaneously. The optimization problem is formulated as a MILP that minimizes the total cost of the supply chain and yields the locations and capacities of GTL refineries as well as feedstock and product connections to their respective locations. This MILP is solved for the nationwide supply chain, six United States regions, and for the state of Texas. Two modes of fuel production are applied, one in which there is no restriction in the fuel ratios and one in which minimum portions of liquid fuels are imposed. The resulting average costs of fuel production for the case studies suggest that GTL supply chains can be highly profitable in the United States. The nationwide case studies have BEOP values ranging from $65.06/bbl-$80.25/bbl and the regional case studies have BEOP values from $65.60/bbl-$90.97/bbl, except for the Western region at $106.27/bbl-$123.09/bbl. Results show that the Southwest and Central regions of the United States are the most profitable regions for GTL supply chains, and that a state GTL supply chain can have a high economic performance in the state of Texas. 62

ACS Paragon Plus Environment

Page 62 of 79

Page 63 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Acknowledgments The authors acknowledge partial financial support from the National Science Foundation (NSF EFRI-0937706; CBET-1158849).

Notation Index s = State Index c = County Index f = Feedstock Index p = Product Index t = Plant Size Index l = Plant Location Index m = Transportation Mode Index r = Port Index q = Carbon Conversion Level Index

Sets S = U.S. States C = U.S. Counties F G = Natural Gas Feedstocks F = Feedstocks P = Products (i.e., gasoline, diesel, kerosene) CF = Feedstock-County Pairs CP = Product-County Pairs CSQ = Carbon Sequestration Sites CW = Locations of Water Resources LF = Candidate Plant Locations LS = Candidate Locations with Solar Plant Capacities

63

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

LW = Candidate Locations with Wind Plant Capacities M = Modes of Transportation (i.e., truck, rail, pipeline, barge) MF = Feedstock-Mode Pairs MP = Product Mode Pairs G MPipe = Natural Gas Pipelines P = Product Pipelines MPipe

R = 50 U.S. ports with high liquid fuels capacity T = Plant Sizes (i.e., 1 kBD, 5 kBD, 10 kBD, 50 kBD, 200 kBD) Q = Fuel Product Ratios (i.e., commensurate with US demands, maximized diesel, maximized kerosene production, unrestricted) EFL = Enumerated Facility Locations FL = Candidate Facility Locations FT = Feasible Feedstock Flow Quadruplets PT = Feasible Product Flow Quadruplets

64

ACS Paragon Plus Environment

Page 64 of 79

Page 65 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Set Definitions

f ∈ F = FG FC = {( f ) | f ∈ F G } MF = {( f , m) ∈ F G × pipeline} MP = {(p, m) ∈ P × {pipeline,truck,barge}} EFL = {( f ,t, l, q) ∈ FC × T × L × Q} FT = {( f , c, l, m) | (p, c) ∈ CF , l ∈ L, FE ( f , m) ∈ MF ,Cost FE f ,c,l,m ≤ MaxCost f }

PT = {(p, l, c, m) | (p, c) ∈ CP , l ∈ L, PE (p, m) ∈ MP ,Cost p,c,l,m ≤ MaxCost pPE }

CF = {( f , c) | ( f , c) ∈ F ×C, FA f ,c > 0} CP = {(p, c) | (p, c) ∈ P ×C, DM p,c > 0} CW = {c | c ∈ C,WAc > 0} CSQ = {c | c ∈ C, SCAPc > 0} LS = {l | l ∈ LF ,CapSl > 0} LW = {l | l ∈ LF ,CapW l > 0} Parameters N = Maximum number of GTL plants built in the United States Ntmax = Maximum number of GTL plants for size t Ntmin = Minimum number of GTL plants for size t LC f ,t,q = GTL levelized investment cost for feed ( f ), size t, and fuel ratio q FRGf,t,q = GTL natural gas requirement for feed ( f ), size t, and fuel ratio q ER f ,t,q = GTL electricity requirement for feed ( f ), size t, and fuel ratio q EPf ,t,q = GTL electricity produced for feed ( f ), size t, and fuel ratio q FW f ,t,q = GTL freshwater requirement for feed ( f ), size t, and fuel ratio q

65

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

SQ f ,t,q = GTL CO2 sequestration flow for feed ( f ), size t, and fuel ratio q FA f ,c = Availability of feedstock f in county c DM p,c = Demand of product p in county c PR p,t,q = Amount of liquid fuel product p for a plant of size t and fuel ratio q WAc = Water availability in location c SCAPc = CO2 sequestration capacity in location c MaxCost FE f = Maximum delivered cost of feedstock f MaxCost pPE = Maximum delivered cost of product p CapSl = Maximum capacity of solar generated electricity at location l CapW l = Maximum capacity of wind generated electricity at location l CapG = Maximum grid electricity usage for the country Cost Ff,c = Cost per unit mass of feedstock f at county c Cost FT f ,c,l,m = Cost per unit mass flow to transport feedstock f from county c to facility l using mode m PT Cost p,l,c,m = Cost per unit mass flow to transport product p from facility l to county c using mode

m CostEl,G = Total investment cost per unit energy of grid electricity CostEl,P = Total profit per unit energy of produced electricity CostEl,S = Total investment cost per unit energy of solar electricity CostEl,W = Total investment cost per unit energy of wind electricity CostcW P = Cost of water purchase at location c per unit flow rate W T = Cost of water transportation by pipeline from source c to facility l Costc,l CO2 ,T Costc,l = Cost of CO2 transportation by pipeline from facility l to sequestration site c CO2 ,In j

Costc

= Cost of CO2 injection at sequestration cite c

Continuous Variables CostlI = Levelized investment cost of facility l FR f ,l = Amount of feedstock f required at facility l EllT = Total electricity required at facility l EllP = Total electricity produced at facility l 66

ACS Paragon Plus Environment

Page 66 of 79

Page 67 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

EllS = Solar electricity required at facility l EllW = Wind electricity required at facility l EllG = Grid electricity required at facility l W Fl = Freshwater requirement for facility l SFl = Sequestered CO2 amount from facility l wc,l = Freshwater flow from source c to facility l seqc,l = CO2 sequestration flow from facility l to injection site c x f ,c,l,m = Flow of feedstock f from county c to facility l using transportation mode m z p,l,c,m = Flow of product p from facility l to county c using transportation mode m

Binary Variables y f ,t,l,q = GTL plant binary at location l with feed ( f ), size t, and fuel ratio q

67

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References [1] Energy Information Administration, Annual Energy Outlook 2013 with Projections to 2035. Document Number: DOE/EIA-0383(2012), http://www.eta.doe.gov/oiaf/aeo/, 2011. [2] Energy Information Administration, Monthly Energy Review – November 2012. Document Number: DOE-EIA-0035(2012/11), http://www.eia.gov/totalenergy/data/ monthly/pdf/mer.pdf, 2012. [3] Floudas, C. A.; Elia, J. A.; Baliban, R. C. Hybrid and Single Feedstock Energy Processes for Liquid Transportation Fuels: A Critical Review. Comp. Chem. Eng. 2012, 41, 24–51. [4] Agrawal, R.; Singh, N. R.; Ribeiro, F. H.; Delgass, W. N. Sustainable fuel for the transportation sector. PNAS 2007, 104, 4828–4833. [5] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Toward Novel Biomass, Coal, and Natural Gas Processes for Satisfying Current Transportation Fuel Demands, 1: Process Alternatives, Gasification Modeling, Process Simulation, and Economic Analysis. Ind. Eng. Chem. Res. 2010, 49, 7343–7370. [6] Elia, J. A.; Baliban, R. C.; Floudas, C. A. Toward Novel Biomass, Coal, and Natural Gas Processes for Satisfying Current Transportation Fuel Demands, 2: Simultaneous Heat and Power Integration. Ind. Eng. Chem. Res. 2010, 49, 7371–7388. [7] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Optimization Framework for the Simultaneous Process Synthesis, Heat and Power Integration of a Thermochemical Hybrid Biomass, Coal, and Natural Gas Facility. Comp. Chem. Eng. 2011, 35, 1647–1690. [8] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Simultaneous Process Synthesis, Heat, Power, and Water Integration of Thermochemical Hybrid Biomass, Coal, and Natural Gas Facilities. Comp. Chem. Eng. 2012, 37, 297–327. [9] Water Science and Technology Board, Water Implications of Biofuels Production in the United States, 2008. 68

ACS Paragon Plus Environment

Page 68 of 79

Page 69 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[10] Lynd, L. R.; Larson, E.; Greene, N.; Laser, M.; Sheehan, J.; Dale, B. E.; McLaughlin, S.; Wang, M. The role of biomass in America’s energy future: framing the analysis. Biofuels, Bioprod. Biorefin. 2009, 3, 113–123. [11] National Academy of Sciences and National Academy of Engineering and National Research Council, Liquid Transportation Fuels from Coal and Biomass: Technological Status, Costs, and Environmental Issues. Prepublication, Washington, D. C., EPA, 2009. [12] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Biomass to Liquid Transportation Fuels (BTL) Systems: Process Synthesis and Global Optimization Framework. Energy Environ. Sci. 2013, 6, 267–287. [13] Baliban, R. C.; Elia, J. A.; Floudas, C. A.; Gurau, B.; Weingarten, M. B.; Klotz, S. D. Hardwood Biomass to Gasoline, Diesel, and Jet Fuel: 1. Process Synthesis and Global Optimization of a Thermochemical Refinery. Energy & Fuels 2013, In press. doi:10.1021/ef302003f. [14] Department of Energy, Biomass as Feedstock for a Bioenergy and Bioproducts Industry: The Technical Feasibility of a Billion-Ton Annual Supply. Document Number: DOE/GO102005-2135, http://www1.eere.energy.gov/biomass/publications.html, 2005. [15] National Research Council, Water Implications of Biofuels Production in the United States, http://www.nap/edu/catalog/12039.html, 2008. [16] de Fraiture, C.; Giordano, M.; Liao, Y. Biofuels and implications for agricultural water use: blue impacts of green energy. Water Policy 2008, 10, 67–81. [17] Kreutz, T. G.; Larson, E. D.; Liu, G.; Williams, R. H. Fischer-Tropsch Fuels from Coal and Biomass, Proceedings of the 25th International Pittsburg Coal Conference, 2008. [18] de Klerk, A. Fischer-Tropsch Refining; Wiley-VCH Verlag & Co. KGaA, 2011. [19] Sasol, for

an

Sasol commences the front-end engineering and design (FEED) phase integrated

gas-to-liquids

and

ethane

cracker

sasollouisianaprojects.com/, 2012.

69

ACS Paragon Plus Environment

complex,

http://www.

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[20] Peng, X. D.; Wang, A. W.; Toseland, B. A.; Tij, P. J. A. Single-step syngas-to-dimethyl ether processes for optimal productivity, minimal emissions, and natural gas-derived syngas. Ind. Eng. Chem. Res. 1999, 38, 4381–4388. [21] Sudiro, M.; Bertucco, A. Synthetic Fuels by a Limited CO2 Emission Process Which Uses Both Fossil and Solar Energy. Energy and Fuels 2007, 21, 3668–3675. [22] Cao, Y.; Gao, Z.; Jin, J.; Zhou, H.; Cohron, M.; Zhao, H.; Liu, H.; Pan, W. Synthesis Gas Production with an Adjustable H2 /CO Ratio through the Coal Gasification Process: Effects of Coal Ranks and Methane Addition. Energy and Fuels 2008, 22, 1720–1730. [23] Zhou, L.; Hu, S.; Li, Y.; Zhou, Q. Study on co-feed and co-production system based on coal and natural gas for producing DME and electricity. Chem. Eng. J. 2008, 136, 31–40. [24] Sudiro, M.; Bertucco, A. Production of synthetic gasoline and diesel fuel by alternative processes using natural gas and coal: Process simulation and optimization. Energy 2009, 34, 2206–2214. [25] Zhou, L.; Hu, S.; Chen, D.; Li, Y.; Zhu, B.; Jin, Y. Study on systems based on coal and natural gas for producing dimethyl ether. Ind. Eng. Chem. Res. 2009, 48, 4101–4108. [26] Adams II, T. A.; Barton, P. I. Combining coal gasification and natural gas reforming for efficient polygeneration. Fuel Proc. Tech. 2011, 92, 639–655. [27] Adams II, T. A.; Barton, P. I. Combining coal gasification, natural gas reforming, and solid oxide fuel cells for efficient polygeneration with CO2 capture and sequestration. Fuel. Proc. Tech. 2011, 92, 2105–2115. [28] Li, Z.; Liu, P.; He, F.; Wang, M.; Pistikopoulos, E. N. Simulation and exergoeconomic analysis of a dual-gas sourced polygeneration process with integrated methanol/DME/DMC catalytic synthesis. Comp. Chem. Eng. 2011, 35, 1857–1862. [29] Borgwardt, R. H. Biomass and natural gas as co-feedstocks for production of fuel for fuelcell vehicles. Biomass and Bioenergy 1997, 12, 333–345.

70

ACS Paragon Plus Environment

Page 70 of 79

Page 71 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[30] Dong, Y.; Steinberg, M. Hynol–an economical process for methanol production from biomass and natural gas with reduced CO2 emission. Int. J. Hydrogen Energy 1997, 22, 971–977. [31] Li, H.; Hong, H.; Jin, H.; Cai, R. Analysis of a feasible polygeneration system for power and methanol production taking natural gas and biomass as materials. Appl. Energy 2010, 87, 2846–2853. [32] Liu, G.; Williams, R. H.; Larson, E. D.; Kreutz, T. G. Design/economics of low-carbon power generation from natural gas and biomass with synthetic fuels co-production. Energy Procedia 2011, 4, 1989–1996. [33] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Biomass and Natural Gas to Liquid Transportation Fuels: Process Synthesis, Global Optimization, and Topology Analysis. Ind. Eng. Chem. Res. 2013, 52, 3381–3406. [34] Baliban, R. C.; Elia, J. A.; Misener, R.; Floudas, C. A. Global Optimization of a MINLP Process Synthesis Model for Thermochemical Based Conversion of Hybrid Coal, Biomass, and Natural Gas to Liquid Fuels. Comp. Chem. Eng. 2012, 42, 64–86. [35] Baliban, R. C.; Elia, J. A.; Weekman, V. W.; Floudas, C. A. Process synthesis of hybrid coal, biomass, and natural gas to liquids via Fischer-Tropsch synthesis, ZSM-5 catalytic conversion, methanol synthesis, methanol-to-gasoline, and methanol-to-olefins/distillate technologies. Comp. Chem. Eng. 2012, 47, 29–56. [36] Dong, L.; Wei, S.; Tan, S.; Zhang, H. GTL or LNG: Which is the best way to monetize ”stranded” natural gas? Petroleum Science 2008, 5, 388–394. [37] National Energy Technology Laboratory, Oilfield Flare Gas Electricity Systems (OFFGASES) Project, DOE Award No.: DE-FC26-02NT15444, 2008. [38] Wood, D. A.; Nwaoha, C.; Towler, B. F. Gas-to-liquids (GTL): A review of an industry offering several routes for monetizing natural gas. Journal of Natural Gas Science and Engineering 2012, 9, 196–208.

71

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[39] Khalilpour, R.; Karimi, I. A. Evaluation of utilization alternatives for stranded natural gas. Energy 2012, 40, 317–328. [40] Elia, J. A.; Baliban, R. C.; Xiao, C. A., X. Floudas Optimal Energy Supply Network Determination and Life Cycle Analysis for Hybrid Coal, Biomass, and Natural Gas to Liquid (CBGTL) Plants Using Carbon-based Hydrogen Production. Comp. Chem. Eng. 2011, 35, 1399–1430. [41] Elia, J. A.; Baliban, R. C.; Floudas, C. A. Nationwide Supply Chain Analysis for Hybrid Energy Processes with Significant CO2 Emissions Reduction. AIChE J. 2012, 58, 2142– 2154. [42] Parker, N.; Tittmann, P.; Hart, Q.; Nelson, R.; Skog, K.; Schmidt, A.; Gray, E.; Jenkins, B. Development of a biorefinery optimized biofuel supply curve for the Western United States. Biomass and Bioenergy 2010, 34, 1597–1607. [43] Gebreslassie, B. H.; Yao, Y.; You, F. Design under uncertainty of hydrocarbon biorefinery supply chains: Multiobjective stochastic programming models, decomposition algorithm, and a comparison between CVaR and downside risk. AIChE Journal 2012, 58, 2155–2179. [44] You, F.; Wang, B. Life Cycle Optimization of Biomass-to-Liquid Supply Chains with Distributed-Centralized Processing Networks. Ind. Eng. Chem. Res. 2011, 50, 10102– 10127. [45] You, F.; Tao, L.; Graziano, D. J.; Snyder, S. W. Optimal design of sustainable cellulosic biofuel supply chains: Multiobjective optimization coupled with life cycle assessment and input-output analysis. AIChE Journal 2012, 58, 1157–1180. [46] Marvin, W. A.; Schmidt, L. D.; Benjaafar, S.; Tiffany, D. G.; Daoutidis, P. Economic optimization of a lignocellulosic biomass-to-ethanol supply chain. Chem. Eng. Sci. 2012, 67, 68–79. [47] Marvin, W. A.; Schmidt, L. D.; Daoutidis, P. Biorefinery Location and Technology Selection Through Supply Chain Optimization. Ind. Eng. Chem. Res. 2013, 52, 3192–3208.

72

ACS Paragon Plus Environment

Page 72 of 79

Page 73 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[48] Daoutidis, P.; Marvin, W. A.; Rangarajan, S.; Torres, A. I. Engineering Biomass Conversion Processes: A Systems Perspective. AIChE Journal 2013, 59, 3–18. [49] An, H.; Wilhelm, W. E.; Searcy, S. W. Biofuel and petroleum-based fuel supply chain research: A literature review. Biomass and Bioenergy 2011, 35, 3763–3774. [50] Awudu, I.; Zhang, J. Uncertainties and sustainability concepts in biofuel supply chain management: A review. Renewable & Sustainable Energy Reviews 2012, 16, 1359–1368. [51] Kim, J.; Realff, M. J.; Lee, J. H.; Whittaker, C.; Furtner, L. Design of biomass processing network for biofuel production using an MILP model. Biomass and Bioenergy 2011, 35, 853–871. [52] Kim, J.; Realff, M. J.; Lee, J. H. Optimal design and global sensitivity analysis of biomass supply chain networks for biofuels under uncertainty. Comp. Chem. Eng. 2011, 35, 1738– 1751. [53] Liu, P.; Whitaker, A.; Pistikopoulos, E. N.; Li, Z. A mixed-integer programming approach to strategic planning of chemical centers: A case study in the UK. Comp. Chem. Eng. 2011, 35, 1359–1373. [54] Papapostolou, C.; Kondili, E.; Kaldellis, J. K. Development and implementation of an optimisation model for biofuels supply chain. Energy 2011, 36, 6019–6026. [55] Sharma, P.; Sarker, B. R.; Romagnoli, J. A. A decision support tool for strategic planning of sustainable biorefineries. Comp. Chem. Eng. 2011, 35, 1767–1781. [56] Andersen, F.; Iturmendi, F.; Espinosa, S.; Diaz, M. S. Optimal design and planning of biodiesel supply chain with land competition. Comp. Chem. Eng. 2012, 47, 170–182. [57] Walther, G.; Schatka, A.; Spengler, T. S. Design of regional production networks for second generation synthetic bio-fuel - A case study in Northern Germany. European Journal of Operational Research 2012, 218, 280–292. [58] Hamedi, M.; Farahani, R. Z.; Husseini, M. M.; Esmaeilian, G. R. A distribution planning model for natural gas supply chain: A case study. Energy Policy 2009, 37, 799–812. 73

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[59] Selot, A.; Kuok, L. K.; Robinson, M.; Mason, T. L.; Barton, P. I. A short-term operational planning model for natural gas production systems. AIChE Journal 2008, 54, 495–515. [60] Contesse, L.; Ferrer, J. C.; Maturana, S. A mixed-integer programming model for gas purchase and transportation. Annals of Operations Research 2005, 139, 39–63. ¨ [61] Ozelkan, E. C.; D’Ambrosio, A.; Teng, S. G. Optimizing liquefied natural gas terminal design for effective supply-chain operations. International Journal of Production Economics 2008, 111, 529–542. [62] dos Santos, S. P.; Leal, J. E.; Oliveira, F. The development of a natural gas transportation logistics management system. Energy Policy 2011, 39, 4774–4784. [63] Baliban, R. C.; Elia, J. A.; Floudas, C. A. Novel natural gas to liquids (GTL) technologies: Process synthesis and global optimization strategies. AIChE J. 2013, 59, 505–531. [64] National Energy Technology Laboratory, Quality Guidelines for Energy System Studies, 2004. [65] Energy Information Administration, Natural Gas Processing: The Crucial Link Between Natural Gas Production and Its Transportation to Market, http://www.eia.gov/pub/ oil_gas/natural_gas/feature_articles/2006/ngprocess/ngprocess.pdf, 2006. [66] Fox, J. M.; T.-P., C.; Degen, B. D. Direct Methane Conversion Process Evaluations, YEAR = 1988, HOWPUBLISHED = Contract No. DE-AC22-87PC79814. [67] Gradassi, M. J.; Green, N. W. Economics of natural gas conversion processes. Fuel Proc. Technology 1995, 42, 65–83. [68] Iandoli, C. L.; Kjelstrup, S. Exergy analysis of a GTL process based on low-temperature slurry F-T reactor technology with a cobalt catalyst. Energy & Fuels 2007, 21, 2317–2324. [69] Gao, L.; Li, H.; Chen, B.; Jin, H.; Lin, R.; Hong, H. Proposal of a natural gas-based polygeneration system for power and methanol production. Energy 2008, 33, 206–212. [70] Hao, X.; Djatmiko, M. E.; Xu, Y.; Wang, Y.; Chang, J.; Li, Y. Simulation analysis of a gas-to-liquid process using Aspen Plus. Chem. Eng. Technol. 2008, 31, 188–196. 74

ACS Paragon Plus Environment

Page 74 of 79

Page 75 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[71] Lee, C. J.; Lim, Y.; Kim, H. S.; Han, C. Optimal gas-to-liquid product selection from natural gas under uncertain price scenarios. Ind. Eng. Chem. Res. 2009, 48, 794–800. [72] Kim, Y. H.; Jun, K. W.; Joo, H.; Han, C.; Song, I. K. A simulation study on gas-to-liquid (natural gas to Fischer-Tropsch synthetic fuel) process optimization. Chem. Eng. J. 2009, 155, 427–432. [73] Bao, B.; El-Halwagi, M. M.; Elbashir, N. O. Simulation, integration, and economic analysis of gas-to-liquid process. Fuel Proc. Technol. 2010, 91, 703–713. [74] Dillerop, C.; van der Berg, H.; van der Ham, A. G. J. Novel syngas production techniques for GTL-FT synthesis of gasoline using reverse flow catalytic membrane reactors. Ind. Eng. Chem. Res. 2010, 49, 12529–12537. [75] Ha, K. S.; Bae, J. W.; Woo, K. J.; Jun, K. W. Efficient utilization of greenhouse gas in a gas-to-liquids process combined with carbon dioxide reforming of methane. Environ. Sci. Technol. 2010, 44, 1412–1417. [76] Heimel, S.; Lowe, C. Technology comparison of CO2 capture for a gas-to-liquids plant. Energy Procedia 2009, 1, 4039–4046. [77] Bin, C.; Hingguang, J.; Lin, G. System study on natural gas-based polygeneration system of DME and electricity. Int. J. Energy Res. 2008, 32, 722–734. [78] Hall, K. R. A new gas to liquids (GTL) or gas to ethylene (GTE) technology. Catal. Today 2005, 106, 243–246. [79] Suzuki, S.; Sasaki, T.; Kojima, T. New process development of natural gas conversion technology to liquid fuels via OCM reaction. Energy & Fuels 1996, 10, 531–536. [80] Horstman, D.; Abata, D.; Keith, J.; Oberto, L. Feasibility study of an on-board natural gas to dimethyl ether reactor for dimethyl ether preinjection and enjanced ignition. J. of Eng. for Gas Turbines and Power 2005, 127, 909–917. [81] Erturk, M. Economic analysis of unconventional liquid fuel sources. Renew. Sustain. Energy Reviews 2011, 15, 2766–2771. 75

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[82] Vliet, O.; Faaij, A.; Turkenburg, W. Fischer-Tropsch diesel production in a well-wheel perspective: A carbon, energy flow and cost analysis. Energy Conversion and Management 2009, 50, 855–876. [83] Baliban, R. C.; Elia, J. A.; Floudas, C. A.; Xiao, X.; Zhang, Z.; Li, J.; Cao, H.; Ma, J.; Qiao, Y.; Hu, X. Thermochemical Conversion of Duckweed Biomass to Gasoline, Diesel, and Jet Fuel: Process Synthesis and Global Optimization. Ind. Eng. Chem. Res. 2013, In press. doi:10.1021/ie3034703. [84] Martin, M.; Grossmann, I. E. Process Optimization of FT-Diesel Production from Lignocellulosic Switchgrass. Ind. Eng. Chem. Res. 2011, 50, 13485–13499. [85] Ellepola, J.; Thijssen, N.; Grievink, J.; Baak, G.; Avhale, A.; van Schijndel, J. Development of a synthesis tool for Gas-To-Liquid complexes. Comp. Chem. Eng. 2012, 42, 2–14. [86] Yue, D.; Kim, M. A.; You, F. Design of Sustainable Product Systems and Supply Chains with Life Cycle Optimization Based on Functional Unit: General Modelling Framework, Mixed-Integer Nonlinear Programming Algorithms and Case Study on Hydrocarbon Biofuels. ACS Sustainable Chem. Eng. 2013, 1, 1003–1014. [87] Wang, B.; Gebreslassie, B. H.; You, F. Sustainable design and synthesis of hydrocarbon biorefinery via gasification pathway: Integrated life cycle assessment and technoeconomic analysis with multiobjective superstructure optimization. Comp. Chem. Eng. 2013, 52, 55– 76. [88] Gebreslassie, B. H.; Waymire, R.; You, F. Sustainable Design and Synthesis of Algae-Based Biorefinery for Simultaneous Hydrocarbon Biofuel Production and Carbon Sequestration. AIChE Journal 2013, 59, 1599–1621. [89] Gebreslassie, B. H.; Slivinsky, M.; Wang, B.; You, F. Life cycle optimization for sustainable design and operations of hydrocarbon biorefinery via fast pyrolysis, hydrotreating and hydrocracking. Comp. Chem. Eng. 2013, 50, 71–91. [90] Floudas, C. A. Deterministic Global Optimization: Theory, Methods and Applications; Kluwer Academic Publishers: Dordrecht, Netherlands, 2000. 76

ACS Paragon Plus Environment

Page 76 of 79

Page 77 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[91] Tawarmalani, M.; Sahinidis, N. V. Convexification and Global Optimization in Continuous and Mixed-Integer Nonlinear Programming: Theory, Applications, Software, and Applications; Kluwer Academic Publishers: Norwell, MA, 2002. [92] Floudas, C. A.; Pardalos, P. M. State of the Art in Global Optimization: Computational Methods and Applications. J. Global Optim. 1995, 7, 113. [93] Floudas, C. A.; Gounaris, C. E. A review of recent advances in global optimization. J. Global Optim. 2009, 45, 3–38. [94] Floudas, C. A.; Akrotirianakis, I. G.; Caratzoulas, S.; Meyer, C. A.; Kallrath, J. Global optimization in the 21st century: Advances and challenges. Comp. Chem. Eng. 2005, 29, 1185–1202. [95] Gounaris, C. E.; Misener, R.; Floudas, C. A. Computational Comparison of PiecewiseLinear Relaxations for Pooling Problems. Ind. Eng. Chem. Res. 2009, 48, 5742–5766. [96] Misener, R.; Gounaris, C. E.; Floudas, C. A. Global Optimization of Gas Lifting Operations: A Comparative Study of Piecewise Linear Formulations. Ind. Eng. Chem. Res. 2009, 48, 6098–6104. [97] Misener, R.; Floudas, C. A. Global Optimization of Large-Scale Generalized Pooling Problems: Quadratically Constrained MINLP Models. Ind. Eng. Chem. Res. 2010, 49, 5424– 5438. [98] Misener, C. E., R. Gounaris; Floudas, C. A. Mathematical modeling and global optimization of large-scale extended pooling problems with the (EPA) complex emissions constraints. Comp. Chem. Eng. 2010, 34, 1432–1456. [99] Misener, R.; Floudas, C. A. Advances for the Pooling Problem: Modeling, Global Optimization, and Computational Studies. Appl. Comput. Math. 2009, 8, 3–22. [100] Misener, R.; Thompson, J. P.; Floudas, C. A. APOGEE: Global optimization of standard, generalized, and extended pooling problems via linear and logarithmic partitioning schemes. Operations Research 2011, 35, 876–892. 77

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[101] Meyer, C. A.; Floudas, C. A. Global Optimization of a Combinatorially Complex Generalized Pooling Problem. AIChE Journal 2006, 52, 1027–1037. [102] Duran, M. A.; Grossmann, I. E. Simultaneous Optimization and Heat Integration of Chemical Processes. AIChE Journal 1986, 32, 123–138. [103] Karuppiah, R.; Grossmann, I. E. Global Optimization for the Synthesis of Integrated Water Systems in Chemical Processes. Computers and Chemical Engineering 2006, 30, 650–673. [104] Ahmetovic, E.; Grossmann, I. E. Optimization of Energy and Water Consumption in CornBased Ethanol Plants. Ind. Eng. Chem. Res. 2010, 49, 7972–7982. [105] Grossmann, I. E.; Mart´ın, M. Energy and Water Optimization in Biofuel Plants. Chinese Journal of Chemical Engineering 2010, 18, 914–922. [106] Ahmetovic, E.; Grossmann, I. E. Global superstructure optimization for the design of integrated process water networks. AIChE J. 2010, 57, 434–457. [107] Energy Information Administration, Annual Energy Outlook 2011 with Projections to 2035. Document Number: DOE/EIA-0383(2011), http://www.eta.doe.gov/oiaf/aeo/, 2011. [108] Energy Information Administration, Refinery Capacity Report, http://www.eia.gov/ petroleum/refinerycapacity/, 2012. [109] World Port Source, Ports in United States, http://www.worldportsource.com, 2010. [110] Energy Information Administration,

Company Level Imports Historical (2008),

http://www.eia.doe.gov/oil\_gas/petroleum/data\_publications/company\ _level\_imports/cli\_historical.html, 2008. [111] Distances.com, World Ports Distances Calculator, http://www.distances.com, 2010. [112] Searcy, E.; Flynn, P.; Ghafoori, E.; Kumar, A. The Relative Cost of Biomass Energy Transport. Applied Biochemistry and Biotechnology 2007, 136-140, 639–652.

78

ACS Paragon Plus Environment

Page 78 of 79

Page 79 of 79

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

[113] Energy Information Administration, Annual Energy Outlook 2010, http://www.eia.doe. gov/oiaf/aeo/, 2010. [114] National Renewable Energy Laboratory, National Renewable Energy Laboratory. National Solar Radiation Database 1991-2005 Update: User’s Manual. Contract No. DE-AC36-99GO10337. [115] National Renewable Energy Laboratory, Wind Integration Datasets, 2010. [116] Kenny, J. F.; Barber, N. L.; Hutson, S. S.; Linsey, K. S.; Lovelace, J. K.; Maupin, M. A. Estimated Use of Water in the United States in 2005, U.S. Geological Survey, 2009. [117] National Energy Technology Laboratory, NETL 2010 Carbon Sequestration Atlas of the United States and Canada (Third Edition, Atlas III), 2010. [118] Ogden, J. M. Conceptial Design of Optimized Fossil Energy Systems with Capture and Sequestration of Carbon Dioxide, Institute of Transportation Studies, University of California, Davis, Research Report UCD-ITS-RR-04-34, 2004. [119] Argonne National Laboratory, GREET 1.8b, The Greenhouse Gases, Regulated Emisssions, and Energy Use in Transportation (GREET) Model, , 2007, release September 2008. [120] National Energy Technology Laboratory, Cost and Performance Baseline for Fossil Energy Plants. Volume 1: Bituminous Coal and Natural Gas to Electricity Final Report. Document Number: DOE/NETL-2007/1281, http://www.netl.doe.gov/energy-analyses/ baseline_studies.html, 2007. [121] Energy Information Administration, Regional Definitions Map, http://www.eia.doe. gov/pub/oil_gas/natural_gas/analysis_publications/ngpipeline/regional_ def.html, 2010.

79

ACS Paragon Plus Environment