NetLogo城市热岛模型构建求助：解决温度指数增长问题

问题背景 #

正在构建课程作业用的城市热岛NetLogo模型，非专业编程人员。原本尝试用公式ΔT = Q / (m·c)模拟热量吸收与释放，因实现困难做了简化。核心需求：

• 城市斑块（urban?为真的灰色斑块）比乡村斑块吸收更多热量（已通过乡村反照率乘数实现）
• 城市斑块储存热量，当周边城市空气温度低于建筑温度（昼夜循环中发生）时缓慢释放，起到降温缓冲作用
  目前遇到的问题：功能无法正常运行，城市斑块的ptemp（建筑温度）与held-heat（储存热量）出现指数增长，需要将数值限制在合理水平。

原代码核心问题分析 #

1. 白天热量累加逻辑循环依赖：用邻居ptemp的平均值计算热量盈余，再叠加到held-heat，随后ptemp又直接等于c-u-air-temp + held-heat，导致每一轮tick都会让ptemp和held-heat互相驱动增长，引发指数飙升。
2. 夜间热量释放公式失控：((0.5 + temp-diff) * temp-diff)会随温差增大呈现二次增长，释放量快速失控，同时还将该值叠加到ptemp，进一步加剧温度异常。
3. 空气温度扩散逻辑无物理依据：diffuse-into-air中的((ptemp * c-u-air-temp) / base-air-temp)公式会导致空气温度异常放大，不符合热交换规律。
4. 缺少数值上限约束：未对ptemp和held-heat设置合理上限，无法模拟实际建筑材料的热容与最高温度限制。

修正后的代码 #

; 新增全局变量，定义温度和储热的合理上限（可根据实际场景调整）
globals [
  max-ptemp ; 建筑最高温度
  max-held-heat ; 最大储热量
  urban-heat-capacity ; 城市斑块热容系数
]

to setup ; 需在setup过程中初始化新增变量
  ; 其他原有初始化代码...
  set max-ptemp 40 ; 示例：建筑最高40℃
  set max-held-heat 10 ; 示例：最大储热量对应10℃温差
  set urban-heat-capacity 0.2 ; 城市斑块热容，控制储热速率
end

to update-day
  set day-timer day-timer + 1

  ; 昼夜切换逻辑保留
  if day? [
    daily-solar-input
    if day-timer >= day-length [
      set day-timer 0
      set day? false
    ]
  ]
  if not day? [
    nightly-temp-loss
    if day-timer >= night-length [
      set day-timer 0
      set day? true
      set days-passed days-passed + 1
      set solar-input solar-input + warming-trend
      
      ask workers [
        if hsp-timer > 1 [
          set hsp-timer hsp-timer - 1
        ]
      ]
    ]
  ]
end

to daily-solar-input
  set c-air-temp c-air-temp + solar-input
  set c-u-air-temp c-u-air-temp + solar-input ; 城市无反照率
  set c-r-air-temp c-r-air-temp + (solar-input * 0.8) ; 乡村反照率0.8
end

to nightly-temp-loss
  set c-air-temp c-air-temp - 0.04
  set c-u-air-temp c-u-air-temp - 0.04
end

to update-urban-air-temp
  if day? [
    ask patches with [urban?] [
      ; 白天：基于太阳输入和热容计算储热，避免循环依赖
      let heat-gain solar-input * urban-heat-capacity
      ; 储热不超过上限
      set held-heat min list max-held-heat max list 0 (held-heat + heat-gain)
      ; 建筑温度 = 城市空气温度 + 储热，不超过最高温度
      set ptemp min list max-ptemp (c-u-air-temp + held-heat)
    ]
  ]
   
  if not day? [ 
    ask patches with [urban?] [
      let temp-diff ptemp - c-u-air-temp
      ; 仅当建筑温度高于空气温度时释放热量
      if temp-diff > 0 [
        ; 线性释放：基于储热和温差的比例，控制释放速率
        let heat-release held-heat * 0.1 * (temp-diff / max-ptemp)
        set held-heat max list 0 (held-heat - heat-release)
        ; 建筑温度降低，同时提升空气温度
        set ptemp ptemp - heat-release
        set c-u-air-temp c-u-air-temp + heat-release
      ]
    ]
  ]
end

; 移除原diffuse-into-air，将热交换逻辑整合到夜间释放中

to update-rural-air-temp
  if day? [
    ask patches with [not urban?] [
      set ptemp c-r-air-temp
    ]
  ]
  
  if not day? [
    set c-r-air-temp c-air-temp
    ask patches with [not urban?] [
      set ptemp c-r-air-temp
    ]
  ]
end

关键修改说明 #

1. 新增约束变量：通过max-ptemp、max-held-heat和urban-heat-capacity直接限制数值范围，模拟实际建筑的热物理特性。
2. 重构白天储热逻辑：改用太阳输入乘以热容系数计算储热，避免依赖邻居ptemp的循环累加，从根源上阻止指数增长。
3. 修正夜间释放逻辑：采用线性比例释放热量，释放速率与当前储热和温差挂钩，确保数值变化平缓可控，同时直接联动空气温度变化，符合热交换规律。
4. 简化热扩散逻辑：将原diffuse-into-air的功能整合到夜间释放过程中，避免冗余计算和异常公式导致的温度飙升。

内容的提问来源于stack exchange，提问作者Noah