1887

Abstract

Summary

New Sequential Fully Implicit (SFI) methods for compositional flow simulation have been recently investigated. These SFI schemes decouple the fully coupled problem into separate pressure and transport problems and have convergence properties comparable to the Fully Implicit (FI) method. The pressure system is a parabolic problem with fixed overall-compositions and the transport system is a hyperbolic problem with fixed pressure and total-velocity. We discuss some aspects of how to design optimal SFI schemes for compositional flow with general equation-of-states by localizing the computations. The different systems are solved sequentially and the Fully Implicit solution is recovered by controlling the a-posteriori splitting errors due to the choice of decoupling.

When the parabolic and the hyperbolic operators are separated, it is possible to design nonlinear domain decomposition schemes taking the advantage of the specific properties of each operator. Usually, for reservoir simulation models, most of the reservoir is converged with SFI methods in one outer-iteration. However, in some localized regions with strong coupling between the pressure and the compositions, the SFI algorithms may need several outer-iterations. Here we propose a domain decomposition method based on a predictor-corrector strategy. As a first step, the nonlinear parabolic pressure equation is solved on the whole domain with the Multiscale Restriction-Smooth Basis (MsRSB) method used as a linear domain decomposition solver. In a second step the compositions system is solved. At the end of this first outer-iteration, most of the reservoir is converged. Based on a-posteriori splitting-errors of the SFI scheme in volume and velocity, we define local regions where additional global outer-iterations would be required in the conventional SFI scheme. We then fix Dirichlet boundary conditions for the pressure and the compositions and solve local problems in these non-converged regions. After convergence of these smaller nonlinear problems, if the boundary conditions are changed by the updated regions, the global pressure problem is revisited. An additional post-processing of local transport iterations makes sure mass is conserved everywhere. The resulting algorithm converges to the same solution as the FI solver, with all simultaneous updates to composition and pressure in localized regions.

We demonstrate the robustness of this nonlinear domain decomposition algorithm across a wide parameter range. Realistic compositional models with gas and water injection are presented and discussed.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802128
2018-09-03
2024-03-28
Loading full text...

Full text loading...

References

  1. Acs, G., Doleschall, S. and Farkas, E.
    [1985] General Purpose Compositional Model. Society of Petroleum Engineers Journal, 25(04), 543–553.
    [Google Scholar]
  2. Coats, K.
    [2000] A Note on IMPES and Some IMPES-Based Simulation Models. SPE Journal, 5(03), 245–251.
    [Google Scholar]
  3. [2003] IMPES Stability: Selection of Stable Timesteps. SPE Journal, 8(02), 181–187.
    [Google Scholar]
  4. Hajibeygi, H. and Tchelepi, H.
    [2014] Compositional Multiscale Finite-Volume Formulation. SPE Journal, 19(02), 316–326.
    [Google Scholar]
  5. Hou, T. and Wu, X.H.
    [1997] A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media.Journal of Computational Physics, 134(1), 169–189.
    [Google Scholar]
  6. Jenny, P., Lee, S. and Tchelepi, H.
    [2003] Multi-scale finite-volume method for elliptic problems in subsurface flow simulation.Journal of Computational Physics, 187(1), 47–67.
    [Google Scholar]
  7. [2006] Adaptive fully implicit multi-scale finite-volume method for multi-phase flow and transport in heterogeneous porous media. Journal of Computational Physics, 217(2), 627–641.
    [Google Scholar]
  8. Killough, J. and Kossack, C.
    [1987] Fifth comparative solution project: evaluation of miscible flood simulators. In: SPE Symposium on Reservoir Simulation. Society of Petroleum Engineers.
    [Google Scholar]
  9. Krogstad, S., Lie, K.A., Møyner, O., Nilsen, H.M., Raynaud, X. and Skaflestad, B.
    [2015] MRST-AD–an open-source framework for rapid prototyping and evaluation of reservoir simulation problems. In: SPE reservoir simulation symposium. Society of Petroleum Engineers.
    [Google Scholar]
  10. Lee, S., Wolfsteiner, C. and Tchelepi, H.
    [2008] Multiscale finite-volume formulation for multiphase flow in porous media: black oil formulation of compressible, three-phase flow with gravity. Computational Geosciences, 12(3), 351–366.
    [Google Scholar]
  11. Mallison, B.T., Gerritsen, M.G., Jessen, K. and Orr, F.M.
    [2005] High order upwind schemes for two-phase, multicomponent flow. SPE Journal, 10(03), 297–311.
    [Google Scholar]
  12. Manzocchi, T., Carter, J.N., Skorstad, A., Fjellvoll, B., Stephen, K.D., Howell, J.A., Matthews, J.D., Walsh, J.J., Nepveu, M., Bos, C., Cole, J., Egberts, P., Flint, S., Hern, C., Holden, L., Hovland, H., Jackson, H., Kolbjørnsen, O., MacDonald, A., Nell, P.A.R., Onyeagoro, K., Strand, J., Syversveen, A.R., Tchistiakov, A., Yang, C., Yielding, G. and Zimmerman, R.W.
    [2008] Sensitivity of the impact of geological uncertainty on production from faulted and unfaulted shallow-marine oil reservoirs: objectives and methods. Petroleum Geoscience, 14(1), 3–15.
    [Google Scholar]
  13. Moncorgé, A., Tchelepi, H. and Jenny, P.
    [2017] Modified sequential fully implicit scheme for compositional flow simulation.Journal of Computational Physics, 337, 98–115.
    [Google Scholar]
  14. [2018] Sequential Fully Implicit Formulation for Compositional Simulation using Natural Variables.Journal of Computational Physics (To appear).
    [Google Scholar]
  15. Møyner, O. and Lie, K.A.
    [2016a] A Multiscale Restriction-Smoothed Basis Method for Compressible Black-Oil Models. SPE Journal, 21(06), 2079–2096.
    [Google Scholar]
  16. [2016b] A multiscale restriction-smoothed basis method for high contrast porous media represented on unstructured grids.Journal of Computational Physics, 304, 46–71.
    [Google Scholar]
  17. Møyner, O. and Tchelepi, H.
    [2018] A Mass-Conservative Sequential Implicit Multiscale Method for Isothermal Equation of State Compositional Problems.SPE Journal (To appear).
    [Google Scholar]
  18. MRST
    [2017b] The MATLAB Reservoir Simulation Toolbox. www.sintef.no/MRST.
  19. Trangenstein, J. and Bell, J.
    [1989] Mathematical Structure of Compositional Reservoir Simulation. SIAM Journal on Scientific and Statistical Computing, 10(5), 817–845.
    [Google Scholar]
  20. Voskov, D.
    [2012] An extended natural variable formulation for compositional simulation based on tie-line parameterization. Transport in porous media, 92(3), 541–557.
    [Google Scholar]
  21. Watts, J.
    [1986] A Compositional Formulation of the Pressure and Saturation Equations. SPE Reservoir Engineering, 1(03), 243–252.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802128
Loading
/content/papers/10.3997/2214-4609.201802128
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error