Premium
Revisiting the Analytical Solutions of Heat Transport in Fractured Reservoirs Using a Generalized Multirate Memory Function
Author(s) -
Zhou Quanlin,
Oldenburg Curtis M.,
Rutqvist Jonny
Publication year - 2019
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2018wr024150
Subject(s) - laplace transform , aquifer , mechanics , convective heat transfer , heat transfer , matrix (chemical analysis) , heat exchanger , boundary value problem , flow (mathematics) , geology , mathematics , mathematical analysis , materials science , thermodynamics , physics , geotechnical engineering , groundwater , composite material
Numerous analytical solutions have been developed for modeling thermal perturbations to underground formations caused by deep‐well injection of fluids. Each solution has been derived for a specific boundary value problem and a simplified flow network with one set of parallel fractures. In this paper, new generalized solutions G * ( x , s ) are developed using (existing) global transfer functions G 0 * x s and a new memory function g * ( s ) , where x and s are the space and Laplace variable. The memory function represents the solutions of conductive heat exchange between fractures and matrix blocks and between fractured aquifers and unfractured aquitards. The memory function is developed to account for multirate exchange induced by different shapes, sizes, properties, and volumetric fractions of matrix blocks bounded by multiple sets of orthogonal fractures with different spacing. The global transfer functions represent the fundamental solutions to convective, convective‐conductive, and convective‐dispersive heat transport in fractures (or aquifers) without exchange and are available for various (1‐D linear, 1‐D radial, 2‐D dipole, and single‐well injection‐withdrawal) flow fields. The new solutions with exchange are developed using G * x s = B * s G 0 * x s 1 + ϑg * s, thereby greatly simplifying solution development in a novel way, where ϑ and B * ( s ) are a fracture‐matrix scaling factor and the boundary condition function. The new solutions are applied to several example problems, showing that heat transport in fractured aquifers is significantly impacted by (1) thermal dispersion in fractures that is rarely considered, (2) multirate heat exchange with a wide range of size and anisotropy of rectangular matrix blocks, and (3) heat exchange between aquifers and aquitards.