地源热泵报价单上最贵的一行是打井，而决定这行数字的那个"每米能取多少瓦"，几乎所有人都当成常数在用——规范给 40~70 W/m，取个中间值，用总负荷一除，孔深就出来了。但把这个数翻译成物理的时候会碰到一个绕不开的问题：岩土不流动，热量只能靠导热往外补，而导热是有时间尺度的。刚开机那一小时，钻孔壁附近的土是满的，随便取；连续取一个冬天，热扰动已经扩散到 4.8 米外，补不上来了。本文用 Kelvin 无限线源解把这条时间轴完整算了一遍，基准取 100 ㎡ 住宅、热指标 50 W/㎡、峰值负荷 5 kW、COP 4.0（从地下取热 3 750 W），竖直单 U 的 De32 PE 管（外径 32 mm、内径 26.2 mm）、钻孔直径 150 mm、岩土导热系数 λs = 1.8 W/(m·K)、体积热容 ρc = 2.4×10⁶ J/(m³·K)、初始地温 14.5 ℃。算下来有四件事跟通行做法不一致：第一，单位孔深取热量根本不是常数，它从开机 1 小时的 86.3 W/m 一路掉到采暖季末的 32.3 W/m，差 2.67 倍；第二，钻孔里那点灌浆和对流热阻 Rb 不是人们以为的小头，短时工况下它能占到总热阻的 53.6%；第三，把管内流量从 0.30 m³/h 提到 0.60 m³/h 跨过 Re = 2300，取热能力涨 15.6%，而泵功只多花 6.86 W，收益是代价的 74 倍——偏偏工程上最常用的 0.3 m³/h 正好卡在层流区；第四，人人担心的"冷堆积"其实被高估了，10 年累积温降只有 2.32 K，同一个冬天的季末温降却有 6.50 K，前者只有后者的三分之一。
一、先把井看成什么：一个时间常数为十年的蓄热库，不是换热器
竖直地埋管换热器的物理图像是：一根单位长度发热功率为 q′（W/m）的线热源埋在无限大介质里，热量沿径向向外扩散。这个经典问题的解是 Kelvin 线源解，钻孔壁处相对初始地温的温升写成热阻形式是：
Rg(t) = E₁(u) / (4πλs)， u = rb² / (4αt)
其中 E₁ 是指数积分，α = λs/(ρc) = 1.8 / 2.4×10⁶ = 7.5×10⁻⁷ ㎡/s 是岩土热扩散率，rb = 0.075 m 是钻孔半径。当时间足够长（u < 0.01，也就是 t > 约 1 天）时，E₁(u) ≈ ln(4αt/rb²) − 0.5772，于是：
Rg(t) ≈ [ln(4αt/rb²) − γ] / (4πλs)
关键就在这个 ln t 上。时间是乘在对数里的，1 小时到 1 天涨了 24 倍，热阻只从 0.0237 涨到 0.1448 m·K/W；1 天到 1 个月再涨 30 倍，热阻是 0.2942；一个采暖季（90 天）0.3428；10 年才 0.5065。整条曲线是压扁的对数，意味着"取多久"这件事比"取多猛"更决定成败。这也顺带说明一件事：1 小时那个点上对数近似误差高达 −86%，短时脉冲必须用完整 E₁，不能套那个简化式。
二、86 W/m 和 32 W/m 是同一口井
给定允许的最低循环液温度 −2 ℃（25% 乙二醇的防冻下限），初始地温 14.5 ℃，可用温差 ΔT = 16.5 K。单位孔深取热量就是 q′ = ΔT / (Rb + Rg)。取普通膨润土灌浆（λg = 1.5 W/(m·K)）、流量 0.30 m³/h 的层流工况，Rb = 0.1675 m·K/W，逐时刻算下来：
| 持续时间 | E₁(u) | Rg (m·K/W) | Rb 占比 | q′ (W/m) |
|---|---|---|---|---|
| 1 小时 | 0.535 | 0.0237 | 87.6% | 86.3 |
| 1 天 | 3.275 | 0.1448 | 53.6% | 52.8 |
| 10 天 | 5.558 | 0.2457 | 40.5% | 39.9 |
| 1 个月 | 6.655 | 0.2942 | 36.3% | 35.7 |
| 采暖季末 90 天 | 7.753 | 0.3428 | 32.8% | 32.3 |
| 1 年 | 9.153 | 0.4047 | 29.3% | 28.8 |
| 10 年 | 11.456 | 0.5065 | 24.9% | 24.5 |
规范里那个 40~70 W/m 的推荐区间，恰好落在"1 天到 1 个月"这一档——它是间歇运行、低负荷率工况下的折算值，不是连续满负荷的保证值。真按它去选孔深，100 ㎡ 的房子会打成什么样？峰值取热 3 750 W 除以 86.3 W/m 只要 43.4 m，除以 32.3 W/m 却要 116.0 m——用短时脉冲的取热能力选孔深，会少打 62.5%，整整差 72.5 米。井打浅了不是"效果差一点"，是第一个冬天过到一半地温就崩了。
三、Rb 不是小头：灌浆和管内对流能占到一半
工程上有个习惯说法："岩土才是控制热阻，钻孔里那点东西无所谓。"上表第三列已经否掉了这句话的一半——1 天工况下 Rb 占 53.6%，即使到采暖季末也还有 32.8%。把 Rb 拆开看三项：
Rb = 1/(2πdih) + ln(do/di)/(4πλp) + ln(rb/(√2·do))/(2πλg)
依次是管内对流（两根管并联）、PE 管壁导热、灌浆导热。第二项只由管材和壁厚决定，是死的：ln(32/26.2)/(4π×0.42) = 0.0379 m·K/W。能动的是另外两项，而第一项的弹性大得离谱。
四、0.30 m³/h 是层流：多花 6.86 瓦泵功换 505 瓦取热
De32 管的内截面积只有 539.1 mm²，25% 乙二醇在 0 ℃ 的运动粘度 ν = 3.0×10⁻⁶ ㎡/s（Pr = 25.35）。算一下 Re = 2300 的临界点：v = 2300×3.0×10⁻⁶/0.0262 = 0.263 m/s，对应流量 Q = 0.263 × 5.391×10⁻⁴ × 3600 = 0.511 m³/h。而工程上给单 U 配的流量通常是 0.3~0.4 m³/h——整个常用区间全在层流区。
层流的 Nu 是常数 4.36，h = 4.36×0.48/0.0262 = 79.9 W/(㎡·K)；一旦跨过 2300，Dittus-Boelter 给出 Nu = 0.023Re0.8Pr0.4，Q = 0.60 m³/h 时 Re = 2700、Nu = 46.6、h = 853.7 W/(㎡·K)，翻了 10.7 倍。对应的对流热阻从 0.0760 掉到 0.0071 m·K/W，整个 Rb 从 0.1675 降到 0.0986。
| 流量 m³/h | 流速 m/s | Re | 流态 | h W/(㎡·K) | Rb | q′ W/m | Δp kPa | 泵功 W |
|---|---|---|---|---|---|---|---|---|
| 0.30 | 0.155 | 1 350 | 层流 | 79.9 | 0.1675 | 32.33 | 4.50 | 1.07 |
| 0.40 | 0.206 | 1 800 | 层流 | 79.9 | 0.1675 | 32.33 | 6.00 | 1.90 |
| 0.511 | 0.263 | 2 299 | 层流 | 79.9 | 0.1675 | 32.33 | 7.66 | 3.11 |
| 0.60 | 0.309 | 2 700 | 湍流 | 853.7 | 0.0986 | 37.38 | 16.65 | 7.93 |
| 0.80 | 0.412 | 3 600 | 湍流 | 1 074.7 | 0.0971 | 37.51 | 27.55 | 17.49 |
注意这张表的形状：层流区里加流量一点用都没有——0.30 到 0.511 把流量加了 70%，q′ 一动不动，因为 Nu 被钉在 4.36 上，只有泵功在白涨。跨过临界点之后才突然开闸。而 0.60 到 0.80 又基本平了，说明 0.6 m³/h 就是性价比拐点：相比 0.30，q′ 从 32.33 涨到 37.38（+15.6%），100 米孔多出 505 W 取热能力，泵功从 1.07 W 涨到 7.93 W，多花 6.86 W。换热收益是泵功代价的 74 倍。
顺便说一下灌浆。把 λg 从 1.5 提到 2.5 W/(m·K)（导热型灌浆），灌浆热阻从 0.0536 降到 0.0322，Rb 从 0.0986 降到 0.0772，q′ 只从 37.38 涨到 39.29，+5.1%。相比之下，把流量拉过临界点给的是 +15.6%。同样一笔预算，花在保证湍流上比花在贵灌浆上划算三倍。
五、100 平米的账：116 米，而且季末地温会掉到 8 度
100 ㎡、50 W/㎡，峰值热负荷 5 000 W，COP 4.0 意味着其中 3 750 W 来自地下。按采暖季末的 q′ = 32.33 W/m 设计，孔深 = 3 750/32.33 = 116.0 m（一口 116 米井，或两口 60 米井）。如果做了湍流 + 导热灌浆，q′ = 39.29 W/m，孔深可以压到 95.4 m——省下 20.6 米，约合六分之一的打井费。
但孔深定下来之后还有一件事要校核：这一口井在冬天结束时会把自己周围的土抽到多冷？全年供热 9 900 kWh（150 天、负荷率 0.55），取热 7 425 kWh，分摊到 116 米上，冬季取热强度 qw′ = 17.78 W/m。用线源解算 150 天：
ΔT = 17.78 × [ln(4×7.5×10⁻⁷×1.296×10⁷/0.075²) − 0.5772] / (4π×1.8) = 17.78 × 0.3654 = 6.50 K
地温从 14.5 ℃ 掉到 8.0 ℃。这个跌幅直接打在 COP 上：蒸发温度按 5 K 端差从 9.5 ℃ 掉到 3.0 ℃，冷凝侧 45+3 = 48 ℃，卡诺效率折 0.5 的压缩机效率系数，COP 从 4.17 掉到 3.57，跌 14.4%。也就是说，同一个系统，11 月开机时的能效和 2 月底的能效能差出一成半，这不是机器坏了，是井的深度决定了它必然如此。想把这个跌幅压到 3 K 以内，孔深得加到约 250 米——成本翻倍换 7% 的季平均能效，划不划算要看电价，但至少这是个能算清楚的选择。
六、冷堆积被高估了：10 年 2.3 度，不如一个冬天掉的 6.5 度
"地源热泵用十年会把地下抽冷"是最常见也最模糊的担忧。把它算清楚需要一个前提：夏季制冷往地下排多少热。100 ㎡ 住宅夏季冷负荷按 35 W/㎡、60 天、负荷率 0.45，年制冷量 2 268 kWh，按 COP 4.5 折算排入地下 2 772 kWh。于是年取热 7 425 kWh、年排热 2 772 kWh，年净取热 4 653 kWh（取/排 = 2.68），折成年均强度 qnet′ = 4.579 W/m。
把这个净强度当连续热流用线源解算长期趋势：
| 孔深 | 冬季取热强度 | 年净强度 | 单采暖季末温降 | 10 年累积 | 20 年累积 | 季末/10年 |
|---|---|---|---|---|---|---|
| 116 m（按季末设计） | 17.78 W/m | 4.579 W/m | 6.50 K | 2.32 K | 2.46 K | 2.80× |
| 200 m | 10.31 W/m | 2.656 W/m | 3.77 K | 1.35 K | 1.43 K | 2.80× |
| 400 m（常见过度配置） | 5.16 W/m | 1.328 W/m | 1.88 K | 0.67 K | 0.71 K | 2.80× |
三个结论。第一，季节波动永远是长期累积的 2.80 倍，这个比值与孔深无关（对数律的平移不变性）——决定系统冬天难不难受的是单季深度，不是十年累积。第二，10 年只降 2.32 K，20 年 2.46 K，第二个十年只多降了 0.14 K（6%），对数律把长期影响压得极扁。第三，冷堆积的量级完全由孔深决定：116 米是 2.32 K，400 米只有 0.67 K，差 3.5 倍。所以"冷堆积"本质上不是独立现象，而是"井打得不够长"的另一种说法——把孔深按季末工况选对，它自己就消失了。
七、孔距 5 米够不够：10 年后邻居贡献 27%
多口井之间会不会互相抢热？同一套线源解换个半径就能回答。热扰动（穿透）半径 rpen = √(4αt)：1 天 0.51 m、1 个月 2.79 m、一个采暖季 4.83 m、1 年 9.73 m、10 年 30.76 m。也就是说 5 米孔距下，第一个冬天结束时邻孔还基本感觉不到，但十年尺度上整个孔群早就在共享同一份热库。
定量看（设计强度 qw′ = 17.78 W/m、qnet′ = 4.579 W/m，孔距 B）：
| 孔距 B | 单季末邻孔贡献 | 10 年邻孔贡献 | 占自身孔壁温降比例（10 年） |
|---|---|---|---|
| 3.0 m | 0.868 K | 0.827 K | 35.7% |
| 4.0 m | 0.537 K | 0.712 K | 30.7% |
| 5.0 m | 0.328 K | 0.624 K | 26.9% |
| 6.0 m | 0.196 K | 0.553 K | 23.8% |
| 8.0 m | 0.063 K | 0.442 K | 19.1% |
反直觉的地方在于：孔距从 5 米拉开到 8 米，10 年干扰只从 26.9% 降到 19.1%，拉了 60% 的距离只换来 7.8 个百分点的改善——还是对数的锅。而四口井两两相邻时，单个邻居贡献 0.624 K，三个邻居合计 1.872 K，已经接近自身孔壁温降（2.319 K）的 81%。所以孔距该怎么定？如果只按单季表现看，5 米绰绰有余（干扰只有自身的 5%）；如果按十年寿命看，5 米带来的额外温降约 1.9 K，相当于把设计孔深偷偷减掉了两成。真正划算的做法不是拉开孔距，而是把省下的场地费直接加到孔深上。
八、结论清单
- q′ 不是常数，是 ln t 的函数：同一口井 1 小时能取 86.3 W/m，采暖季末只有 32.3 W/m，10 年尺度 24.5 W/m。规范 40~70 W/m 对应 1 天~1 月的间歇工况。
- 用短时取热能力选孔深会少打 62.5%：100 ㎡ 按峰值算是 43.4 m，按季末算是 116.0 m，差 72.5 米。
- Rb 能占到总热阻的一半：1 天工况 53.6%，即使季末也有 32.8%；钻孔内部不是可忽略项。
- 0.30 m³/h 在层流区，加了流量也没用：Re = 2300 的临界流量是 0.511 m³/h，跨过去 h 从 79.9 跳到 853.7 W/(㎡·K)，q′ +15.6%。
- 湍流比贵灌浆划算三倍：提流量给 +15.6%，λg 从 1.5 提到 2.5 只给 +5.1%。
- 收益/代价 74 倍：0.30→0.60 m³/h，泵功多花 6.86 W，100 米孔多 505 W 取热能力。
- 季末地温掉 6.50 K，COP 跟着跌 14.4%：14.5 ℃ → 8.0 ℃，能效 4.17 → 3.57，这是孔深决定的必然结果，不是故障。
- 冷堆积只有 2.32 K/10 年，是季末波动的 1/2.8；第二个十年只多降 0.14 K。它本质是"井不够长"的同义词。
- 孔距 5 米在单季尺度绰绰有余，在十年尺度等于偷减两成孔深；拉开到 8 米只换回 7.8 个百分点。
- 对数律是这套系统的总纲：时间、孔距、孔深全都走 ln，所以"加倍投入"换来的永远是固定增量，而不是比例增量。
九、被漏掉的一项：钻孔壁的接触热阻能吃掉半个灌浆热阻
前面把灌浆当成均匀介质算，默认它和孔壁严丝合缝。实际下管灌浆时，PE 管在孔里靠管卡撑着、灌浆从中心管注入自下而上返浆，如果遇到缩径、塌孔或者浆液泌水收缩，孔壁与浆体之间会留下薄薄一层空隙。这层东西的热阻不在设计文件里，却实打实串在 Rb 上。按面积热阻估算：
Rc = δ / (λair · π · db)
取 δ = 0.5 mm 的空气隙、λair = 0.026 W/(m·K)、db = 0.15 m：Rc = 0.0005 / (0.026 × π × 0.15) = 0.0408 m·K/W。而前面算过的 λg = 1.5 灌浆热阻本身才 0.0536 m·K/W——半毫米空气缝就把灌浆段的热阻抬高了 76%。把它加进 Rb，0.0986 变成 0.1394，季末 q′ 从 37.38 掉回 33.5 W/m，前面靠提流量挣来的 15.6% 有七成被这层缝吃回去。
这也是为什么导热灌浆（λg = 2.5）在现场有时"没效果"：它的 0.0322 m·K/W 本来就小，一旦串上 0.0408 的接触热阻，总灌浆段变成 0.0730，比普通灌浆加接触热阻的 0.0944 只好了 23%，远不是实验室里说的 40%。灌浆真正的技术含量在"灌得实"而不是"灌得导热"——压力注浆、泌水率控制、孔径与管束的匹配，这些施工参数对结果的影响量级和材料导热系数是同一档的。
本文不算的部分
本文只算了单孔径向导热这一条路径，以下几件事没有纳入，它们各自还能单独成篇：其一，地下水渗流（热弥散）的对流补热——本文假设岩土是纯导热介质，实际富水地层中渗流会把热扰动带走，等效于把 λs 抬高，但抬高多少取决于流速与弥散度，是另一个量纲；其二，水平埋管与竖埋的对比，这里只算了竖直孔；其三，Eskilson g 函数与本文 E₁ 解的差异——g 函数考虑了钻孔有限长度和地表边界，长时间尺度上比无限线源解偏低约 10%~20%，本文的 10 年数字偏保守；其四，机组侧的部分负荷特性与变频压缩机的 COP 曲线，本文只用了固定 COP 4.0 做折算，而季末地温下降恰好与部分负荷率叠加，两者的耦合效应没有被拆开。
