Three-dimensional modeling of heat transport in deep hydrothermal reservoirs

Abstract

The need for alternative energy increases steadily-especially due to the climate change and the limited availability of fossil fuels. Geothermal power uses the intrinsic heat which is stored in the accessible part of the Earth's crust. Its importance among the renewable energy resources originates from the almost unlimited energy supply of the Earth and its independence from external influences such as seasonal or even daily climatic variability. Nevertheless, there are risks which have to be assessed. From a mathematical point of view there are four building blocks of the characterization of deep geothermal systems: seismic exploration, gravimetry, modeling transport processes, and modeling the stress field. In particular, local depletion poses a significant risk during the industrial utilization of geothermal reservoirs. In order to reduce this risk, reliable techniques to predict the heat transport and the production temperature are required. To this end, a three-dim ensional model to simulate the heat transport in hydrothermal systems is developed which is based on a transient advection-diffusion-equation for a two-phase porous medium. The existence, uniqueness, and continuity of the weak solution of the resulting initial boundary value problem is verified. For the numerical realization, a Galerkin scheme is introduced on the basis of scalar kernels. Its convergence is proven. Furthermore, exemplary applications of this method are investigated for the biharmonic kernel

Mathematische Methoden in der Geothermie

“…Die damit verbundene Theorie und numerischen Lösungsverfahren sind z.B. in [5,16,18,23,44,68,69,76] zu finden.…”

“…Dabei entspricht das lokal kompakte Wavelet ψ τ j der Differenz der Skalierungsfunktionen φ τ j +1 und φ τ j , d.h. Im Folgenden werden die für hydrothermale Systeme in [68,69] und für petrothermale Systeme in [59] vorgestellten Methoden zur Simulation des Wärmetransports näher betrachtet.…”

“…[15,85]). Sie erfüllen die nötigen Bedingungen des Galerkin-Verfahrens (siehe [68,69] für Details). Die Ergebnisse eines numerischen Tests für den Kubus (−1, 1) 3 im Falle eines diffusionsdominanten Problems (mit Péclet-Zahl Pe = 0.01), einer konstanten Anfangsbedingung von T 0 = 393.15 K und einer (Punkt-)Injektion im Ursprung mit einer Injektionstemperatur, die (zeitlich) konstant 50K unter der Anfangsbedingung liegt, sind in Abb.…”

“…In dieser Arbeit werden -basierend auf den Betrachtungen in [3,28,30,33,35,37,45,46,59,68,69] -zuerst kurz der geothermische Hintergrund und insbesondere die relevanten Reservoirtypen dargestellt. Im Anschluss werden für diese Reservoirtypen die den einzelnen Säulen zugrundeliegenden Differential-und Integralgleichungen beschrieben.…”

“…1) heraus. Modelle und (numerische) Lösungsansätze zu den beschriebenen Problemen im Bereich der Transportprozesse findet man unter anderem in[21,43,46,53,54,59,68,69,71] und der darin aufgeführten Literatur.…”