Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.


Motivation

The Temperature Flat Source Solution @model is one of the fundamental solutions of temperature diffusion equations modelling the temperature conduction in linear direction (see Fig. 1).

...


Subsurface Temperature Profile around Lateral Flow makes adjustments to Geothermal Temperature Profile 

LaTeX Math Inline
bodyT_G(z)
 to account for the lateral reservoir flow with a constant temperature (see Fig. 1 and Fig. 2).


Image Added

Image Added

Fig. 1. Sample Subsurface Temperature Profile around a 

LaTeX Math Inline
bodyh_f
 height lateral flow at depth 
LaTeX Math Inline
body

...

z_f

...

with temperature 

LaTeX Math Inline
bodyT_

...

Image Removed

f

Fig.

1
around Lateral Flow

 around two lateral flows with temperature 

LaTeX Math Inline
body--uriencoded--T_%7Bf1%7D
 and 
LaTeX Math Inline
body--uriencoded--T_%7Bf2%7D


Outputs

...

LaTeX Math Inline
bodyT_be(t, z)

Temperature Subsurface temperature distribution


Inputs

...

LaTeX Math Inline
bodyt

Time lapse after the temperature step from 

LaTeX Math Inline
bodyT_

b

e(z=0) =0
  up to 
LaTeX Math Inline
bodyT_

b

e(z=0) =T_f

LaTeX Math Inline
bodyz

Spatial coordinate along the transversal direction to constant temperature 

LaTeX Math Inline
bodyT_

b

e(z)= T_f
plane 
LaTeX Math Inline
bodyz=0

LaTeX Math Inline
bodyz_f

TVDss of the top of the lateral flow unit

LaTeX Math Inline
bodyh_f

True vertical thickness of the the lateral flow unit

LaTeX Math Inline
bodyT_f

Boundary temperature at 

LaTeX Math Inline
bodyz=0

LaTeX Math Inline
bodya_e

Thermal diffusivity of the surroundings

LaTeX Math Inline
bodyT_G(z)

Geothermal Temperature Profile


Equations

...

Driving equationInitial conditions Boundary conditions


LaTeX Math Block
anchor1
alignmentleft
\frac{\partial T_be}{\partial t} = a^2 a_e^2 \, \Delta T_be = a^2a_e^2 \, \frac{\partial^2 T_be}{\partial z^2}



LaTeX Math Block
anchor1
alignmentleft
T_be(t=0, z) = T_G(z)



LaTeX Math Block
anchor7685E
alignmentleft
T_be(t, z=0 z_f \leq z \leq z_f + h_f) = T_f = {\rm const}


LaTeX Math Block
anchor7685E
alignmentleft
T_be(t, z \rightarrow \infty) = T_G(z)


...


LaTeX Math Block
anchor1
alignmentleft
\mbox{if} \, z < z_f \; \Longrightarrow \;T_be(t,z) = T_f + (T_G(z) - T_f) \cdot \mbox{erf} \left( \frac{2z_f-z}{\sqrt{\pi}} \int_0^{z/\sqrt{4at}} e^{-\xi^2} d\xi4 a_e t}} \right)



LaTeX Math Block
anchor1
alignmentleft
\mbox{if} \, z_f \leq z \leq z_f + h_f  \; \Longrightarrow \; T_e(t,z) = T_f



LaTeX Math Block
anchor1
alignmentleft
\mbox{if} \, z > z_f + h_f  \; \Longrightarrow \; T_e(t,z) = T_f + (T_G(z) - T_f) \cdot \mbox{erf} \left( \frac{z-z_f-h_f}{\sqrt{4 a_e t}} \right)


where

LaTeX Math Inline
body--uriencoded--\mbox%7Berf%7D(x)

Error function


See Also

...

Geology / Geothermal Temperature Field / Geothermal Temperature Profile

...