|
/ a# l& k5 A* T. ~6 h
6 J7 [" _7 P6 @. @# A2 V 文/毛小平
) O' t1 `! o- w' i8 |1 q 8 F- W2 s+ H, u
大地热流的本质是什么？
! |" u$ f0 ]- M 这个问题似乎大家都知道，问此问题也似乎被认为在搞笑。不就是地球深部向地表传热的快慢吗？还用问，多此一问。遇到高热导率的地层，传热快，热流就高；遇低热导率地层，传热慢，热流就低。因为大地热流公式长这样：Q=K*△t/△z，与热导率K成正比啊，与地温梯度△t/△z成正比啊，难道是成反比不成？
6 N+ T9 w4 w, b% ?% t: L: l; `# n8 ] 比如郎旭娟等（2016）计算得到的扎仓寺地热深钻 ZR1孔热流测量的结果，如表1所示，感觉热导率高的层段，热流就高，得出不同深度段热流数值范围是9.9-85.3mW/m2，结果很分散，最小值与最大值相差8倍以上。当然，作者也考虑了对流因素。
P7 L. a9 @ t$ g5 ]$ Z- s 表1 扎仓寺地热深钻ZR1孔热流测量及相关参数 # @* e1 M2 q# C+ B f3 J5 c& I
3 p: ^, ^* m. H2 }
汪缉安等(1985)研究了辽河断陷热流，得出长古1井Ed, Es1, Es3,Es4分别为0.598HFU, 0.891 HFU, 1.721 HFU, 1.952 HFU。由深到浅逐渐变低，到东营组Ed比最大值小4倍，几乎等于“断流”了，即热不流动了。 : O% K, z( P2 ~& J3 z* s- ?
王一波等(2019)发表在《地球物理学报》的文章中，得出了一些高质量热流数据。赶脚和测井曲线完全一样？变化幅度大，如下图所示。 $ p8 [3 R* a- U" \1 q
4 W; P5 v! y3 _& h" B2 T 图1 热流的垂向变化 % [& ~; B; O' t
这样曲里拐弯的变化正常吗？真是这样的变化吗？大地热流的本质是什么？ * g5 c+ a, n7 b! ` ` Y, a- @
不正常，不是大错，而是特错。地球物理学报可是大刊啊？也会错？
3 b% q% [9 ?( q9 V( r 我认为是酱紫的：
& J: b0 i9 f* b$ } 在地层接近水平或物性横向差异不太大的地区，在不考虑对流的情况下，大地热流是一个常数。
3 u/ v1 h& o+ p0 N! q+ O( Y, k/ t 理由如下，按一维来推理，深度为Y方向，t为时间，T为温度,K为热导率，A为放射性生热，P为密度，C为热容。对目前的勘探尺度内，放射性A几乎可以忽略，即可令A=0，在温度场进入稳态后，温度T不随时间变化，这是下面公式右端项为0：
+ W' v0 W6 p R% K image005.png 0 H4 w/ Q- E- r, S6 |
于是，上述公式就变成了： / l+ B. m# j" Y, B+ @ H1 D
image007.png 8 j* w9 k+ e4 L) H
令A=0，则该方程的解是：k∂T/∂y=C(常数），这个常数就是热流。为了理解热流不是可变的，以下图为例。设模型垂向分为两段，分别为1000米，上段热导率为1W/mk，下段2w/mk，顶部恒温20度，底部恒温80度，中间界面初始温度为50度。即为30度/公里的正常地温梯度，这是初始时刻。上段的热流Q1=30mw/m2，下段为60mW/m2。显然，随着时间的推移，中间的50度是不可能维持的。因为上述公式中，∂(k∂T/∂y)/∂y可写成∂(Q)/∂y, Q就是热流，代入上式，表示热流梯度应该为0，方程才平衡，如果有高热流，必然后向低热流地区传导热。最后进入稳态后，中间的温度将变成60度，这时，上、下两段的热流均为40mW/m2。按前述实测数据，底部热导率高，热流就一定大？ 8 N( y1 u$ O( j, {
所以，在不考虑对流、放射性生热等条件时，热流在不同深度可变是不对的，是特错！即热流在垂向上是一个常数：大地热流常数，不同地区有不同的常数。可以这样讲，目前所定义的大地热流，其本质相当于是一个单位面积内，单位时间内由地球内部向外散发的热量，是一个热通量，虽然有时间概念。即，以地质历史尺度衡量，在垂向上，几乎是一个常数。
2 N3 W6 Y5 }$ E0 w7 e1 u
' s6 A+ t `. O8 E6 }" d0 K# ^ （a) 初始温度及热 (b)稳态后的温度与热流
# ]8 ^5 m! Z. f: I 图2 热流常数的形成示意图 8 p* p. p! x, C- I4 k/ q& d
就这么一个小公式，居然出了屁漏，问题出在哪儿？是本文只考虑了热传导，没考虑放射性生热及对流？这种计算得到的类似于测井曲线的大地热流，让热传导方程的脸往哪儿搁？我在想，即使考虑了放射性生热及热对流，也很难得出这种忽高忽低的热流，这在理论上是说不通的，单因素分析能明晰规律。事实上，放射性生热，就近地表10公里范围内，其贡献是不大的，如刘绍文(2017)研究了塔里木盆地热流分配，8公里以内的沉积物贡献只有9mW/m2，且应是递增的，和上面计算的热流起伏动且正负40-80mW/m2，不是一个量级啊。地下介质能对流的储集层少，大量情况是不对流的，能对流的情况从地温曲线上是能看出来的，如赵平和多吉(2003)所展示的西藏羊八井最高温度的钻孔ZK4001和ZK4002，在800米以内是碎裂花岗岩，而之下则是原生的，未破坏的花岗岩，无孔隙，无水，但温度就不增加。 z7 Q: I) E3 b: Z7 T3 Q2 v
( d; F" A( D" w: g% C% X: B
图3 西藏羊八井附近多口钻井测温曲线
: a2 p5 u6 u" [5 o8 p/ v2 i ^ 温度测量没问题，很精确，那问题可能出在热导率测量的代表性问题，它是核心，“世界没有变化，是你的心乱了”。同一套地层多个样品的热导率测量结果的平均值不代表该地层的综合热导率，而是偏向于小值；热导率还有各向异性，顺层和垂直于层理其热导率不同，同时，地热盖层一般为第四系，其热导率误差极大，含水量、孔隙率。为此，提出以下两个概念（好的不学，我也喜欢造概念了，杯具）：
& m5 j5 @, s! q) |# z 大地热流常数：指某一地区由地球深部向浅层的传热通量，对沉积盆地而言，一般某地垂向上只需要求取一个数值即可。某一地层的大地热流常数，取决于垂向上最低热导率，如果为0，即，完全的热阻层，则此地热流为0。如图3所示。
$ ]* |) [3 u, F) _ 综合热导率：某一地层的热导率，并非各样品热导率的平均值Ks，样品热导率和样品的层理方向还有关系。公式如下所示Za,Zb分别为地层顶、底界的深度，K(z)为热导率。 / F! D& _& l: ~( ]* [! _+ r- y7 J
# Z3 u+ f9 N; M6 G2 c/ n c5 h5 l 如果横向热传导非均匀性较强，则热传导会产生一些偏转，有学者称热折射，导致垂向上热流产生一些变化，如下图所示，左侧有一个热阻层，热导率为0。 F1 y( I# `8 L( @! y; H1 S
: A9 C* q5 D" t) F9 G+ d
图4 具有高热阻层下的热流
3 |3 R( j$ P1 e% v9 d (a)温度场℃；(b)热流场mW/m2 " r* j5 A) o9 c
下图为某地区地质模型，红字为热导率，之下二图分别为垂向热流和水平方向的热流，可以看出，虽然由于横向物性非均匀性，局部热传导方向有所偏转，但总的来说，垂向热导率基本为一常数。
* N) V; h C" o ) j8 M/ g3 M+ x' h) t# e" q$ y3 H
图5 某断陷盆地热流强度变化剖面图 ) A0 N7 H, i# e* i
(a)温度剖面；(b)垂向热流；(c)水平方向热流
# f" F$ j7 D N) c 此观点为2020年11月28日在武汉召开的“中国深部地热论坛”所汇报的内容之一，为目前对大地热流的误区二。实际有三个误区：
7 ~. I7 n# H- Z; l 误区一，大量地热田，在浅层出现了高地温梯度，它实际对应为低大地热流。当然也有高地温梯度的确对应了高热流，如共和盆地、海口。
* o( M0 r4 m9 D- a1 V0 c 误区二，如上，就是本博文，大地热流垂向可变，甚至相隔几米就会有变化，这是错误的。
$ W' n9 W9 b4 n) S" p+ y9 S 误区三，洋壳的热流一定高于陆壳，这可能是最严重的错误，在武汉会议上我汇报了，但未引起太多关注，且听下回分解。该错误被板块理论等地球动力学引用，并作为重要证据。 # T/ A* a% i4 q
参考文献
' Z; ~& f% a! Q2 P0 X; i 龚建洛,张金功,惠涛,黄传卿,张林晔,孙志刚,陈晓军. 2014. 沉积岩平行层面与垂直层面方向热导率与孔隙连通性之间的关系.吉林大学学报:地球科学版, 44(6): 1789-1797. % J2 T# D6 P$ n7 L' k9 I' n
郝春艳, 刘绍文, 王华玉, 王良书. 2014. 全球大地热流研究进展[J]. 地质科学, 49(003):754-770. + I5 R: ~. E" q s/ R# M" R
姜光政,高堋,饶松,张林友, 唐晓音, 黄方, 赵平, 庞忠和, 何丽娟, 胡圣标, 汪集旸.2016.中国大陆地区大地热流数据汇编(第四版).地球物理学报,59(8):2892-2910.
U3 \) W6 i. v, x# Z8 g 金春爽,付修根,陈文彬乔德武, 葛佳, 祝有海, 卢振权, 王平康. 2019. 藏北地区钻孔大地热流值测量.地球物理学报, 62(8): 3095-3105. 9 B4 k& ^9 i& k0 h( h
郎旭娟, 刘峰, 刘志明, 蔺文静, 邢林啸, 王贵玲. 2016. 青海省贵德盆地大地热流研究[J]. 地质科技情报, 035(003):227-232. 3 U; d$ z* `: g2 T0 O
雷晓东, 胡圣标, 李娟, 姜光政, 杨全合, 李巧灵. 2018. 北京平原区西北部大地热流与深部地温分布特征[J]. 地球物理学报, 61(09):227-240. . S0 J. i; d. X P) l Y
李林果, 李百祥. 2017. 从青海共和—贵德盆地与山地地温场特征探讨热源机制和地热系统[J]. 物探与化探,41(1): 29-34. 7 U; ?* K2 T0 S* p! W
李宗星, 高俊, 郑策,刘成林, 马寅生, 赵为永. 2015. 柴达木盆地现今大地热流与晚古生代以来构造-热演化[J]. 地球物理学报, 58(010):3687-3705. + M! d/ s. w7 x: S
饶松,胡圣标,朱传庆,唐晓音,李卫卫,汪集旸. 2013. 准噶尔盆地大地热流特征与岩石圈热结构.地球物理学报, 56(8):2760-2770. V( C- V5 s8 j. p1 g
沈显杰, 邓肖粤, 康文华,等. 1985. 热流值的校正和藏南热流测量中的湖底温度波动校正[J]. 地球物理学报, 28(z1):70-92.
" G$ H% n0 o. p! T 沈显杰, 张文仁. 1989. 西藏热流数据最新报道[J]. 科学通报, 5:373-376. # x& y" D+ j: p
沈显杰. 1983. 西藏高原的热流测量[J].科学通报, 14: 876-877.
{ y! e6 E; T5 Y. {7 }+ z 汪缉安, 汪集晹, 王永玲, 张忠义. 1985. 辽河断陷地热基本特征[J]. 石油与天然气地质, 04:347-358.
3 H) L! ~( g+ K0 u" C! I/ g" p5 G 汪集旸,胡圣标,庞忠和,何丽娟, 赵平,朱传庆, 饶松, 唐晓音,孔彦龙, 罗璐, 李卫卫. 2012. 中国大陆干热岩地热资源潜力评估[J]．科技导报, 30(32)：25-31. ) h4 {& A* O5 y1 G: | A% z4 h
王君珂, 朱喜, 刘彦广, 张庆莲, 牛飞, 马峰. 2020. 献县地热田地温场特征及控热因素研究[J]. 能源与环保, 42(1): 113-120. q; d0 V/ g1 Z0 ?6 `. U9 p# d8 j
王一波, 胡圣标, 聂栋刚, 张克松, 姜光政, 王朱亭. 2019. 郯庐断裂带是热异常带吗:来自断裂带南段热流的约束[J]. 地球物理学报, 062(008): 3078-3094.
& L' E, f8 p3 }7 @& T. |% Q 肖琳,李晓昭,胡增辉, 赵晓豹, 赵高峰,陈立根. 2009.土体内在结构对热导率计算模型的影响研究. 地质论评， 55（4）：598-605.
. h d7 }) _+ I% n" m' u 袁彩萍, 徐思煌. 西藏伦坡拉盆地地温场特征及烃源岩热演化史[J]. 石油实验地质, 2000, 22(2):156-160.
' D$ |; K- p# l" [) L n- Y 张超, 胡圣标, 宋荣彩, 左银辉, 姜光政, 雷玉德, 张盛生, 王朱亭. 2020. 共和盆地干热岩地热资源的成因机制:来自岩石放射性生热率的约束.地球物理学报, 63(7): 2697-2709.
' \6 Q$ q: f6 h8 V9 w' c1 S 赵平, 多吉, 谢鄂军, 金建. 2003. 中国典型高温热田热水的锶同位素研究. 岩石学报, 19(3): 569-576. 0 [5 z" K% Z2 X2 T+ h3 h
左银辉, 邱楠生, 邓已寻, 饶松, 徐深谋, 李建国. 查干凹陷大地热流[J]. 地球物理学报, 2013, 56(9):3038-3050. ) ]/ i' I6 c/ N1 t- r
+ F7 e& o2 r, G7 k/ n: p- V6 Y+ s- H
; j4 a- e2 x/ X" ?7 J% Q0 ~
9 y2 }$ `( y. B7 l+ O% x2 a: S; x. ~0 j6 J; f
| 
