Multi-Objective Optimization of CO 2 Sequestration in Heterogeneous Saline Aquifers under Geological Uncertainty

: This paper presents a Pareto-based multi-objective optimization for operating CO 2 sequestration with a multi-well system under geological uncertainty; the optimal well allocation, i.e., the optimal allocation of CO 2 rates at injection wells, is obtained when there is minimum operation pressure as well as maximum sequestration efﬁciency. The distance-based generalized sensitivity analysis evaluates the inﬂuence of geological uncertainty on the amount of CO 2 sequestration through four injection wells at 3D heterogeneous saline aquifers. The spatial properties signiﬁcantly inﬂuencing the trapping volume, in descending order of inﬂuence, are mean sandstone porosity, mean sandstone permeability, shale volume ratio, and the Dykstra–Parsons coefﬁcient of permeability. This conﬁrms the importance of storable capacity and heterogeneity in quantitatively analyzing the trapping mechanisms. Multi-objective optimization involves the use of two aquifer models relevant to heterogeneity; one is highly heterogeneous and the other is less so. The optimal well allocations converge to non-dominated solutions and result in a large injection through one speciﬁc well, which generates the wide spread of a highly mobile CO 2 plume. As the aquifer becomes heterogeneous with a large shale volume and a high Dykstra–Parsons coefﬁcient, the trapping performances of the combined structural and residual sequestration plateau relatively early. The results discuss the effects of spatial heterogeneity on achieving CO 2 geological storage, and they provide an operation strategy including multi-objective optimization. porosity, merical beginning of then horizontally. properties near step analyses for CO heterogeneous, able to reduce the plume mobility, while unable to maintain the low injection pressure. The results demonstrate the importance of considering geolog- ical uncertainty in the performance of CO 2 sequestration in a more realistic manner. Asymmetric parameter interactions with spatiotemporal variables remain a chal- lenge for the analysis of quantitatively supercritical CO 2 flow under geological uncer- tainty. Conceptualization, and and sequestration, in descending order of inﬂuence: mean sandstone porosity, mean sandstone permeability, shale volume ratio, and the Dykstra–Parsons coefﬁcient. The optimal well allocation suggested a large amount of injection through one well, and it induced a wide distribution of CO 2 plume while maintaining low injection pressures. Numerical simulations reliably explained the transport of the highly mobile CO 2 plume in brine-saturated media; the plume moved upward at the beginning of injection, reached the impermeable cap rock, and then spread horizontally. Obtaining a detailed description of the properties at a vertical section near an injector is an essential step toward achieving quantitative analyses for CO 2 sequestration. As the aquifer became heterogeneous, the impermeable shale was able to reduce the plume mobility, while it was unable to maintain the low injection pressure. The results demonstrate the importance of considering geological uncertainty in analyzing the performance of CO 2 sequestration in a more realistic manner. Asymmetric parameter interactions with spatiotemporal variables remain a challenge for the analysis of quantitatively supercritical CO 2 ﬂow under geological uncertainty.


