Professional Documents
Culture Documents
Drainage System
Drainage System
Giovanni Badino
Dip. Fisica Generale – Università di Torino, Associazione La Venta
E-mail: [email protected]
Republished from : Acta Carsologica 2005, 34 (2), 277-316
Abstract
The paper presents an analysis of the interaction between the geothermal flux and the water or air- deep drainage
networks. The problem of geothermal power intercepted by deep structures and, in general, the temperature field
calculations, is converted to classical thermo-engineering problems in terms of shape factors. It is shown that the fluid flow
in a conduit perturbs the whole deep rock temperature field until the geothermal flux of a large area is focalised onto the
conduit. It is shown that either small water masses flowing into a mountain are able to perturb the rock temperature up to
the surface, on sizes that do not depend on water mass dimension, but on its depth, and then on enormous volumes. The
introduction of the “geothermal cross section” of an underground drainage structure allows us to improve the classical
formula of minimum provenance depth of geothermal water. Enlarging factors are applied to the classical estimation in
dependence of the ratio between the actual average discharge and the critical discharge Qc, which depends on the conduit
geothermal cross section. The geothermal “umbra cones” created in the overlying rock by deep underground structures are
described.
It is shown that the geothermal flux can play a significant role in the underground drainage phenomenology.
1. The geothermal flux and the energy The flux is some 60 kW per square kilometre,
contents of rocks which corresponds to a total release of 3×1013
W on the whole planet. For comparison the
An introduction to the geothermal energy flux energy flux received from the Sun is 1.7×1017
W, therefore the geothermal flux is around 5000
It is widely known that the rocks below us times smaller than the main energy source for
have temperatures that increase with depth. The Earth. It cannot play a role in the free
reason is that the internal part of Earth is hot atmosphere phenomenology, but we are going
and the surface cold; there are then two “heat to see that in the case of underground
sources” (but in this work I have adopted the atmospheres it can and it does.
suggestion of Bohren (Bohren, 1998), avoiding
use of the word “heat”), and the thermal energy At first, this appears not to be true in the case
flows between them with the rules given by the of cave atmospheres, that are really quite cold
usual conduction equations. (essentially at the yearly average external
temperature, -from hereafter Tave ; Badino,
Table 1, adapted from (Lee, 1966) gives 2004) while the atmosphere of mines can be
typical values, widely variable, of geothermal very hot (Badino and Forti, 2005). Actually the
flux, estimated by measures of deep first aim of this work was to show that the
underground temperature gradients. geothermal flux could not play any role in the
The world average (Verhoogen, 1956) is deep karst microclimates and genesis, because it
is shielded by deep drainage conduits: This is
Fgt = 0.06 W m −2 exactly the contrary of what we are going to
show...
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.2
TABLE 1
Typical values of geothermal flux
Regions Mean geothermal flux Fgt Number of measures
[mW m-2] or [kW km-2]
Continents
Africa 36-61 15
America 25-150 85
Australia 35-160 65
Asia 22-150 60
Europe 26-140 60
Oceans
Atlantic 3.4-250 250
Indian 5.9-220 250
Pacific 2.9-340 600
Arctic 33-62 20
Geothermal regions
Larderello (I) 450
Oora Kei Korake (NZ) 4200
Matsukawa (J) 630
Mines, caves and tunnels The Gran Sasso rock is essentially limestone,
The first point to discuss is the difference and drillings have shown a cave layer 550 m
between the temperature of artificial and natural above the tunnel altitude, and a general water
cavities in order to clarify the common table extended up to that level. The infiltration
confusion between “thermal flux” and waters at the surface are essentially at Tave, and
“temperature”. cross the whole mountain in nearly adiabatic
conditions, which means that they are only very
It is useful to compare two Underground slowly heated along the fall. Therefore in the Mt
Neutrino Observatories, in Mont Blanc Blanc’s depth there are essentially “mine”
(between France and Italy) and in Gran Sasso waters in thermal equilibrium with hot, deep
(in central Italy), which are assembled in halls rocks, whilst the deep Gran Sasso waters are
in motorway tunnels. The depth of the first, essentially meteoric waters, in equilibrium with
dismissed in 2001, was some 1800 m below the the atmosphere.
surface, at an altitude of 1300 m asl, whilst the
depth of the second is around 1050 m at an It is useful to discuss a little more the internal
altitude of 1000 m asl. Their temperatures are water heating in karst. The reason for water
nevertheless completely different, in the first it temperature increase during underground fall is
is around 32 °C, in the second at 6 °C, the two the gravitational energy which is converted in
unsuitable for working. Therefore it is necessary thermal (a process that gives a water adiabatic
to act in the opposite sense, in the first to cool, lapse rate –2.34 °C km-1) and, in the case of
in the second to warm the experimental halls. flowing in vadose conditions, also to thermal
exchanges with moist air, always characterised
The reason for the two different temperatures by a different adiabatic lapse rate, around -5 °C
is the different rock permeability. km-1. The actual caves’ lapse rates are between
The Mont Blanc rock is mainly granite, the these two extremes (Badino, 2000; Luetscher
waters met by the tunnel were essentially fossil and Jeannin, 2004), a fact that has huge
waters, the meteoric water circulation being consequences on the caves energetic balances,
quite epidermal (up to depth of 100-200 m which nevertheless are outside this work aim;
below the surface), with some exceptions we return briefly to it in the next chapter, but a
localised along large major rock discontinuities. discussion can be found in (Badino, 2005). We
concentrate here on the fact that the infiltrations
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.3
below the surface, has been encountered. The cannot be extremely dangerous. Otherwise, if
Author says that “discrepancies are due to during excavation appears a tendency to a
intensive methane infiltration through the temperature reduction, it is a sure signal of an
Eocene clay shale”, which looks quite strange. approaching water stream in direct contact with
Nevertheless it is quite common to meet high the surface, that at these depths has obviously
temperatures connected with local thermal water enormous pressure, which is able to create
infiltrations, often connected with hyperkarstic extremely dangerous situations, also because it
phenomena (Badino and Forti, 2005). Here we is surely associated to dramatic rock
deal only with a “standard” situation to analyse discontinuities. This is the reason why a lot of
the geothermal impact on our usual karst. work is made in the field of rock temperature
estimations in deep tunnels (Koenigsberg, 1906;
Convective mountains Goy et al., 1996; Badino, unpubl. 2005).
The graphics (Fig. 5) show in short the above Underground high temperatures are
discussed data (temperature vs. depth) where connected with good thermal insulations, which
each type of data point indicator describes a means that a tunnel, or a mine, can be excavated
different tunnel. It is possible to see two across rocks that are very hot (high temperature)
completely different behaviours, the “hot” because, i) they have almost no contact with the
mountains without internal water fluxes, with surface and ii) they have then acquired
positive temperature gradient, and the “karstic” equilibrium with the geothermal flux. Then the
mountains with slowly negative or zero relative rock insulation has allowed to a so
gradients. small thermal flux to heat up to high
temperatures enormous quantities of matter.
As we have seen above, it is possible to
consider mountains with caves as good thermal
conductors and then in general (let us forget for
a while the geothermalism) they are in thermal
equilibrium with the external atmosphere.
As larger are the caves, as smaller it is the
impedance for water crossing the whole
mountain, then deep water circulation is
possible and the whole structure is crossed,
which is a big difference from a mountain
without caves where only water circulation is
Fig. 5. Rock temperature versus depth in the large
quite epidermal.
alpine tunnels. Another consequence to be pointed out is that
if the mountain is highly permeable to water
fluxes the thermal energy transfers inside it are
In fact these mountains are in thermal absolutely dominated by fluid motion, then the
contact, and in equilibrium, with the local pure conductive terms into the rock (“heat”
climate, that they follow with a delay that
diffusion) become negligible.
depends on the mountain depth (Badino, 2004).
It is possible to add that low-impedance water
This thermal contact surface-underground (or air) transfers inside a mountain can be
due to water vein, is the base for the traditional considered a special case of thermal convective
and fundamental “remote sensing” during tunnel movements, and then the karstic mountains are
excavations. Continuous rock temperature examples of convective thermal contact with the
measurements are performed during work atmosphere, whilst the impermeable rocks have
(Guichonnet, 1967); a regular temperature only a conductive contact with it, with a very
increase with the surface distance is a signal of
poor efficiency.
compact rock. Occasional water veins are in
thermal equilibrium with geothermal flux,
The hot mountains energy contents
which is possible only if these water reservoirs
are relatively small and without a hydro- Let us estimate the energy contents of a hot
geological connection on long ranges, and then rock prism, with a surface A and an altitude of
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.5
H. Its temperature at the beginning is T0=Tave, power, but it has worked from a very long time,
local average temperature. the rock is a very efficient thermal insulator and
the final result is an efficient energy storage.
The geothermal flux will heat the rock until
the temperature gradient in it becomes able to
evacuate the geothermal power on the surface. Rock heating time
This equilibrium condition (“stationarity”, i.e. To estimate correctly the time needed for
no time dependence) is obtained when the heating, it is necessary to take into account that
gradient is the rock prism is in contact with others all
around; the calculation would have to include
∂T
Fgt = K R these in the estimation, reconstructing the whole
∂z temperature field and its dependence on time.
The problem is complex but we can do a crude
Where we have used partial derivative estimation of the prism heating time scale
because in general the temperature is function considering it as thermally isolated from its
also of time. This is exactly the condition that surroundings (which is equivalent to assume a
allows the measure of geothermal flux from flat surface and uniform, uni-dimensional
temperature gradient data, assuming a steady- heating), then
state condition.
In this case we can assume a rock thermal Q Fgt Z 2 1
∆t heat = = A C R ρ R =
conductivity KR=2.5 Wm-1 K-1, typical for AFgt K R 2 AFgt
granite, whilst the limestone conductivity is
some 10% less. Then the temperature gradient C ρ Z2
at the equilibrium is
= R R
KR 2
∂T Fgt 0.06
= = = 0.024 K m −1 The last term can be rewritten in terms of
∂z K R 2 .5 thermal diffusivity coefficient defined as
Which is 24 °C per kilometre. It is easy to KR
calculate the energy necessary to heat at this aR =
constant temperature gradient a rock prism of CR ρ R
surface A from depth Z up to the surface; its
total available energy is In the case of rock
For comparison, a 20 kTon atomic bomb Let us recall the classical thermal fields
yields 1014 J. solution of a homogeneous thermal conductor to
a sinusoidal and to a sudden (step) temperature
It is important to note that this energy, that a change.
power plant of 1 GW produces in two years, has
been released to our prism rock by the small In the first case a thermal wave propagates
geothermal heater, only 60 kW. It is a small inside, fading exponentially (Badino, 2004;
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.6
The last equation suggests that a step thermal temperature of the surface, but when its
wave is able to reach a depth x in a time t or, temperature has attained the “stationary
vice versa, after a time t the thermal disturbance temperature field” geothermal gradient seen in
has reached the depth x, and the relation the previous chapter. Nevertheless it makes no
between the two quantities is given by the sense to try to perform exact calculations, which
argument of Error Function. More precisely it is in any case work with unrealistic forms of
possible to show (for example, classically mountain.
(Laidler and Meiser, 1995) that We can conclude that the heating and the
cooling time-scale up to the equilibrium
x 2 ≈ 2aRt configuration are almost the same, and they do
not depend on the temperature drop, then at
Then the penetration in rock up to depth Z of
depth Z in metres we have
a cool wave requires a time-scale ∆tcool
Z2
Z2 ∆t eq = ≈ ∆t cool ≈ ∆t heat [1.1]
∆t cool = 2a R
2 aR
In the case of compact rock we can assume
The estimation is heavily approximated, that the equilibration time scale in years is given
because we are not looking for the complete by
cooling of the mountain, but for the equilibrium
temperature field formation inside it, what is ∆t eq = 0.01 Z 2 [y]
attained not when the whole rock is at the same
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.7
This estimation is very important for our complex, but it is better to postpone a detailed
discussion. In the next chapters we are going to discussion to a future work.
consider water and geothermal fluxes that have We can spend some words about the exact
attained stationary conditions. It is then obvious value of T0 and its relation with Tave. Really Tave
that if the water flow has begun from a time that depends on the altitude, and on average
is much smaller than that of equilibration time, decreases of 6-6.5°C per kilometre, as described
it is not possible to assume that the system has with the International Standard Atmosphere.
attained a general equilibrium configuration, Also at a first approach it looks better to think
what in fact happens in the artificial that the rock assumes not the average yearly
excavations. This is the reason why the mines temperature of the atmosphere Tave, but the
are hot. average temperature of waters at the infiltration
point, that is quite lower because the rain waters
in alpine karst are generally associated with
2. The basic approach to the problem colder periods (but in tropics with warmer
periods); really, many other corrections are
The problem of infiltration temperature necessary to estimate the local rock
temperatures (Badino, 2004).
The classical Desio’s formula (Celico, 1986) Another corrective term, already cited above,
that gives the minimum depth attained by appears during the underground flow because
geothermal water is based on the assumption the temperature increase of underground waters
that underground there exists a first layer along their travel is different from outside,
(heterotermic) in thermal contact with the where the waters follow the ISA mean lapse rate
atmosphere with essentially the temperature of (6 °C per kilometre of fall). Underground, in
local inflowing fluids, which have a temperature adiabatic conditions the theoretical value of
T0 quite precisely (but in general, a little lower, water temperature increase is 2.34°C per
(Badino, 2004)) equal to the local yearly kilometre of fall, but the experimental values
temperature Tave. Below this first layer the rock into the caves, where energy exchanges with the
temperature increases in conductive regime air are possible, are around -3 and -4°C km-1.
(homotermic region). Really, we can assume that the water
The heterotermic layer is also called “active temperature that has infiltrated at altitude H
layer” (US Bureau of Mines, 1996), mainly [km], when it arrives at level 0 is some (3-4)×H
because the seasonal variations can create icing. [°C] hotter, then sensibly colder than the
It is necessary to note, nevertheless, that corresponding infiltrating waters at that altitude.
some problems do exist in the definition of In detailed calculations these effects, that
homotermic layer (Schoeller, 1962; Celico, create a difference between the actual cave
1986), because it is considered the layer where temperature T0 and the local yearly temperature
temperature does not depend on atmospheric at its same altitude Tave, have to be taken into
temperature variations, not the layer where account, but in our discussion they are
temperature is equal to the local average of the completely negligible.
atmosphere. It is then defined on the basis of its
temperature stability, not on the basis of its The energy release to groundwater
thermal contact, if with atmosphere or with deep
rocks. Really in rock with deep aquifers we There is an obvious approach to consider the
meet a heterotermic layer with seasonal role of geothermal flux.
fluctuations (some score of metres), a layer In the upper Earth surface layers the
above the aquifer (included) at T0 (Luetscher, geothermal energy is essentially intercepted by
2004; Badino, 2005), a thin layer of thermal water that releases it to the atmosphere when it
contact deep rocks-aquifer, where a relatively goes out from springs.
sudden temperature increase, probably The energy that comes onto a large surface A
dependent of aquifer permeability, is possible is obviously FgtA, and it is very regular in time.
(Goy, 1996; Badino, unpubl. 2005). We would Let us suppose that it is absorbed by a mass M:
call it the “geothermal exchange layer”. Below Its temperature increase rate is then given by the
this region, a regular temperature increase in condition
deep rocks is found. The scenario is then more
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.8
In alpine karsts P* is some 1000 mm a-1 and assumptions true for similar “discrete” systems?
then the water average temperature increase due The general answer is that no, they are not
to geothermal energy is some 0.5°C. generally correct.
It is a small term, very difficult to measure, The water flows along definite branches, that
and really it has never been measured. But it is do not cover a large surface and with a general
not always so small, a temperature increase of shape that is far from regular. The complete
5°C between the caves and the springs has been shielding assumption is not reasonable for
measured in Cuatro Ciénegas, a figure that karstic drainage.
corresponds to an infiltration (P*-Pout)=100 mm
a-1 in this desert region of Coahuila, Mexico
(AA. VV., 2004). 3. The underground temperature field
with a drainage network
Therefore it would theoretically be possible
to estimate the average infiltrating flow
measuring this temperature increase, but it is a The problem
measure quite difficult to do with some The scenario described above (a regular, flat,
accuracy, because it is the difference of two diffuse water table) can sometimes be correct
uncertain parameters, and other processes surely but in general it gives completely wrong results
occur along the water rock crossing. in the internal rock temperature field
estimations.
The plane watertable
To study a more real model it is necessary to
We have used these calculations in a previous estimate the energy interception made by a
work (Badino, 1995) to explain why the karstic system (a thermodynamical sink) that is buried
mountains are so cold, therefore excluding a in a semi-infinite medium where a thermal
significant role of geothermal flux in energy flux is flowing from infinite.
underground climate definition, because the
Let us consider the problem details.
energy flow from Earth depths is easily
shielded. This very natural approach to refuse, We have a semi-infinite rock volume in
in general, any role to the geothermal flux, it which a thermal flux Fgt is coming from below.
was not original, it is a quite traditional point of It creates a temperature gradient given by:
view.
∂ T Fgt
Bögli (Bogli, 1980) estimates reasonably 0.2 =
°C “to prevent the karstified zone above from ∂ z Kr
being geo-thermally heated up”, that is a very Where we have assumed the depth z as
interesting idea that we shall meet again. In a positive downwards. The temperature field at
very interesting and complete paper Mathey depth H below the surface (or, better, below the
(Mathey, 1974) estimates a maximum of 0.55 heterotermic layer) it is therefore given by
°C. More recently Jeannin et al. (1997) estimate
the specific discharge of a karstic spring (the ∂ T F
equivalent of infiltration) between 30 and 3000 T (H ) = T0 + H = T0 + gt H
mm a-1. The first figure seems too low (almost ∂z Kr
three times dryer than the North Mexico
We have previously seen that T0 is essentially
deserts), but in this paper it is used an energy
the local Tave. It is very natural to suppose that
flux too low by a factor thousand, to finish to
the geothermal energy intercepted by a deep
say that the temperature increase, that has
structure is that given in previous chapters, that
physical dimension [°C], is “less than 0.1 °C/a”,
where really the 30 mm a-1 case would be is (geothermal flux)×(structure area), therefore
warmed by some 15 °C. W = AFgt [3.1]
Let us return to this estimation. We have
obtained the average water temperature increase This means that the energy intercepted by a
and we could now begin to calculate the “cave” is proportional to its surface in the
fluctuations dependence on discharge and so on. direction of energy flux. It is a very natural
But we are dealing with conduits, not with plane assumption, but it is false.
watertables. Are these calculations and
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.10
The fact that, up to now, has not been the stationarity assumption. Whatever initial
considered is this: If we bury a system able to temperature field condition will converge to
intercept and to evacuate elsewhere the energy, asymptotic values which are solution of Laplace
the whole temperature field in the rock is equation, but this convergence requires time.
altered and then the flux itself changes. During this time, which is of the order of
equilibration time scale introduced above [Eq.
The geothermal field with a cave 1.1], the difference between the real field (in
transient phase) and the asympotic one
The problem of temperature field calculation (stationary phase) can be important. If the
in this configuration has to be solved with these equilibration time is comparable with the typical
boundary conditions: changes of boundary conditions (global climatic
1) The geothermal flux from infinite is changes, infiltration of hot waters or so on), the
constant, Fgt; system can never be considered in a true
2) The temperature on the surface is T0 stationary phase, and the equation given by this
constant everywhere; assumption has to be considered heavily
approximated. This limit of stationary
3) The cave temperature is T0, the same as on
assumption gives strong uncertainty in the
the surface (we have sufficiently discussed the
temperature fields estimations for the new deep
limits of this assumption);
alpine tunnels, which have very long
4) At the infinite the temperature field is not equilibration times (Badino, unpubl. 2005), and
disturbed by the cave existence. affects also our next considerations.
These conditions imply a flat external
topography and assume that the infiltrated water
heating from surface to the cave is negligible.
We said that if the cave has an area A, it is
natural to assume that the absorbed geothermal
energy is Fgt×A [Eq 3.1]. Is this correct? Let us
firstly discuss it qualitatively, drawing the
isothermal surfaces in the rock. The figure (Fig.
8) shows a reasonable situation that respects the
boundary conditions. It is possible to see two
things:
1) The isothermal surfaces have a tendency to
converge, then to be focussed, onto the cave;
2) They are “compressed” around the cave.
If we remember that the thermal flux flows
along the maximum T variation (i.e. along the
grad (T), which means perpendicularly to the
isothermal surfaces) and that its value is Fig. 8. Qualitative view of the stationary thermal
field due to the interaction between the geothermal
proportional to the gradient of T, we have that
flux and a conduit with strong drainage (system S).
the two features are equivalent to say that: i) the The isothermal surfaces are affected in a wide
cave focus on itself the geothermal flux and, ii) region, and the cave temperature is T0.
in the rock surrounding the cave the geothermal
flux (and the geothermal gradient) is much more
intense than the natural one...
The general solution
So, the assumption that gave us the [Eq. 3.1]
is surely wrong. Therefore the problem of calculation of
stationary temperature fields is very difficult to
But is it possible to calculate the correct
solve just with the easy boundary conditions
value? Quite surprisingly it seems that nobody
given in the previous chapter. What to expect,
has studied this important problem.
then, when we will have to assume finite energy
Before we look for the solution, we have to transfer rate inside the cave or situations in
make some note about the enormous weight of which the temperature of the cave itself is
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.11
determined by the geothermal flux? And what where H is the cave depth from the surface. We
monstrous form may assume the solution in a then assume that the cave in S” has exactly the
non-stationary case, if we want to consider, for temperature of the rock at depth H in the system
instance, the cooling of a mountain during the S’ plus T0.
karst creation inside it? Let us consider now the three temperature
Nevertheless there is a simple way to lead fields. They are the solution of general Fourier
this problem to typical situations of engineering equation (Isachenko, 1969)
thermal exchange.
1 ∂T
Let us show how, considering three different ∇ ⋅ ∇T ( x, y , z ) =
systems, S, S’ and S”: a ∂t
1) The main system just now described, the That in our case, not time dependent, it
real case, which we call from hereafter S (Fig. reduces to the Laplace equation
8);
2) The system composed only of semi-infinite ∇ 2T = 0
undisturbed rock, without caves and with
It means that the T fields behave like a huge
external temperature equal to 0, which we call
class of phenomena for which the sum of the
from hereafter S’ (Fig. 9);
three spatial second derivatives is zero. The
3) The more complex system S”, (Fig. 10), functions that satisfy these conditions are called
composed by a cave at a particular temperature “Harmonic Functions”, and are among the most
T” buried in a semi-infinite rock, that releases important and studied functions in Physics
energy to the surface at temperature T0. Also the (Carslaw and Jaeger, 1959; Bejan, 1993;
rock at the infinite is assumed to have Balcerzak and Raynor, 1961; Nashchokin,
temperature T0. In this last ideal system S” there 1979). By the way, it would be possible to solve
is no geothermal flux at all. our fields using solutions given for different and
We are then ready to do the final step, well-studied problems like, for instance, the
assuming that T” in S” has a value given by: electric field due to particular charge
distributions, but we can do better for our
Fgt purpose.
T " = T0 + H [3.2]
Kr
Fig. 9. Qualitative view of the stationary thermal Fig. 10. Qualitative view of the stationary thermal
field due to the interaction between the geothermal field between a “hot conduit” at temperature T” and
flux and a conduit with no drainage (system S’). The the surface at T0 (system S”). The isothermal
isothermal surfaces are unaffected, and the cave surfaces are finite, closed and contain the conduit.
temperature is higher than T0. Heat diffusion problems like this one are quite usual
in thermal exchanges engineering.
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.12
TABLE 2
Semi-infinite medium with 4πR T1
isothermal surface and isothermal SF =
spherical cavity of radius R at depth 1 − R / 2H
H
H
T2
R
Semi-infinite medium with 4πR T1
isothermal surface and isothermal SF =
disc of radius R parallel to the π / 2 R − tan −1 (R / 2 H ) T2
surface at depth H H
R
Semi-infinite medium with 2πL T1
isothermal surface and isothermal SF =
cylindrical cavity of length L of cosh −1 (H / R ) T2
radius R, parallel to the surface at
H
R
depth H
L
Semi-infinite medium with 2π H T1
isothermal surface and isothermal SF =
cylindrical hole of radius R drilled to ln (2 H / R ) T2
a depth H normal to the surface.
H
R
Semi-infinite medium with 2π W T1
isothermal surface and isothermal SF =
plate (width W, long L, H>>W) ln (4 L / W ) T2
parallel to the surface at depth H H
W L
If we know the shape factors of two large that each temperature field is not affected
independent systems S1 and S2, say, one by the other. Only in this latter case the shape
composed by two cylinders and the other of a factor of the two conduits is the sum of the
sphere and a cylinder, can we consider a third shape factor of two single conduits, but in
system S3 (in this case, two cylinders and a general it is not so.
sphere) as composed by some “sum” of the two
firsts, and consider that its shape factor is given
4. The interception of geothermal flux by
by some “sum” of the two?
caves
Unfortunately the answer is: No, we cannot.
The temperature fields are linear and then the S3 The geothermal cross sections of caves
temperature field can really be calculated from
the S1 and S2, but it changes completely the It is now possible to apply the previous
equipotential surfaces on which the integration results to the problem of interaction between
is performed to “average” the flux in the caves and geothermal flux. It is stated above
integral [Eq. 3.5]. It is therefore necessary to that, to satisfy the boundary conditions, the T”
recalculate these surfaces and to repeat the temperature has to be
integration that will give a result that has no
Fgt
direct connection with the integrations of S1 and T " = T0 + H
S2 fields. Kr
This means that, for instance, the knowledge
As a consequence, the equation that gives the
of the shape factor of a conduit buried in a semi-
power intercepted
infinite medium tells us almost nothing about
the shape factor of two parallel conduits in the W = K r S F ∆T
same medium, unless their distance would be so
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.15
2π L 1900
SF = = = 250 m
−1 H 7.6
cosh
r
With Fgt=2×106 J m2a-1 we have
TABLE 3
Conduit radius r Shape factor SF Intercepted power Average Flux on
[m] [m] [J a-1] conduit [W m-2]
0.01 165 1.7×1011 900
0.1 210 2.1×1011 110
0.5 250 2.5×1011 25
2 310 3.1×1011 8
5 360 3.6×1011 3.8
We have then that the maximum energy water is minimum and the water flows out quite
extraction efficiency is obtained if the exit hot, at T”.
temperature T is almost equal to T0. Let us define the “critical fluid flux” Qc that
We have the opposite limit, if the warming is divides these two scenarios in a usually
so large that the water temperature T becomes idealised way. We look for a water flux Qc that
T”, the temperature of undisturbed rock. In this enters at temperature T0 and flows out at T” in
case the conduit becomes completely stationary conditions. The enthalpy subtracted
“transparent” to the geothermal flux, the to the system is
temperature field assumes a regular geothermal
gradient not affected by the cave presence and,
dE = C w Qc dt (T "−T0 )
as a consequence, just a little energy is If we admit that the system parameters do not
intercepted. The water into the conduit is then a depend on time, this enthalpy deficit must be
hot “mine water” and the classical, “wrong” given by the incoming geothermal flux Wdt.
solution [Eq. 3.1] becomes correct. Then using [Eq. 4.1] we have
The real cases are intermediate between these
two extremes, because in first approximation the Fgt SF H dt = C w Qc dt (T "−T0 )
water really warms but, as a consequence, its
capability to intercept geothermal energy is But T” is given by
reduced, because the temperature difference
Fgt
between the rock and the water becomes T " = T0 + H
smaller. Really the scenario is still more KR
complex, because the water temperature
increases along its path underground: It enters And we obtain
cool, very efficient in geothermal energy
[kg s−1 ]
dM K R SF [4.2]
focussing, but as long as it warms downstream Qc = =
dt Cw
its capability to intercept the geothermal flux
decrease. The non-linearity of SF forbids correct Then the critical flux is simply the conduit
analytical solutions, but we can make some shape factor “scaled” by a dimensional term
other step. (rock conductivity divided by the thermal
capacity of flowing fluid). In terms of volume
The critical shielding discharge flux
It is possible to estimate the water heating
[m 3s − 1 ]
dV Qc K R SF
along a deep conduit and its final temperature as = =
a function of conduit parameters. dt ρ w C w ρ w
We have just seen that two extreme scenarios From another point of view, Qc can be
are possible. If the fluid flux is very large the considered the critical flux below which it is
rock temperature field is completely changed, possible to consider that the rock temperature
the geothermal energy flow interception is field is undisturbed. Or, from still another point
maximum and the water flows out at T0. At the of view, we can be sure that a water flux much
other extreme, if the water flow is very small, larger than Qc perturbs the rock.
the rock temperature field is completely
undisturbed, the geothermal energy flow to the
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.17
The solution [Eq. 4.2], which is surprisingly We are going to improve this point of view
simple in comparison with the bloody analytical and discover that Qc has another, still more
difficulties of the problem, is able to give also important, meaning.
the velocity that water must have to satisfy the
equation. If the conduit radius is r, then Geothermal power intercepted
1 K S 1 We can calculate the effective cave
v = Qc = R F temperature T at the equilibrium and solve the
π r ρw
2
Cw π r 2 ρ w
inverse problem, the estimation of flowing
In the case of a water draining conduit in depth of hot spring waters.
limestone we have The final system temperature must lie
between T0 (near it for high water discharges)
2.3
Qc,water = SF = 5.5 ×10 −4 SF [kg s−1 ] and T” (near it for low discharges). Let us call
4.2 ×10 3 [4.3] WM the maximum power that it can be
subtracted by our system
If the conduit drains air, the critical flux
becomes W M = Fgt SF H = 0.06 SF H [W]
[4.4]
2.3
S = 2.3×10 −4 SF [kg s−1 ]
Which is valuable for very large water flows
Qc,air =
10 3 F and outflowing temperature around T0.
If we call T the real (unknown) conduit
The shape factor is in general comparable
temperature at the equilibrium, the residual
with the conduit length, usually some 103-104
outgoing upward flux is not zero, because it
m, therefore the water flux able to create
“filters” an energy given by
changes in the rock temperature fields is in
general quite small. Neither the air flux Wup = K R S F (T − T0 )
requested to cool the rock it is too large, in
absolute, but its small thermal capacity and its The difference between WM and Wup is the
very small density cause a large volume flow net flux entering inside the cave from below. If
request. the conduit is at temperature T, the energy
It is nevertheless necessary to remember that conservation states
temperature field changes are only possible if
the fluid fluxes have had sufficient time (more Win = WM − Wup = Fgt S F H − K R S F (T − T0 )
than ∆teq defined in [Eq. 1.1]) to converge at the
stationary (equilibrium) state. But from [Eq. 3.2] we have
The critical flux Qc for air or water are Fgt H
extremely important for another reason: They KR =
are the air or water fluxes able to shadow the T "−T0
upper rock from the geothermal flux, forcing its
Then
temperature near to the average yearly
temperature of external atmosphere, Tave. T − T0 T "−T
Finally, it is important to note that this work Win = Fgt SF H 1− = WM
T "−T0 T "−T0 [4.5]
was originally a chapter of an underground
climate physics book. Its purpose was to
If T=T0 the equation reduces to Win=WM
estimate when we could consider a cave as
(which describes the “system S” situation), and
“shielded” from the upward geothermal energy
if T=T” the term Win vanishes, as expected.
flux. If the flux is much larger than Qc we are
sure that the mountain parts (and caves) above Therefore, the geothermal flux intercepted by
the conduit are shielded; if the flux is well the cave is reduced as long as its final
below the critical value the conduit is temperature T increases: The cave is becoming
“transparent” to the geothermal energy and we “transparent”, and this equation describes its
have to include also its contribution to analyse “fading” inside the temperature field.
the underground climate in the rock above the
conduit.
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.18
Temperature and deposited power versus the ratio between this difference and the
discharge theoretical, maximum difference T”-T0. Then
Usually we have very small possibilities to T "+qT0 1
establish a natural conduit discharge. We deal T − T0 = − T0 = (T "−T0 )
with actual discharge Q, and we can only try to 1+ q 1+ q
understand how this water (or air) flow is
And calling the “excess temperature ratio” of
affected by the geothermal flux.
groundwater, that is the amount of actual
Let us then return to the application of the heating in comparison with the maximum
First Principle [Eq. 4.5]; if we call Q the total attainable, we have
actual water flux in the system, it gives T − T0 1
rT = = [4.10]
T "−T T"−T0 1+ q
WM = C w Q (T − T0 ) [4.6]
T "−T0 With the assumption [Eq. 4.9] we can rewrite
Where we have assumed that the whole the [Eq. 4.5]
entering energy flux goes to warm the water
flux. This statement is true if we admit that the T "−T q
Win = WM = WM [4.11]
thermal energy is released on the whole system, T "−T0 1 + q
i.e. it is true if the temperature T, which is
actually the output water temperature, can be We have then two fundamental equations, [Eq.
used to describe the whole system, also in its 4.10] and [Eq. 4.11], which connect the internal
further parts. It is a strong, and in general false, drainage Q to the outgoing temperature and to
hypothesis, but it is better to assume it as true the intercepted geothermal flux. The two
and only afterwards have a look on what graphics (Fig. 12) describe the behaviour of the
happens in more real situations. out-flowing water temperature T and of water
With the previous results and trivial absorbed energy as a function of discharge,
calculations, we have obviously measured in function of our nice
scale-discharge Qc (it can be adapted to air flow
WM T "−T with trivial changes).
Q= =
C w (T − T0 ) T "−T0 We have previously discussed the Qc as the
[4.7] “shielding flux” and WM as the “maximal
WM T "−T T "−T intercepted flux”. Now we see that they are
= = Qc
C w (T "−T0 ) T − T0 T − T0 mainly the natural scales of fluid flow and of
geothermal power flux intercepted, exactly as
happens with the Similarity Numbers, always a
This important equation relates the actual flux ratio between a parameter and a scale-value for
Q to the output temperature T in terms of the it.
critical shielding flux Qc and the surrounding
temperatures. It is very easy to solve it to obtain
Q
T "+T0
Qc T "+ qT0
T= = [4.8]
1+
Q 1+ q
Qc
Q
q= [4.9]
Qc
Fig. 12. Water temperature increase and intercepted
The excess temperature above the “external power by a deep conduit versus water discharge.
average” T0 is then T-T0, but its natural scale is
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.19
In general a typical deep circuit has a very It is possible to perform a last step,
large size, many kilometres. If we call L* its considering a long conduit L, along which the
length expressed in kilometres we have, for the water is heated. The term T is now the
circuit described above temperature at the length x, in a section dx with
shape factor sF. The thermal flow and the
Qc ≈ 0.5 L * temperature increase in that section is then
given by [Eq. 4.5]
And finally
dW = dWM − dWup =
2Q
H ≈ H 0 1 + [4.16] Fgt H
L* = Fgt s F H − s F (T − T0 ) = C w Q dT
T "−T0
It is a simple formula that estimates the
drainage depth as a function of discharge. And then
As an example, let us consider a spring with Fgt S F H T "−T sF
an average discharge Q=20 kg s-1, and a = dT
temperature 20°C above the local average. The CwQ T "−T0 S F
depth estimation [4.12] gives H0=830 m.
This is an equation that could solve the
If the hydro-geological context suggests a problem, if we would able to integrate the left
circuit length L*=20 km, the critical discharge part, but unfortunately this is not possible.
is Qc=11 kg s-1, and then the q term is almost
equal to 2. Our formula estimates then a depth We have written sF and not dSF (as would be
H of 2.5 km for the circuit, much more than natural) because it is not possible to pass from
supposed… the equation that gives SF as a function of L, to
the contribution of a part dx of L to SF. At the
With this value we can return to the end of the third chapter, we have noted that SF
geological context and, if we have other does not linearly depend on each dx part,
information (like the circulation time), to more because it comes from an average on the whole
accurate estimations of the ratio between the space and system, and it is not possible to
depth and the radius of a deep drainage system. consider it as the result of an integration on
A better estimation of provenance depth with some dx.
[Eq. 4.13] and [Eq. 4.14] is then possible.
For instance, the contribution of the dx at one
It is nevertheless better to remember that we conduit edge covers a cone above it, and the
are working with the assumption that the system local sF is like that given by a small sphere,
has attained stationary condition; the above whereas the dx in the middle of conduit gives a
formula is then correct for water fluxes that very small contribution.
persist from very long times, much longer than
We can nevertheless integrate the last
∆teq of [Eq. 1.1].
equation assuming the false approximations
sF=dSF. It makes no analytical sense but it
Temperature changes into the system probably introduces an error smaller than the
It is possible to perform a last “calibration”. assumption of uniform system temperature.
We have already noted that we are assuming Then for a conduit buried in a semi-infinite
that it is possible to define univocally a system medium at depth H
temperature T, but this is not always true. It is
Fgt S F H T "−T dx
possible to speak univocally of “system = dT
temperature” if each part of the system is CwQ T "−T0 L
uniformly heated by the energy flux, for
instance if the head and the tail of water flux are Using Qc becomes
mixed (for instance, when the water enters in a
spherical deposit). Nevertheless this is not the 1 Fgt H T"−T dx Qc dx
= (T "−T ) = dT
usual situation because in a real conduit the q K r T "−T0 L Q L
water enters with a temperature T0 and flows
warming up to the final temperature T. And integrating on x from 0 to L and on T
from T0 to T we have
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.21
During deep flow the underground waters are carbonate dissolution releases ions that are in
warmed from their starting temperature T0 up to part the same already present in water.
the final T, that has a theoretical maximum at This complex system forces to find the
T”, as shown above. solution of many different equations describing
What is the typical temperature increase? We equilibrium kinetics, charge and mass
have seen that a conduit is able to shield the conservation. General solution charts are given
geothermal flux like a plane watertable does, in (Ford and Williams, 1989); they show the
then the scale temperature increase is around the saturation values at the equilibrium for various
value given by [Eq. 2.1] initial partial carbon dioxide pressures. It results
that in open systems (with release of carbon
0.06 500
∆Tgt = = [°C] dioxide excesses) the warming of a calcite
4.2 ×10 (P − Pout ) P *
3
saturated water gives, without exceptions, a
super-saturation and then provokes a calcite, or
Therefore, in real cases it is in the range aragonite, deposition.
between 0.2 and 3 °C, a temperature drop that
In a closed system this behaviour changes in
water gains during flow between the cave
a complex way. Generally a super-saturation is
bottom and the springs. Now we know that this
produced, but if the initial CO2 partial pressure
temperature change happens also along
is below 200 Pa (0.002 atm) and the temperature
conduits, not only in the “watertable”, and that
is below 30 °C, a calcite under-saturation
the power release is concentrated on the conduit
appears as result of water heating, as larger as
surface walls. How does this warming affect the
colder is the water.
water chemistry?
The typical carbon dioxide partial pressure in
We can outline the chemical behaviour of
free atmosphere is 3.5×10-4 atm, then at 10 °C
water saturated of calcium carbonate entering in
the calcite equivalent content at the saturation is
flooded conduits.
around 12 mg l-1. A water temperature increase
It is well known that the carbonate of 1 °C result in a saturation value of 0.02-0.04
dissolution in water is very complex (Snoeyink, mg larger, i.e. with a flux of 1 m3s-1 it gives
1980), (Ford and Williams, 1989), (Dreybrodt, around 103 kg of dissolved rock per year.
2000) because the equilibrium state results from
It is a small figure but it suggests that further
the combined equilibriums of different,
studies are necessary to a more complete
interconnected reactions, which depends on
understanding of saturation conditions as a
temperature, local pH and the presence of other
function of temperature, of chemically complex
dissolved salts with common ions.
waters in a closed system.
In the simplest case, the first equilibrium
In any case indirect evidences of
reaction gives the amount of dissolved carbon
effectiveness of speleogenetic processes
dioxide, for which in usual conditions the Henry
induced by geothermal heating in phreatic
Law holds, stating that the dissolved gas
conduits can be found, because if these
decreases with temperature and it is
processes are possible, they have to affect the
proportional to its partial pressure above the
network morphologies:
water surface. Therefore, its quantity depends
also on the kinetics of gas transport until the 1) the geothermal energy is released only in
surface, if it does exist. the lowest conduit walls, then the dissolving
characteristics have to depend on the rock
The other reactions, which involve only water
surface orientation;
and carbon dioxide, are the dissociation of
carbonic acid in calcium bicarbonate and H+, the 2) a deep conduit shadows completely the
dissociation of bicarbonate and the equilibrium upper rock, then the formation of a conduit that
H+ and OH- in water. These dissociations tend cross the rock above another is hampered, and
to increase with temperature thanks to the this affects the whole drainage conduit structure.
increase of available energy. Similar processes can probably play a part
The last main reaction describes the also in the deep drainage network formation in
equilibrium between the calcium carbonate and glaciers (Badino, 2002), but either ice or
water enriched with carbon dioxide. The limestone, a lot of work has still to be done for a
better understanding of geothermal role in karst.
Giovanni Badino / Speleogenesis and Evolution of Karst Aquifers, 2005, 3 (2), p.24
Tiefenstufe durch Berge und Taler. Eclog. Ozisik M. and Necati. 1993. Heat
Geol. Helv. IX (1). Conduction,.Wiley InterScience.
Laidler K. and Meiser J. 1995 Physical Schoeller H. 1962. Les eaux souterraines.
Chemistry. Houghton-Mifflin Co. Masson.
Lismonde B. 2002. Aérologie des Systèmes Snoeyink V.and Jenkins D. 1980. Water
Karstiques. CDS Isére. Chemistry. John Wiley & Sons.
Luetscher M. and Jeannin P.-Y. 2004. Szechy K. 1973. The Art of Tunnelling,
Temperature distribution in karst systems: Akademiai Kiado, Budapest.
the role of air and water fluxes. U.S. Bureau of Mines. 1996. Dictionary of
Speleogenesis and Evolution of Karst Mining, Mineral, and Related Terms, CD-
Aquifers 2 (2), (from Terra Nova 16):, 344– ROM.
350
Verhoogen J. 1956. Temperatures within the
Nashchokin V. 1979. Engineering Earth. In Physics and Chemistry of the
Thermodynamics and Heat Transfer. Mir, Earth, Pergamon Press, I.
573 pp.