Simulation of foam in Enhanced Oil Recovery (EOR) and Carbon Capture and Storage (CCS) applications

Foam is an agglomeration of gas bubbles separated from each other by thin liquid films. Foam injection has been proved to be an efficient way to increase sweep efficiency in gas-injection enhanced-oil-recovery (EOR) processes or enlarge the storage space for trapping of CO2 in aquifers. Numerical simulation, an important tool developed by combining physics, mathematics, and computer programming, provides an efficient way to understand the complex fluid flow in subsurface reservoirs with applications to the evaluation of hydrocarbon recovery, energy efficiency, performance analysis, and various optimization problems. However, simulation of foam-related displacement in reservoirs is still an expensive process for conventional simulation due to the strongly nonlinear physics. In this dissertation, a newly developed numerical simulator, called Delft Advanced Research Terra Simulator (DARTS), is introduced to investigate the complicated flow behavior in porous media with the presence of foam.Petroleum Engineerin

Hybrid Berth Allocation for Bulk Ports with Unavailability and Stock Level Constraints

Berth allocation is fundamental to port-related operations in maritime shipping. Port managers have to deal with the increasing demands either by expanding the terminals or by improving efficiency to maintain competitiveness. Port expansion is a long-term project, and it requires much capital investment. Thus, the question of how to enhance the efficiency of berth allocation has received much research interest. Research on the Berth Allocation Problem (BAP) in container ports is quite advanced. However, only limited research focuses on BAP in bulk ports, although some similarities exist. Contributing to Operations Research approaches on the BAP, this paper develops a hybrid BAP mixed-integer optimization model dedicated to bulk ports. In addition to considering the handling characteristics of bulk ports, we also incorporate more practical factors such as unavailability and stock levels. The objective of the proposed model is to minimize the demurrage fee for all vessels under consideration of unavailability and stock constraints. We use the commercial software CPLEX to obtain the optimal solutions for a set of distinct instances, explicitly considering the situation of multiple cargo types on one vessel, which provides a better fit for the loading or discharging operations in real-world bulk ports. This is the first study to our knowledge that dedicates itself to the BAP in bulk ports and considers unavailability and stock constraints simultaneously. Our solutions can provide timely and effective decision support to bulk port managers.Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Transport Engineering and Logistic

Operator-based linearization approach for modeling of multiphase flow with buoyancy and capillarity

Numerical simulation of coupled multiphase multicomponent flow and transport in porous media is a crucial tool for understanding and forecasting of complex industrial applications related to the subsurface. The discretized governing equations are highly nonlinear and usually need to be solved with Newton's method, which corresponds with high computational cost and complexity. With the presence of capillary and gravity forces, the nonlinearity of the problem is amplified even further, which usually leads to a higher numerical cost. A recently proposed operator-based linearization (OBL) approach effectively improves the performance of complex physical modeling by transforming the discretized nonlinear conservation equations into a quasilinear form according to state-dependent operators. These operators are approximated by means of a discrete representation on a uniform mesh in physical parameter space. Continuous representation is achieved through the multilinear interpolation. This approach provides a unique framework for the multifidelity representation of physics in general-purpose simulation. The applicability of the OBL approach was demonstrated for various energy subsurface applications with multiphase flow of mass and heat in the presence of buoyancy and diffusive forces. In this work, the OBL approach is extended for multiphase multicomponent systems with capillarity. Through the comparisons with a legacy commercial simulator using a set of benchmark tests, we demonstrate that the extended OBL scheme significantly improves the computational efficiency with the controlled accuracy of approximation and converges to the results of the conventional continuous approach with an increased resolution of parametrization.Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Petroleum Engineerin

Accurate Modeling of Near-Wellbore Effects Induced by Supercritical CO2 Injection

During injection of supercritical CO2 into saline aquifers or depleted gas reservoirs, the complex interaction of CO2 and impurities with reservoir fluids plays a very important role and can significantly alternate the injectivity. Brine evaporation into the CO2-rich phase can lead to salt precipitation which will reduce the effective permeability of the porous rock. A tangible cooling of the near-wellbore region due to the Joule-Thomson effect can lead to hydrate formation which will reduce injectivity even more. Complex phase behavior of supercritical CO2 with brine and hydrocarbon components in highly heterogeneous porous media accompanied by all these phenomena will strongly affect pressure distribution which is in turn related to mechanical risks.In this work, we present a unified simulation framework for modelling near-wellbore effects induced by supercritical CO2 injection developed in the Delft Advanced Research Terra Simulator (DARTS) platform. This framework uses the Operator-Based Linearization (OBL) technique for incorporating all complex physical phenomena in a fully coupled fully implicit manner. A general multicomponent multiphase flash based on a combination of classic cubic equations of state (e.g., Peng-Robinson) for hydrocarbon/CO2-rich phases and an activity model for the aqueous phase is implemented. Hydrate phase behavior is modelled using a modified Van der Waals-Platteeuw hydrate equation of state. Formation dry-out and salt precipitation are incorporated by using the Element Balance approach coupled with thermodynamics. Thermophysical property correlations relevant to the thermodynamic conditions of interest are implemented and validated against lab experiments.We demonstrate that all important physical phenomena, such as the Joule-Thomson effect, hydrate formation and salt precipitation can be effectively captured by the OBL approach. We use several existing numerical benchmarks to validate the accuracy of the developed framework in the dynamic representation of all these effects. The interplay between these complex phenomena and reservoir heterogeneity is demonstrated in an unstructured heterogeneous near-wellbore reservoir model.Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Reservoir Engineerin

