Synonym: Geothermal Temperature Profile @model = Constant Areal Geothermal Temperature Profile @model
Motivation
...
In many subsurface applications which require the knowledge of subsurface temperature distributions the land area of the study is small enough to consider the subsurface heat flux
...
...
along-hole geothermal temperature profile
--uriencoded--%7B \bf j%7D(x,y,z) = \%7B j_x, \, j_y, \, j_z \%7D |
|
and Thermal Conductivity LaTeX Math Inline |
---|
body | --uriencoded--\lambda_e(%7B\bf r%7D) |
---|
|
to be homogeneous across location area:
LaTeX Math Block |
---|
| {\bf j}({\bf r}) ={\bf j}(x,y,z)={ \bf j}(z) |
| LaTeX Math Block |
---|
| \lambda_e({\bf r}) =\lambda_e(z) |
|
where
is true vertical depth.Since the heat flux is conservative (see LaTeX Math Block Reference |
---|
anchor | rot_j |
---|
page | Geothermal Temperature Field @model |
---|
|
) then it immediately implies that:
LaTeX Math Block |
---|
|
{\bf j}=\{ j_x = {\rm const}, \, j_y = {\rm const} , \, j_z(z) \} |
Further admitting that a surface temperature over the study area is constant:
one can see that lateral components of the heat flux are vanishing: LaTeX Math Block |
---|
|
{ \bf j}(x,y,z) = \{ j_x = 0, \, j_y = 0 , \, j_z(z) \} |
Normally there are no heat sources within a subsurface volume under study other than upward Earth's Heat Flux which means that true vertical component
LaTeX Math Inline |
---|
body | j_z(z) = j_z = \rm const |
---|
|
is constant along true vertical direction. It varies across the Earth but local value is usually well known.This simplifies the procedure of modelling a Geothermal Temperature Field
LaTeX Math Inline |
---|
body | --uriencoded--%7B \bf j%7D(x,y,z) = \%7B 0, \, 0 , \, j_z \%7D |
---|
|
to modelling of a Constant Areal Geothermal Temperature Profile along a given wellbore trajectory.Outputs
Inputs
| Local Calendar Time | | Annual average surface temperature variation based on weather reports |
| | | Period of annual temperature variation cycle: LaTeX Math Inline |
---|
body | --uriencoded--A_T = 1 \, %7B\rm year%7D |
---|
|
|
| True vertical component of regional Earth's Heat Flux | | Time shift of annual highest temperature with respect to January 1 |
| Local annual average surface temperature |
t | Astronomic time | | Annual average surface temperature variation l | wellbore trajectory Measured Depth with reference to Earth's surface ( | Daily average surface temperature variation based on weather reports |
|
l=0 | )A annual daily temperature variation cycle: |
T year%7D(lWellbore trajectory. True Vertical Depth Sub-Sea (TVDss)A moment annual daily highest temperature with respect to |
January 1j_z | True vertical component of regional Earth's heat flux (usually --uriencoded--\delta T_%7B\rm cut%7D |
|
| Temperature measurement threshold (usually |
j_z = 25 \div 55 \ mW/m%5e2)\delta T_%7B\rm cut%7D = 0.01 \, %7B\rm °C%7D |
| ) |
where
T_D | Daily average surface temperature variation based on weather reports | | Annual average surface temperature based on weather reportsD_T | Period of daily temperature variation cycle: --uriencoded--A_D = 1 \, %7B\rm day%7D | TVDss of the Earth's surface in a given location. In case the Earth's surface is at sea level then |
aeAverage Thermal diffusivity of the soil between Earth's surface and neutral layer
Assumptions
\delta t_D | Time moment of daily highest temperature with respect to Midnight 00:00--uriencoded--%7B \bf j%7D(x,y,z) = \%7B 0, \, 0 , \, j_z = %7B\rm const%7D \%7D |
|
| |
---|
r(l)Rock Thermal Conductivity along-hole profile | Assumptions
Equations
...
Equations
...
LaTeX Math Block |
---|
| T_G(t, z) = T_s + \ |
| frac | j_ | }{\lambda | e | z- | _s) dz + T_Y(t, z) + T_D(t, z) |
| | =\frac{d T_g}{d | z} | r |
| AY(t,z) = \delta T_A \, \exp \left[ \, {(z_s-z}) \sqrt{\frac{\pi}{a_ |
| e{en} \, A_T}} \, \right] \, \cos \left[ \, 2 \pi \frac{t - \delta t_A}{A_T} + (z_s -z) \sqrt {\frac{\pi}{a_ |
| e |
LaTeX Math Block |
---|
| T_D(t,z) = \delta T_D \, \exp \left[ \, {(z_s-z}) \sqrt{\frac{\pi}{a_ |
| e | A | DT}} \, \right] \, \cos \left[ \, 2 \pi |
| \frac{t - \delta t_D}{D_T} + (z_s -z) \sqrt {\frac{\pi}{a_e \, D_T}} \, \right]
where
...
...
TVDss of the Earth's surface in a given location. In case the Earth's surface is at sea level then
\frac{t - \delta t_D}{D_T} + (z_s -z) \sqrt {\frac{\pi}{a_{en} \, D_T}} \, \right] |
|
Neutral Layer |
---|
LaTeX Math Block |
---|
| z_n = z_s + H_n |
| LaTeX Math Block |
---|
| H_n = \sqrt{\frac{a_{en} \, A_T }{\pi}} \, \ln \frac{\delta T_A }{\delta T_{\rm cut} } |
|
See Also
...
Geology / Geothermal Temperature Field
See Also
Geology / Geothermal Temperature Profile
[ Geothermal Temperature Field @model ] [ Geothermal Temperature Gradient ]
[ Neutral Temperature Layer @model ]
References
...
Kasuda, T., and Archenbach, P.R. "Earth Temperature and Thermal Diffusivity at Selected Stations in the United States", ASHRAE Transactions, Vol. 71, Part 1, 1965.
GeothermalTemperatureProfile.xlsx
...