Upscaling of Modeling of Thermal Dispersion in Stratified Geothermal Formations
J. Tang (United Arab Emirates University)
Yang Wang (Qingdao University of Technology)
WR Rossen (TU Delft - Atmospheric Remote Sensing)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
Upscaling of geothermal properties is necessary given the computational cost of numerical simulations. Nevertheless, accurate upscaling of thermo-physical properties of layers combined in simulation grid blocks has been a long-standing challenge. In stratified porous media, non-uniform velocity between layers combined with transverse thermal conduction across layers causes spreading of the thermal front: thermal Taylor dispersion. Neither effect of heterogeneity is accounted for in conventional upscaling. Based on thermal Taylor dispersion, we develop a new upscaling technique for simulation of geothermal processes in stratified formations. In particular, we derive a model for effective longitudinal thermal diffusivity in the direction of flow, αeff, to represent this phenomenon in two-layer media. αeff, accounting for differences in velocity and transverse thermal conduction, is much greater than the thermal diffusivity of the rock itself, leading to a remarkably larger effective dispersion. We define a dimensionless number, NTC, a ratio of times for longitudinal convection to transverse conduction, as an indicator transverse thermal equilibration of the system during cold-water injection. Both NTC and αeff equations are verified by a match to numerical solutions for convection/conduction in two-layer systems. We find that for NTC > 5, thermal dispersion in the system behaves as a single layer with αeff This suggests a two-layer medium satisfying NTC > 5 can be combined into a single layer with an effective longitudinal thermal diffusivity αeff. Compared with conventional approaches by averaging, the αeff model provides more accurate upscaling of thermal diffusivity and thus more-accurate prediction of cooling-front breakthrough. In stratified geothermal reservoirs with a sequence of layers, upscaling can be conducted in stages, e.g. combining two layers satisfying the NTC criterion in each stage. The application of the new technique to upscaling geothermal well-log data will be presented in a companion paper.