CHEMICAL ENGINEERING TRANSACTIONS VOL. 56, 2017 A publication of The Italian Association of Chemical Engineering Online at www.aidic.it/cet Guest Editors: Jiří Jaromír Klemeš, Peng Yen Liew, Wai Shin Ho, Jeng Shiun Lim Copyright © 2017, AIDIC Servizi S.r.l., ISBN 978-88-95608-47-1; ISSN 2283-9216 Selection of CO2 Utilization Options in Carbon Capture, Utilization and Storage (CCUS) Systems using Analytic Hierarchy Process- Data Envelopment Analysis (AHP-DEA) Approach John Frederick Tapiaa, Michael Angelo Promentillaa, Ming Lang Tsengb, Raymond Tan*,a aChemical Engineering Department, De La Salle University-Manila, Philippines bDepartment of Business Administration, Lunghwa University of Science and Technology, Taiwan raymond.tan@dlsu.edu.ph Carbon capture, utilization and storage (CCUS) is one of the most important technologies for reducing greenhouse gas emissions into the atmosphere. Carbon Dioxide (CO2) utilization enables the use of CO2 emissions as input for processes to gain additional revenue. Options for CO2 utilization include CO2-enhanced oil recovery (CO2 - EOR) and CO2-enhanced coal methane (CO2 - ECBM) recovery. These techniques involve injection of CO2 into a geological reservoir enabling the increased recovery of oil (CO2 - EOR) and gas (CO2 - ECBM and CO2 - EGR) and storing CO2 emissions into the ground (geological sequestration) simultaneously. Integrating these CO2 utilization operations into a large-scale CCUS system requires selection of oil and gas reservoirs to develop an efficient CCUS infrastructure. In this study, a site screening framework based on the analytic hierarchy process (AHP) and data envelopment analysis (DEA) approaches is developed to select reservoirs for CO2 utilization operations. AHP-based approach is used to aggregate evaluation of qualitative data (reservoir’s structural integrity and injection well security) to be integrated into DEA approach in determining site efficiencies. A case study is presented to illustrate the framework. 1. Introduction Majority of the worldwide CO2 emissions amounting to 32.3 Gt comes from combustion of fossil fuels (IEA, 2015). Based on a business-as-usual case, it is expected that these emissions will increase by 37 % from its current level by 2040 (IEA, 2012). One of the technologies that contribute to the decrease in greenhouse gas (GHG) emissions is carbon capture and storage (CCS). It involves capturing CO2 from flue gas and injecting it to storage reservoirs (Pires et al., 2011). CCS enables the reduction of CO2 emission into the atmosphere and the use of fossil fuel for energy production. Despite its potential for mitigating the effects of climate change, it is expected to incur substantial additional cost for capture, transportation and injection of CO2 (De Coninck, 2014). Thus, options for productive use of CO2 are needed to reduce the cost associated with CO2 emissions reduction. Such options include enhanced oil recovery (EOR), enhanced coal bed methane (ECBM) and enhanced gas recovery (EGR) from shale gas (Wei et al., 2015).These options enable the recovery of valuable product (e.g. oil and natural gas) by injecting CO2 into the reservoirs. Moreover, some of injected CO2 is stored into the reservoir. CCS systems can be integrated with these options to form a large carbon capture, utilization and storage (CCUS) system. However, properly selecting and screening of candidate options for CO2 utilization is important to avoid investing into inefficient and expensive operations and to maximize for the available CO2 source. Mathematical models have been proposed for designing systems of CO2 sources and storage options. Such models include a continuous-time (Tan et al., 2012) which is later improved (Lee and Chen, 2012) and discrete- time (Tan et al., 2013) scheduling approach to match CO2 sources and sinks with different times of availability. A fuzzy mixed integer linear program (MILP) has been proposed for balancing risks due to uncertainties in DOI: 10.3303/CET1756079 Please cite this article as: Tapia J.F., Promentilla M.A., Tseng M.L., Tan R.R., 2017, Selection of co2 utilization options in carbon capture, utilization and storage (ccus) systems using analytic hierarchy process-data envelopment analysis (ahp-dea) approach, Chemical Engineering Transactions, 56, 469-474 DOI:10.3303/CET1756079 469 storage reservoir capacities and injectivity limits (Tapia and Tan, 2014). On the other hand, adjusting a CCS system due to sudden availability of new information has been developed using a two - step optimization approach (Tapia and Tan, 2015). A unified CCS source - sink model is also proposed to address both temporal issue and power generation make-up in CCS (Lee et al., 2014). For EOR, several approaches have been developed to design EOR systems. Such studies include allocation and scheduling of EOR operations as a strip packing problem (Tapia et al., 2015). On the other hand, Mohd Nawi et al. (2015) proposed a pinch analysis approach for CCUS systems. These studies, however, requires that the sources and sinks are already qualified for the design. Selection of sources and sinks has been only addressed for CCS operations (Promentilla et al., 2013) and has been limited to geological sequestration options only. In this study, an analytic hierarchy process - data envelopment analysis-based (AHP - DEA) framework is developed for screening of candidate options for CO2 utilization. The presence of both qualitative expert judgment and quantitative assessment of these CO2 utilization options needs both AHP and DEA to properly assess which options are to be selected. In this study, the strength of DEA-based evaluation is maximized in quantitative data while AHP - based technique is utilized for qualitative data. AHP - based pairwise comparison is used to make qualitative expert judgment into an output data for DEA - based approach. On the other hand, DEA - based approach determines the efficiencies of all options and screens out those which are inefficient. The rest of the paper is organized as follows: Section 2 presents the problem statement while Section 3 elaborates the AHP - DEA framework. Section 4 presents a case study for illustration and lastly, Section 5 gives the conclusions and future works. 2. Problem Statement In addressing the screening of option for CO2 utilization, the formal problem statement was given below. This defines the system that needs to be solved by the AHP - DEA framework:  The system consists of n different options for CO2 utilization.  Each options are evaluated based on the following criteria: o Distance from CO2 source - in this criterion, each reservoir associated with these options is evaluated based on the nearest CO2 source. o Minimum Flow Rate Requirement - the minimum flow rate to start CO2 flooding (in Mt CO2 / y). o Injectivity Limit - the maximum flow rate to maintain the structural integrity of the reservoir (in Mt CO2 / y). o Operating Life - length of operation (in y) for a specific option o Product Yield - amount of product yield per CO2 injected (Mtoe / Mt CO2). o Product Value - price of commodity (e.g. oil, gas etc.) recovered (M$/ Mtoe). o Sequestration Parameter- amount of CO2 stored per CO2 injected. o Reservoir Capacity - total CO2 that can be stored to the reservoir. o Well Security - refers to the security of CO2 from escaping in CO2 wells. o Structural Integrity - refers to the security of CO2 from escaping from the geological formation. This also includes risk of CO2 from escaping to a nearby groundwater source etc.  For well security and structural integrity, an AHP - based pairwise comparison approach will be used to determine the weights of each option for these criteria. The method is best applied to measure expert judgement by pairwise comparison.  The efficiency of each CO2 utilization options will be evaluated using the Charnes - Cooper - Rhodes (CCR) model for DEA. The input criteria that will be considered are those preferred when the scores is higher (i.e. minimum pipeline distance and flow rate requirement). On the other hand, the rest are considered the output criteria.  The qualified options for CO2 utilization are those with efficiency equal to 1. 3. AHP - DEA Framework The AHP - DEA framework that will be used for this study is illustrated in Figure 1. The quantitative data should be evaluated based on measuring tools and seismic survey. On the other hand, the qualitative criteria presented are evaluated based on a pairwise comparison approach to quantify the judgments made by experts. Quantifying expert judgment is done using a 9 - point scale proposed by Saaty (2003).The pairwise comparison matrix is aggregated using the eigenvector method to determine the weights of each alternative for both well security and structural integrity criteria. When all data were expressed quantitatively, the CO2 utilization options were evaluated using CCR model used for DEA. This model calculates the Pareto efficiency of each decision making unit (DMU) based on an aggregated input and aggregated output. The objective of the CCR model is to maximize the efficiency of one 470 option subject to having an aggregated input of 1 and all other options have efficiency equal to 1 or less. The model is solved based on the number of option to be evaluated and the efficiency is calculated. Note that more than one option can have efficiency equal to one. Therefore, in this study, options with efficiency equal to 1 will be selected for CCUS systems. Figure 1: AHP - DEA Framework for Selection of CO2 Utilization Options (Methodology based from Veni et al., 2012) 4. Case Study To illustrate the methodology used, a hypothetical case study is made which consists of 9 CO2 utilization options. These are composed of four EOR operations with depleted oil reservoir as storage medium, three ECBM operations with coal beds as storage medium and two EGR with shale reservoir as storage medium. The output data is shown in Table 1 while the input data is shown in Table 2. The pairwise comparison matrices for both well security and structural integrity are shown in Tables 3 and 4. These were results from an expert judgment in which two options are evaluated at a time. For instance, EOR1 is 7 times better than that of EOR1 in terms of well security, thus a score of 7 is placed in the pairwise comparison matrix at row 1, column 3. The reverse is calculated by obtaining the reciprocal and placing it in row 3, column 1. This is done for both well security and structural integrity. 471 Table 1: Output Data for DEA Approach DMU Operating Life (y) Reservoir CO2 Capacity (Mt) Injectivity Limit (Mt/y) Product Yield (Mtoe/Mt CO2) Product Value (M$/Mtoe) Sequestration Parameter EOR1 25 100 10 4.5 580 0.35 EOR2 15 200 12 7.5 580 0.45 EOR3 20 150 14 8.9 580 0.67 EOR4 15 150 13 2.3 580 0.50 ECBM1 14 50 5.6 4.5 120 0.89 ECBM2 10 300 6.7 3.5 120 0.88 ECBM3 20 35 7.8 8.6 120 0.76 ShaleGas1 10 65 8.0 2.3 120 0.67 ShaleGas2 8 35 7.6 1.5 120 0.75 Table 2: Input Data for DEA Approach DMU Distance from CO2 Source (km) Flow Rate Required (Mt/y) EOR1 150 3.4 EOR2 120 5.8 EOR3 110 3.4 EOR4 90 7.6 ECBM1 100 5.6 ECBM2 24 2.4 ECBM3 15 4.4 ShaleGas1 50 2.3 ShaleGas2 45 7.8 Table 3: Pairwise Comparison Matrix for Well Security Criteria DMU EOR1 EOR2 EOR3 EOR4 ECBM1 ECBM2 ECBM3 Shale Gas1 Shale Gas2 EOR1 1.000 0.333 7.000 0.111 0.143 0.167 1.000 0.111 0.333 EOR2 3.000 1.000 9.000 0.556 0.333 1.000 4.000 0.500 2.000 EOR3 0.143 0.111 1.000 0.111 0.111 0.143 0.167 0.111 0.111 EOR4 9.000 1.800 9.000 1.000 2.000 2.000 8.000 1.000 3.000 ECBM1 7.000 3.000 9.000 0.500 1.000 2.000 5.000 1.000 3.000 ECBM2 6.000 1.000 7.000 0.500 0.500 1.000 4.000 3.000 2.000 ECBM3 1.000 0.250 6.000 0.125 0.200 0.250 1.000 0.167 0.333 ShaleGas1 9.000 2.000 9.000 1.000 1.000 0.333 6.000 1.000 3.000 ShaleGas2 3.000 0.500 9.000 0.333 0.333 0.500 3.000 0.333 1.000 Table 4: Pairwise Comparison Matrix for Structural Integrity Criteria DMU EOR1 EOR2 EOR3 EOR4 ECBM1 ECBM2 ECBM3 Shale Gas1 Shale Gas2 EOR1 1.000 0.167 0.111 0.500 0.111 1.000 0.111 0.500 0.111 EOR2 6.000 1.000 2.000 3.000 1.000 5.000 0.333 3.000 0.333 EOR3 9.000 0.500 1.000 5.000 1.000 7.000 0.500 5.000 0.500 EOR4 2.000 0.333 0.200 1.000 0.200 2.000 0.111 1.000 0.111 ECBM1 9.000 1.000 1.000 5.000 1.000 6.000 0.333 4.500 0.500 ECBM2 1.000 0.200 0.143 0.500 0.167 1.000 0.111 7.000 0.111 ECBM3 9.000 3.000 2.000 9.000 3.000 9.000 1.000 7.000 1.000 ShaleGas1 2.000 0.333 0.200 1.000 0.222 0.143 0.143 1.000 0.111 ShaleGas2 9.000 3.000 2.000 9.000 2.000 9.000 1.000 9.000 1.000 472 The weights of each option for both criteria are evaluated using the eigenvector method. Table 5 shows the scores for each option in terms of both well security and structural integrity. These scores are then treated as output for the CCR model which in turn solves the efficiency of each option. Table 6 shows the solution of the CCR model. Note that each calculation for the efficiency has negligible computational time. Table 5: Eigenvector Values for Well Security and Structural Integrity Criteria DMU Well Security Structural Integrity EOR1 0.03234 0.01948 EOR2 0.10582 0.11839 EOR3 0.01337 0.13312 EOR4 0.22269 0.03272 ECBM1 0.18811 0.12956 ECBM2 0.15877 0.03854 ECBM3 0.03381 0.25467 ShaleGas1 0.17171 0.02658 ShaleGas2 0.07339 0.24694 Table 6: Efficiency Calculations based from CCR Model DMU Efficiency Efficient? EOR1 1.000 YES EOR2 0.872 NO EOR3 1.000 YES EOR4 1.000 YES ECBM1 0.734 NO ECBM2 1.000 YES ECBM3 1.000 YES ShaleGas1 1.000 YES ShaleGas2 0.604 NO Based from the calculations, EOR2, ECBM1 and ShaleGas2 are screened out since they have efficiencies less than one. On the other hand, options with efficiency equal to one are selected for next design stage based on the criteria given. The options with efficiency equal to one dominate the three options with efficiency less than one. For this example, an efficient CCUS system can be made using three options of EOR, two options and ECBM and one EGR operation. More detailed CCUS design can then be made from these options. 5. Conclusions An AHP - DEA approach was developed for application in selecting CO2 utilization options in CCUS systems. This makes use of the pairwise comparison from AHP to measure qualitative judgement such as well security and structural integrity. The basis for selecting utilizations options is the DEA efficiencies calculated from CCR model. The strength of DEA - based approach in this study is to select multiple options from the candidate alternatives without determining which criterion is better than the other. On the other hand, AHP - based pairwise comparison enables the conversion of expert qualitative judgment into DEA-usable quantitative data. Future work includes extension of the method to address uncertainties in both expert judgment and available data. Also, the method will be extended to include CO2 sources as one of the factors that influences the selection of CO2 utilization and storage options, Acknowledgments J.F.D. Tapia would like to acknowledge the Philippine Department of Science and Technology (DOST) through its Engineering Research and Development for Technology (ERDT) Program and the Commission on Higher Education- Philippine Higher Education Research Network (CHED-PHERNet) for financial support. Reference De Coninck, H., Benson, S.M., 2014, Carbon Dioxide Capture and Storage: Issues and Prospects, Annu. Rev. Environ. Resour. 39, 243–270. International Energy Agency (IEA), 2012, World energy outlook, International Energy Agency, Paris, France. 473 International Energy Agency (IEA), 2015, World energy investment outlook, International Energy Agency, Paris, France. Veni K.K., Pugazhendhi S., Pugazhendhi S., 2012, Development of Decision Making Model Using Integrated AHP and DEA for Vendor Selection, Procedia Eng. 38, 3700–3708. Lee J.-Y., Chen C.-L., 2012, Comments on “continuous-time optimization model for source-sink matching in carbon capture and storage systems”, Ind. Eng. Chem. Res. 51, 11590–11591. Lee J.-Y., Tan R. R., Chen C.-L., 2014, A unified model for the deployment of carbon capture and storage, Appl. Ener. 121, 140-148. Mohd Nawi W.N.R., Wan Alwi S.R., Manan Z.A., Klemeš J.J., 2015, A new algebraic pinch analysis tool for optimising CO2 capture, utilisation and storage, Chemical Engineering Transactions 45, 265–270. Pires J.C.M., Martins FG., Alvim-Ferraz M.C.M., Simões M., 2011, Recent developments on carbon capture and storage: An overview, Chem. Eng. Res. Des. 89, 1446–1460. Promentilla M.A.B., Tapia J.F.D., Arcilla C.A., Dugos N.P., Gaspillo P.D., Roces S.A., Tan R.R., 2013, Interdependent ranking of sources and sinks in CCS systems using the analytic network process, Environ. Modell. Softw. 50, 21-24. Saaty T.L., 2003, Decision-making with the AHP: Why is the principal eigenvector necessary, Eur. J. Oper. Res. 145(1), 85-91. Tan R.R., Aviso K.B., Bandyopadhyay S., Ng D.K.S., 2012, Continuous-Time Optimization Model for Source– Sink Matching in Carbon Capture and Storage Systems, Ind. Eng. Chem. Res. 51, 10015–10020. Tan R.R., Aviso K.B., Bandyopadhyay S., Ng D.K.S., 2013, Optimal source-sink matching in carbon capture and storage systems with time, injection rate, and capacity constraints, Environ. Prog. Sustain. Energy 32, 411–416. Tapia J.F.D., Lee J.-Y., Ooi R.E.H., Foo D.C.Y., Tan R.R., 2015, Design and scheduling of CO2 enhanced oil recovery with geological sequestration operations as a strip packing problem, Chemical Engineering Transactions 45, 1615-1620. Tapia J.F.D., Tan R.R., 2014, Fuzzy optimization of multi-period carbon capture and storage systems with parametric uncertainties, Process Saf. Environ. Prot. 92, 545–554. Tapia J.F.D., Tan R.R., 2015, Optimal revamp of multi-region carbon capture and storage (CCS) systems by two-step linear optimization, Energy Syst. 6, 269-289. Wei N., Fang Z., Bai B., Li Q., Liu S., Jia Y., 2015, Regional resource distribution of onshore carbon geological utilization in China, J. CO2 Util. 11, 20–30. 474