reverse_osmosis_model

RO膜污染监控与CIP预测系统

系统简介

基于压差预测数据分析膜污染速率，自动计算最优CIP清洗时机。

核心原理：

• 通过滑动窗口计算污染速率 k 值（压差上升速度）
• 识别 k 值持续上升的时间段，找出污染最严重的时机
• 综合一段、二段数据，给出最优CIP建议

主要功能：

• 支持 RO1-RO4 四个机组同时分析
• 提供 4 种决策策略（最早/最晚/加权/严重度）
• 自动生成日志、图表、HTML报告
• 支持 API 回调推送决策结果

快速开始

1. 基本运行

from main_simple import main

# 使用默认参数（策略3-加权平均，当前时间，所有机组）
result_df = main()

# 指定策略和时间
result_df = main(strategy=3, start_date='2025-08-26 00:00:00')

# 只分析特定机组
result_df = main(strategy=1, unit_filter='RO1')

2. 自动监控模式

# 自动持续监控，到期触发预测
python smart_monitor.py

监控系统会：

• 检查 config.json 中的 CIP 历史时间
• 到期自动运行预测并更新下次时间
• 循环监控，无需人工干预

3. 策略说明

策略	说明	适用场景
1	最早时机 - 取一段/二段中较早的	保守维护
2	最晚时机 - 取一段/二段中较晚的	最大化运行时间
3	加权平均 - 根据k值加权（推荐）	平衡维护和运行
4	污染严重 - 选k值最大的段	基于实际污染程度

配置文件 (config.json)

{
  "api": {
    "base_url": "http://example.com",
    "callback_endpoint": "/api/cip/callback",
    "jwt_token": "your_token"
  },
  "scada": {
    "project_id": "12345"
  },
  "cip_history": {
    "RO1": "2025-08-05 10:00:00",
    "RO2": "2025-07-20 15:00:00",
    "RO3": "2025-08-01 08:00:00",
    "RO4": "2025-07-15 12:00:00"
  }
}

说明：

• cip_history：记录各机组最近一次CIP时间，用于计算预测周期
• api：可选，配置后会自动推送决策结果到回调接口
• scada.project_id：项目标识，用于API回调

算法原理

核心逻辑

1. 获取预测数据 → 从污染预测模型获取未来压差数据
2. 计算k值 → 7天滑动窗口线性回归，得到污染速率序列
3. 识别上升段 → 找出k值连续上升的时间段
4. 应用约束 → 至少距离起点30天（避免过早CIP）
5. 找最优点 → 在有效时间段内选k值最大的点
6. 策略综合 → 对一段、二段按策略综合决策

参数调整

核心参数在 OptimalCIPPredictor 类中：

• window_days=7：k值计算窗口，增大平滑噪声，减小提高敏感度
• min_continuous_rising=3：k值连续上升最小点数，增大避免误报
• min_delay_days=30：最小延迟天数，根据实际CIP周期调整

输出示例

控制台输出

分析机组 RO1
  预测天数: 90天
  找到RO1压差列: 2个
  C.M.RO1_DB@DPT_1_pred: 2025-10-15 18:00:00 (第71天, k=0.000234)
  C.M.RO1_DB@DPT_2_pred: 2025-10-18 12:00:00 (第75天, k=0.000198)
  RO1最优CIP时机: 2025-10-16 08:00:00

==================================================
分析完成
==================================================
机组类型           CIP时机              策略说明
   RO1  2025-10-16 08:00:00  加权平均策略 (第72天)
   RO2  2025-11-04 06:00:00  加权平均策略 (第91天)
   RO3  2025-10-28 14:00:00  加权平均策略 (第84天)
   RO4              None      无有效CIP时机

日志与图表

系统自动生成完整的分析记录，保存在 analysis_logs/ 目录：

目录结构

analysis_logs/
├── CIP_Analysis_20251015_095138.log              # 主日志
├── data/                                          # 数据文件
│   ├── CIP_Analysis_20251015_095138_prediction_data.csv    # 原始预测数据
│   ├── CIP_Analysis_20251015_095138_RO1_pressure_data.csv  # RO1压差数据
│   ├── CIP_Analysis_20251015_095138_RO2_pressure_data.csv
│   └── CIP_Analysis_20251015_095138_final_results.csv      # 最终结果
├── plots/                                         # 图表
│   ├── CIP_Analysis_20251015_095138_pressure_trends.png    # 压差趋势图
│   └── CIP_Analysis_20251015_095138_unit_comparison.png    # 机组对比图
└── reports/                                       # 报告
    ├── CIP_Analysis_20251015_095138_analysis_data.json     # 完整分析数据
    └── CIP_Analysis_20251015_095138_report.html            # HTML报告

主要输出

1. 压差趋势图 (pressure_trends.png)

• 4个子图，分别展示RO1-RO4的压差变化
• 一段、二段压差曲线对比
• 自动截取到各机组预测周期

2. 机组对比图 (unit_comparison.png)

• 左图：各机组预测周期对比（柱状图）
• 右图：CIP建议时机分布（距起点天数）

3. HTML报告 (report.html)

• 会话信息、输入参数、分析结果
• 带样式的表格展示
• 适合分享和存档

4. 主日志 (.log)

• 详细记录每个机组的分析过程
• k值计算、时间段识别、策略选择
• API调用结果、耗时统计

系统集成

API回调

配置了API的话，系统会自动推送决策结果：

推送格式：

{
  "list": [
    {
      "type": "RO1",
      "project_id": "12345",
      "ctime": "2025-10-16 08:00:00"
    }
  ]
}

失败会自动重试3次，间隔10秒。

依赖模块

• fouling_model_0922.predict.Predictor：压差预测模型
• cip.run_this.main：计算各机组预测周期
• logging_system.CIPAnalysisLogger：日志和图表生成

---

注意事项

1. 预测周期不够：确保预测数据至少覆盖30天以上
2. 频繁无结果：可以降低 min_continuous_rising 或 min_delay_days
3. 结果不合理：检查预测模型质量，或尝试不同策略
4. 查看详细过程：打开 analysis_logs/ 目录下的日志文件