Felix Feldmann, Oddbjørn Nødland, Jan Sagen, Børre Antonsen, Terje Sira, Jan Ludvig Vinningland, Robert Moe, Aksel Hiorth
{"title":"IORSim: A Mathematical Workflow for Field-Scale Geochemistry Simulations in Porous Media","authors":"Felix Feldmann, Oddbjørn Nødland, Jan Sagen, Børre Antonsen, Terje Sira, Jan Ludvig Vinningland, Robert Moe, Aksel Hiorth","doi":"10.1007/s11242-024-02094-9","DOIUrl":null,"url":null,"abstract":"<h3 data-test=\"abstract-sub-heading\">Abstract</h3><p>Reservoir modeling consists of two key components: the reproduction of the historical performance and the prediction of the future reservoir performance. Industry-standard reservoir simulators must run fast on enormous and possibly unstructured grids while yet guaranteeing a reasonable representation of physical and chemical processes. However, computational demands limit simulators in capturing involved physical and geochemical mechanisms, especially when chemical reactions interfere with reservoir flow. This paper presents a mathematical workflow, implemented in <i>IORSim</i>, that makes it possible to add geochemical calculations to porous media flow simulators without access to the source code of the original host simulator. An industry-standard reservoir simulator calculates velocity fields of the fluid phases (e.g., water, oil, and gas), while IORSim calculates the transport and reaction of geochemical components. Depending on the simulation mode, the geochemical solver estimates updated relative and/or capillary pressure curves to modify the global fluid flow. As one of the key innovations of the coupling mechanism, IORSim uses a sorting algorithm to permute the grid cells along flow directions. Instead of solving an over-dimensionalized global matrix calling a Newton–Raphson solver, the geochemical software tool treats the species balance as a set of local nonlinear problems. Moreover, IORSim applies basis swapping and splay tree techniques to accelerate geochemical computations in complex full-field reservoir models. The presented work introduces the mathematical IORSim concept, verifies the chemical species advection, and demonstrates the IORSim computation efficiency. After validating the geochemical solver against reference software, IORSim is used to investigate the impact of seawater injection on the NCS Ekofisk reservoir chemistry.</p><h3 data-test=\"abstract-sub-heading\">Article Highlights</h3><ul>\n<li>\n<p>The IORSim sorting algorithm decouples the nonlinear geochemical reaction calculations into recurring one-dimensional problems to assure numerical stability and computation efficiency. To the best of our knowledge, this work presents the mathematical concept, implementation, and application of topological sorting for the first time on (industry) field-scale problems.</p>\n</li>\n<li>\n<p>IORSim combines topological sorting with basis swapping and splay trees to significantly reduce computation times. Moreover, a high-speed forward simulation mode was developed to allow the post-advection of chemical components to visualize species distribution, water chemistry, and mineral interactions. If the geochemical reactions interfere with the fluid flow, the IORSim backward mode uses relative permeability curves to update the global fluid flow at each time step.</p>\n</li>\n<li>\n<p>We validate the implemented topological scheme on a reservoir grid, show the computation efficiency, and compare the impact of explicit, implicit, and grid refinement on numerical dispersion.</p>\n</li>\n<li>\n<p>The decoupled flow simulator and geochemical reaction calculations allow seamless integration of full-field reservoir models that contain complex geological structures, a large number of wells, and long production histories.</p>\n</li>\n<li>\n<p>The computation capabilities of IORSim are demonstrated by simulating and reproducing the impact of seawater injection in the southern segment of the giant Ekofisk field (more than 50 years of injection and production history). IORSim shows that seawater injection changed the Ekofisk mineralogy and impacted the produced water chemistry. In the investigated Ekofisk case, seawater promoted calcite dissolution and led to the precipitation of magnesite and anhydrite. Moreover, surface complexation modeling revealed that sulfate is adsorbed on the calcite surface.</p>\n</li>\n</ul>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":null,"pages":null},"PeriodicalIF":2.7000,"publicationDate":"2024-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Transport in Porous Media","FirstCategoryId":"5","ListUrlMain":"https://doi.org/10.1007/s11242-024-02094-9","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"ENGINEERING, CHEMICAL","Score":null,"Total":0}
引用次数: 0
Abstract
Reservoir modeling consists of two key components: the reproduction of the historical performance and the prediction of the future reservoir performance. Industry-standard reservoir simulators must run fast on enormous and possibly unstructured grids while yet guaranteeing a reasonable representation of physical and chemical processes. However, computational demands limit simulators in capturing involved physical and geochemical mechanisms, especially when chemical reactions interfere with reservoir flow. This paper presents a mathematical workflow, implemented in IORSim, that makes it possible to add geochemical calculations to porous media flow simulators without access to the source code of the original host simulator. An industry-standard reservoir simulator calculates velocity fields of the fluid phases (e.g., water, oil, and gas), while IORSim calculates the transport and reaction of geochemical components. Depending on the simulation mode, the geochemical solver estimates updated relative and/or capillary pressure curves to modify the global fluid flow. As one of the key innovations of the coupling mechanism, IORSim uses a sorting algorithm to permute the grid cells along flow directions. Instead of solving an over-dimensionalized global matrix calling a Newton–Raphson solver, the geochemical software tool treats the species balance as a set of local nonlinear problems. Moreover, IORSim applies basis swapping and splay tree techniques to accelerate geochemical computations in complex full-field reservoir models. The presented work introduces the mathematical IORSim concept, verifies the chemical species advection, and demonstrates the IORSim computation efficiency. After validating the geochemical solver against reference software, IORSim is used to investigate the impact of seawater injection on the NCS Ekofisk reservoir chemistry.
Article Highlights
The IORSim sorting algorithm decouples the nonlinear geochemical reaction calculations into recurring one-dimensional problems to assure numerical stability and computation efficiency. To the best of our knowledge, this work presents the mathematical concept, implementation, and application of topological sorting for the first time on (industry) field-scale problems.
IORSim combines topological sorting with basis swapping and splay trees to significantly reduce computation times. Moreover, a high-speed forward simulation mode was developed to allow the post-advection of chemical components to visualize species distribution, water chemistry, and mineral interactions. If the geochemical reactions interfere with the fluid flow, the IORSim backward mode uses relative permeability curves to update the global fluid flow at each time step.
We validate the implemented topological scheme on a reservoir grid, show the computation efficiency, and compare the impact of explicit, implicit, and grid refinement on numerical dispersion.
The decoupled flow simulator and geochemical reaction calculations allow seamless integration of full-field reservoir models that contain complex geological structures, a large number of wells, and long production histories.
The computation capabilities of IORSim are demonstrated by simulating and reproducing the impact of seawater injection in the southern segment of the giant Ekofisk field (more than 50 years of injection and production history). IORSim shows that seawater injection changed the Ekofisk mineralogy and impacted the produced water chemistry. In the investigated Ekofisk case, seawater promoted calcite dissolution and led to the precipitation of magnesite and anhydrite. Moreover, surface complexation modeling revealed that sulfate is adsorbed on the calcite surface.
期刊介绍：
-Publishes original research on physical, chemical, and biological aspects of transport in porous media-
Papers on porous media research may originate in various areas of physics, chemistry, biology, natural or materials science, and engineering (chemical, civil, agricultural, petroleum, environmental, electrical, and mechanical engineering)-
Emphasizes theory, (numerical) modelling, laboratory work, and non-routine applications-
Publishes work of a fundamental nature, of interest to a wide readership, that provides novel insight into porous media processes-
Expanded in 2007 from 12 to 15 issues per year.
Transport in Porous Media publishes original research on physical and chemical aspects of transport phenomena in rigid and deformable porous media. These phenomena, occurring in single and multiphase flow in porous domains, can be governed by extensive quantities such as mass of a fluid phase, mass of component of a phase, momentum, or energy. Moreover, porous medium deformations can be induced by the transport phenomena, by chemical and electro-chemical activities such as swelling, or by external loading through forces and displacements. These porous media phenomena may be studied by researchers from various areas of physics, chemistry, biology, natural or materials science, and engineering (chemical, civil, agricultural, petroleum, environmental, electrical, and mechanical engineering).
京公网安备 11010802042870号
京ICP备2023020795号-1
分享
求助内容：
应助结果提醒方式：
扫码关注我们
