数学预警模型 · 第二版技术文档
版本: v2(第二版) 算法名称: 传感器标定对齐 + 快速滑动基准预警模型 日期: 2026-08-14
一、v2 相对 v1 的改动
| 项目 | v1 | v2 | 原因 |
|---|---|---|---|
| 预警等级命名 | 1=三级 / 3=一级(一级最严重) | 0=正常 / 1=一级 / 2=二级 / 3=三级(三级最严重) | 等级体系规范:数字越大越严重,最高为三级 |
| 泥石流权重 | 0.40S+0.30Q+0.20dQ+0.10ΔS | 0.30S+0.35Q+0.20dQ+0.15ΔS | 土壤易饱和顶分,降权重;流量(山洪信号)升最高权重 |
| Q_MAX(满槽流量) | 30 L/min | 3 L/min | 对齐 YF-S401 实测量程(实测最大 ~1.7) |
| DQ_MAX(变化率上限) | 5 L/min/s | 0.5 L/min/s | 实测秒级跳变 ~0.3 |
| 尖峰回退阈值 | 2.0 | 0.5 | 同上,让流量数据真正参与推理 |
| 滑动基准冷启动 | 2 小时 | 15 分钟 | 更快自适应当地气候 |
| 山洪独立路径 F | 0.6Q̄+0.4dQ̄ | 保留 | 流量暴涨独立触发,不受土壤门槛限制 |
| LCD 显示 | 无 | 风速/土湿度/水流速/预警等级 | 现场直观显示 |
| 流速脉冲诊断 | 单沿计数 | 上升/下降沿分开计数 + D8 电平 | 快速定位接线/沿方向问题 |
二、算法结构(v2)
最终等级 = max(泥石流 D, 山洪 F, 风速, 流量尖峰, 突变, 土壤饱和) ← 0~3
等级定义(数字越大越严重):
| 等级 | 含义 | 触发示例 |
|---|---|---|
| 0 | 正常 | 各因子均低于阈值 |
| 1 | 一级预警(最轻) | D/F ≥ 0.45,或风速 ≥25 km/h,或突变 |
| 2 | 二级预警 | D/F ≥ 0.55,或风速 ≥40 km/h,或尖峰 ≥P95 |
| 3 | 三级预警(最严重) | D/F ≥ 0.75,或风速 ≥60 km/h |
1. 泥石流路径
D = 0.30 × S̄ + 0.35 × Q̄ + 0.20 × dQ̄ + 0.15 × ΔS̄
- 土壤 <60%(未饱和)时跳过本路径(山洪/风速不受限)
- 特殊规则:土壤超当地 P90 + 流量上涨 + 5 分钟持续上升 → 直接三级;土壤超当地 P90 + 流量 ≥70% 满槽 → 二级;土壤 10 分钟涨幅 ≥20% → 至少一级
2. 山洪路径(独立指数)
F = 0.60 × Q̄ + 0.40 × dQ̄
- 不受土壤门槛限制,流量暴涨即触发
- 曼宁公式 70% 分量的流量近似替代(本方案不采用超声波水位传感器——水位与流速功能重叠,流速为核心)
3. 风速路径
≥25 km/h → 一级;≥40 km/h → 二级;≥60 km/h → 三级
4. 流量尖峰路径
秒级变化率超当地历史 P95 → 二级;≥0.8 归一化 → 一级
5. 突变路径(秒级响应)
土壤 5 秒内涨 ≥20% / 风速 60 秒内涨 ≥12 km/h → 立即一级
6. 土壤饱和路径(v2 新增)
土壤 ≥85% 或超当地 P90 → 至少一级(土体饱和本身就是危险,不依赖流量;持续生效,直到土壤回落)
7. 滞回保持(v2 新增,防抖)
- 升级:立即生效并保持
- 降级:需超过保持时间才生效 —— 一级 10 秒 / 二级 15 秒 / 三级 20 秒
- 防止单点抖动/瞬时回落造成的等级闪烁
三、滑动基准(v2 调整)
- 窗口 24h 每分钟采样,把当前值映射到当地 [P10,P95] 区间 → 0~1
- 冷启动缩短至 15 分钟:数据攒 15 分钟即启用,之前的 2 小时等待期砍掉
- 数据不足时回退绝对值:Qn=flow/3、dQn=变化率/0.5(对齐传感器量程)
四、LCD 显示(v2 新增)
- 硬件:I2C LCD 1602(SCL→A5,SDA→A4,地址 0x27)
- 显示内容(ASCII,每秒刷新):
W12.3 S45% F0.50 ← 风速 km/h · 土湿度 % · 水流速 L/min
ALARM:3 L3 ← 预警等级 0~3(0=正常,3=三级最严重)
- 风速与预警等级由服务端(server.py)通过串口下发(
D<风速>,<等级>协议),土壤与流速由 Arduino 本地读取
五、优势
- 流量数据真正参与推理:标定参数对齐传感器量程后,Q̄/dQ̄ 从"压扁的蚊子腿"变成有效因子,山洪路径响应真实水流变化
- 滑动基准快速启用:15 分钟冷启动,插电一刻钟后自适应当地气候
- 等级体系规范:0~3 数字越大越严重,LCD/网页/文档口径统一(三级=最高)
- 现场可视化:LCD 实时显示四路核心数据,不依赖电脑/网页即可观测
- 诊断能力:脉冲双向沿计数 + 电平输出,接线/信号问题 30 秒定位
- 保留 v1 全部优点:多路径独立、简单可解释、冷启动可用
六、缺陷
- 流量计量程仍窄:YF-S401 实测 0.05~3 L/min,真实洪水流量远超量程——标定参数是为当前测试环境(水泵+管道)服务的,换环境需重新标定
- 风速数据依赖外部电脑(192.168.1.200:8001),该服务离线时风速路径失效(LCD 风速显示 0)
- 滑动基准 15 分钟仍可能不够稳:突发天气下前 15 分钟基准质量差,回退绝对值精度有限
- 突变阈值仍是经验值(土壤 5s+20%、风速 60s+12km/h),未做统计标定
- LCD 仅 ASCII:无法显示中文,现场人员需熟悉 0~3 数字含义
- 等级硬阈值:0.45/0.55/0.75 阶梯跳变,无置信度/持续性判定(单点抖动可能误触发)
七、特点
- 实时性:每秒重算 + LCD 每秒刷新 + 网页 0.5 秒刷新
- 双端显示:网页(完整指标)+ LCD(现场核心指标)
- 参数可标定:Q_MAX/DQ_MAX/尖峰阈值来自实测,非拍脑袋
- 模块化:5 条路径独立,后续可单独迭代
八、版本规划
- v2.x 候选:突变阈值按历史分布标定(替代经验值);
预警加持续性判定(连续 N 秒超阈值才触发,防抖)—— 已实现:滞回保持(一级 10s / 二级 15s / 三级 20s) - v3 候选:风速接入本机(传感器迁移到综合板,摆脱外部电脑依赖);LCD 显示山洪 F/泥石流 D 数值