Page tree

Versions Compared

Key

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

@wikipedia 

Motivation


In some specific subsurface applications which require the knowledge of subsurface temperature distributions the assumption of the Constant Areal Geothermal Temperature Profile is not valid and the problem requires a proper 3D modelling solution.


Outputs

LaTeX Math Inline
body--uriencoded--T_G(t, %7B\bf r%7D)

LaTeX Math Inline
body--uriencoded--G_T(%7B\bf r%7D)

Geothermal Temperature Gradient

LaTeX Math Inline
bodyH_n

Neutral Temperature Layer (NTL)


Inputs

LaTeX Math Inline
bodyt

Local Calendar Time

LaTeX Math Inline
body

\delta T_AAnnual average surface temperature variation based on weather reports

--uriencoded-- %7B\bf r%7D

Position vector

LaTeX Math Inline
bodyz(l)

LaTeX Math Inline
bodyA_T

Period of annual temperature variation cycle:

LaTeX Math Inline
body--uriencoded--

A_T = 1 \,

%7B\

rm year%7D LaTeX Math Inlinebody--uriencoded--%7B\

bf j%7D(

%7B\bf r%7D

x,y, z = z_%7Bref%7D)

Earth's Heat Flux at some reference depth 

LaTeX Math Inline
body

\delta t_ATime shift of annual highest temperature with respect to January 1

--uriencoded--z = z_%7Bref%7D
 as function of 
LaTeX Math Inline
body

T_s

Local annual average surface temperature based on weather reports

LaTeX Math Inline
body\delta T_D

Daily average surface temperature variation based on weather reports

(x, \, y)

LaTeX Math Inline
body--uriencoded--

a_%7Ben%7DLocal average Thermal diffusivity of the soil between Earth's surface and NTL

q(%7B\bf r%7D)

Volumetric density of heat sources distributed throughout the subsurface rock volume

LaTeX Math Inline
body

D_TPeriod of daily temperature variation cycle:

T_s(t, x, y)

Surface temperature based on weather reports

LaTeX Math Inline
body--uriencoded--

A_D = 1 \, %7B

\

rm day%7D LaTeX Math Inlinebody--uriencoded--\

lambda_e(%7B\bf r%7D)

Subsurface Thermal Conductivity profile as function of position vector

LaTeX Math Inline
body

\delta t_DTime shift of daily highest temperature with respect to Midnight 00:00 LaTeX Math Inlinebody

--uriencoded--

\delta T_%7B\rm cut%7DTemperature measurement threshold (usually
LaTeX Math Inline
body--uriencoded--\delta T_%7B\rm cut%7D = 0.01 \, %7B\rm °C%7D

a_%7Be%7D(%7B\bf r%7D)

Subsurface Thermal diffusivity profile as function of position vector

where

LaTeX Math Inline
bodyl

Measured Depth of wellbore trajectory with reference to Earth's surface (

LaTeX Math Inline
bodyl=0
)

LaTeX Math Inline
bodyz_s = z(l=0)

TVDss of the Earth's surface in a given location. In case the Earth's surface is at sea level then 

LaTeX Math Inline
bodyz_s = 0


Assumptions




Equations


LaTeX Math Block
anchorT_z
alignmentleft
T_G(t, {\bf r}) = T_{GS}({\bf r})  + T_Y(t, z) + T_D(t, z)
LaTeX Math Block
anchorG_T
alignmentleft
G_T({\bf r}) = \frac{j_z}{\lambda_r({\bf r})}
LaTeX Math Block
anchor1
alignmentleft
\nabla T_{GS} = \lambda^{-1}
\rho_e \, c_e \frac{\partial T_G}{\partial t} + \nabla \left( \lambda_e \nabla T_G \right) = q({\bf r})
\cdot {\bf j}



LaTeX Math Block
anchor
rot_jalignmentleft LaTeX Math Blockanchor
\nabla  \times {\bf j} = 0
1
alignmentleft
 T_
{GS}(
G(t, x, y, z = z_s) = T_s
LaTeX Math Block
anchorT_z
alignmentleft
T_Y
(t,
z) = \delta T_A \, \exp \left[ \, {(z_s-z}) \sqrt{\frac{\pi}{a_{en} \, A_T}} \, \right] \, \cos \left[ \, 2 \pi \frac{t - \delta t_A}{A_T} + (z_s -z) \sqrt {\frac{\pi}{a_{en} \, A_T}} \, \right]
 x, y)




LaTeX Math Block
anchor
T_z
1
alignmentleft
T_D(t,z) = \delta
\Big[ \lambda_e \nabla T_
D
G \
, \exp \left[ \, {(z_s-z}) \sqrt{\frac{\pi}{a_{en} \, D_T}} \, \right] \, \cos \left[ \, 2 \pi \frac{t - \delta t_D}{D_T} + (z_s -z) \sqrt {\frac{\pi}{a_{en} \, D_T}} \, \right]
Big]_{z=z_{ref}} = {\bf j}(x,y, z = z_{ref})



LaTeX Math Block
anchorG_T
alignmentleft
G_T({\bf r}) = \frac{j_z({\bf r})}{\lambda_e({\bf r})
Neutral Layer
LaTeX Math Block
anchorz_N
alignmentleft
z_n = z_s + H_n
LaTeX Math Block
anchorH_N
alignmentleft
H_n = \sqrt{\frac{a_{en} \, A_T }{\pi}} \, \ln \frac{\delta T_A }{\delta T_{\rm cut}
}



See Also

...

Geology / Geothermal Temperature Field

Constant Areal Geothermal Temperature Profile @model ] [ Geothermal Temperature Gradient ]

...