Fully compositional multi-scale reservoir simulation of various CO 2 sequestration mechanisms

Voskov, Denis; Henley, Heath; Lucia, Ángelo

doi:10.1016/j.compchemeng.2016.09.021

Cited by 27 publications

(8 citation statements)

References 53 publications

(57 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The proposed model is implemented in Automatic Differentiation General Purpose Research Simulator (AD-GPRS) (Voskov et al, 2017;Rin, 2017;Garipov et al, 2018). The coupled system of equations is spatially discretized employing a finite-volume scheme with a two-point flux approximation and a first-order upwind scheme for advection terms.…”

Section: Methodsmentioning

confidence: 99%

“…In this work, we provide a detailed description of numerical framework based on DBS and present a comparison between Darcy and DBS models based on numerical implementation in the new version of Automatic Differentiation General Purpose Research Simulator (AD-GPRS) (AD-GPRS, 2018;Voskov et al, 2017;Rin et al, 2017;Rin, 2017;Garipov et al, 2018). In the first section, we will briefly describe governing equations for the proposed framework.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Robust And Accurate Formulation For Modeling Of Acid Stimulation

Tomin¹,

Voskov²

2018

Proceedings

View full text Add to dashboard Cite

Accurate representation of processes associated with energy extraction from subsurface formations often requires models which account for chemical interactions between different species in the presence of multiphase flow. In this study, we focus on modeling of acid stimulation in the near-well region. For the chemical processes which include a dissolution of rock material, an issue arises with the predictive representation of flow. Taking into account the spatial scale of discretization, some of simulation control volumes can have values of porosity close to 1, which makes an application of Darcy's law inconsistent and requires employing a true momentum equation such as the Darcy-Brinkman-Stokes (DBS) equation. The DBS equation automatically switches the description between Darcy equation in control volumes with low porosity and Stokes equation in grid blocks with high porosity. For chemical reactions, we propose a local nonlinear solution technique that allows solving the balance of solid species separately yet retaining the full coupling with rest of the equations. Finally, we study the impact of multiphase flow. The DBS approach is not well established for multiphase flow description. Therefore we employ a hybrid approach, where we assume that the single-phase DBS flow and the multiphase Darcy flow occur in separate regions. We test the accuracy and performance of both approaches on realistic models of practical interest.

show abstract

Section: Methodsmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Robust And Accurate Formulation For Modeling Of Acid Stimulation

Tomin¹,

Voskov²

2018

Proceedings

View full text Add to dashboard Cite

show abstract

“…The description of kerogen pyrolysis and subsequent oil cracking kinetic reaction model is adapted from Braun and Burnham, 1992. The reaction module is fully coupled with the mass and energy balance of the system using an extension of natural formulation for a solid phase treatment (Voskov et al, 2016). Fluid properties and phase equilibrium are modeled using Peng-Robinson equation of state (EOS) (Peng and Robinson, 1976).…”

Section: Work Description and Research Objectivesmentioning

confidence: 99%

Hierarchical Coarsening of Simulation Model for In-Situ Upgrading Process

Fucinos

Voskov

Burnham

2017

Day 1 Mon, February 20, 2017

View full text Add to dashboard Cite

Faculty of Civil Engineering and Geosciences Department of Geosciences and Engineering Master of ScienceHierarchical coarsening of simulation model for in-situ upgrading process by Raul Fucinos M.Oil shales are sedimentary rocks containing organic matter in the form of kerogen which accounts for more than 5 trillion barrels of oil in place according to Birol, 2010; therefore, oil shales represents a plausible solution for the constantly increasing demand for hydrocarbons. Oil shale production is currently done using two different techniques, surface retorting and in-situ retorting. The last one being the focus of this study. During this process, the sedimentary rock containing the kerogen is brought into a high-temperature environment with oxygen deficit. At this stage, the organic matter is subject to a thermo-chemical decomposition that finally releases the hydrocarbon in liquid and gas forms. This process is also known as pyrolysis. During this process, solid and fluid components experience compositional and physical changes, which requires complex chemical models represented by multiple species and several governing relations.In this work, we first developed a numerical solver for closed systems with simple kinetics models. This initial work allowed us to analyze the dynamic behavior of each component during the chemical decomposition of the kerogen and its impact in the porosity of the system. Then, we described an accurate base model for chemical decomposition of kerogen. This model was then implemented in our inhouse simulator ADGPRS. The model is based on the most recent understanding of pyrolysis process, and it incorporates coupling of chemical kinetics to heat and mass transport. Due to the high number of species, variations of porosity as consequence of the transformation of solid species into fluid products and complex multi-scale structure of porous media, the simulation performance of the high-fidelity model is limited. Therefore, in the second step of this work, we introduce a hierarchy of coarser models to improve the run-time of forward solution without significant reduction in accuracy. We applied coarsening in time, space, and chemical representation, and quantify errors introduced at each coarsening level. In conclusion, we provided recommendations for large-scale modeling of in-situ upgrading process. v

show abstract

“…In addition, coupling with chemical reactions requires a combination of thermodynamic and chemical equilibria [17,24]. This can significantly increase the cost of phase-behavior computations in compositional simulation [34].…”

Section: Introductionmentioning

confidence: 99%

Optimization of CO2 injection using multi-scale reconstruction of composition transport

Chen

Voskov

2019

Comput Geosci

View full text Add to dashboard Cite

The current situation with green gas emission requires the development of low-carbon energy solutions. However, a significant part of the modern energy industry still relies on fossil fuels. To combine these two contradictory targets, we investigate a strategy based on a combination of CO 2 sequestration with enhanced oil recovery (EOR) in the hydrocarbon reservoirs. In such technology, the development of miscibility is the most attractive strategy from both technological and economic aspects. Modeling of this process involves solving complex nonlinear problem describing compositional flow and transport in highly heterogeneous porous media. An accurate capture of the miscibility development usually requires an extensive number of components to be present in the compositional problem which makes simulation run-time prohibitive for optimization. Here, we apply a multi-scale reconstructing of compositional transport to the optimization of CO 2 injection. In this approach, a prolongation operator, based on the parametrization of injection and production tie-lines, is constructed following the fractional flow theory. This operator is tabulated as a function of pressure and pseudo-composition which then is used in the operator-based linearization (OBL) framework for simulation. As a result, a pseudo two-component solution of the multidimensional problem will match the position of trailing and leading shocks of the original problem which helps to accurately predict phase distribution. The reconstructed multicomponent solution can be used then as an effective proxy-model mimicking the behavior of the original multicomponent system. Next, we use this proxy-model in the optimization procedure which helps to improve the performance of the process several fold. An additional benefit of the proposed methodology is based on the fact that important technological features of CO 2 injection process can be captured with lower degrees of freedom which makes the optimization solution more feasible.

show abstract

Fully compositional multi-scale reservoir simulation of various CO 2 sequestration mechanisms

Cited by 27 publications

References 53 publications

Robust And Accurate Formulation For Modeling Of Acid Stimulation

Robust And Accurate Formulation For Modeling Of Acid Stimulation

Hierarchical Coarsening of Simulation Model for In-Situ Upgrading Process

Optimization of CO2 injection using multi-scale reconstruction of composition transport

Contact Info

Product

Resources

About