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
...
...
...
wellbore trajectory Measured Depth with reference to Earth's surface (
)...
...
Period of annual temperature variation cycle:
LaTeX Math Inline |
---|
body | --uriencoded--A_T = 1 \, %7B\rm year%7D |
---|
|
...
...
...
--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) \} |
...
...
True vertical component of regional Earth's heat flux (usually
LaTeX Math Inline |
---|
body | --uriencoded--j_z = 25 \div 55 \ mW/m%5e2 |
---|
|
)...
...
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 based on weather reports | | Daily average surface temperature variation based on weather reports |
LaTeX Math Inline |
---|
body | --uriencoded--a_%7Ben%7D |
---|
|
| Local average Thermal diffusivity of the soil between Earth's surface and NTL | | Period of daily temperature variation cycle: LaTeX Math Inline |
---|
body | --uriencoded--A_D = 1 \, %7B\rm day%7D |
---|
|
|
| Subsurface Thermal Conductivity profile as function of TVDss |
Annual average surface temperature based on weather reports | | Period of daily temperature variation cycle: LaTeX Math Inline |
---|
body | --uriencoded--A_D = 1 \, %7B\rm day%7D |
---|
|
| | Average Thermal diffusivity of the soil between Earth's surface and neutral layer moment shift of daily highest temperature with respect to Midnight |
00:0000:00 |
|
| LaTeX Math Inline |
---|
body | --uriencoded--\delta T_%7B\rm cut%7D |
---|
|
| Temperature measurement threshold (usually |
\lambda_r(l) | Rock Thermal Conductivity along-hole profile | Assumptions
Equations
...
LaTeX Math Block |
---|
|
z_N = z_s + H_N |
LaTeX Math Block |
---|
|
H_N = \sqrt{\frac{a_e \, \delta_T }{\pi}} \, \ln \frac{T_A }{\delta T_{cut} } |
...
--uriencoded--\delta T_%7B\rm cut%7D = 0.01 \, %7B\rm °C%7D |
| ) |
where
| |
| TVDss of the Earth's surface in a given location. In case the Earth's surface is at sea level then |
Assumptions
LaTeX Math Inline |
---|
body | --uriencoded--%7B \bf j%7D(x,y,z) = \%7B 0, \, 0 , \, j_z = %7B\rm const%7D \%7D |
---|
|
| LaTeX Math Inline |
---|
body | \lambda_e(x,y,z) = \lambda_e(z) |
---|
|
|
---|
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
...