本发明公开了一种三维地温场动态数值模拟方法、设备及介质，将目标区域进行空间离散化和时间离散化，在空间上，结合背景场控制方程和异常场控制方程迭代计算空间域总场温度；在时间上，结合显式差分格式的地温场递推公式可以由当前时间节点的空间‑波数域异常场温度得到下一时间节点的空间‑波数域异常场温度，进而迭代计算下一时间节点的空间域总场温度；以此类推得到所有时间节点的空间域总场温度。实现了三维地温场动态数值模拟，能够真实反映地温场的动态变化，将背景场控制方程、异常场控制方程与地温场递推技术结合，保证精细模拟温度场动态数值的同时，降低了计算及存储成本，提高了计算效率。