Simulation-Based Optimization of Microbial Enhanced Oil Recovery with a Model Integrating Temperature, Pressure, and Salinity Effects

: The microbial enhanced oil recovery (MEOR) method is an eco-friendly and economical alternative technology. The technology involves a variety of uncertainties, and its success depends on controlling microbial growth and metabolism. Though a few numerical studies have been carried out to reduce the uncertainties, no attempt has been made to consider temperature, pressure, and salinity in an integrated manner. In this study, a new modeling method incorporating these environmental impacts was proposed, and MEOR analysis was performed. As a result, accurate modeling was possible to prevent overestimating the performance of MEOR. In addition, oil recovery was maximized through sensitivity analysis and optimization based on an integrative model. Finally, applying MEOR to an actual reservoir model showed a 7% increase in oil recovery compared to waterﬂooding. This result proved the practical applicability of the method.


Introduction
Among various methods used to increase the productivity of the reservoir, chemical flooding, the gas injection method, and thermal methods are the most widely used processes for enhanced oil recovery (EOR) [1][2][3][4][5][6][7]. In addition to these methods, there is an alternative technology called microbial enhanced oil recovery (MEOR). The MEOR method is a technology that uses microbial metabolism and products to improve oil production [8,9]. This technology is eco-friendly in that the materials used are biodegradable, and it is economical in that the cost to produce them is low. The MEOR process changes the petrophysical and/or petrochemical properties of the reservoir system depending on the microorganisms and the microbial materials [10]. Among the microorganisms used in MEOR, Leuconostoc mesenteroiedes is a microbial species that produces a biopolymer called dextran. It is used to block the high permeability zone, which decreases productivity due to early water breakthrough after waterflooding. Dextran improves productivity by the bypass effect for injected water by selective plugging in the reservoir.
Since the 1950s, MEOR has been applied in many countries in various ways, among which selective plugging has also been used as a major mechanism to increase oil productivity. The North Burbank Unit in Oklahoma [11] and North Blowhorn Creek Oil Unit in Alabama [12], USA, showed increased oil production through selective plugging. In Canada, plugging by Leuconostoc was applied to the MEOR [13], and a similar method was attempted for the fractured reservoir [14]. Dutch researchers also showed the result of increasing oil productivity and improving water-oil ratio through the selective plugging effect by Betacoccus dextranicus [15]. In addition, MEORs using selective plugging have been applied in various fields in China, UK, Saudi Arabia, etc. [16]. However, for the result of applying MEOR, one out of ten was evaluated as ineffective [17] because proper application failed due to uncertainty in the method.
Various environmental factors affect the growth of microorganisms [10,[18][19][20]. Temperature has the greatest influence on the growth of microbes [21][22][23][24][25][26][27][28][29], and there are a range of possible growth temperatures and an optimum growth temperature depending on the microbial species. The effect of temperature on the rate of microbial growth is similar to the activity of an enzyme. The growth rate increases as the temperature increases to a certain point where protein denaturation occurs. The growth rate begins to decrease from the moment the temperature exceeds that temperature [30]. Pressure affects microbial survival [31][32][33][34], and higher pressure conditions adversely affect microbial survival. The destruction of microbial cells depends on the level of pressure and the time of exposure to the pressure [35]. Salinity is also known to affect microbial growth [36][37][38][39][40][41][42][43]. Salt concentration has a great influence on enzymatic catalysis [44]. In general, the salinity range of the reservoir is known to be 100 mg/L to 300 g/L or more [45]. According to Maudgalya et al. [46], successful cases of MEOR have been reported in a reservoir where the temperature is below 200 • F and the salinity is above 1000 ppm.
A few numerical studies on selective plugging by bacteria have been conducted. Early models described microbial growth based on the Monod equation, and oil recovery was calculated using a simple fluid flow model [47][48][49]. Delshad et al. [50] depicted MEOR using UTCHEM, and Stewart and Kim [51] simulated the selective plugging effect using a biofilm evolution-removal model. Vilcáez et al. [52] and Surasani et al. [53] studied Leconostoc mesenteroides growth and plugging phenomena using a model called CrunchFlow. Recently, selective plugging simulation and oil productivity analysis have been performed using CMG STARS [54,55]. These previous studies have not been able to simultaneously reflect the effects of temperature, pressure, and salinity on the growth of microorganisms, and have not performed optimization considering these environmental factors.
Most of the earliest studies were limited to describing microbial metabolism in simple transport models. They were not suitable for analyzing oil productivity for complex reservoir conditions in filed scale. Only a few recent studies have attempted to analyze MEOR effects in the reservoir scale. However, studies that have evaluated and optimized considering the subsurface environment are not available yet. A typical reservoir environment is high temperature, high pressure, and high salinity, which are fatal to the growth of microorganisms. If these factors are not considered, it is impossible to accurately evaluate the MEOR performance. Therefore, this study attempted to perform an MEOR simulation considering three environmental variables at the same time.
The first objective of this study was to develop a microbial growth model that comprehensively considers the effects of temperature, pressure, and salinity through the Arrhenius equation. Experimental results were used to verify the proposed model, and comparisons with previously developed models were made. Oil recovery was compared between the models with and without reflecting the environmental impact to accurately predict the MEOR result. Optimization was performed to maximize oil recovery. This process employed a response surface methodology, and the injection scenarios were set as design parameters. In addition, the sensitivity of each parameter to oil recovery was analyzed. Finally, the MEOR results for the actual heterogeneous reservoir were shown to represent the applicability of the technology.

Fluid Flow in Permeable Media
In this study, the multiphase and multicomponent flow in permeable media is calculated through the conservation equations as follows [56]: Here, subscript i is the component, subscript j is the phase, subscript s is the stationary phase, N P is the number of phases, N C is the number of components, ϕ is the porosity, S is the saturation, m is the mass fraction, → u is the superficial velocity, → → K is the dispersion tensor, and r is the kinetic reaction rate. The first term on the left side is the accumulation, the second term is the flux, and the right side is the source term.
Under spatially discretized conditions, CMG STARS used in this study calculates the conservation equation of flowing component i as follows: Here, subscripts w, o, g, k are water phase, oil phase, gas phase, and layer k, respectively. The w i , x i , y i are mole fractions of component i in water phase, oil phase, and gas phase, respectively. n f is the number of neighboring regions or grid block faces, K is the transmissibility, Φ is the potential, D ji is the component dispersibility, and q jk is the phase rate in layer k.

Stoichiometric Equations
The model bacteria used in this study were Leuconostoc mesenterioides. They generate a polysaccharidic biopolymer known as dextran under sucrose-rich conditions. Their metabolic mechanisms are well known, and modelling in this study was performed based on experimental results [57,58]. Reactions that happen during microbial metabolism were divided into four steps: (1) microbial growth, (2) sucrose hydrolysis, (3) dextran production, and (4) microbial decay.
The stoichiometric equation for microbial growth is described as follows [52]: Here, C 12 H 22 O 11 is sucrose, C 5 H 7 O 2 N is the model bacteria, C 6 H 12 O 6 is fructose, CH 3 CHOHCOO − is lactate, CH 3 COO − is acetate, C 6 H 14 O 6 is mannitol, and CH 3 CH 2 OH is ethanol.
The sucrose is hydrolyzed into glucose and fructose by microbial enzymatic reactions or is converted into dextran by the linking of glucosidic chains. The hydrolysis reaction is as follows: and the dextran formation is as follows [59]: The molecular weight of dextran measured in this study was 10,053 lb/mole, and based on this, n was assumed to be 6.2.
To account for the microbial decay reaction, we followed Bültemeier et al. [60], who assumed that the dead microbes did not participate in the chemical reactions but were changed to water as follows:

Microbial Reaction Rate
In this study, the rate of microbial reactions was described using the Arrhenius equation. It is mainly used in chemical reactions instead of the Monod equation, which is widely used in environmental engineering. By adding a division factor to the general Arrhenius equation, the role of sucrose as a limiting factor is described as follows: Here, r is reaction rate, F Freq is the frequency factor, F div is the division factor, E a is the activation energy, R is the gas constant, T is temperature, n c is the number of components, c i is the concentration of component i, c sucrose is the sucrose concentration, and A and B are constants.
The microbial decay rate is expressed as a function related to the frequency factor and the microbial population as follows: where r decay is the reaction rate for bacterial decay.

Temperature
The cardinal temperature model proposed by Rosso et al. [61] is expressed as follows: Here, r temp is the reaction rate as a function of temperature, r opt is optimum reaction rate, T opt is the optimum temperature for microbial growth, T max is the upper limit for growth temperature, and T min is the lower limit for growth temperature. This temperature model is described by modifying the exponential term of the Arrhenius equation as follows [62]: where A j is the constant of the j-th step, E a,j is the activation energy of the j-th step, and T j is the temperature of the j-th step. The final form of the equation reflecting the temperature effect is given as:

Pressure
Basak et al. [35] developed a decay model of L. mesenteroides as a function of pressure. The population of microbes followed first-order kinetics during the pressurehold time as follows: where N o is the initial number of microbes, N is the number of surviving microbes, α is the rate constant, and t p is the pressure-hold time. The time taken to reduce the number of microbes to one tenth is defined as a decimal reduction time (D), and Equation (13) is as D = 2.303/k. The relationship between pressure and decimal time is as follows: where D j is the decimal time of the j-th step, p j is the pressure of the j-th step, and z p is the negative reciprocal slope of log D vs. p. Experiments showed that log D and p have a linear relationship. Through this relationship, α in Equation (13) can be estimated and microbial destruction due to pressure can be calculated. In this study, the frequency factor in Equation (7) is determined as a constant by pressure as follows: r decay,p = F Freq,p c bacteria (15) where r decay,p is the reaction rate as affected by pressure, and F Freq,p is the frequency factor for each pressure step. F Freq,p is then used as a matching parameter for history match, and the value for each pressure step is derived.

Salinity
Leroi et al. [63] presented the reaction rate as a function of NaCl concentration based on the cardinal model as follows: where r NaCl is the reaction rate as a function of NaCl concentration, c NaCl is NaCl concentration, c NaCl opt is the optimum NaCl concentration for microbial growth, and c NaCl max is the upper limit of NaCl concentration for microbial growth. Since the optimal NaCl concentration of the model bacteria was close to zero [64][65][66], Equation (16) can be defined as follows: r NaCl = r opt − r opt c NaCl c NaCl max 2 (18) In this study, the growth inhibition by NaCl is described as the backward reaction of Equation (2), and the rate is expressed as follows: where r NaCl,backward is the backward reaction rate by NaCl, and n is a constant.

Permeability Reduction Model
The generated dextran is insoluble, and it is adsorbed or trapped in the pores, which lowers the permeability of the reservoir. In this study, dextran was assumed to be a solid phase, and the change in pore volume caused by the produced dextran affected the change in permeability as follows: where ϕ is the changed porosity, ϕ o is the initial porosity, c s is the concentration of the solid phase in pore space, ρ s is the density of the solid phase, k is the changed permeability, k o is the initial permeability, and k mul is the multiplier factor. c s /ρ s represents the volume of the solid phase in the pore space. When the volume occupied by the solid phase in the pores increases, the porosity decreases, which leads to a decrease in permeability.

History Matching Error
The history matching error means the relative difference between the measured data and simulation results for each objective function. The matching error can be estimated as follows [62]: Here, E is the history matching error, t and t elap are time and elapsed time, respectively, Y s t is the simulated result, Y m t is the measured result, and s is the normalization scale. In this study, the normalization scale is the maximum of the following quantities: where ∆Y m is the measured maximum change, and M err is the measurement error. The value of M err means that if the simulated result is between (measured value − M err ) and (measured value + M err ), the result of this matching is satisfactory. That is, M err is half the absolute error range.

Proxy Model for Sensitivity Analysis and Optimization
In this study, an optimization process was performed using the response surface methodology (RSM) to maximize the oil recovery. The RSM is a method to search the correlation between input parameter and response (objective function). The main idea is to use a proxy model to represent the original simulation result. Linear or quadratic form is mainly used to make a proxy model, and the latter is applied in this study. This proxy model is also used to analyze the sensitivity of each parameter to the objective function. A quadratic polynomial method was applied to build a proxy model. The quadratic model is expressed as follows [62,67]: Here, y is the objective function (response), x j expresses the linear effect of parameters, x 2 j are the quadratic effects of the parameters, x i x j indicate the interaction effects of parameters, a 0 is the intercept, a j are the coefficients of linear terms, a jj are the coefficients of quadratic terms, and a ij are the coefficients of interaction terms.
The sensitivity indicates how much a parameter change affects the objective function change. The greater the value of coefficient a (effect estimate), the greater the sensitivity of the parameter to the objective function. At this time, each parameter must have scale invariance. All parameters have an average value of zero and are set to vary from −1 to 1.

Batch and Sandpack Model Simulation
The results of microbial growth, dextran formation, and the subsequent decrease in permeability were taken from previous studies [54,58]. Verification of the microbial growth and dextran production was performed through the batch model shown in Table A1. Verification of the permeability reduction was performed using the sandpack simulation shown in Table A2. The batch simulation was performed under the conditions of 1 atm and 30 • C for a total volume of 1 L. The initial sucrose concentration in the batch model was 15 g/L, and the growth was carried out for 24 h. In the sandpack simulation, sucrose injection and soaking were conducted 20 times for 44 h. The concentration of the injection was 15 g/L, and the flow rate was 0.018 bbl/day. The detailed process of the simulations was referenced from Jeong et al. [54].
The reaction rates of the four types of microbial metabolism are calculated using Equations (6) and (7), and the history matching for the experiment was performed with matching parameters F Freq , A, and B. The values of the variables minimizing the matching error are shown in Table 1. Figure 1a,b describes the microbial growth and dextran production, respectively. The pressure difference measured in the pressure measurements shown in Table A2 is shown in Figure 1c, and the simulation results are matched with the experimental data. The matching parameter in this simulation was k mul in Equation (21), and it was calculated to be 244.8. The history match errors measured in the batch model and sandpack model are 1.6% and 4.6%, respectively.

Verification of Models Reflecting Environmental Factors
The previous batch model was simulated under constant temperature, pressure, and salinity conditions. Comparisons with existing methods were conducted after applying the proposed environmental factor reflection model to this model. The reaction rate calculated by the proposed model was compared with the result calculated by the cardinal model, and the results are depicted in Figure 2a. Microbial growth temperature conditions for calculating Equation (9) were referenced in previous studies [55]. It was shown that  Figure   2b shows the result of calculating the microbial decay as a function of the pressure step using Equations (13) and (14). This graph shows the microbial decay rate according to the pressure used in this study, and the same result is calculated with Equation (15). Microbial growth behavior as a function of NaCl concentration is depicted in Figure 2c. As shown, the salinity effect on bacterial growth was described using Equation (19). Like the existing model, the proposed model was able to derive results that well reflected microbial growth and decay reactions according to temperature, pressure, and salinity.

Verification of Models Reflecting Environmental Factors
The previous batch model was simulated under constant temperature, pressure, and salinity conditions. Comparisons with existing methods were conducted after applying the proposed environmental factor reflection model to this model. The reaction rate calculated by the proposed model was compared with the result calculated by the cardinal model, and the results are depicted in Figure 2a. Microbial growth temperature conditions for calculating Equation (9) were referenced in previous studies [55]. It was shown that f T j in Equation (10) could accurately calculate the γ(T) of the cardinal model. Figure 2b shows the result of calculating the microbial decay as a function of the pressure step using Equations (13) and (14). This graph shows the microbial decay rate according to the pressure used in this study, and the same result is calculated with Equation (15). Microbial growth behavior as a function of NaCl concentration is depicted in Figure 2c. As shown, the salinity effect on bacterial growth was described using Equation (19). Like the existing model, the proposed model was able to derive results that well reflected microbial growth and decay reactions according to temperature, pressure, and salinity.

Optimization of MEOR
The conditions of the reservoir used for the MEOR process are indicated in Table A3. The reservoir had a high probability of a conformance problem because of the high transmittance channels. Production was assumed to continue for three years. After implementing waterflooding during the first year, MEOR treatment was performed for eight weeks, and then waterflooding was implemented again. The model that did not reflect the environmental impact was named the "ideal case", and the model accounting for the environmental impact was named the "real case". Before the optimization, the injection rate in both cases was 50 bbl/day, and the sucrose concentration was 0.005 (mole fraction). In the real case, the temperature of the injected water and the NaCl concentration were set to 75 °F and 1%, respectively. These no-optimization models were used as a base case for each case. As shown in Figure 3, the model considering environmental factors tended to decrease the productivity of oil by 21% compared to the ideal case. Nevertheless, it could be seen that the oil recovery of MEOR in the real case was 25% higher than that of waterflooding. 40 60 80 100 120

Optimization of MEOR
The conditions of the reservoir used for the MEOR process are indicated in Table A3. The reservoir had a high probability of a conformance problem because of the high transmittance channels. Production was assumed to continue for three years. After implementing waterflooding during the first year, MEOR treatment was performed for eight weeks, and then waterflooding was implemented again. The model that did not reflect the environmental impact was named the "ideal case", and the model accounting for the environmental impact was named the "real case". Before the optimization, the injection rate in both cases was 50 bbl/day, and the sucrose concentration was 0.005 (mole fraction). In the real case, the temperature of the injected water and the NaCl concentration were set to 75 • F and 1%, respectively. These no-optimization models were used as a base case for each case. As shown in Figure 3, the model considering environmental factors tended to decrease the productivity of oil by 21% compared to the ideal case. Nevertheless, it could be seen that the oil recovery of MEOR in the real case was 25% higher than that of waterflooding. To maximize oil recovery, the ideal and real case injection scenarios were optimized using a proxy model created by the RSM. The proxy model set the oil recovery as the objective function and calculated design parameters that maximize the oil recovery. The design parameters for optimum injection included sucrose concentration, treatment period, injection rate, temperature, and salinity. The detailed ranges of all parameters and the optimum designs calculated by the proxy model are listed in Table 2. The optimal injection design differs in the two cases, especially in sucrose concentration and injection temperature. In addition, the optimum temperature for injection was measured somewhat lower than the optimum temperature for microbial growth, i.e., 86 °F [55]. The influence of each design parameter on oil recovery could be explained based on the sensitivity analysis. Injection salinity (%) --0-3 0 Figure 4 shows a tornado plot resulting from the sensitivity analysis of the design parameters. The values in figures called effect estimates are the expected change of the oil recovery when the scaled parameters travel from −1 to 1. In other words, the larger these values are, the greater the influence on oil recovery change. The sensitivity analysis showed that there were differences in the priority of parameters affecting oil recovery in both cases. In both cases, the quadratic effect of the treatment period had the greatest effect on recovery. However, the quadratic effect of injection temperature had the second largest effect on recovery in the real case. The maximum oil production predicted through this optimization is shown in Figure 5. The maximum oil recovery for the ideal case was 48.3%, which was a 33% increase over the no-optimized model. The recovery of the optimized To maximize oil recovery, the ideal and real case injection scenarios were optimized using a proxy model created by the RSM. The proxy model set the oil recovery as the objective function and calculated design parameters that maximize the oil recovery. The design parameters for optimum injection included sucrose concentration, treatment period, injection rate, temperature, and salinity. The detailed ranges of all parameters and the optimum designs calculated by the proxy model are listed in Table 2. The optimal injection design differs in the two cases, especially in sucrose concentration and injection temperature. In addition, the optimum temperature for injection was measured somewhat lower than the optimum temperature for microbial growth, i.e., 86 • F [55]. The influence of each design parameter on oil recovery could be explained based on the sensitivity analysis. Injection salinity (%) --0-3 0 Figure 4 shows a tornado plot resulting from the sensitivity analysis of the design parameters. The values in figures called effect estimates are the expected change of the oil recovery when the scaled parameters travel from −1 to 1. In other words, the larger these values are, the greater the influence on oil recovery change. The sensitivity analysis showed that there were differences in the priority of parameters affecting oil recovery in both cases. In both cases, the quadratic effect of the treatment period had the greatest effect on recovery. However, the quadratic effect of injection temperature had the second largest effect on recovery in the real case. The maximum oil production predicted through this optimization is shown in Figure 5. The maximum oil recovery for the ideal case was 48.3%, which was a 33% increase over the no-optimized model. The recovery of the optimized real case was 41.8%, which was a 47% increase over the no-optimized model. These results showed that it was important to consider environmental effects for the accurate prediction and optimization of productivity. The actual proxy models used for optimization are the same as Equations (A1) and (A2). real case was 41.8%, which was a 47% increase over the no-optimized model. These results showed that it was important to consider environmental effects for the accurate prediction and optimization of productivity. The actual proxy models used for optimization are the same as Equations (A1) and (A2).

Field Application
To evaluate the effectiveness of MEOR through the developed model, the process was applied to the highly-heterogeneous actual field model based on logging data from King Island [53]. The reservoir conditions are shown in Table 3. Since the reservoir depth was approximately 5000 ft, the reservoir pressure was assumed to be 2500 psi, which was calculated using a typical pressure gradient of 0.5 psi/ft. The reservoir pressure was also assumed to be 100 °F. Waterflooding and MEOR treatment were applied in the same way as the previous models, and the injection temperature and salinity were 80 °F and 1% NaCl, respectively. Figure 6 represents the influence of the MEOR process on the oil recovery and water-oil ratio (WOR). The recovery factor of MEOR was 36.4%, which was improved by 2.6% compared with waterflooding. WOR was also reduced by the MEOR, which was decreased by 11%. The distribution of dextran produced in the reservoir and the resulting permeability reduction in the highly permeable zone are shown in Figure 7. Most of the injected fluid flowed through the highly permeable layers, and it was found that most of the dextran was generated near the injection well. At this time, even if dextran was formed in the high permeability layers as shown in Figure 7a, the effect of reducing the permeability was not very large, as shown in Figure 7b.   real case was 41.8%, which was a 47% increase over the no-optimized model. These results showed that it was important to consider environmental effects for the accurate prediction and optimization of productivity. The actual proxy models used for optimization are the same as Equations (A1) and (A2).

Field Application
To evaluate the effectiveness of MEOR through the developed model, the process was applied to the highly-heterogeneous actual field model based on logging data from King Island [53]. The reservoir conditions are shown in Table 3. Since the reservoir depth was approximately 5000 ft, the reservoir pressure was assumed to be 2500 psi, which was calculated using a typical pressure gradient of 0.5 psi/ft. The reservoir pressure was also assumed to be 100 °F. Waterflooding and MEOR treatment were applied in the same way as the previous models, and the injection temperature and salinity were 80 °F and 1% NaCl, respectively. Figure 6 represents the influence of the MEOR process on the oil recovery and water-oil ratio (WOR). The recovery factor of MEOR was 36.4%, which was improved by 2.6% compared with waterflooding. WOR was also reduced by the MEOR, which was decreased by 11%. The distribution of dextran produced in the reservoir and the resulting permeability reduction in the highly permeable zone are shown in Figure 7. Most of the injected fluid flowed through the highly permeable layers, and it was found that most of the dextran was generated near the injection well. At this time, even if dextran was formed in the high permeability layers as shown in Figure 7a, the effect of reducing the permeability was not very large, as shown in Figure 7b.

Field Application
To evaluate the effectiveness of MEOR through the developed model, the process was applied to the highly-heterogeneous actual field model based on logging data from King Island [53]. The reservoir conditions are shown in Table 3. Since the reservoir depth was approximately 5000 ft, the reservoir pressure was assumed to be 2500 psi, which was calculated using a typical pressure gradient of 0.5 psi/ft. The reservoir pressure was also assumed to be 100 • F. Waterflooding and MEOR treatment were applied in the same way as the previous models, and the injection temperature and salinity were 80 • F and 1% NaCl, respectively. Figure 6 represents the influence of the MEOR process on the oil recovery and water-oil ratio (WOR). The recovery factor of MEOR was 36.4%, which was improved by 2.6% compared with waterflooding. WOR was also reduced by the MEOR, which was decreased by 11%. The distribution of dextran produced in the reservoir and the resulting permeability reduction in the highly permeable zone are shown in Figure 7. Most of the injected fluid flowed through the highly permeable layers, and it was found that most of the dextran was generated near the injection well. At this time, even if dextran was formed in the high permeability layers as shown in Figure 7a, the effect of reducing the permeability was not very large, as shown in Figure 7b.           Particularly in layered systems, the results of 2D and 3D models will not show a significant difference. However, some differences may occur due to the curvilinear flowline in the 3D model. To confirm this, the extension of the 2D model to the 3D heterogeneous model was shown in Figure 8a. Though the average permeability of each layer was the same as in the previous 2D model, the permeability distributions in the horizontal direction were geostatistically generated to have a Dykstra-Parsons coefficient of 0.4 (Figure 8b). This coefficient can be calculated as follows [68]: where V DP is the Dykstra-Parsons coefficient, k 50 is the mean permeability at 50% probability, and k 84.1 is the permeability at 84.1% probability. This coefficient is the most widely used to describe permeability variability in the petroleum industry. When it is zero, it means a homogeneous reservoir, and the closer to 1, the higher the heterogeneity. Particularly in layered systems, the results of 2D and 3D models will not show a significant difference. However, some differences may occur due to the curvilinear flowline in the 3D model. To confirm this, the extension of the 2D model to the 3D heterogeneous model was shown in Figure 8a. Though the average permeability of each layer was the same as in the previous 2D model, the permeability distributions in the horizontal direction were geostatistically generated to have a Dykstra-Parsons coefficient of 0.4 ( Figure  8b). This coefficient can be calculated as follows [68]: where DP V is the Dykstra-Parsons coefficient, 50 k is the mean permeability at 50% probability, and 1 . 84 k is the permeability at 84.1% probability. This coefficient is the most widely used to describe permeability variability in the petroleum industry. When it is zero, it means a homogeneous reservoir, and the closer to 1, the higher the heterogeneity. The injection conditions were also the same as the previous model, except that the injection rate increased by 9.8 times to 490 bbl/day in consideration of the increase in the reservoir volume. The results are shown in Figure 9. The oil recovery increased 3.2% in MEOR compared to waterflooding (Figure 9a). This result explained that when the reservoir was heterogeneous even in the horizontal direction, the effect of MEOR also increased compared to waterflooding. In addition, it was confirmed that a lot of dextran was generated in the vicinity of the injector in high permeable zones (Figure 9b). There was a slight difference between 2D and 3D results. Since the overall trend of the results was the same, subsequent optimization analysis was conducted on a 2D layered model that was advantageous in terms of computation time. The injection conditions were also the same as the previous model, except that the injection rate increased by 9.8 times to 490 bbl/day in consideration of the increase in the reservoir volume. The results are shown in Figure 9. The oil recovery increased 3.2% in MEOR compared to waterflooding (Figure 9a). This result explained that when the reservoir was heterogeneous even in the horizontal direction, the effect of MEOR also increased compared to waterflooding. In addition, it was confirmed that a lot of dextran was generated in the vicinity of the injector in high permeable zones (Figure 9b). There was a slight difference between 2D and 3D results. Since the overall trend of the results was the same, subsequent optimization analysis was conducted on a 2D layered model that was advantageous in terms of computation time. The optimization process was also applied to maximize oil recovery in this field. The range of design parameters for optimal injection and results are listed in Table 4 (the proxy model used is shown in Equation (A3)). Optimization increased the oil recovery to 37.8%, which was a 7% increase over waterflooding (Figure 10). Figure 11a shows the permeability reduction due to the use of MEOR. The permeability reduction in the high permeability The optimization process was also applied to maximize oil recovery in this field. The range of design parameters for optimal injection and results are listed in Table 4 (the proxy model used is shown in Equation (A3)). Optimization increased the oil recovery to 37.8%, which was a 7% increase over waterflooding ( Figure 10). Figure 11a shows the permeability reduction due to the use of MEOR. The permeability reduction in the high permeability zone was up to 18% in the no-optimization model and up to 41% in the optimized model. The oil saturation distribution of each layer is depicted in Figure 11b. Comparing the waterflooding and optimized models, oil saturation was reduced by up to 16% in some high permeability layers and up to 10% in some low permeability layers. Due to very large permeability contrast, the permeability of the high permeable zones was still larger than that of the low permeable zones even after plugging. This explains why the oil in the low permeable region was not swept sufficiently.

Conclusions
In this study, a new MEOR model incorporating environmental influences was developed to describe microbial metabolism in actual field conditions. An MEOR optimization method was discussed based on this model. The existing microbial growth models for temperature, pressure, and salinity were integrated using the Arrhenius method. To verify the proposed method, the batch model and sandpack model were employed, and

Conclusions
In this study, a new MEOR model incorporating environmental influences was developed to describe microbial metabolism in actual field conditions. An MEOR optimization method was discussed based on this model. The existing microbial growth models for temperature, pressure, and salinity were integrated using the Arrhenius method. To

Conclusions
In this study, a new MEOR model incorporating environmental influences was developed to describe microbial metabolism in actual field conditions. An MEOR optimization method was discussed based on this model. The existing microbial growth models for temperature, pressure, and salinity were integrated using the Arrhenius method. To verify the proposed method, the batch model and sandpack model were employed, and the environmental impacts on microbial growth were successfully accounted for using the proposed method. The importance of these environmental models was that they could accurately analyze MEOR productivity. Based on our analysis of ideal and real cases for the hypothetical reservoir, the productivity of MEOR was overestimated when the environmental impact was not considered. Additionally, a sensitivity analysis showed that the temperature effect had a great influence on the oil recovery in a real case. In particular, it was found that the injection temperature should consider both the reservoir temperature and the optimum temperature for microbial growth. Optimization showed that the optimal design of nutrient injection could further increase oil recovery. The King Island model, which was an actual reservoir, also showed that it was possible to predict the improved oil recovery accurately using the MEOR method that reflected environmental factors. This study showed that the MEOR analysis could be performed more accurately through the proposed method, and optimization for improved recovery efficiency was also essential.    Table A1. Batch model for bacterial growth and dextran generation [54].        Table A3. Description of the initial conditions of the target reservoir [54].

Permeability Distribution
Energies 2021, 14, x FOR PEER REVIEW Table A3. Description of the initial conditions of the target reservoir [54].