Simulation of foam-assisted co2 storage in saline aquifers

Geological storage of CO2 is a crucial emerging technology to reduce anthropogenic greenhouse gas emissions. Due to the buoyant characteristic of injected gas and the complex geology of subsurface reservoirs, most injected CO2 either rapidly migrates to the top of the reservoir or fingers through high-permeability layers due to instability in the convection-dominated displacement. Both of these phenomena reduce the storage capacity of subsurface media. CO2-foam injection is a promising technology for reducing gas mobility and increasing trapping within the swept region in deep brine aquifers. A consistent thermodynamic model based on a combination of a classic cubic equation of state (EOS) for gas components with an activity model for the aqueous phase has been implemented to describe the phase behavior of the CO2-brine system with impurities. This phase-behavior module is combined with representation of foam by an implicit-texture (IT) model with two flow regimes. This combination can accurately capture the complicated dynamics of miscible CO2 foam at various stages of the sequestration process. The Operator-Based Linearization (OBL) approach is applied to reduce the nonlinearity of the CO2-foam problem by transforming the discretized conservation equations into space-dependent and statedependent operators. Surfactant-alternating-gas (SAG) injection is applied to overcome injectivity problems related to pressure build-up in the near-well region. In this study, a 3D large-scale heterogeneous reservoir is used to examine CO2-foam behaviour and its effects on CO2 storage. Simulation studies show foams can reduce gas mobility effectively by trapping gas bubbles and inhibit CO2 from migrating upward in the presence of gravity, which in turn improves remarkably the sweep efficiency and opens the unswept region for CO2 storage. We also study how surfactant injection and forming of foam affect enhanced dissolution of CO2 at various thermodynamic conditions. This work provides a possible strategy to develop robust and efficient CO2 storage technology.Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Petroleum Engineerin

Gravity segregation with CO2 foam in heterogeneous reservoirs

Foam injection is one efficient way to mitigate gravity segregation during CO2 injection into porous media. The effect of gravity segregation on foam propagation in heterogeneous porous media is not yet fully resolved. To assess CO2 foam transport for enhanced oil recovery (EOR) and for CO2 storage processes in heterogeneous reservoirs, an accurate prediction of foam behavior is essential. In this study, we investigate the effect of heterogeneity on gravity segregation in the presence of foam. For nonlinear analysis, we use an extension of an Operator-Based Linearization (OBL) approach proposed recently. The OBL approach helps to reduce the nonlinearity of complex physical problems by transforming the discretized nonlinear conservation equations into a quasi-linear form based on state-dependent physical operators. The state-dependent operators are approximated by discrete representation on a uniform mesh in parameter space. In our study, foam in porous media is described using an implicit-Texture (IT) foam model with two flow regimes. We first validate the numerical accuracy of the foam simulation with OBL by comparing segregation length using the IT foam model with Newtonian rheology to analytical solutions. Next, the foam-model parameters are fit to foam-quality scan data for four sandstone formations ranging in permeability by an order of magnitude using a least-squares optimization approach. We then construct several hypothetical models containing two communicating layers with different permeability and thickness ratios to examine foam s effect on gravity segregation. The numerical results of the segregation length in homogeneous domains show good agreement with analytical solutions, except in a transition zone beneath the override zone which is not included in the analytical model. Through fractional-flow theory, we find that the transition zone is not a numerical artefact, but caused by low gas relative-mobility during the transient displacement process. Permeability affects both the mobility reduction of wet foam in the low-quality regime and the limiting capillary pressure at which foam collapses. Thus the segregation length varies with permeability and foam strength. In two-layer models, the thickness of the top layer plays an important role in the ultimate segregation length. A thin top layer does not affect segregation in the bottom layer, while a thicker top layer dominates the segregation length, with less influence of the bottom layer.Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Petroleum Engineerin

Modeling of Non-Newtonian Polymer Flooding with Adsorption and Retention Using Parametrization Approach

