Algebraic Dynamic Multilevel Method for Fractured Geothermal Reservoir Simulation
- Mousa HosseiniMehr (Delft University of Technology) | Rhadityo Arbarim (Delft University of Technology) | Matteo Cusini (Delft University of Technology) | Cornelis Vuik (Delft University of Technology) | Hadi Hajibeygi (Delft University of Technology)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Simulation Conference, 10-11 April, Galveston, Texas, USA
- Publication Date
- Document Type
- Conference Paper
- 2019. Society of Petroleum Engineers
- 5.9.2 Geothermal Resources, 5.9 Non-Traditional Resources, 5.5 Reservoir Simulation, 0.2 Wellbore Design, 1.2.3 Rock properties, 5 Reservoir Desciption & Dynamics
- Physics-Based Non-linear Solvers, Fractured geothermal reservoirs, Scalable Linear and Nonlinear Solvers, Embedded Discrete Fracture Model, Multilevel Multiscale Method
- 3 in the last 30 days
- 115 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 9.50|
|SPE Non-Member Price:||USD 28.00|
A dynamic multilevel method for fully-coupled simulation of flow and heat transfer in heterogeneous and fractured geothermal reservoirs is presented (FG-ADM). The FG-ADM develops an advanced simulation method which maintains its efficiency when scaled up to field-scale applications, at the same time, it remains accurate in presence of complex fluid physics and heterogeneous rock properties. The embedded discrete fracture model is employed to accurately represent fractures without the necessity of unstructured complex grids. On the fine-scale system, FG-ADM introduces a multi-resolution nested dynamic grid, based on the dynamic time-dependent solution of the heat and mass transport equations. The fully-coupled implicit simulation strategy, in addition to the multilevel multiscale framework, makes FG-ADM to be stable and efficient in presence of strong flow-heat coupling terms. Furthermore, its finite-volume formulation preserves local conservation for both mass and heat fluxes. Multi-level local basis functions for pressure and temperature are introduced, in order to accurately represent the heterogeneous fractured rocks. These basis functions are constructed at the beginning of the simulation, and are reused during the entire dynamic time-dependent simulation. For several heterogeneous test cases with complex fracture networks we show that, by employing only a fraction of the fine-scale grid cells, FG-ADM can accurately represent the complex flow-heat solutions in the fractured subsurface formations.
|File Size||2 MB||Number of Pages||16|
Burnell, J., OSullivan, M., OSullivan, J., Kissling, W., Croucher, A., Pogacnik, J., Pearson, S., Caldwell, G., Ellis, S., Zarrouk, S., and Climo, M. (2015). Geothermal supermodels: the next generation of integrated geophysical, chemical and flow simulation modelling tools. In In Proceedings World Geothermal Congress, pages 19-21.
Cusini, M., Gielisse, R., Groot, H., van Kruijsdijk, C., and Hajibeygi, H. (2018b). Incomplete mixing in porous media: Todd-longstaff upscaling approach versus a dynamic local grid refinement method. Computational Geosciences, 10.1007/s10596-018-9802-0.
Dogru, A. H., Fung, L. S. K., Middya, U., Al-Shaalan, T., and Pita, J. A. (2009). A next-generation parallel reservoir simulator for giant reservoirs. 2009 SPE Reservoir Simulation Symposium, The Woodlands, Texas. 10.2118/119272-MS.
Klemetsdal, O., Moyner, O., and Lie, K.-A. (2018). Use of dynamically adapted basis functions to accelerate multiscale simulation of complex geomodels. ECMOR XVI - 16th European Conference on the Mathematics of Oil Recovery. http://dpi.org/10.3997/2214-4609.201802251.
Singh, G., Leung, W., and Wheeler, M. (2018). Multiscale methods for model order reduction of non-linear multiphase flow problems. Comput Geosci. 10.1007/s10596-018-9798-5.
Wallis, J. R., Kendall, R. P., Little, T. E., and Nolen, J. S. (1985). Constrained residual acceleration of conjugate residual methods. SPE Reservoir Simulation Symposium, 10.2118/13536-MS.