Faruk O. Alpak, Mohamad Jammoul, Mary F. Wheeler, Kachi Onyeagoro
{"title":"Less-Intrusive Consistent Discretization Methods for Reservoir Simulation on Cut-cell Grids – Algorithms, Implementation, and Testing","authors":"Faruk O. Alpak, Mohamad Jammoul, Mary F. Wheeler, Kachi Onyeagoro","doi":"10.1007/s10596-024-10299-6","DOIUrl":null,"url":null,"abstract":"<p>Consistent discretization methods are a natural fit for the novel cut-cell gridding technique for reservoir simulation, which preserves the orthogonality characteristic in the lateral direction. Both uniform (global) and novel hybrid (local) variants of consistent discretization methods are implemented and tested in the vicinity of fault representations on cut-cell grids. Novel consistent discretization methods, which do not require major intrusive changes to the solver structure of industrial-grade reservoir simulators, are investigated in this work. Cell-centered methods such as multi-point flux approximation (MPFA), average multi-point flux approximation (AvgMPFA), and nonlinear two-point flux approximation (NTPFA) methods fit naturally into the framework of existing industrial-grade simulators. Thus, cut-cell compatible variants of AvgMPFA and NTPFA and their novel hybridizations with TPFA are implemented and tested. An implementation of the relatively more computationally expensive MPFA is also made to serve as accuracy reference to AvgMPFA and NTPFA. AvgMPFA and NTPFA multiphase simulation results are compared in terms of accuracy and computational performance against the ones computed with reference MPFA and TPFA methods on a set of synthetic cut-cell grid models of varying complexity including conceptual models and a field-scale model. It is observed that AvgMPFA consistently yields more accurate and computationally efficient simulations than NTPFA on cut-cell grids. Moreover, AvgMPFA-TPFA hybrids run faster than NTPFA-TPFA hybrids when compared on the same problem for the same hybridization strategy. On the other hand, the computational performance of AvgMPFA degrades more rapidly compared to NTPFA with increasing “rings” of orthogonal blocks around cut-cells. Auspiciously, only one or two “rings” of orthogonal blocks around cut cells are sufficient for AvgMPFA to deliver high accuracy.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"66 1","pages":""},"PeriodicalIF":2.1000,"publicationDate":"2024-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational Geosciences","FirstCategoryId":"89","ListUrlMain":"https://doi.org/10.1007/s10596-024-10299-6","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Consistent discretization methods are a natural fit for the novel cut-cell gridding technique for reservoir simulation, which preserves the orthogonality characteristic in the lateral direction. Both uniform (global) and novel hybrid (local) variants of consistent discretization methods are implemented and tested in the vicinity of fault representations on cut-cell grids. Novel consistent discretization methods, which do not require major intrusive changes to the solver structure of industrial-grade reservoir simulators, are investigated in this work. Cell-centered methods such as multi-point flux approximation (MPFA), average multi-point flux approximation (AvgMPFA), and nonlinear two-point flux approximation (NTPFA) methods fit naturally into the framework of existing industrial-grade simulators. Thus, cut-cell compatible variants of AvgMPFA and NTPFA and their novel hybridizations with TPFA are implemented and tested. An implementation of the relatively more computationally expensive MPFA is also made to serve as accuracy reference to AvgMPFA and NTPFA. AvgMPFA and NTPFA multiphase simulation results are compared in terms of accuracy and computational performance against the ones computed with reference MPFA and TPFA methods on a set of synthetic cut-cell grid models of varying complexity including conceptual models and a field-scale model. It is observed that AvgMPFA consistently yields more accurate and computationally efficient simulations than NTPFA on cut-cell grids. Moreover, AvgMPFA-TPFA hybrids run faster than NTPFA-TPFA hybrids when compared on the same problem for the same hybridization strategy. On the other hand, the computational performance of AvgMPFA degrades more rapidly compared to NTPFA with increasing “rings” of orthogonal blocks around cut-cells. Auspiciously, only one or two “rings” of orthogonal blocks around cut cells are sufficient for AvgMPFA to deliver high accuracy.
期刊介绍:
Computational Geosciences publishes high quality papers on mathematical modeling, simulation, numerical analysis, and other computational aspects of the geosciences. In particular the journal is focused on advanced numerical methods for the simulation of subsurface flow and transport, and associated aspects such as discretization, gridding, upscaling, optimization, data assimilation, uncertainty assessment, and high performance parallel and grid computing.
Papers treating similar topics but with applications to other fields in the geosciences, such as geomechanics, geophysics, oceanography, or meteorology, will also be considered.
The journal provides a platform for interaction and multidisciplinary collaboration among diverse scientific groups, from both academia and industry, which share an interest in developing mathematical models and efficient algorithms for solving them, such as mathematicians, engineers, chemists, physicists, and geoscientists.