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
LaTeX Math Inline |
---|
body | --uriencoded--%7B \bf j%7D(x,y,z) = \%7B j_x, \, j_y, \, j_z \%7D |
---|
|
to be independent on areal location: and Thermal Conductivity %7B bf j%7D(x,y,z) = j%7D(z) and vanishing lateral components: LaTeX Math Inline |
---|
body | --uriencoded--%7B \bf j%7D(x,y,z) = \%7B j_x = 0, \, j_y = 0 , \, j_z(z) \%7D |
---|
|
.Normally there are no heat sources within a subsurface volume under study other than 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.This simplifies the procedure of modelling the Geothermal Temperature Field
LaTeX Math Inline |
---|
body | --uriencoded--%7B \bf j%7D(x,y,z) = \%7B 0, \, 0 , \, j_z \%7D |
---|
|
along a given wellbore trajectory.Outputs
...
...
along-hole geothermal temperature profile
...
...
...
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
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 | Annual average surface temperature based on weather reports | | Daily average surface temperature variation based on weather reports | Average Thermal diffusivity of the soil between Earth's surface and neutral layer | | Period of daily temperature variation cycle: LaTeX Math Inline |
---|
body | --uriencoded--A_D = 1 \, %7B\rm day%7D |
---|
|
| Rock Thermal Conductivity along-hole profile | | Time shift of daily highest temperature with respect to Midnight 00:00 | LaTeX Math Inline |
---|
body | --uriencoded--\delta T_%7B\rm cut%7D |
---|
|
| Temperature measurement threshold (usually | Local Calendar Time | | Annual average surface temperature variation based on weather reports |
| | | Period of annual temperature variation cycle: |
\delta T_ cut%7D = 0.001 \, %7B\rm °C%7D) ...
wellbore trajectory Measured Depth with reference to l | surface (l=0 | ) ...
| Time shift of annual highest temperature with respect to January 1 |
|
--uriencoded--%7B \bf j%7D(x,y,z) = \%7B 0, \, 0 , \, j_z = %7B\rm const%7D \%7D | Equations
...
LaTeX Math Block |
---|
|
z_N = z_s + H_N |
LaTeX Math Block |
---|
|
H_N = \sqrt{\frac{a_e \, A_T }{\pi}} \, \ln \frac{\delta T_A }{\delta T_{\rm cut} } |
...
...
...
...
| 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 | | Time shift of daily highest temperature with respect to Midnight 00:00 |
|
| LaTeX Math Inline |
---|
body | --uriencoded--\delta T_%7B\rm cut%7D |
---|
|
| Temperature measurement threshold (usually LaTeX Math Inline |
---|
body | --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{en} \, D_T}} \, \right] \, \cos \left |
| [ \, 2 \pi \frac{t - \delta t_D}{D_T} + (z_s -z) \sqrt {\frac{\pi}{a_e \, D_T}} \, \right]
where
[ \, 2 \pi \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} } |
|
| TVDss of the Earth's surface in a given location. In case the Earth's surface is at sea level then LaTeX Math Inline |
---|
body | z_s = 0
See Also
...
Geology / Geothermal Temperature Field / 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
...