寒潮预警一来,最怕地暖管里的水结冰涨裂。但未来几天的气温不是一个数,而是"最低可能−8℃、也可能−15℃"的波动。与其拍脑袋,不如用蒙特卡洛模拟:把每个不确定因素按各自的概率分布随机抽样,重复跑几千次,直接算出"管内存水结冰"的概率,帮你在"排空防冻"和"再观察"之间做决定。
为什么用蒙特卡洛而不是算一个代表值
地暖管是否冻裂,取决于管内水温有没有跌破 0℃。而管内水温由室外温度曲线、系统流量、管道保温、甚至室内残余热量共同决定,这些输入个个都有不确定性。蒙特卡洛不追求一个"平均天气",而是承认不确定、用随机抽样去覆盖各种可能,再看坏情况出现的频率——这正好契合"风险"的本质。
模拟怎么做:四步
第一步,给每个不确定输入定分布:比如室外最低温服从均值 −10℃、标准差 3℃ 的正态分布;风速、保温质量、流量各取自己的分布。第二步,每次随机抽一组输入,用简化传热模型算出这条天气路径下管内的最低水温。第三步,重复抽样几千次,记录每次的最低水温。第四步,统计"最低水温 < 0℃"出现的次数占比,就是冻裂概率。
一个直觉例子
假设抽 2000 组天气样本,其中 120 组管内水温跌到冰点以下,冻裂概率约 6%,属于可接受、可观察。若寒潮更强,最低温均值降到 −14℃、标准差增大,同样模型下跌破冰点的样本可能涨到 800 组,概率升到 40%——这时就该果断排空或加防冻液,而不是赌预报里的"最高气温"。
它和确定性算法的关系
简单工况可以用解析或数值模型直接算临界温度;但输入多、关系非线性时,蒙特卡洛几乎是唯一可行的"暴力"办法。它和用牛顿迭代法反解供水温度、欧拉法算降温曲线同属数值家族,只是它处理的是随机输入而非确定方程。结论用概率说话,比单看一个预报数字稳得多。