Introduction
A challenging problem in engineering analytics has been the existence of many objectives. Multi-objective optimization attempts to find the optimal trade-offs that are most acceptable to the decision maker among all objective functions [1][2][3]. The Pareto front, i.e., a set of Pareto solutions, illustrates the trade-offs for which algorithms should secure solution diversity as well as make comparative evaluations among the potential solutions [4][5][6][7][8]. Evolutionary multi-objective optimization (EMO) algorithms, e.g., the non-dominated sorting genetic algorithm (NSGA; NSGA-II; NSGA-III), strength Pareto evolutionary algorithm (SPEA; SPEA-II), Pareto envelope-based selection algorithm (PESA; PESA-II), and multi-objective evolutionary algorithm based on decomposition (MOEA/D), have continuously improved fitness assignment and diversity control [1][2][3][4][5][6][7][8][9][10]. NSGA is a well-known scheme designed to preserve non-dominated points in objective space that also has a wide solution-searching capability with a genetic algorithm. Its strength is that it can provide non-dominated trade-offs in the comparison of objective functions, while its weaknesses are the fact that it requires a large amount of computing power and its decreasing convergence with an increasing number of objective functions, i.e., 'the curse of dimensionality' [5][6][7][8][9][10].
Geological uncertainty limits quantitative analyses of subsurface flow [11][12][13][14]. The input uncertainties in geo-modeling are non-linearly interacted; some are heterogeneous (e.g., porosity, permeability, and lithology) while others are functional (e.g., relative permeability, solubility, and capillary pressure) or scenario-based (e.g., the depositional system) [13][14][15][16][17]. The flow responses are spatiotemporal. Sensitivity measures the relationship between the inputs and the responses and, thus, denotes some influential parameters that significantly affect the responses. The stochastic features between the inputs and the responses do not result in the emergence of a certain parameter as a significant factor but instead require multi-way parameter interactions. The distance-based generalized sensitivity analysis (DGSA) has been successfully used to evaluate the significance of scale-variant properties in heterogeneous geo-models [10,[15][16][17][18][19]. Fenwick et al. [15] quantified the interaction of asymmetric parameters such as residual oil saturation, maximum water relative permeability, and training images. Park et al. [16] analyzed spatial uncertainty using kernel principal components and self-organizing maps. CO 2 sequestration into deep saline aquifers, e.g., around 1 km below ground level, is known to have enormous storage potential, but little is known about its geological characterization and project experiences compared to those of depleted oil or gas reservoirs. The significant uncertainties can be categorized into capacity, injectivity, and containment: the capacity determines the storable amount, the injectivity estimates the possible injection rates, and the containment evaluates the risks of leakage [20][21][22][23][24]. Bachu [20] explained the importance of the hydrodynamic and buoyancy forces needed for the CO 2 plume to propagate into the homogeneous saline formation. Kumar et al. [21] discussed the key factors related to the rock and fluid properties, reservoir conditions, and injection strategy as storage mechanisms (structural, residual, solubility, and mineral trapping) in saline aquifers; these include reservoir heterogeneity, depth, permeability, pressure, and temperature. Along with these geological properties, the injection pressure is one of the crucial operating parameters depending on the containment and the security. One reason for this is the well injectivity-defined as the CO 2 injection mass rate per unit pressure differential-which determines the total CO 2 amount injected into the geological formation. On the other hand, the safety constraint is that the injection pressure (the bottom hole pressure; BHP) should be lower than the formation fracture pressure so as to prevent CO 2 leakage or fast movement without trapping [25][26][27][28][29][30]. When CO 2 is injected into a porous medium using multi-injectors, the well allocation, i.e., the allocation of CO 2 amount at each injector, worsens the problem; the well allocation at the heterogeneous aquifers would involve a lot of scenarios, so the optimization process could be subject to the issue of convergence [31][32][33][34][35][36][37][38][39][40][41][42].
This study intends to search for trade-offs between cost and containment efficiency as a multi-objective problem for CO 2 sequestration into heterogeneous saline aquifers; the cost is related to the CO 2 injection pressure used to operate the facilities, while the containment efficiency involves storing as much CO 2 as possible while considering the amount of CO 2 injected. Well allocation constructs variable spaces, while two objective functions (the injection pressure and the sequestration efficiency) constitute the posterior domain. The uncertainty assessment of spatial properties should be made prior to multiobjective optimization because the subsurface aquifers are heterogeneous in nature and their properties are uncertain.

3D Heterogeneous Aquifer Models and Simulation Conditions
A total of 1024 heterogeneous aquifers are constructed geostatistically under spatial uncertainty of permeability and porosity. The spatial distribution of these properties is modeled using Petrel (Schlumberger, Houston, TX, USA). Three dimensional unstructured grid systems are generated to demonstrate both the upward CO 2 movement and the horizontal flow. Total number of unstructured grids is 22,230 (= 38 × 45 × 13). The saline aquifer size is (x, y, z) = (6310 m, 7076 m, 250 m). It is assumed that the impermeable cap rock is located at 840 m true vertical depth. Dykstra-Parsons coefficients (v DP ) measure permeability heterogeneity between 0 and 1 (Equation (1)).
In Equation (1), k 0.5 is the mean value and k 0.841 is the value placed at one standard deviation plus the mean [43,44]. v DP = 0 means that the medium is homogeneous. The higher the Dykstra-Parsons coefficient is, the more heterogeneous the medium is. The ratio of vertical to horizontal permeability is 0.1 and the permeabilities are isotropic horizontally (x and y directions). Table 1 summarizes the input ranges of the spatial properties for constructing 3D heterogeneous aquifers. The stochastic realizations, i.e., geostatistical property modeling, are carried out with the inputs randomly selected ( Table 1). The heterogeneous geo-models are made with massive sandstone intercalated with almost impermeable shale ( Figure  1). Shale is placed at the bottom of aquifer and interbedded in sandstone; the former is to prevent leakage at the early time and the latter is to demonstrate the realistic heterogeneous aquifers. The shale volume ratio is defined as the ratio of shale volume to the rock volume and represents the effects of shale layers; this is selected between 2 and 20%. The shale porosities (ranging from 0.1 and 0.15) are smaller than those of sandstone. The shale permeabilities (ranging from 0.05 and 0.5 millidarcy) represent almost impermeable barriers The standard deviations of shale permeability and porosity are negligible compared to sandstone. The end point of relative permeability at the irreducible water saturation (=0.1) is 0.85, while that at the irreducible CO 2 saturation (=0.2) sets as 0.2.  Figure 1 illustrates an example of heterogeneous aquifer models with four injectors (namely, I1, I2, I3 and I4). Figure 1a depicts the spatial distribution of permeability and Figure 1b describes that of porosity. The injection wells are installed to evenly divide the entire aquifer; (x, y) grid positions are as follows, I1 = (10, 12), I2 = (28, 12), I3 = (10, 34), and I4 = (28,34). The perforation zone in z direction is from 1032.3 m to 1051.53 m, the lower part of the aquifer. The wells are partially perforated, over the less permeable basement rock (the shale layer), to demonstrate the upward movement of CO 2 plume. MRST (MATLAB Reservoir Simulation Toolbox) is used to simulate CO 2 sequestration in saline aquifers [44,45]. It provides various trapping volumes: a volume of combined structural and residual trapping, a volume of residual trapping, a volume of residual trapping with plume, a volume of structural trapping, and a volume of movable plume, respectively. CO 2 properties are assessed using CO2 Lab module in MRST [44]. The reference pore-pressure is 120 bar and the temperature is 45 degrees Celsius at the top of the aquifer (840 m depth) where the formation fracture pressure is 150 bar. CO 2 viscosity is set as 0.051 centipoise (1 centipoise = 10 −3 Pa·s) and its density is 675.324 kg/m 3 at the given reference pressure and temperature. Brine density is 975.86 kg/m 3 .

DGSA for Evaluating the Significance of Spatial Properties
DGSA, based on regionalized sensitivity analysis, defines the cumulative density functions (CDF) distance [15,16]. The distance reflects variational properties of the responses or the dimensionality-reduced responses [16]. After classifying of the output response with a clustering scheme, the class can generate CDF and the distance can be defined as L1 norm between two CDFs (Equation (2)).
In Equation (2),d c,i is the CDF distance (L1 norm distance, d L1 ) between the input CDF for parameter p i and the CDF of the parameter in the class, c.Ĝ(p i ) is an empirical distribution function from sampling.Ĝ(p i |c) is a class-conditioned empirical distribution function after classification. C is the number of class that the modeler pre-assigned and m is the number of parameters. The measure of sensitivity is standardized by the resam-pling quantile of the distance. DGSA introduced multidimensional scaling, characterized the population and the clusters using CDF, and tested the asymmetric interactions of spatiotemporal parameters using the distance between CDFs (L1 norm) [15][16][17][18][19].

NSGA-II for Multi-Objective Optimization Calibrating Well Allocations
Pareto optimality, i.e., the trade-off solutions, is the non-domination where no solution can be improved with respect to any objective without worsening at least another objective [4,5,8]. Multi-objective evolutionary algorithm searches the trade-offs that would give a vector of decision variables for all the objective functions under any given constraints (Equation (3)).
In Equation (3), f (X) is an objective function with the vector (X) of decision variables, e.g., well allocation in this study, and the subscript n is the number of objective functions. NSGA-II algorithm consists of competing individuals with a genetic algorithm, non-dominated sorting, and crowding distance [5,[7][8][9][10]. The genetic algorithm-based operations, e.g., crossover and mutation, form the offspring and then the individuals in the mixture of parent and offspring population are competing with non-dominated sorting and the crowding distance. The non-dominated sorting strategy of NSGA-II proceeds by comparing each individual one by one, and when a non-dominated front is found, it is temporarily stored, and then the next non-dominated front is searched. This non-dominated sorting continues until all individuals are ranked [1,5] (Figure 2). The solutions can evolve to the direction of the real Pareto optimization set, as a larger probability is assigned to the individuals in the lowest dominant level in order to select the superior individuals into the next generation of the population. To maintain the solution diversity, the crowding distance calculates the congestion degree of each individual and other individuals [7][8][9] (Figure 2). Equation (4) defines the k-th individual's crowding distance (δ k ).
In Equation (4), f i (k + 1) and f i (k − 1) indicate i-th objective function values for (k + 1)-th and (k − 1)-th individual, respectively. f max i is the maximum value, and f min i is the minimum for i-th objective function. By selecting the individuals with lower crowding distance value in the same dominant level, the evenness and diversity of the evolutional population can be guaranteed [9,10]. Pareto optimality can be obtained by reaching convergence conditions, or by repeating the process up to the maximum number of evolutions, and, thus, multi-objective optimization is costly because of enormous computational simulations.

Results and Discussion
3.1. Spatial Properties Influencing CO 2 Trapping DGSA parameters are assigned as follows: k-medoid clustering, with six clusters, classifies 1024 aquifer models and the resampling quantile of distance is 0.95. A response for DGSA is the trapping volume of combined structural and residual trapping, representing relatively permanent and safe sequestration. The trapping amount is a cumulative volume during 200 years from the first injection. Six uncertain parameters (C = 6 in Equation (2)) are mean sandstone porosity, mean sandstone permeability, standard deviation of sandstone porosity, standard deviation of sandstone permeability, shale volume ratio, and Dykstra-Parsons coefficient, respectively (m = 6 in Equation (2); Table 1). To remove the effects of well allocation, a constant injection rate, 4000 m 3 /day for each injector is assumed. Total injection rate per day is 16,000 m 3 /day because four injectors are placed (Figure 1a). CO 2 is injected continuously for 30 years and the trapping trend is monitored for 200 years, i.e., an additional 170 years from the end of CO 2 injection. Figure 3 depicts the result of DGSA with spatial parameters; the sensitive parameters are the mean porosity of sandstone (PoroSand), mean permeability of sandstone (PermSand), shale volume ratio (SVR), and Dykstra-Parsons coefficient (VDP) in each row. The vertical line (the standardized sensitivity = 1; the significant level) indicates whether the parameter influences the response. The larger standardized sensitivity means more influence. The sensitive parameters support the importance of heterogeneity and aquifer properties on CO 2 trapping: PoroSand determines the capacity size, PermSand affects CO 2 mobility, SVR and VDP represent the effects of the shale barrier on storage and transport, respectively. The trapping amounts significantly depend on the pore volume of sandstone. The other parameters over the significant level are closely related to CO 2 flow. With increasing SVR and VDP, shale is likely to obstruct CO 2 flow. However, this would be a subject of discussion as to whether the large amount of shale always has a positive effect on the trapping volume. The asymmetric parameter interactions can make this debate more complicated as a parameter can simultaneously influence the different responses. If the operating conditions are included, deriving a reasonable conclusion would become a conundrum.  Table 1. PoroSand, PermSand, SVR, VDP, StdPoro, and StdPerm represent the mean porosity of sandstone, the mean permeability of sandstone, shale volume ratio, Dykstra-Parsons coefficient, the standard deviation of sandstone porosity, and the standard deviation of sandstone permeability, respectively. The Pareto bars are colored according to the percentile values. The horizontal black line represents confidence interval within a parameter that is still accepted as influencing. The vertical line means the significant level (if the standardized sensitivity is greater than 1, it means that the parameter is sensitive to the response).

Multi-Objective Optimization with Well Allocations
Multi-objective optimization requires a lot of simulation runs and, thereby, this work selects two different aquifers based on the DGDA result (Table 2; Figure 4); one is the less heterogeneous (L aquifer; less heterogeneous relevant to the small value of Dykstra-Parsons coefficient) and the other is highly heterogeneous (H aquifer; the high heterogeneity). Figure 4 demonstrates spatial distributions of the key properties influencing the trapping amount: Figure 4a,b show permeability and porosity of L aquifer, while Figure 4c,d illustrate those of H aquifer. Upon comparison of the two aquifers, L aquifer is capable of storing more CO 2 and induces slow fluid transport through sandstones (because L aquifer has the smaller sandstone permeability and the larger sandstone porosity; refer to Table 2).  On the other hand, H aquifer is more ambiguous in terms of fluid transport. It has a large portion of shale formations (negative to the continuity of fluid flow) but more permeable sandstones (favorable to a transporting velocity). The former, i.e., shale beds, would be advantageous for long-term CO 2 geological storage, while a drawback would be having to maintain the lower injection pressure since shale can increase the operating cost to inject CO 2 continuously. The more permeable sandstone has double-sidedness as well; the higher transmissivity contributes towards the spread of supercritical CO 2 widely, to occupy additional pore spaces, but can increase the risk of leakage. Therefore, this paper selects H aquifer with overall mean permeability including both sandstone and shale similar to L aquifer model, notwithstanding that H aquifer is more heterogeneous; the entire mean permeability of L aquifer is 227 millidarcy and it is 256 millidarcy for H aquifer, which occupies less permeable shales and reduces the difference in PermSand from 147.2 millidarcy (=448.7-301.5) to 29 millidarcy (=256-227). Shales distribute sparsely, are intercalated with the massive sandstone, and cannot sequestrate CO 2 plumes completely but delay CO 2 movement. In short, the H aquifer model has more shale, higher PermSand, and smaller pore volume than the L aquifer model. The H aquifer model is more heterogeneous with a higher Dykstra-Parsons coefficient.
The multi-objective problem is set to minimize two objectives (Equation (5)): one ( f 1 ) is to minimize the sum of maximum BHPs (that is, injection pressure at the well bottom) assigned at each injector, and the other ( f 2 ) is related to the trapping volume observed after 200 years. The trapping volume (V Trap ) is the same as the response of DGSA, i.e., the volume of combined structural and residual trapping, to evaluate the actual value of safe sequestration. Because the V Trap would be large, f 2 is defined as the minus value of the trapping volume divided by the total injection amount. To make a minimizing problem, the minus sign is added. Injecting more CO 2 requires the higher BHP at the facilities and negatively affects the trapping efficiency, f 2 (X), because the injection amount (V Inj ) is located at the denominator. Thus, the two objective functions are dependent.
In Equation (5), P w means BHPs observed at a well, w (∀w = I1, I2, I3, I4), which has been monitored for 200 years. argmax(P w ) means the maximum P w that should be kept smaller than the formation fracture pressure ( = 150 bar assumed at the reference depth, 840 m). If argmax(P w ) is greater than the fracture pressure, this scenario is excluded at the optimization process. X consists of daily CO 2 injection rate assigned at each injector; x w ∈ X, x w is the CO 2 injection rate (m 3 /day) at w-th well. The parameters in the genetic algorithm are as follows; crossover is 0.7, mutation rate is 0.02, and the maximum iteration is 50. The daily CO 2 injection rate is fixed as 16,000 m 3 /day and four vertical wells (I1, I2, I3, and I4) split this value; the injection rates, less than 4000 m 3 /day for each well, are assigned randomly to three wells (I2, I3, and I4) and the remaining value is set at I1. This study sets the constraint that the rate allocated at the well should be kept until 30 years. The total volume of CO 2 injected over 30 years (V Inj ; total volume of CO 2 injection) is 175.2 million cubic meters (= 16,000 m 3 /day × 365 days × 30 years). Figure 5 depicts the Pareto solutions with the lowest rank and non-domination: Figure 5a presents Pareto solutions for L aquifer and Figure 5b for H aquifer. In terms of f 1 , two cases are chosen for each aquifer: L-1 means the lowest f 1 and L-2 is located at the opposite side. H-1 and H-2 are selected at the boundaries of f 1 for H aquifer in the same way. The distribution of Pareto solutions show that the L aquifer model can store more CO 2 with smaller BHP compared to the highly heterogeneous H aquifer. The maximum trappable volume of L aquifer (the mechanism of combined structural and residual trapping) is 754,304 m 3 , while H aquifer remains at 582,005 m 3 . The sum of maximum BHPs of L aquifer is 505.37 bar, but that of H aquifer is 506.5 bar ( Table 3). The maximum BHP is assigned at I1 well: 132.55 bar for L aquifer and 134.35 bar for H aquifer. To sum up, these results show that a less heterogeneous aquifer with large pore volume would be positive for the trapping and cost-effective operation.   Table 3 summarizes the results of the optimal well allocations, the maximum BHP at an injection well across 200 years, and the trapping volume. The optimal well allocations suggest the large amount of CO 2 injection through I1 well with a high BHP. CO 2 plumes can move from the I1 well area to the others. As the injection progresses, the pore pressure increases from 120 bar (assigned before the sequestrating operation). The small difference between the BHP and the pore pressure is able to obstruct CO 2 injection at the late period. Severe inter-well interference might occur, disturbing the facility operation. Under the constraint of maintaining the allocated rate, this strategy of well allocation reasonably supports the maximization of the permanent trapping but also the minimization of the sum of the operating pressures.
The CO 2 volume injected into the aquifer is divided in the order of movable plume, residual trapping with plume, residual trapping, structural trapping, and the combined structural and residual trapping in a row. The form of movable plume accounts for 60 to 70% of the total injection. The simulation results are in agreement with a typical feature of supercritical CO 2 plume movement; supercritical CO 2 is much more mobile than the formation water and, thereby, moves faster to the aquifer top than brine flows downward [44]. Figure 6 depicts the height of the CO 2 column over time. At 30 years, the column is highest near I1 well where most of the CO 2 is injected. As time passes, many of the mobile CO 2 plumes move upward and spread horizontally when they reach the impermeable cap rock. Eventually, most of the CO 2 is present near the top of the aquifer. The column near I1 well gradually decreases and the distributed area expands below the cap rock. Figures 7 and 8 show the spatial distributions of CO 2 saturation at the cross sections (A and B line in the x-y plane) for L-1 ( Figure 7) and H-1 (Figure 8). CO 2 saturation is observed at 200 years. CO 2 saturation near the injection wells increases similar to the cone shape, but the area where CO 2 is occupied is limited at the upper zones of the aquifer (Figure 7b). The high mobility of supercritical CO 2 generates an upward flow mainly, even though the vertical permeability is 0.1 of the horizontal permeability. The horizontal spreading phenomenon is possible when CO 2 accumulates below the impermeable cap rock. Notwithstanding that this result explains the characteristics of high-mobile CO 2 flow, it requires the detailed description of heterogeneity along a vertical direction.   A notable result related to the heterogeneity is found at the decrement trend of the movable CO 2 plume ( Figure 9); the percentage of movable plume means the volume ratio of movable plume to total injection. It accounts for most of the CO 2 sequestration, e.g., 73% for L aquifer (estimated at 30 years). The plume has possibilities to turn into the trap, to maintain itself, or to leak out of the boundaries. Thus, the large amount of movable plume is negative on sequestration. From 30 years to 200 years, the percentage at L aquifer decreases from 73% to 64%, while that of H aquifer is reduced from 68% to 61%. Thus estimates that the heterogeneity degree (Dykstra-Parsons coefficient) contributes to make these differences; the highly heterogeneous H aquifer has more obstruction for plume movement so that the shale beds would be positive in reducing the leakage risk. Figure 10 plots the trends of the combined structural and residual trapping of both aquifers. The less heterogeneous L aquifer increases the combined trapping until 200 years, while the increment in trapping volume of the highly heterogeneous H aquifer stops at around 80 years. This result recommends that if a decision maker wants to sequestrate a small volume within a relatively short time, the high heterogeneous aquifer would be favorable with a large amount of shale layers in the sandstone. However, an increment in the operating cost would be needed to maintain the high injection pressure. Figure 9. Comparison of movable plume (the ratio of CO 2 volumes that are movable as the plume forms to the total injection amount) between the L and H aquifer models. The total CO 2 amount injected is 1.752 million cubic meter for each aquifer. The difference made by the optimal conditions at the given aquifer model, e.g., L-1 and L-2 or H-1 and H-2, is negligible.

Conclusions
This paper discussed the multi-objective problem related to CO2 sequestration under geological uncertainty. It searched for optimal well allocations to achieve both minimum bottom hole pressures as well as maximum trapping efficiency. The distance-based generalized sensitivity analysis confirmed the influence of the following spatial properties on The results can be used to identify the properties significantly influencing both operation condition and trapping efficiency. Pareto optimality is obtained with non-domination among the solutions, and it provides a well allocation strategy that satisfies the aforementioned multi-objectives. This method is limited in its ability to search for non-dominated solutions; this study does not assume any tolerance or constraints of objective functions, so, while the small differences in objective values can be used to distinguish between the solutions, the differences are not significant. One solution to this limitation would involve the use of constraint scenarios, e.g., setting the two values as the same value if the difference in BHPs is negligible from the operator's view. The simulation results confirm the typical characteristics of a highly mobile CO 2 plume in a brine-saturated system; it moves upward in the early period, reaches the impermeable cap rock, and then spreads horizontally. Thus, the geological features of the upper part of heterogeneous aquifers govern the plume flow. The trapping amounts are mainly influenced by the pore volume, while CO 2 movement is significantly affected by the aquifer heterogeneity. Assessing uncertainty is an important step in elucidating the characteristics of CO 2 sequestration in heterogeneous media. To achieve more realistic applications, detailed descriptions of the geological as well as fluid properties-e.g., capillary trapping and hysteresis-are essential [30].
The challenges associated with this work are as follows. An individual well has been constrained to an allocated fixed rate during the entire period of the CO 2 injecting process. While this does reduce the enormous number of potential cases, this constraint limits this method to making rational decisions based on continuously updating dynamic data. The objective functions in the optimization problems must be appropriately relevant to result in trade-offs that satisfy the decision makers. The conflicting goals that may arise represent a problem in the proper decision-making process. This paper focuses on the sensitivity of spatial properties related to geological uncertainty. Further sensitivity analyses that include the operating conditions would be able to evaluate whether or not the objective functions have been reliably defined. In addition, this paper does not consider facility cost, so future studies must extend upon this research by classifying operating scenarios in terms of both the injection pressures and the facility cost. Although the shut-in and out of well operations are frequent in field applications, this simulation does not consider operational flexibility. It is crucial to determine the optimum amount of CO 2 sequestration without leakage, but aquifer heterogeneity makes quantitative analyses more difficult [25,26]. The optimal solutions determined in this way may not be the best ones. Depending on the constraints, various strategies can be used to address this issue. However, it is challenging to set the constraints in such a way that demonstrates realistic operations [31]. Multiobjective optimization implements a non-dominated relationship to move the solution fronts to the lowest rank. If a small change occurs in the objective values, then the nondomination is broken, thus leading to the return of unreasonable and inefficient strategies. The local minima problem has been a representative limitation for the optimization of multi-objective functions. To avoid local minima, it is important to both conserve the solution diversity and improve the convergence to the global optimality. The optimization problems must consider both the efficiency of on-site operation and cost-effectiveness. To guarantee rationality in decision making, it is also important to set a good indicator as the objective function. A future study could test the use of efficient surrogate modeling as a replacement for time-consuming numerical simulation to enhance the optimization processes. This paper shows that the spatial heterogeneity affects the CO 2 sequestration performance. However, asymmetric parameter interaction with scale-variant features represents another challenge in analyzing multi-phase flow in heterogeneous media.

Conclusions
This paper discussed the multi-objective problem related to CO 2 sequestration under geological uncertainty. It searched for optimal well allocations to achieve both minimum bottom hole pressures as well as maximum trapping efficiency. The distance-based generalized sensitivity analysis confirmed the influence of the following spatial properties on CO 2 sequestration, in descending order of influence: mean sandstone porosity, mean sandstone permeability, shale volume ratio, and the Dykstra-Parsons coefficient. The optimal well allocation suggested a large amount of injection through one well, and it induced a wide distribution of CO 2 plume while maintaining low injection pressures. Numerical simulations reliably explained the transport of the highly mobile CO 2 plume in brine-saturated media; the plume moved upward at the beginning of injection, reached the impermeable cap rock, and then spread horizontally. Obtaining a detailed description of the properties at a vertical section near an injector is an essential step toward achieving quantitative analyses for CO 2 sequestration. As the aquifer became heterogeneous, the impermeable shale was able to reduce the plume mobility, while it was unable to maintain the low injection pressure. The results demonstrate the importance of considering geological uncertainty in analyzing the performance of CO 2 sequestration in a more realistic manner. Asymmetric parameter interactions with spatiotemporal variables remain a challenge for the analysis of quantitatively supercritical CO 2 flow under geological uncertainty.