Polymer flooding is one efficient EOR technology by overcoming non-uniform and unstable displacement caused by water injection. Polymer flooding in reservoirs is a complicated process that involves strongly nonlinear physics, e.g., non-Newtonian rheology in porous media with retention and adsorption. In the presence of multi-scale heterogeneity, high-fidelity simulations are usually required to capture such nonlinear behavior, which is a time-consuming process for conventional reservoir modelling.In this study, we extend an advanced linearization strategy, called the Operator-Based Linearization (OBL) approach, to simulate non-Newtonian polymer flooding with retention and adsorption mechanisms using the fully implicit method. A velocity-dependent viscosity multiplier compliments the operator form of governing equations to represent the non-Newtonian rheology of the high-molecular-compound polymer. The retention of polymer, reducing the porosity, is represented by a Langmuir-type adsorption model. Several simplified models have been used for validation of the developed numerical framework. The numerical results show good agreement with both the analytical solutions and the coreflood experimental data though some negligible discrepancies can be observed in simulation results.A highly resolved near-well model is used to test the performance of polymer flooding in realistic reservoir conditions. Both shear-thinning and thickening regimes, depending on the injection velocity and polymer concentration, are recognized in the near-wellbore zone. The injected polymer concentration and brine salinity significantly affect the shear viscosity, and consequently, polymer injectivity. Polymer retention and adsorption have a substantial effect on the rate of polymer propagation through porous media. Overall, polymer flooding shows its advantages to mitigate water fingering in field-scale operations and improves the ultimate sweep of the reservoir. However, optimal injectivity is one essential factor that affects the performance of polymer flooding. The computational superiority of the proposed model allows us to optimize the parameters of polymer flooding in realistic reservoirs and operational settings.Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Reservoir Engineerin

Collaborative Berth Allocation with Row Generation Methods for the Core and Nucleolus

Transport Engineering and Logistic

Simulation of foam enhanced-oil-recovery processes using operator-based linearization approach

Foam injection is a promising enhanced-oil-recovery (EOR) technology that significantly improves the sweep efficiency of gas injection. Simulation of foam/oil displacement in reservoirs is an expensive process for conventional simulation because of the strongly nonlinear physics, such as multiphase flow and transport with oil/foam interactions. In this work, an operator-based linearization (OBL) approach, combined with the representation of foam by an implicit-texture (IT) model with two flow regimes, is extended for the simulation of the foam EOR process. The OBL approach improves the efficiency of the highly nonlinear foam-simulation problem by transforming the discretized nonlinear conservation equations into a quasilinear form using state-dependent operators. The state-dependent operators are approximated by discrete representation on a uniform mesh in parameter space. The numerical-simulation results are validated by using three-phase fractional-flow theory for foam/oil flow. Starting with an initial guess depending on the fitting of steady-state experimental data with oil, the OBL foam model is regressed to experimental observations using a gradient-optimization technique. A series of numerical validation studies is performed to investigate the accuracy of the proposed approach. The numerical model shows good agreement with analytical solutions at different conditions and with different foam parameters. With finer grids, the resolution of the simulation is better, but at the cost of more expensive computations. The foam-quality scan is accurately fitted to steady-state experimental data, except in the low-quality regime. In this regime, the used IT foam model cannot capture the upward-tilting pressure gradient (or apparent viscosity) contours. 1D and 3D simulation results clearly demonstrate two stages of foam propagation from inlet to outlet, as seen in the computed-tomography (CT) coreflood experiments: weak foam displaces most of the oil, followed by a propagation of stronger foam at lower oil saturation. OBL is a direct method to reduce nonlinearity in complex physical problems, which can significantly improve computational performance. Taking its accuracy and efficiency into account, the data-driven OBL-based approach could serve as a platform for efficient numerical upscaling to field-scale applications.Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.Petroleum Engineerin

A Collaborative Berth Planning Approach for Disruption Recovery

Traditionally, terminal operators create an initial berthing plan before the arrival of incoming vessels. This plan involves decisions on when and where to load or discharge containers for the calling vessels. However, disruptive unforeseen events (i.e., arrival delays, equipment breakdowns, tides, or extreme weather) interfere with the implementation of this initial plan. For terminals, berths and quay cranes are both crucial resources, and their capacity limits the efficiency of port operations. Thus, one way to minimize the adverse effects caused by disruption is to ally different terminals to share berthing resources. In some challenging situations, terminal operators also need to consider the extensive transshipment connections between feeder and mother vessels. Therefore, in this work, we investigate a collaborative variant of the berth allocation recovery problem which focuses on the collaboration among terminals and transshipment connections between vessels. We propose a mixed-integer programming model to (re)-optimize the initial berth and quay crane allocation plan and develop a Squeaky Wheel Optimization metaheuristic to find near-optimal solutions for large-scale instances. The results from the performed computational experiments, considering multiple scenarios with disruptive events, show consistent improvements of up to 40% for the suggested collaborative strategy (in terms of costs for the terminal operators).Transport Engineering and Logistic