EGU General Assembly 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Poroelastic aspects in geothermics

Bianca Kretz1, Willi Freeden2, and Volker Michel1
Bianca Kretz et al.
  • 1Geomathematics Group, University of Siegen, Germany
  • 2University of Kaiserslautern, Germany

The aspect of poroelasticity is anywhere interesting where a solid material and a fluid come into play and have an effect on each other. This is the case in many applications and we want to focus on geothermics. It is useful to consider this aspect since the replacement of the water in the reservoir below the Earth's surface has an effect on the sorrounding material and vice versa. The underlying physical processes can be described by partial differential equations, called the quasistatic equations of poroelasticity (QEP). From a mathematical point of view, we have a set of three (for two space and one time dimension) partial differential equations with the unknowns u (displacement) and p (pore pressure) depending on the space and the time.

Our aim is to do a decomposition of the data given for u and p in order that we can see underlying structures in the different decomposition scales that cannot be seen in the whole data.
For this process, we need the fundamental solution tensor of the QEP (cf. [1],[5]).
That means we assume that we have given data for u and p (they can be obtained for example by a method of fundamental solutions, cf. [1]) and want to investigate a post-processing method to these data. Here we follow the basic approaches for the Laplace-, Helmholtz- and d'Alembert equation (cf. [2],[4],[6]) and the  Cauchy-Navier equation as a tensor-valued ansatz (cf. [3]). That means we want to modify our elements of the fundamental solution tensor in such a way that we smooth the singularity concerning a parameter set τ=(τxt). 
With the help of these modified functions, we construct scaling functions which have to fulfil the properties of an approximate identity.
They are convolved with the given data to extract more details of u and p.


[1] M. Augustin: A method of fundamental solutions in poroelasticity to model the stress field in geothermal reservoirs, PhD Thesis, University of Kaiserslautern, 2015, Birkhäuser, New York, 2015.
[2] C. Blick, Multiscale potential methods in geothermal research: decorrelation reflected post-processing and locally based inversion, PhD Thesis, Geomathematics Group, Department of Mathematics, University of Kaiserslautern, 2015.
[3] C. Blick, S. Eberle, Multiscale density decorrelation by Cauchy-Navier wavelets, Int. J. Geomath. 10, 2019, article 24.
[4] C. Blick, W. Freeden, H. Nutz: Feature extraction of geological signatures by multiscale gravimetry. Int. J. Geomath. 8: 57-83, 2017.
[5] A.H.D. Cheng and E. Detournay: On singular integral equations and fundamental solutions of poroelasticity. Int. J. Solid. Struct. 35, 4521-4555, 1998.
[6] W. Freeden, C. Blick: Signal decorrelation by means of multiscale methods, World of Mining, 65(5):304-317, 2013.

How to cite: Kretz, B., Freeden, W., and Michel, V.: Poroelastic aspects in geothermics, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-7100,, 2020

Comments on the presentation

AC: Author Comment | CC: Community Comment | Report abuse

Presentation version 1 – uploaded on 30 Apr 2020 , no comments