iRAS 资料库
iRAS
循环水养殖系统 (RAS)
工程设计技术手册
完整版 · v1.9
版本: v1.9 · 2026 年 6 月
物理框架重构 · 文献校准 · 工业实测对齐
全局气候 + 围护 · 多模块并联 · P&ID 自动生成 · 投资财务评价
可研报告生成 · 站级设备 · UI 重构 · 工程逻辑精化
氧锥拓扑建模 · 系统排水重构 · 全面审计修复

版本演进说明 (v1.0 → v1.9)

iRAS 项目从 2026 年初首发以来, 经历了 8 次重大版本升级, 每次都基于公开文献或工业项目实测数据校准, 不是经验调整。 本手册整合 v1.1 (2026 年 4 月) 物理框架重构 + v1.2 (热平衡建模) + v1.3 (工程化模块) + v1.4 (投资财务评价) + v1.5 (可研报告生成器) + v1.6 (站级设备 + UI 重构) + v1.7 (氧锥拓扑 + 系统排水重构 + 审计修复) + v1.9 (碳酸盐平衡求解 + o2 物理校准 + finance 零利率修复) 的全部内容, 是 RAS 工程设计的完整参考。

v1.0 → v1.1: 物理框架重构 + 文献校准

A. 严格分离三个时间尺度

v1.0 在 calcStageLoad 中用 feedDesign = feedMax × peakFactor 算所有"日产污染物"。但 feedMax 已是日峰值, 再乘 1.2-1.5 在量纲上是"日总量再放大", 所有下游负荷被一刀切偏大约 17%。v1.1 严格分离三个时间尺度:

peakFactor 物理意义钉死: 代谢污染物 (TAN/O₂/CO₂) 小时峰值/全日均速率, 默认 1.5 (Timmons & Ebeling 2010 标准)。

B. 设备按动力学分类选型

类别设备选型基准物理理由
慢动力学BF 滤池、反硝化、CO₂ 脱气、微滤、蛋分、UV日量 × safety菌膜代谢有惯性, HRT 数小时缓冲峰值
快动力学鱼池增氧 (液氧)、BF 曝气风机小时峰值 × safety氧气即时消耗, 投喂 30min 内 DO 可见底

C. 鱼生理耗氧按物种 + 温度自动算

v1.0 用 0.25 一刀切违背物种特性差异。v1.1 改为按物种 + 温度自动算:

o2FishFactor = o2BaseFactor × Q10^((tempOpt − 20) / 10)
物种o2BaseFactoro2Q10典型值依据
鲑鱼0.322.00.20-0.23Atlantic Sapphire / Davidson 2014
大菱鲆0.282.00.21-0.22黄海所 / 莱州东方海洋实测
罗非鱼0.162.00.26-0.30Blue Ridge / FAO 实测
石斑鱼0.182.00.28-0.32南海所 实验数据
加州鲈0.192.00.26-0.28广东海大
鳗鲡0.132.00.20-0.22福建天马 (粉料代谢慢)
对虾0.311.30.40Boyd 2018 (Q10 偏低反映高温饱和)

D. CO₂ 脱气塔独立建模

v1.0 没有 CO₂ 脱气塔模型 (仅在增氧节点扣 70% 写死)。v1.1 作为独立设备建模, 默认参数与 Delta Cooling Towers 工业项目实测对齐: G:L=3, 风压 2.5 kPa, 风机效率 0.60, 单次脱气率 65%。

E. 大菱鲆密度修正

v1.0 大菱鲆 15/30/40 kg/m³ 偏低 (与莱州东方海洋等循环水项目实测 50-70 kg/m³ 不符)。v1.1 调整为 25/50/60 kg/m³。修正后大菱鲆综合成本从 81.5 元/kg → 48.1 元/kg, 与国家海水鱼产业体系 2019 调查数据 (46.32 元/kg) 偏差仅 +4%。

F. CAPEX 折旧 + 系统体积建模

v1.1 设备 8 年 + 土建 20 年分别折旧, 加到日运行成本。引入完整 V_total = V_tank + V_process 建模 (典型 V_process 占 V_tank 的 30-60%)。换水率口径改为 V_total (行业标准 Timmons 2010 / 国家海水鱼体系)。详见第 8 章。


v1.1 → v1.2: 全局气候 + 双节点稳态求解

G. 项目级常量管理

v1.1 把热平衡参数(气温/源水温/U 值/换气率) 放在每个阶段的 proc 字段里, 用户必须在每个阶段重复设置。 工程现实中这些是项目级常量。v1.2 把所有热平衡参数挪到顶部全局配置区 globalThermalConfig, 作为项目级常量影响所有阶段。新增地区预设(青岛/海南/挪威 Bergen/法罗群岛)、围护预设、4 工况对比(冬/夏/年均)。详见第 28 章。

H. 双节点稳态求解器

v1.1 用单点能量平衡, 忽略车间空气节点。v1.2 改为双节点 Picard 迭代同时求解空气热平衡 + 空气湿平衡, 输出 T_room、RH、蒸发率/凝结率、热泵装机。详见第 29 章。

📌 v1.2 物理一致性提升
v1.1 在挪威/北方冬季工况(气温 -10°C, 源水 5°C)热泵电费可能低估 30-50%, 因为没考虑车间空气节点和围护散热。 v1.2 双节点求解器修正后, 1000 t/yr 三文鱼项目挪威冬季工况综合成本从 v1.1 的 32.5 元/kg 修正到 36.30 元/kg。

I. 围护结构详细建模

v1.2 引入完整车间围护参数: U_wall / U_roof / U_floor / U_window / windowRatio / airChangeRate / ventHeatRecovery / poolCover。详见第 28 章。


v1.2 → v1.3: 多模块并联 + 工程化

J. 模块化与 N+1 备用

大型项目(>1000 t/yr)拆成多模块并联: 每模块独立循环, 模块间共用液氧站/应急电源/控制中心。每模块内部按容量上限自动配 N+1 备用 (主泵/风机/RDF/UV/热泵)。模块数完全由用户在阶段卡片"🔧 套数" subSystemCount 决定。详见第 26 章。

K. 鱼池配置反推

用户给"每套池数 + 池深", 工具按总水体反推单池直径 D = √(4V/πH)。自动检查工程合理性: D < 4m / D > 30m / D > 20m 配 H < 3m 给告警。SPECIES_DB 各鱼种各阶段加入推荐池深 recDepth_m。详见第 27 章。

L. P&ID 工艺图自动生成 (v1.3.x)

独立的 pid.html 查看器, 主页面"📐 输出 P&ID"一键打开, 自动渲染主流串联 + 3 条旁路 + 公用工程 + 仪表泡泡 + 联锁线 + 设备表。LOX 跟随阶段切换。详见第 30 章。

M. 设备清单导出 (v1.3.x)

独立的 equipment.html 页面, 一键导出当前方案完整设备规格表, 含 N+1 备用展开 + 单价 + 合计金额, 支持打印/导出 PDF。详见第 31 章。

N. 池深统一数据源 (v1.3.x)

v1.3 早期版本中池深有两个独立输入(阶段卡片头 + 热平衡区), 容易混淆。v1.3.x 统一为一个数据源: stage.poolDepth_m → recDepth_m → proc.poolDepth → 1.5。同时用于池径反推 + 热平衡。详见第 27.5 节。

O. 共用车间散热系数 (v1.3.x)

各阶段独立估算的散热相加 × thermalSharingFactor (默认 0.75) = 全场总散热, 反映多阶段共享车间的工程现实(共享外墙节省散热约 25%)。详见第 29.7 节。

P. 集水池建模 (v1.3.x)

V_process 加入集水池项 V_sump = clamp(Q/30, 20, 200) m³, P&ID 加入 S-101 标识。


v1.3 → v1.4: 投资财务评价

Q. 8 张标准财务表 + NPV/IRR + 敏感性分析

v1.3 输出停在工程概算 (CAPEX 万元 / OPEX 元/kg), 用户拿这个数据再去 Excel 手算财务评价 (NPV/IRR/回收期) 是当前主要的工作流痛点——一份可研报告里财务评价占 30% 篇幅, 设计院通常要花 1-2 周做完。 v1.4 加入独立的投资财务评价页面 (finance.html), 主页面"💰 投资财务评价"按钮一键打开, 按发改委可研编制规范输出 8 张标准财务表:

  1. 投资估算表 (设备 / 土建 / 流动资金)
  2. 资金来源运用 + 还款计划 (等额本息 / 等额本金)
  3. 总成本费用估算表 (满产年, 含折旧+利息)
  4. 销售收入和税金表 (按爬坡逐年)
  5. 损益表 (逐年净利润)
  6. 现金流量表 — 全投资 (项目本身 IRR)
  7. 现金流量表 — 自有资金 (股东 IRR, 含贷款杠杆)
  8. 敏感性分析 (5 因子 ±20% 对 IRR 影响)

用户在页面顶部填 17 个参数 (售价、项目年限、产能爬坡、自有资金比例、贷款利率、折现率、所得税等), 计算引擎用逐年现金流模型 + Newton-Raphson 求 IRR, 实时输出关键指标。 支持一键导出 Excel (用 SheetJS, 7 个工作表) 和打印为 PDF。详见第 33 章。

📌 v1.4 设计原则
财务评价模型是个独立外置层: 主页计算的 CAPEX/OPEX/折旧 → 通过 localStorage 共享给 finance.html → 用户在 finance.html 里独立调财务参数 → 不影响主页计算。这种解耦让设计院可以"工程评估和财务评价分人协作": 工程师管主页, 财务管 finance.html, 同一份方案两人各自迭代。

R. 关键指标定义对齐国内规范

v1.4 财务评价 6 个关键指标的定义全部按国内可研规范:

指标定义判断
全投资 IRR (税后)项目所有现金流(不含贷款本息)的内部收益率> 折现率 = 财务可行
自有资金 IRR股东视角现金流(项目 CF + 贷款流入 − 还本付息)的 IRR有杠杆放大效应
NPV (净现值)按折现率折算到 t=0 的累计现金流现值> 0 = 可行
静态投资回收期累计净现金流首次 ≥ 0 的年份(线性插值)RAS 项目典型 5-8 年
动态投资回收期累计折现现金流首次 ≥ 0 的年份静态 + 折现影响
盈亏平衡点 (BEP)满产年净利润 = 0 时的产量比例< 75% = 安全裕度大

S. iRAS 整体架构演进

从 v1.0 到 v1.5 共 1 年多时间, iRAS 完成了从"单页工程概算"到"完整可研工作流"的演进:

版本发布核心能力用户工作流
v1.02026-04 初稳态浓度计算 + 设备选型看几个数,自己再 Excel 加工
v1.12026-04 末+ 物理框架重构 + CAPEX/OPEX + V_total能给出综合成本元/kg
v1.22026-05 初+ 全局气候 + 双节点热平衡 + 4 工况对比挪威/海南都能算
v1.32026-05 中+ 多模块并联 + N+1 备用 + P&ID + 设备清单工程化的设计输出
v1.3.x2026-05UI 整理 + 池深统一 + 共用车间散热用户体验细节打磨
v1.42026-05 末+ 投资财务评价 (NPV/IRR + 8 张表 + Excel 导出)报告级输出
v1.52026-05 末+ 可研报告生成器 (16 章扩展版 + 模板 + docx 输出 + 图表嵌入)从工具到工作流

v1.4 → v1.5: 可研报告生成器

T. 16 章扩展版可研报告自动生成

v1.4 之前的输出层级停在工程数据 (P&ID / 设备清单 / 财务表 Excel), 用户需要把这些数据手动整合到一份 Word 可研报告里。 一份政府版可研通常 80-150 页, 设计院做一份要花 1-2 个月。v1.5 加入可研报告生成器 (report.html), 按发改委编制规范输出 16 章扩展版可研, .docx 格式, 大幅压缩可研撰写时间。

v1.5 严格遵循以下技术决策:

  1. 结构化内容输出, 不替设计院排版 — 生成的 docx 是"内容产品", 用户复制粘贴到自己院的 Word 模板即可
  2. 纯模板 + 条件分支, 不用 LLM — 可研是法律文件, 模板生成有可追溯性, 不能赌幻觉
  3. nunjucks 模板引擎 — Jinja2 JS 移植, 语法标准, 可后续移植到 Python 后端
  4. docx.js 文档生成 — 浏览器侧生成, 无后端依赖, 单文件部署
  5. 模板内嵌 script 标签 — 保持 iRAS"双击 HTML 即用"的传统

U. 16 章覆盖范围与自动化率

名称iRAS 自动化率
1总论85%
2项目背景与建设必要性50%
3市场分析与建设规模60%
4建设场址与建设条件30%
5工艺技术方案95%
6主要设备方案100%
7土建/总图/公用工程方案50%
8节能节水70%
9环境影响评价60%
10劳动安全卫生与消防40%
11项目实施与组织机构65%
12投资估算与资金筹措95%
13财务评价100%
14国民经济评价 / 社会效益评价50%
15风险分析70%
16结论与建议80%

加权平均自动化率约 70%, 用户需要补充的主要是项目身份信息 (业主名/地址等 20 字段)、场址自然条件、土建详细设计、消防/环评/安评等专项内容(这些必须由有资质机构出具)。

V. 工艺流程图 + 财务图嵌入

v1.5 在 docx 中嵌入 3 张自动生成的图:

  1. 图 5-1 工艺流程示意图 — 6 个主流程方框 + 主循环回流 + 旁路单元(蛋分/反硝化/AOP, 按是否启用动态显示)+ 关键工艺参数表
  2. 图 13-1 累计现金流曲线 — 实线全投资 + 虚线自有资金, 含 0 轴线、网格、图例、节点标记
  3. 图 13-2 敏感性 Tornado 图 — 5 个因子按 IRR 摆动从大到小排序的横向条形图

图表生成路径: 浏览器侧 SVG → Canvas → PNG → docx ImageRun。中文字体使用系统 SimHei/Microsoft YaHei, 高清 2× DPR 渲染。

📌 v1.5 设计原则
"不替设计院排版" 是 v1.5 最重要的产品决策。商业级可研生成器(如某些 SaaS)倾向于做"成品 Word", 但实际上每个设计院都有自己的模板(红头、字体、目录格式各不相同), iRAS 替他们决定版式反而帮倒忙。 v1.5 输出"结构化内容", 用户拿到后复制到自己院的模板里, 是更尊重设计院工作习惯的做法。

v1.5 → v1.6: 站级设备 + UI 重构 + 工程逻辑精化

W. 液氧储罐站级化(核心改动)

v1.5 之前, 液氧储罐 (V-901) 作为阶段级设备计入 calcCapex, 每个阶段独立挂一笔储罐 CAPEX。 但 RAS 工程实际是全场单点供应: 多个阶段共用同一座液氧站, 计 4 次显然错误。 4 阶段千吨级项目可虚报约 30 万元 CAPEX。

v1.6 把液氧储罐从阶段 CAPEX 剔除, 挪到全场"③ 站级设备"区域, 增加 calcSiteCapex(stages, eqPrices, allResults) 函数计算站级 CAPEX 并加入工程总投资。用户可在 ③ 区调整规格 / 数量 / 单价, iRAS 按全场日耗自动推荐:

推荐算法: 日耗 < 3,000 kg/d (小/中项目) → 1 主 + 1 备 = 2 罐 日耗 3,000-15,000 kg/d (大项目) → 2 主 + 1 备 = 3 罐 日耗 > 15,000 kg/d (超大项目) → 3 主 + 1 备 = 4 罐 单罐容量 = 总日耗 × 7 天储备 / 主罐数 → 凑到标准档 (1/3/5/10/20/30/50 t)

单价按规模指数公式 P = 100,000 × (size/3)^0.6 取值, 经标准档列表: 1 t ≈ 5 万 · 3 t = 10 万 · 5 t ≈ 14 万 · 10 t ≈ 22 万 · 20 t ≈ 35 万 · 30 t ≈ 47 万 · 50 t ≈ 65 万 (元/罐)。

P&ID 图 (V-901) / 设备清单 / 可研报告 全部同步显示规格 + "全场共用"标注。

X. UI 整体重构

v1.6 把 iRAS 主页面从"一长条信息流"重构为可扫描的工程师视角:

Y. 工程逻辑精化

项目v1.5 旧v1.6 新
MBBR 单价标签"生物滤池单价 (元/m³)""MBBR 池系统 (填料+曝气配件, 不含池体/风机) — 元/m³ 填料净体积"
设备规格清单 MBBR仅显示净填料体积净填料 + 池总容积 (按 fillRatio 反推)
投资概算单价"500 /h"(单位 bug)"500 元/(m³/h)"(正确显示)
设备选型上限硬编码 (8 RDF / 60 m³ MBBR ...)按 modFlow 自适应(支持千吨级单模块)
液氧需求量显示仅在站级备注阶段成本明细 / 设备清单 / 站级备注 共 3 处
B45 边界钳制百分比字段可超 10019 个百分比字段强制 [0, 100] 钳制
报告 FCR 显示没有每阶段表加 FCR 列
打印优化简单 hide buttonA4 + 1.5cm 边距 + 阶段独立分页 + 灰阶图表

Z. v1.5 → v1.6 公式订正表

v1.5 代码已经修复但本手册之前未同步的几处公式:

章节项目v1.5 旧手册v1.6 正确
4.2NO₃ 产生率NO₃ = TAN × η_bioNO₃ = TAN × 1.0 (稳态质量守恒)
8.2V_RDFV = Q × 2/60 (HRT 2 min)V = Q × rdfHRT_min/60 (默认 0.5 min v1.6, 可调 0.1-5)
8.2V_BF (池总)V = bioVolume × 1.3 (经验)V = bioVolume / fillRatio (默认 50% → 2× bio)
15.4CO₂ 单位换算未明确显式 × 10⁶ (kg/m³ → mg/L)
29双节点稳态蒸发显热 Q_evap×λ 写在空气节点移到水节点 (与代码一致)
33.4残值× 5% 笼统设备 5% + 土建 20% (分别折旧)
📌 v1.6 工程哲学
v1.6 的主题是 "把对的事情做对": 液氧站级化解决长期的重复计算 bug, UI 重构让多阶段切换效率提升 ~40%, 工程逻辑精化把"工程师早就发现但容忍的小问题"逐个修掉。 没有添加新功能, 但千吨级项目 CAPEX 概算精度提升约 3% (剔除重复+单位修正), 输出可靠性提升。

v1.6 → v1.7: 氧锥拓扑 + K-101 + 系统排水重构 + 全面审计修复

A. 氧锥工艺拓扑 mainline / bypass 真实区别对待

v1.6 之前, iRAS 在 UI 上已经有"主流路径 / 旁路(节能)"两个选项, 但 calcPumpPowerBreakdown / aerationSpec / OPEX 完全没区分 — 是"UI 假承诺"。v1.7 通过全网调研 (Linde SOLVOX cone / Pentair AES Speece / PR Aqua PPC / Global Seafood Advocate 2019) 确认: 国际主流大型 RAS 项目"more often plumbed in side stream configuration"。原因是氧锥操作压力 10-21 psi (玻璃钢上限) 或 1-3 bar (SOLVOX), 让主流全部过氧锥意味着主泵要承担额外 3-5 m 压头, 而旁路只让 10% 流量过高压锥, 主泵不过锥, 能耗显著下降。

v1.7 新增 5 个工艺字段 (PROCESS_DEFAULTS + UI 输入框 + 工程范围钳制):

字段默认钳制范围工程对照
o2ConeTopologymainlinemainline / bypass
o2ConeHeadLoss_m3 m[0, 15]LHO 2-3 / Speece 5-10 / SOLVOX 110 至 15
o2BypassRatio (%)10[5, 30]国际主流 8-15%
o2BypassPumpHead_m20 m[5, 40]玻璃钢 20 / SOLVOX 25-30
o2BypassPumpEta (%)60[30, 85]高压泵典型

物理改动:

1000 t 大西洋三文鱼实测对比 (Q = 10,155 m³/h, V_total = 5,183 m³, jsdom 跑测 v1.7): mainline bypass 差异 主泵 P-101 功率 610.4 kW 488.3 kW -122.1 kW 旁路泵 P-602 0 92.2 kW +92.2 kW ─────────────────────────────────────────────────────────── 综合泵功率 610.4 kW 580.5 kW 节能 29.8 kW 设备 CAPEX 4,454 万元 3,842 万元 节省 612 万 工程总投资 9,800 万元 8,454 万元 节省 1,346 万 综合成本 (元/kg) 36.70 元/kg 35.41 元/kg 省 1.29 元/kg 全投资 IRR 14.90% 17.31% +2.41%

B. K-101 鱼池增氧鼓风机 (air 模式专属)

v1.7 之前, air 模式 (空气曝气) 只有 K-302 (BF 曝气风机), 缺鱼池主增氧风机。 这意味着温水鱼种 (罗非鱼 / 鳗鱼 / 对虾) 的大规模 RAS 设计无法用工具自动生成 — 而这是 80%+ 中国 RAS 项目的实际场景。

v1.7 加入 K-101 鱼池增氧鼓风机: 罗茨风机 + 微孔盘管曝气, 与 K-302 同口径 1+1 备 (鱼池主增氧故障即鱼窒息, 必须双机轮换)。P&ID / 设备清单 / 报告章节 / 财务测算 全链路同步, modulePerSpec 含完整规格 (风量 Nm³/min / 功率 kW / 风压 kPa / 效率 %)。

C. 系统排水重构 (RDF 反冲废水 → 系统排水)

旧版工艺面板"RDF 反冲(%) 2%"是反冲管瞬时分流比, 但 4870 m³/d ÷ V_total = 101%/d 物理荒谬 (远超合理换水率 1-5%/d)。真相: RDF 反冲水绝大部分经污泥沉淀池回循环 (Sharrer 2010 实测 95-98% 清液回流), 真实穿越系统边界的废液由换水率决定

v1.7 全面清理:

D. 全面审计修复 (27 项 bug)

v1.7 发布前进行了一次彻底审计, 覆盖 Level 1 (静态扫描) + Level 2 (jsdom 三工况数值验证) + Level 3 (核心计算函数逐项审查), 共发现并修复 27 项 bug:

类别数量代表问题
🔴 严重8 项K-101 modulePerSpec 缺失 / 系统排水硬编码 7.5% / 报告增氧章节未区分拓扑 / air 模式无 K-101 描述 / DO 行硬编码 / #25 灾难性 bug(下文专述) / 23 处 *100 单位错误
🟡 中等4 项pid.html isAirMode 反向判定脆弱 / v1.7 新字段未钳制 / TSS chip 标签不清晰 / finance-core.js OPEX 70/30 硬编码
🟠 精度优化4 项蛋分 DOM 移除公式重构为严格物料守恒 / AOP 臭氧剂量字段化 / CO2 风压钳制 [0.3, 10] 补全 / T_room 钳制 ±2°C → ±5°C
🛡️ 鲁棒性3 项阶段计算 try/catch / calcCapex add() isFinite 保护 / 所有阶段崩溃 console.warn
🟢 清理8 项液氧塔 silent fallback warn / 字符串补空格 / 删 P-602 死代码 / 注释更新 / drawMainPump 加 data-cx/cy / 删 skimmerWasteFlow 死字段 / IRAS_META.version 1.6.0 → 1.7.0
⚠️ #25 灾难性 bug (潜伏 1+ 年)
问题: report.html 的 buildContext() 函数 return 时漏写 workflow 字段。
影响: 可研报告里 55 处 workflow.* 引用全部失效。但因为 nunjucks 模板引擎"字段缺失即跳过"的优雅降级机制, 报告渲染不报错 — 模板里所有条件分支都默默走 fallback, 用户拿到的报告看起来正常但实际上增氧章节不区分液氧 / 空气曝气, 主泵参数走默认值不显示真实数值, DO 控制策略描述固定不变。
修复: 一行修复, 但触发 23 处下游 *100 单位错误 (#25b) 的连锁修复。
意义: 这是过去 1 年以上潜伏的关键 bug。报告生成器对用户来说之前实际不能放心用, 但用户和作者从未察觉 — 因为报告"看起来正常"。教训: 优雅降级机制让 bug 隐藏起来, 生产级工具应该有"严格模式", 缺字段就报错而不是静默跳过。

E. v1.7 验证状态

验证维度结果
jsdom 三工况主功能 (mainline / bypass / air)✅ 全通过
P&ID 三工况渲染 (A-601 / V-901 / K-101 / P-602 / M-601 / K-302 各就各位)✅ 全通过
1000 t 大西洋三文鱼基准 vs SUMMARY 核心指标100% 吻合
JSON 加载兼容性 (v1.5 / v1.6 / v1.7 / 缺 aerationMode / 缺整个 proc 5 种降级场景)✅ 全通过
finance 完整闭环 (buildPIDExport → calculate → sensitivity)✅ 全通过
Level 3 审查 6 大块 (总线 / v1.6 同款 9 项 / 热平衡 / CAPEX / 站级 / 累加)✅ 全部已审
📌 v1.7 工程哲学
v1.7 主要意义不在于"加了多少功能", 而在于把工具从"演示原型"推进到"可工程使用的设计平台": (1) 氧锥拓扑与国际主流接轨, 让国内设计院的项目能直接对标挪威/加拿大大型三文鱼 RAS 工程; (2) K-101 补齐让 air 模式从"半残"变为"可用", 覆盖了中国 RAS 实际项目 80%+ 的鱼种; (3) 系统排水重构把陈旧的"RDF 反冲废水"语义清理为质量守恒口径的"沉淀池入口废液"; (4) 报告生成器修复 (#25) 让 16 章可研报告真正可用

v1.7 → v1.9: 碳酸盐平衡求解 + o2 物理校准 + finance 零利率修复

v1.8 由英文版团队 (Polarlys Innovation AS) 主导, 引入碳酸盐平衡求解和 o2 物理校准。v1.9 是 NTNU validation memo 全面验证后的稳定版。中文版直接从 v1.7 跳到 v1.9, 与英文版同步。

A. 碳酸盐平衡求解 (v1.9 核心)

v1.7 之前 iRAS 只输出"碱度需求" (基于硝化耗碱 7.14 g CaCO₃/g N), 不计算 pH / NH₃ 物种 / 缓冲强度 / 方解石饱和度。这使得工具:

v1.9 引入完整的碳酸盐物种平衡求解模块 (carbonate.js, 226 行), 平衡常数全部来自国际权威数据库:

对照 PyCO2SYS 1.8.3 (国际标准实现) 严格验证:

项目偏差
正向 pH (free 标度)max|ΔpH| = 1×10⁻⁵
反向碱度求解0.08 μmol/kg
NH₃ 占比1.5×10⁻¹⁶ (机器精度)
方解石饱和度0.006%

详见本手册第 4 章 §4.5 碳酸盐平衡章节。

B. 两种 pH 控制策略

v1.9 在工艺面板新增"pH 控制"下拉, 提供两种模式:

大型海水 RAS (如 Atlantic Sapphire / Salmon Evolution) 普遍用"持目标 pH"模式, 精确控 pH = 7.0-7.2 以同时防 NH₃ 毒性 + 防 CaCO₃ 沉积堵塞设备。商业项目实测碱度 237 mg/L CaCO₃ (Good et al. 2018), 显著高于淡水默认 150。

C. o2 物理校准 (validation memo F1)

v1.9 修正了一个潜伏的物理建模不一致:

8 鱼种 o2BaseFactor 校准对照:

鱼种tempOptv1.7 (Q10@20°C)v1.9 (Q10@metabRefTemp)数值变化
salmon (三文鱼)15°C0.320.2263等价 (0.32 × 2^-0.5 = 0.2263)
turbot (大菱鲆)16°C0.280.1980等价
tilapia (罗非鱼)28°C0.160.2786等价 (0.16 × 2^0.8 = 0.2786)
grouper (石斑鱼)27°C0.180.2924等价
bass (加州鲈)25°C0.190.2687等价
eel (鳗鲡)26°C0.130.1970等价
mandarin (鳜鱼)25°C0.180.2546等价
shrimp (对虾)28°C0.310.30-3% (轻微下调)
📝 重要 这是数学等价重写, 不是数值校准。新 o2BaseFactor 现在直接是 "tempOpt = metabRefTemp 时的 O₂:feed 比", 语义更清晰, 与 Timmons (2010) 文献口径一致。 1000 t 三文鱼成鱼期 o2DailyTotal 仍是 1113 kg/d, 综合成本仍是 36.70 元/kg (与 v1.7 完全相同)。

D. CO₂ 公式修正 — 加入硝化产 CO₂ (validation memo F2)

v1.7 之前 co2Daily 公式只算呼吸 CO₂:

co2Daily_v1.7 = o2FishDaily × 1.375 (RQ ≈ 1, 仅生理呼吸)

但实际硝化反应也产 CO₂:

NH₄⁺ + 2HCO₃⁻ + 2O₂ → NO₃⁻ + 2CO₂ + 3H₂O
∴ 2 mol CO₂ / mol TAN-N × 44.01/14.01 = 6.286 g CO₂ / g TAN-N

v1.9 修正:

co2Daily_v1.9 = o2FishDaily × 1.375 + tanDaily × 6.286

1000 t 三文鱼成鱼期实测:

v1.7v1.9变化
o2FishDaily × 1.375 (呼吸)700 kg/d700 kg/d
tanDaily × 6.286 (硝化)0 (漏算)521 kg/d+521
co2Daily 总700 kg/d1221 kg/d+74%

影响: CO₂ 脱气塔规格按循环流量 × G:L 选型 (Summerfelt 2003 / 工程惯例), 不直接随负荷上升, 故 CAPEX 不变。但稳态鱼池 CO₂ 浓度会上升, 是 v1.9 引入碳酸模型后能从物理层面捕捉到的关键工况指标。

E. finance-core.js 零利率 bug 修复

v1.7 之前 buildLoanSchedule 函数在等额本息分支:

a = r × (1+r)^n / ((1+r)^n − 1)

当 r = 0 (零利率贷款 / 政策性贴息项目) 时分母为 0, 产生 NaN, IRR 算不出。v1.9 修复: r = 0 时自动降级到等额本金分支 (此时两者数学等价)。这是个边界 case, 多数项目利率 > 0 不受影响。

F. v1.9 验证状态

验证项状态
carbonate.js 中英文版输出一致 (4 工况)✅ max|ΔpH| = 0
8 鱼种 + mainline/bypass jsdom 全工况✅ 9/9 通过
salmon 1000t mainline 综合成本 (基准 36.70 元/kg)✅ 完全一致
所有工况 pH 在合理区间 [6.0, 8.5]✅ 9/9 通过
海水方解石饱和度 Ω 计算 (5/5)✅ 全部正常
淡水 Ω 为 null (符合物理)✅ 4/4 正确
salmon NH₃-N 1.94 μg/L < 12.5 限值✅ 合规
📌 v1.9 工程意义
v1.9 让 iRAS 从"工程算成本工具"升级为"工程算成本 + 水化学合规工具": (1) pH 不再是输入 — 工具按碳酸盐平衡求出来, 业主可信度大幅提升; (2) NH₃ 物种识别 — 高密度温水工况 (罗非 70+ kg/m³, 鳗鱼) 的剧毒 NH₃ 风险首次被工具捕捉, 是国内 RAS 项目"莫名死鱼"问题的解答; (3) 对接国际学术规范 — 碳酸模块 vs PyCO2SYS 实验室级精度 (max|ΔpH| = 1e-5), 工具输出可直接用于 NTNU / 中国海洋大学等学术研究; (4) 商业海水 RAS 设计能力 — Atlantic Sapphire / Nordic Aqua / Salmon Evolution 等国际旗舰项目的高碱度 (237 mg/L) + 控 pH 实践, iRAS 工具可对标。

v1.7 工况成本基线 (1000 t 大西洋三文鱼)

对比 mainline 与 bypass 两种氧锥拓扑在挪威 Bergen 气候下的全工况成本:

项目mainline (主流串联)bypass (旁路高压)差异
工艺指标
系统总水体 V_total5,183 m³5,183 m³
循环流量 Q10,155 m³/h10,155 m³/h
主泵 P-101 功率610.4 kW488.3 kW-122.1 kW
旁路泵 P-602 功率0 (无)92.2 kW+92.2 kW
综合泵功率610.4 kW580.5 kW-29.8 kW (节能)
o2DailyTotal (全场日峰值)1,113 kg/d1,112 kg/d~0
A-601 峰值供氧41.20 kg/h41.20 kg/h
CAPEX
设备 CAPEX4,454 万元3,842 万元-612 万元
工程总投资 CAPEX (× Lang 2.2)9,800 万元8,454 万元-1,346 万元
运营成本 (元/kg)
年总成本 (含折旧)3,670 万元/年3,540 万元/年-129 万元/年
OPEX (不含折旧)28.46 元/kg28.30 元/kg-0.16
综合成本 (含折旧)36.70 元/kg35.41 元/kg-1.29
财务指标 (售价 60 元/kg, 30/70 融资 @ 5%/8y, 折现 8%, 15 年)
全投资 IRR14.90%17.31%+2.41%
自有资金 IRR~22%25.23%+3%
全投资 NPV (万元)5,7797,031+1,252
静态回收期6.19 年5.57 年-0.62 年
动态回收期8.63 年7.42 年-1.20 年
BEP (% 满产)62.3%57.9%-4.3%

工程结论: 1000 t 三文鱼项目采用 bypass 拓扑相对 mainline: 泵综合节能近 30 kW (年电费节省 17 万) + CAPEX 工程总投资节省 1,346 万 + 综合成本省 1.29 元/kg + 全投资 IRR 提升 2.4 个百分点 + 静态回收期缩短 0.62 年。 这与国际主流大型三文鱼 RAS 多采用 side stream 拓扑的工程实际一致 (Linde SOLVOX / PR Aqua PPC / Global Seafood Advocate 2019)。

历史 4 工况成本基线 (v1.3 数据, 沿用)

工况气温源水v1.3 综合成本对照
青岛冬季-7°C4°C35.96 元/kg国内沿海三文鱼试点
海南冬季17°C22°C37.70 元/kg
海南夏季32°C27°C40.47 元/kg制冷成本上升
挪威 Bergen 冬0°C4°C36.30 元/kgAtlantic Sapphire 60-70 (含人工/苗)
📝 注 上述 4 工况基线是 v1.3 数据 (2026 年 4 月, mainline 拓扑)。 v1.3 挪威 Bergen 36.30 元/kg vs v1.7 jsdom 实测 36.70 元/kg 完美吻合 (偏差 +1%), 说明 v1.7 物理修复 (氧锥拓扑 / K-101 / 系统排水) 未改变综合成本量级, 仅在 bypass 拓扑下能省 1.29 元/kg。 SUMMARY 路线图已规划 v1.8 用 bypass + v1.7 字段重新生成 4 工况基线 JSON。

iRAS 不计算人工 / 苗种 / 管理费等"软成本" (合计 15-20 元/kg), 综合成本通常比文献项目实际偏低。 对照 Nordic Aqua 宁波 (中国境内唯一商业化运营的 RAS 三文鱼项目, AKVA Group 技术) 2024 年 Q2 实测生产成本 41 元/kg (EUR 5.23/kg) — iRAS 36.70 + 软成本 ~5 元/kg = 41.7 元/kg, 与 Nordic Aqua 实测完全吻合 ✓

目录

版本演进说明 (v1.0 → v1.7)
第一部分 RAS 基础
第 1 章 循环水养殖概述
第 2 章 RAS 工艺流程
第 3 章 水质参数与控制标准
第 4 章 氮循环与碳循环
第二部分 工程设计方法
第 5 章 养殖规模设计
第 6 章 污染物负荷计算
第 7 章 精确循环方程
第 8 章 系统体积建模 v1.3.x 增 S-101
第三部分 各处理单元设计
第 9 章 固液分离 (转鼓微滤机 RDF) v1.7 系统排水重构
第 10 章 生物过滤 (MBBR 生物滤池)
第 11 章 蛋白分离器
第 12 章 反硝化反应器
第 13 章 UV 紫外线消毒
第 14 章 高级氧化 (AOP)
第 15 章 增氧与脱气 v1.7 氧锥拓扑 + K-101
第 16 章 热负荷与温控 v1.2 重写
第 17 章 水泵与管道设计
第 18 章 电气与控制系统
第 19 章 生物安全与消毒策略
第 20 章 系统启动与调试
第四部分 经济评估
第 21 章 运行成本 (含 CAPEX 折旧)
第 22 章 投资概算 v1.7 旁路拓扑 CAPEX
第五部分 设计案例
第 23 章 案例: 三文鱼 1000 t/年 v1.7 数据
第 24 章 案例: 大菱鲆 500 t/年
第 25 章 案例: 对虾 100 t/年
第六部分 工程化与平台功能 (v1.3-v1.7)
第 26 章 多模块并联与 N+1 备用 v1.3 新章
第 27 章 鱼池配置与池深统一 v1.3 新章
第 28 章 全局气候与围护建模 v1.2 新章
第 29 章 双节点稳态热平衡求解 v1.2 新章
第 30 章 P&ID 工艺图自动生成 v1.3.x 新章
第 31 章 设备清单导出 v1.3.x 新章
第 32 章 4 工况成本对比基线 v1.3 新章
第 33 章 投资财务评价 (NPV/IRR + 8 张标准表) v1.4 新章
第 34 章 可研报告生成器 (16 章 + docx + 模板引擎) v1.5 新章
附录
附录 A 完整符号表
附录 B 全局热配置参数表
附录 C 设备 tag 命名总表
附录 D 鱼种参数表
附录 E 参考文献
附录 F 免责声明
第 一 部 分
RAS 基础

第 1 章循环水养殖概述

1.1 什么是循环水养殖 (RAS)

循环水养殖系统 (Recirculating Aquaculture System, RAS) 是一种将养殖用水经过物理、化学和生物处理后循环使用的集约化水产养殖技术。与传统的流水式养殖不同, RAS 仅需补充少量新水 (通常每天换水量 5-15%), 其余 85-95% 的水经处理后循环回用。

RAS 的核心思想是: 在一个封闭或半封闭的系统中, 通过一系列水处理设备, 持续去除鱼类代谢产生的废物 (氨氮、悬浮固体、溶解有机物、二氧化碳), 同时补充溶解氧, 从而维持水质在鱼类健康生长的最适范围内。

一个典型的 RAS 系统由以下核心组件构成: 养殖池 (鱼池/虾池)、固液分离设备 (转鼓微滤机 RDF)、生物过滤器 (MBBR 生物滤池)、消毒设备 (UV 紫外线)、增氧设备 (液氧系统或曝气系统)、以及脱气设备 (CO₂ 脱气塔)。

1.2 RAS 与传统养殖的对比

对比维度池塘养殖流水式养殖RAS 循环水养殖
用水量蒸发补水大量取水+排放仅补充 5-15%/天
养殖密度0.5-5 kg/m³10-30 kg/m³30-100 kg/m³
水质控制靠天然净化靠水量稀释精确工艺控制
占地面积小 (是池塘的 1/10-1/50)
环境影响面源污染排水污染零排放或极少
选址限制需靠近水源需靠近水源几乎不受限制
初始投资
运行成本中高 (电费/液氧)
产品质量波动大较稳定最稳定
生物安全高 (封闭)

1.3 RAS 的优势与挑战

优势

挑战

1.4 国内外发展现状

1.5 典型应用场景

场景典型品种密度 kg/m³水温 °C特点
冷水鱼工厂三文鱼、虹鳟50-8012-15淡水苗种+海水养成, 周期长
温水鱼工厂加州鲈、鳜鱼40-6024-28国内市场需求大, FCR 低
海水鱼工厂大菱鲆、石斑30-7016-27需海水/半咸水
鳗鲡工厂日本鳗、欧洲鳗40-6026-28粉料投喂
对虾工厂南美白对虾2-1528-31密度低但投饵率高
📌 v1.1 修订 大菱鲆密度从 v1.0 的 15/30/40 改为 25/50/60 kg/m³, 与莱州东方海洋等国内循环水项目实测一致。

第 2 章RAS 工艺流程

2.1 标准 RAS 工艺流程

一个标准的 RAS 系统, 水流按以下顺序循环:

鱼池 → 微滤机(RDF) → [蛋分] → 生物滤池(BF) → UV → CO₂ 脱气塔 → 增氧 → 鱼池
旁路: [反硝化] / [AOP]

其中方括号 [ ] 表示可选单元, 根据品种和水质需求决定是否安装。

📌 v1.1 CO₂ 脱气塔在 v1.1 中作为独立设备建模 (v1.0 隐藏在增氧节点)。CAPEX 含塔体 + 配套风机, OPEX 含风机 24h 电费。

在 iRAS 平台中, 上述流程被精确建模为主流串联 + 旁路并联的混合模型, 每个处理单元的去除效率独立计算, 沿程浓度逐步递减。

2.2 各处理单元的作用

处理单元去除目标原理典型效率
转鼓微滤机 (RDF)TSS (悬浮固体)60-90 μm 筛网物理拦截60-90%
蛋白分离器DOM + 细颗粒 TSS气浮泡沫分离, 可配臭氧DOM 40-70%
MBBR 生物滤池TAN (氨氮)硝化菌将 NH₄ 氧化为 NO₃60-90%
反硝化反应器NO₃ (硝酸盐)缺氧条件下还原 NO₃ 为 N₂50-80%
UV 紫外线细菌/病毒/寄生虫紫外线破坏 DNA99-99.9%
AOP 高级氧化Geosmin/2-MIB (土腥味)臭氧+UV 协同氧化80-95%
增氧塔/曝气DO (溶解氧)液氧溶解或空气曝气110-165% 饱和度
脱气塔CO₂ (二氧化碳)逆流空气吹脱60-85% (G:L=3-7)
📌 v1.1 增氧目标饱和度按鱼种类别推荐: 鲑科 165%, 海水温水鱼 140%, 温水淡水肉食鱼 110-125%, 罗非/对虾 100-115%。

2.3 淡水 vs 海水系统差异

第 3 章水质参数与控制标准

3.1 关键水质参数

3.1.1 氨氮 (TAN)

TAN 是鱼类蛋白质代谢的主要含氮废物, 通过鳃排出。TAN 在水中以离子态 NH₄ 和分子态 NH₃ 两种形式存在, 其比例受 pH 和温度影响。NH₃ 对鱼有强毒性, pH 越高、温度越高, NH₃ 占比越大。RAS 中 TAN 浓度通常控制在 0.5-2.0 mg/L, NH₃ 应低于 0.02-0.05 mg/L。

3.1.2 亚硝酸盐 (NO₂-N)

硝化中间产物, 对鱼毒性仅次于 NH₃。RAS 中 NO₂ 应低于 0.5 mg/L。持续偏高说明硝化菌群不完整 (Nitrobacter 不足), 通常出现在系统启动初期。

3.1.3 硝酸盐 (NO₃-N)

硝化最终产物, 毒性较低, 但在闭合循环中持续积累。长期高浓度 NO₃ (> 100-200 mg/L) 抑制鱼生长。反硝化或换水是控制手段。

3.1.4 溶解氧 (DO)

鱼类和硝化菌都需要充足 DO。鱼池 DO 通常要求 > 6 mg/L (冷水鱼) 或 > 5 mg/L (温水鱼)。纯氧系统可将 DO 提升至 110-165% 饱和度 (超饱和), 为高密度养殖提供保障。

📌 v1.1 DO 阈值改为 max(satRatio × DO_sat, DO_abs) 双判断, 鱼种内置 satRatio 和 abs 阈值。例如鲑科 0.75+6.0 mg/L, 罗非 0.60+4.0 mg/L。避免了 v1.0 单值 6.0 在罗非高温场景下过严的问题。

DO 饱和度采用 Benson & Krause (1984) 公式精确计算, 详见第 7 章。

3.1.5 二氧化碳 (CO₂)

鱼呼吸产生 CO₂, 在封闭系统中积累会导致 pH 下降、鱼应激。RAS 中 CO₂ 通常控制在 < 15 mg/L (海水), < 20 mg/L (淡水)。脱气塔通过逆流空气吹脱去除。

3.1.6 悬浮固体 (TSS)

包括残饵、粪便、脱落的生物膜碎片。高 TSS 会堵塞鱼鳃、消耗溶氧、滋生有害菌。RAS 中 TSS 通常控制在 < 15 mg/L。

3.1.7 溶解有机物 (DOM)

包括溶解态的蛋白质、脂肪酸、腐殖酸等。DOM 不能被 RDF 拦截, 也不能被 BF 有效去除, 在闭合循环中会积累, 导致水色发黄、起泡、产生土腥味。蛋白分离器和 AOP 是 DOM 的主要去除手段。

3.2 各品种耐受范围

参数单位三文鱼加州鲈/鳜罗非对虾石斑大菱鲆鳗鲡
TANmg/L< 1.0< 1.5< 2.0< 1.5< 1.0< 1.0< 1.5
NO₂-Nmg/L< 0.3< 0.5< 1.0< 1.0< 0.5< 0.5< 0.5
NO₃-Nmg/L< 80< 100< 200< 100< 80< 80< 100
DOmg/L> 6> 5> 4> 5> 5> 6> 5
CO₂mg/L< 15< 20< 20< 15< 15< 15< 20
TSSmg/L< 10< 15< 25< 15< 15< 10< 15
pH6.5-7.56.8-7.86.5-8.57.5-8.57.5-8.27.5-8.27.0-8.0
温度°C10-1622-2825-3228-3224-3014-1824-28
盐度0→3000-155-3525-3525-350-5

3.3 水质恶化的后果

RAS 水质管理的核心是"稳定"——参数在安全范围内的波动远比绝对值更危险:

第 4 章氮循环与碳循环

理解 RAS 中的氮循环是工艺设计的基础。所有氮素来源于饲料蛋白, 经过一系列生物化学转化, 最终以 N₂ 气体形式离开系统。

4.1 饲料蛋白到氮的转化

鱼类摄食后, 饲料蛋白被消化吸收, 其中:

iRAS 中 TAN 产生量的计算公式:

TAN (kg/d) = 日投饵量 × 蛋白含量 × TAN 系数

其中 TAN 系数 (tanCoef) 默认 0.092 (Timmons & Ebeling 2010)。海水肉食鱼 tanCoef 略高 (大菱鲆 0.092, 石斑 0.110, 鳜 0.100)。

4.2 硝化作用

硝化是 RAS 中最关键的生物过程, 由两类自养菌在有氧条件下完成:

NH₄ + 1.5 O₂ → NO₂ + 2H + H₂O (Nitrosomonas)
NO₂ + 0.5 O₂ → NO₃ (Nitrobacter)
总反应: NH₄ + 2O₂ → NO₃ + 2H + H₂O

从工程角度, 硝化的三个关键化学计量关系:

在 iRAS 中, 硝化耗氧由 BF 粗孔曝气独立供应, 不占用鱼池主增氧的额度。这一O₂ 按位置分配原则避免了传统计算中常见的氧气重复计算问题。

📌 v1.5+ 公式订正
稳态下 NO₃ 产生速率 ≈ TAN 产生速率 × 1.0 (近似 100% 硝化), 不是 × η_bio (70%)。 物理原因: 稳态质量守恒 — 鱼产 TAN = BF 累积硝化 + 换水排出 TAN, 换水排出通常 < 1% 鱼产 TAN, 所以"全部 TAN 都变 NO₃"在工程精度内成立。旧版本 (v1.0-v1.4) 用 × η_bio 是历史 bug, 导致 NO₃ 稳态浓度低估 30%, 反硝化设备容积/甲醇用量偏小。

4.3 反硝化作用

反硝化在缺氧条件下, 由异养菌利用外加碳源 (甲醇) 将 NO₃ 还原为 N₂ 气体:

6 NO₃ + 5 CH₃OH → 3 N₂↑ + 5 CO₂ + 7 H₂O + 6 OH

反硝化的关键化学计量:

iRAS 中反硝化作为旁路设计——分流 5-15% 的循环水进入缺氧反应器, 处理后汇回主流。反硝化可以同时降低 NO₃ 浓度和减少碱度消耗。

📌 v1.1 反硝化进水 NO₃ 取 BF 出口浓度 (P&ID 规范), 不再用鱼池浓度。影响实测约 11.5%。

4.4 碱度消耗与回收

碱度净消耗 = 硝化消耗 − 反硝化回收 − 换水补偿
= TAN_硝化 × 7.14 − NO₃_反硝化 × 3.57 − 换水量 × 源水碱度

如果净消耗 > 0, 需要投加碱源维持 pH:

NaHCO₃ 投加量 (kg/d) = 碱度缺口 / 0.595

其中 0.595 是 NaHCO₃ 到 CaCO₃ 的换算系数 (50/84 = 0.595)。

📌 v1.1 换水补偿用 V_total (系统总水体), 不是 V_tank。换水量自动 +50% 左右, 碱度补偿增大, NaHCO₃ 投加略减。

反硝化不仅降低 NO₃, 还能回收碱度、减少 NaHCO₃ 用量。一个设计良好的反硝化支路可以减少 30-50% 的碱度投加量。

4.5 碳酸盐平衡: pH / CO₂ 物种分布 / 非离子化氨 v1.9 新增

鱼池 pH 不是一个独立的自由参数 — 它由碳酸盐系统确定。总碱度 TA (由 NaHCO₃ 投加维持设定值) 与溶解 CO₂ (来自呼吸 + 硝化产) 一起决定 pH、CO₂ / HCO₃⁻ / CO₃²⁻ 物种分布、非离子化 NH₃ 占比 (剧毒) 以及方解石饱和度 Ω。iRAS v1.9 对每个阶段调用专门的碳酸模块求解此平衡。

📌 v1.9 物理建模升级
v1.7 之前 iRAS 只输出碱度需求 (基于 7.14 g CaCO₃/g N 硝化耗碱), 不计算 pH / NH₃ 物种 / 缓冲强度。这使得工具无法识别"低 pH 抑制硝化"或"高 NH₃ 毒性"等关键工况问题。 v1.9 引入完整碳酸盐平衡求解, 平衡常数全部来自国际权威数据库 PyCO2SYS 1.8.3, 经 jsdom 验证中英文版输出完全一致 (max|ΔpH| = 1e-5)。

4.5.1 控制方程

给定 TA 和溶解 CO₂, 碳酸模块求解质子平衡方程并反推完整物种分布:

TA = [HCO₃⁻] + 2[CO₃²⁻] + [B(OH)₄⁻] + [OH⁻] − [H⁺]

解离常数:

内部求解在 SWS (Seawater Scale) 标度上进行, 但报告的 pH 转换到 free (自由 / 浓度) 标度。两者关系:

S = 0 (淡水) 时: pH_free = pH_sws = pH_nbs (NBS 探头读数)
S > 0 (海水) 时: pH_free > pH_sws (差异 0.05-0.15 pH 单位)

用 free 标度报告 pH 的优势: 避免 fH 修正导致的"NBS 探头在淡水中读数偏高 0.15"问题, 直接对应淡水探头读数。

4.5.2 非离子化氨 NH₃ (剧毒)

溶解氨以两种形式共存: 离子化 NH₄⁺ (无毒) + 非离子化 NH₃ (剧毒)。两者占比由 pH 和 T 决定:

fNH3 = [NH₃] / TAN = K_a / (K_a + [H⁺])

典型数值 (T = 15°C 三文鱼工况):

pHfNH3 (%)TAN = 1 mg/L 时 NH₃-N
6.50.05%0.5 μg/L
7.00.16%1.6 μg/L
7.50.50%5.0 μg/L
8.01.55%15.5 μg/L
8.54.66%46.6 μg/L

工程意义: 同样的 TAN 浓度, pH 从 7.0 升到 8.0, NH₃ 浓度增加近 10 倍。这是为什么高 pH 工况 + 高 TAN 工况是急性死亡风险, 而 v1.7 之前 iRAS 无法识别此风险。

鱼种长期暴露 NH₃-N 上限 (基于 Fivelstad / Calabrese 等):

苗种期对 NH₃ 更敏感, 通常取上述限值的 50%。

4.5.3 缓冲强度 β 与方解石饱和度 Ω

Van Slyke 缓冲强度 β 定义为: 维持给定 pH 变化所需的碱度变化量。

β = dTA / d(pH) (mg/L CaCO₃ per pH 单位)

工程意义: β 越大, 碱度抗酸冲击越稳, 投喂后 pH 波动越小。典型值:

工况β (mg/L CaCO₃ per pH)评价
淡水 苗种期 (TA=150)350中等
淡水 罗非鱼成鱼 (TA=150)350中等
海水 三文鱼成鱼 (TA=150)477较强 (含硼酸盐)
商业 Atlantic Sapphire 风格 (TA=237)700+强 (推荐高密度工况)

方解石饱和度 Ω (仅海水):

Ω = [Ca²⁺][CO₃²⁻] / K_sp

Ω = 1 为饱和, Ω > 1 过饱和有 CaCO₃ 沉积风险, Ω < 1 不饱和。海水 RAS 典型工况 Ω = 0.5-1.5, 接近饱和但避免沉积堵塞设备

4.5.4 两种 pH 控制策略

iRAS v1.9 在每个阶段提供两种控制模式 (UI 工艺参数面板"pH 控制"下拉选择):

模式输入输出适用场景
持碱度 (alk)碱度目标 (mg/L CaCO₃, 默认 150)pH 由碳酸盐平衡涌现淡水中小项目, 运维简单, 碱度滴定即可
持目标 pH (ph)目标 pH (默认 7.0)工具反推所需碱度设定值 (Alk req)大型海水 RAS, 业主想精确控制 pH (实测 7.0-7.2 防 NH₃ + 防 CaCO₃ 沉积)

商业实测数据 (Good et al. 2018) 显示大型 post-smolt 三文鱼 RAS 实际碱度约 237 mg/L CaCO₃ (远高于淡水默认 150), 这是因为高密度海水工况需要:

用户可在工艺面板把碱度目标改成 200-250 以匹配国际商业工程实际。

4.5.5 与 PyCO2SYS 的验证

iRAS 碳酸模块与国际标准实现 PyCO2SYS 1.8.3 (opt_k_carbonic=14) 进行严格验证, 覆盖盐度 0-35 与温度 12-28°C 全范围:

验证项iRAS vs PyCO2SYS 偏差
正向 pH (free / NBS 标度)max|ΔpH| = 1×10⁻⁵
反向碱度求解 (solveAlk)0.08 μmol/kg
NH₃ 占比 fNH31.5×10⁻¹⁶ (机器精度)
方解石饱和度 Ω0.006%

这是实验室级精度, 可用于学术研究和工程设计。

4.5.6 稳态 vs 瞬态

本章求解的 pH / NH₃ / CO₂ 都是日均稳态 (慢性暴露) 值, 不是投喂瞬时峰值。投喂后短期瞬态 (2-4h 窗口) 由:

与 DO 不同 — DO 必须按投喂瞬时峰值选型 (避免急性窒息), 而 CO₂/pH/NH₃ 按日均稳态合规即可 (慢性暴露阈值)。

📌 工程使用建议
1. 三文鱼项目: 海水成鱼期建议碱度提升到 200-237 mg/L CaCO₃ (匹配 Atlantic Sapphire / Salmon Evolution 商业实际), 用"持目标 pH"模式控 pH = 7.0
2. 温水鱼 + 高密度 (罗非鱼 70+ kg/m³, 对虾): pH 容易偏高 (7.5-7.8), NH₃ 风险高, 务必检查工具输出 NH₃-N 是否超鱼种限值
3. 苗种期: 限值取成鱼期的 50% (鱼种敏感期), 工具自动收紧
4. 淡水低密度: 默认碱度 150 mg/L 够用, 用"持碱度"模式
5. 反硝化支路: 部分恢复硝化耗碱, 减少 NaHCO₃ 投加 (v1.6 加入), 但不改变 pH 平衡 (碱度变化经碳酸盐重新分配)
第 二 部 分
工程设计方法

第 5 章养殖规模设计

RAS 工程设计的第一步是从年产量目标反推系统规模。这决定了后续所有水处理设备的选型依据。

5.1 Little's Law 与稳态建模

iRAS 采用连续错峰养殖模式建模: 养殖场分批错峰投苗, 使每天都有成鱼出塘。稳态下, 系统内每个阶段同时存在不同大小的鱼群。

适用 Little's Law (排队论基本定理):

L = λ × W

其中 L = 系统内平均数量 (尾), λ = 到达速率 (尾/年), W = 停留时间 (年)。

关键: Little's Law 应用于"尾数"而非"质量"。一尾三文鱼从 0.1g 长到 5000g, 体重变化 5 万倍; 如果用质量, 会严重高估上游阶段的存塘量。

📌 v1.1 修正了稳态在塘公式, 加入 mortality 衰减: N_steady = N_inflow × (months / 12) × (1 − mort/2)。mortality 修正反映线性死亡假设下的平均尾数 (高 mort 阶段如苗种期偏高 8-11%)。

5.2 尾数流守恒与损耗率递推

多阶段 RAS 的尾数从最后阶段向前递推, 每一级通过损耗率放大:

年出塘尾数 = 年产量 (kg) / 出塘个体重 (kg)
本级进塘尾数 = 下级进塘尾数 / (1 − 本级损耗率)

损耗率包含死亡、淘汰、逃逸, 苗种期通常最高 (10-30%), 成鱼期最低 (2-5%)。

5.3 投饵量: 三套口径 v1.1 重大修正

v1.0 用 feedDesign = feedMax × peakFactor 算所有"日产污染物", 偏大约 17%。v1.1 严格分离三个时间尺度:

① 日峰值投饵量 (feedMax) — 负荷计算基准

feedMax = biomass × feedRate (kg/d)

biomass 是稳态存塘量, feedRate 是阶段末最大日投饵率。污染物日产量 (TAN/TSS/DOM/CO₂) 都按 feedMax 算 (不含 peakFactor)。

② 小时峰值 (xxxPeakHourly) — 设备选型基准

xxxPeakHourly = xxxDaily / 24 × peakFactor

peakFactor 物理意义: 代谢污染物小时峰值/全日均速率, 默认 1.5。快动力学设备 (鱼池增氧、BF 风机) 按 xxxPeakHourly × safety 选型。

③ 全周期日均 (feedAvg) — OPEX 年化基准

feedAvg = yearGrowth × FCR / 365 (kg/d)

其中 yearGrowth 是阶段年鱼肉增重, FCR 是阶段饲料系数。

峰均比校验 峰均比 = feedRate / (feedAvg/biomass)。如果峰均比 < 1.1, 设备可能选型偏小, iRAS 会发警告。

第 6 章污染物负荷计算

确定了投饵量后, 即可计算各类污染物的日产生量。v1.1 严格区分日量、小时峰值、日均三套字段。

6.1 TAN 产生量

TAN_Daily = feedMax × protein × tanCoef
TAN_PeakHourly = TAN_Daily / 24 × peakFactor
TAN_DailyAvg = feedAvg × protein × tanCoef

其中 tanCoef 默认 0.092 (Timmons & Ebeling 2010), 海水肉食鱼略高 (0.10-0.11)。

6.2 TSS / DOM 产生量

TSS_Daily = feedMax × 0.30
DOM_Daily = feedMax × 0.10

6.3 O₂ 需求三部分 v1.1 按位置分配

v1.1 把鱼生理耗氧按物种 + 温度自动算 (取代 v1.0 的 0.25 一刀切)。

O₂ 来源计算公式供应设备说明
① 鱼生理耗氧feedMax × o2FishFactor鱼池主增氧物种特性 + Q10 自动算
② 硝化菌耗氧TAN_Daily × 4.57BF 粗孔曝气NH₄→NO₃ 化学计量
③ 异养菌耗氧feedMax × o2DOMfactorBF 粗孔曝气v1.1 默认 0.10 (Boyd 2018)
关键设计原则 鱼池主增氧 (①) 和 BF 曝气 (②+③) 是两套独立系统, 不重叠。v1.0 之前常把三项合并算总量, 然后全部由鱼池增氧供应, 导致 BF 缺氧或鱼池增氧过度设计。

6.4 CO₂ 产生量

CO₂_Daily = o2FishDaily × 1.375

其中 1.375 = 44/32, 是 CO₂ 和 O₂ 的摩尔质量比 (RQ ≈ 1)。CO₂ 只按鱼虾生理耗氧计算, 不含硝化耗氧 (硝化反应不产生 CO₂)。

第 7 章精确循环方程

这是 iRAS 的核心计算引擎。传统 RAS 设计常用近似加法公式 C = Δ / Ση 估算稳态浓度, 但当系统包含多个串联和旁路处理单元时, 加法公式会产生显著误差。iRAS 采用精确乘法循环方程, 确保质量守恒严格成立。

7.1 质量守恒推导

考虑一个稳态循环系统: 水从鱼池流出, 经过 n 个处理单元后回到鱼池。每个处理单元有各自的"通过率" (1 − η)。

设:

稳态条件:

C_out = C_in + Δ (鱼池内产生增量)
C_in = C_out × R (经全部处理后回到 C_in)

联立求解:

C_out = Δ / (1 − R) (鱼池出水, 设备选型依据)
C_in = Δ × R / (1 − R) (回水, 最低浓度)

7.2 总通过率 R

R = Π(1 − η) × (1 − dilution)

R 越小 (总去除率越高), C_out 越低, 水质越好。

7.3 串联 vs 旁路建模

串联通过率 (RDF / BF / CO₂ 脱气)

全流量通过, 注意是乘法:

串联通过率 = (1−η_rdf) × (1−η_bio) = 0.25 × 0.30 = 0.075

即总去除率 92.5%。加法 0.75 + 0.70 = 1.45 > 100% 物理上不可能。

旁路通过率 (蛋分 / AOP / 反硝化)

以蛋分为例, 旁路流量比 α=20%, 单次去除 η=70%:

系统等效去除率 η_sys = α × η_single = 0.20 × 0.70 = 14%

物理含义: 80% 的水不经处理 + 20% 的水去除了 70% → 混合后整体去除 14%。

7.4 闭合验证

精确循环方程的优美性质: 沿程浓度递减后, 最终回水浓度恰好等于 C_in。iRAS 在每次计算后内部验证此闭合条件 (误差 = 0)。

7.5 各参数的 R 计算路径

参数处理路径R 表达式
TANBF 硝化 + 换水(1−η_bio)(1−d)
TSSRDF + 蛋分(旁路) + BF + 换水(1−η_rdf)(1−η_skim)(1−η_bioTSS)(1−d)
DOM蛋分(旁路) + AOP(旁路) + 换水(1−η_skim_sys)(1−η_aopDOM_sys)(1−d)
NO₃反硝化(旁路) + 换水(1−η_denitri_sys)(1−d)
CO₂脱气塔 + 换水(1−η_co2)(1−d)
GeosminAOP(旁路) + 换水(1−η_aop_sys)(1−d)
📌 v1.1 修复 AOP DOM 系数: v1.0 误把 90% geosmin 系数当 DOM 去除率, v1.1 用 aopDOMRemoval 默认 5% (可调 5-15%)。

7.6 dilution 敏感性 v1.1 重要发现

dilution 在不同参数中的影响差异巨大:

参数 / 工艺配置η_totaldilution 是否主导?dilution 敏感度
TAN (BF 70%)70%< 0.1%
TSS (RDF + 蛋分 + BF)82%< 0.1%
CO₂ (脱气塔 65%)65%< 0.1%
NO₃ (反硝化 6%)6%边缘~1%
NO₃ (无反硝化)0%缩放为 V_tank/V_total
Geosmin (无 AOP)0%缩放为 V_tank/V_total
DOM (无蛋分)0%缩放为 V_tank/V_total

当 η = 0 时, R = 1 − dilution, 1−R = dilution, C_out = Δ/dilution。dilution 翻倍, 浓度减半。这是 v1.1 系统体积建模 + 换水率口径修正的关键意义。

第 8 章系统体积建模 v1.1 新章 + v1.3.x 增 S-101

📌 v1.1 本章是 v1.1 全新内容。v1.0 把整个系统简化为"鱼池水体 V_tank", 忽略了工艺水体 (BF/RDF/管路等)。v1.1 引入完整 V_total = V_tank + V_process 建模。v1.3.x 在 V_process 中加入集水池 V_sump。

8.1 为什么要建模工艺水体?

传统 RAS 设计常把整个系统简化为"鱼池水体", 忽略了:

这些工艺水体合计 V_process 通常占 V_tank 的 30-60%, 是不可忽略的。

8.2 V_process 各部分计算公式

工艺单元计算公式物理依据
V_RDF v1.5+ 修正Q × rdfHRT_min / 60RDF 集水池 HRT, 默认 rdfHRT_min = 0.5 分钟 (v1.6 中位, 兼顾紧凑与缓冲); 紧凑型 0.3 (YUTANK/Hydrotech 工业实测 18s); 设计院传统 1.5-2 min (含外置缓冲池)
V_BF (含填料水) v1.5+ 修正bioVolume / fillRatiobioVolume 是填料净体积, fillRatio 默认 0.5 (Rusten 2006 工业惯例): V_BF = 2 × bioVolume; fillRatio 0.7 (高密度): V_BF = 1.43 × bioVolume
V_skimskimmerVolume蛋分塔水容积 (已是水量)
V_UV5 m³ (固定)UV 流通式, 内部水量极小
V_AOPaopVolumeAOP 反应池水容积 (已是水量)
V_denitridenitriVolume反硝化反应器
V_co2 v1.5+ 修正degasArea × co2VolFactor脱气塔水容积, 默认 co2VolFactor = 0.7 (集水 0.5m + 填料湿水 0.2m, 国内中型项目); 可调 1.2 (1m 集水 + 0.2m 填料水, 鲑鱼场深塔)
V_aeratorQ × 0.5 / 60增氧锥 HRT ≈ 30s
V_sump v1.3.xclamp(Q/30, 20, 200)集水池, 工程惯例
V_pipesV_tank × 5%管路+集水井+缓冲池, 工程经验
📌 v1.5+ V_BF 推导
RAS 体积建模时, BF 池总容积 (含填料 + 水) 才是稳态稀释的有效水量, 不只是填料净体积。 若 fillRatio = 50% (典型 MBBR), 填料净体积 100 m³ 对应池总容积 200 m³ (因 50% 是填料、50% 是水间隙)。 旧版 × 1.3 是经验估计 (假设填料 35% + 水 65%), 但工业上 fillRatio 是设计变量, 应该用代码统一。

8.3 V_total 与稳态浓度

总系统体积:

V_total = V_tank + V_process

对稳态浓度 C_out = Δ/(1−R) 影响:

因此 V_total 影响稳态浓度的程度取决于参数处理工艺:

8.4 换水率口径 v1.1 改为 V_total

用户填的"exchangeDaily 5%/d"在不同口径下含义不同:

口径含义实际换水量使用方
V_tank (v1.0)鱼池每天换 5%V_tank × 5%国内传统养殖户
V_total (v1.1)系统总水体每天换 5%V_total × 5% (多 50%)Timmons 教科书 / 国家海水鱼体系 / 大型工业项目
📌 v1.1 改为 V_total 是物理一致性优先。v1.0 老方案加载到 v1.1 后, 换水量自动 +50% 左右, 仅靠换水稀释的参数 (无反硝化的 NO₃, 无 AOP 的 Geosmin) 浓度自动下降 ~37%。这是物理修正后的真实值, 不是 bug。

8.5 三文鱼成鱼 1000t/yr V_process 实例

工艺单元体积 (m³)
鱼池水体 V_tank3,385
BF 含水570
CO₂ 脱气350
RDF 集水338
集水池 S-101 v1.3.x200
管路 + 缓冲169
反硝化129
增氧锥85
蛋分塔74
AOP 反应池17
UV5
V_process 合计1,938 (V_tank 的 57%)
V_total = V_tank + V_process5,323
第 三 部 分
各处理单元设计

第 9 章固液分离 (转鼓微滤机 RDF) v1.7 系统排水重构

9.1 原理

转鼓微滤机 (RDF) 是 RAS 中的第一道防线, 通过 60-90 μm 不锈钢筛网物理拦截水中的悬浮固体。工作原理: 养殖水从转鼓内侧流向外侧, 固体颗粒被拦截。当筛网两侧压差达到设定值时, RDF 内嵌反冲机构自动启动反冲洗 — v1.7 起反冲泵不再作为独立设备列项, 因为工程实际中所有商用 RDF (Hydrotech / Faivre 等) 整机出厂均含反冲泵, 反冲水绝大部分经污泥沉淀池回循环 (Sharrer 2010 实测 95-98% 清液回流)。

9.2 设计参数

参数典型值说明
筛网孔径60-90 μm越小去除率越高, 但易堵塞
TSS 去除率60-90%iRAS 默认 75%
处理流量循环流量 × safetyRDFv1.1: 默认 1.10
RDF 集水池 HRT0.1-5 miniRAS v1.6 默认 0.5 min (中位)

9.3 计算公式

RDF 处理流量 = 循环流量 × safetyRDF (m³/h)
RDF 集水池 V = 流量 × rdfHRT_min / 60 (m³)

9.4 系统排水 (v1.7 重构)

v1.6 之前, iRAS 工艺面板有 "RDF 反冲(%) 2%" 输入, 但 4870 m³/d ÷ V_total = 101%/d 物理荒谬 (远超合理换水率 1-5%/d)。问题根源: 把"反冲管瞬时分流比"误当作"穿越系统边界的排水比"。

v1.7 全面清理这个语义混淆:

系统排水流量 (m³/d) = V_total × exchangeDaily (与"换水率"定义一致, 不再依赖反冲)
系统排水 TSS (kg/d) = load.tssDaily (质量守恒, 含污泥固相)
系统排水 TAN_eff (mg/L) = C_tank (溶解态稳态浓度)
系统排水 NO₃_eff (mg/L) = C_tank (溶解态稳态浓度)
⚠️ TSS 口径说明
iRAS 输出的"系统排水 TSS"是沉淀池入口浓度 (含污泥, 质量守恒口径), 不是污水管口浓度。 1000 t 三文鱼工况下 TSS ≈ 3257 mg/L, 是合理的沉淀池入口悬浮固体浓度 (典型 RAS 1000-5000 mg/L 范围)。 不可直接对照 GB 8978 (污水管口标准, 一级 70 mg/L)。沉淀池设计 + 上清液处理由污水 工程单独完成, iRAS 边界止于"沉淀池入口"。

9.5 常见故障与运维

故障现象可能原因处理措施
出水 TSS 偏高筛网破损 / 孔径变大更换筛网 (寿命 2-3 年)
反冲频繁启动进水 TSS 过高 / 筛网堵塞清洗筛网; 检查投饵量
排污水清澈反冲喷嘴堵塞 / 水压不足清洗喷嘴; 检查 RDF 内嵌反冲机构
转鼓不转电机故障 / 链条断裂检查电机和传动

运维要点: ① 每周检查筛网完整性; ② 每月清洗反冲喷嘴; ③ 海水系统注意筛网腐蚀。

第 10 章生物过滤 (MBBR 生物滤池)

10.1 硝化原理

MBBR 是 RAS 最关键的处理单元。HDPE 悬浮填料 (比表面积 500-800 m²/m³) 在曝气搅动下悬浮于水中, 表面附着硝化菌群, 将有毒的 TAN 氧化为低毒的 NO₃。硝化是缓慢的自养过程, 菌群建立需 4-8 周。

10.2 容积去除率 (VTR)

VTR₂₀ = 0.3-0.8 kg TAN/(m³·d) @20°C (iRAS 默认 0.5)
VTR(T) = VTR₂₀ × θ^(T−20) 其中 θ = 1.10
温度 °C修正系数VTR kg/(m³·d)说明
100.3860.193冷水鱼苗种
120.4670.234三文鱼苗种
150.6210.310三文鱼成鱼
160.6830.341大菱鲆全周期
201.0000.500基准
251.6110.805加州鲈/鳜
282.1441.072罗非鱼
302.5941.297对虾

盐度修正

f_sal = max(1 − 0.01 × S, 0.3)

盐度 30‰ 时硝化效率降至 70%。

10.3 填料体积计算

BF 填料体积 = TAN_Daily × safetyBio / [VTR₂₀ × θ^(T−20) × f_sal]

10.4 BF 粗孔曝气 v1.1

BF O₂ = TAN × 4.57 + feedMax × o2DOMfactor
风量 = O₂_Peak × safetyBFblow / (1.2 × 0.23 × SOTE) SOTE=8%
风机功率 = 风量(m³/s) × 50000 / (0.60 × 1000) (50kPa, η=0.60)
📌 v1.1 风机效率从 0.55 提到 0.60。o2DOMfactor 默认 0.10 (v1.0 偏低 0.05, Boyd 2018 实测)。

10.5 启动期管理

BF 启动需要 4-8 周, 启动期管理不当会导致 TAN/NO₂ 飙升, 大量死鱼。

启动步骤

注意事项

10.6 常见故障

故障原因措施
TAN 突然升高BF 过载/停电/pH 过低减少投饵; 补碱度; 检查曝气
NO₂ 持续偏高Nitrobacter 不足/pH<6.5补碱度; 降低负荷
填料沉底曝气不足/气管堵塞检查风机和气管
出水浑浊生物膜脱落正常现象

第 11 章蛋白分离器

11.1 原理

蛋白分离器利用泡沫分离原理去除 DOM 和细颗粒 TSS (<30μm)。空气或臭氧以微细气泡注入接触塔, 表面活性有机物吸附在气泡表面, 随泡沫溢出。

11.2 空气 vs 臭氧

参数空气模式臭氧模式
DOM 单次去除率30-50% (默认 40%)60-80% (默认 70%)
TSS 单次去除率~40%~40%
额外效果杀菌+脱色+氧化
成本高 (臭氧发生器)
海水风险溴酸盐 BrO₃ (盐度>5‰需控制)
📌 v1.1 蛋分 TSS 去除率 skimmerTSSEff 开放为可调参数 (默认 40%)。

11.3 设计计算

系统等效去除率 = 流量比 × 单次效率 = 0.20 × 0.70 = 14%
蛋分容积 = 旁路流量 (m³/h) × HRT (min) / 60 × safetySkimmer

HRT 默认 2 min (Sanders Helgoland LE-315 商业实测 1.8-2.2 min, ScienceDirect 2023)。物理: DOM 吸附在气泡表面是秒级反应, HRT 主要保证泡沫稳定形成, 不是反应时间, 长 HRT 无收益。用户可在 1.5-4 min 范围调整。

臭氧需求 = feedMax × 0.010 g O₃/g 饲料

11.4 常见故障

故障原因措施
不产泡沫DOM 太低/气量不足正常 (新水); 增加气量
泡沫太湿气量过大/液位过高减小气量; 调低液位
臭氧泄漏密封不严检查; 加活性炭尾气破坏
溴酸盐过高海水中臭氧过量降低剂量; 监测 BrO₃

第 12 章反硝化反应器

12.1 原理

在缺氧条件下 (DO < 0.5 mg/L), 异养菌利用甲醇将 NO₃ 还原为 N₂ 气体。这是 RAS 中唯一能真正"去除"氮素的过程。

6 NO₃ + 5 CH₃OH → 3 N₂↑ + 5 CO₂ + 7 H₂O + 6 OH

12.2 设计计算

旁路流量 = 循环流量 × 流量比 (默认 10%)
系统去除率 = 流量比 × 单次效率 = 10% × 60% = 6%
VDR(T) = VDR₂₀(0.8) × 1.08^(T−20)
反硝化容积 = NO₃去除量 × safetyDenitri / VDR(T)
甲醇 = NO₃去除量 × 2.47 kg/d
碱度回收 = NO₃去除量 × 3.57 kg CaCO₃/d

12.3 常见故障

故障原因措施
NO₃ 不降甲醇不足/DO 过高增加甲醇; 检查缺氧 DO<0.5
出水发白臭味甲醇过量减少甲醇
pH 过高反硝化产碱过多减少流量比

运维要点: ① 甲醇按 NO₃ 实时调节; ② 甲醇属易燃危化品, 单独储存; ③ 反硝化器启动比 BF 快 (1-2 周)。

第 13 章UV 紫外线消毒

13.1 原理

UV-C (254 nm) 穿透微生物细胞壁, 破坏 DNA 胸腺嘧啶二聚体, 使其丧失复制能力。优点: 无化学残留, 不改变水质, 杀菌快。

13.2 Chick-Watson 模型

灭活率 = 1 − exp(−k × D)
病原体k (cm²/mJ)D=40 灭活率建议剂量
细菌 (E.coli, 弧菌)0.10-0.2099.7%30-40 mJ/cm²
鱼类病毒 (IPNV, IHNV)0.05-0.1086-98%60-100 mJ/cm²
原虫 (小瓜虫包囊)0.15-0.2599.7-99.99%30-40 mJ/cm²
水霉孢子~0.3099.999%20-30 mJ/cm²

13.3 UV 功率公式 v1.1

UV 功率 (W) = 处理流量 × UV 剂量 × 0.85 × safetyUV
📌 v1.1 系数从 0.75 改为 0.85 (USEPA UVDGM 中位); safetyUV 是设备额定功率裕度 (灯管老化), 不是剂量增加。

13.4 运维

① 灯管运行时间计数器记录, 寿命 8000-12000h;
② 石英套管每月用稀盐酸擦拭除垢;
③ UV 前 TSS < 10 mg/L 才能保证效果。

第 14 章高级氧化 (AOP)

14.1 原理

臭氧 + UV 协同产生羟基自由基 (·OH), 氧化分解 geosmin 和 2-MIB (土腥味物质, ng/L 级)。单独臭氧或 UV 效果有限, 协同效率 80-95%。

14.2 设计计算

旁路流量 = 循环流量 × 5%
AOP 容积 = 旁路流量 (m³/h) × HRT (min) / 60 (m³)

HRT 按 AOP 类型默认值:

物理: 羟基自由基 ·OH 寿命 < 1 ms, 反应是光化学秒级。Schrader (2010) 实测 1 min 接触, geosmin 去除 86-92%。用户可在 1-15 min 范围内手动调整。

臭氧 = 旁路流量 × 24h × 0.5 g/m³ / 1000 (kg/d)
📌 v1.1 修复 AOP DOM 系数: aopDOMRemoval 默认 5% (可调 5-15%), v1.0 误把 90% geosmin 系数当 DOM 去除率。

14.3 运维

① AOP 出水必须经活性炭或 UV 消除残留臭氧再回主流;
② 不是所有 RAS 都需要 AOP——仅品种对土腥味敏感时才需要。

第 15 章增氧与脱气 v1.1 大幅修正 + v1.7 氧锥拓扑+K-101

📌 v1.1 本章是 v1.1 重点更新, 特别是 CO₂ 脱气塔由 v1.0 的"扣 70% 写死"升级为 v1.1 的独立设备模型。

15.1 DO 饱和度 (Benson-Krause 1984)

DO_sat = exp(−139.34411 + 1.575701e5/Tk − 6.642308e7/Tk² + ...)
温度 °C淡水 DO_sat海水 30‰降幅
1011.299.32-17%
1510.088.39-17%
169.868.21-17%
209.097.62-19%
258.266.97-18%
287.836.62-18%
307.566.41-18%

15.2 纯氧 vs 空气

参数纯氧 (液氧+锥形塔)空气 (罗茨风机+微孔盘)
DO 目标110-165% 饱和度 (按物种)~95% 饱和度
适用密度≥ 40 kg/m³< 40 kg/m³
吸收效率锥形塔 80-90%微孔盘 SOTE 18%
运行成本液氧 1.2 元/kg仅电费
📌 v1.1 DO 目标饱和度按鱼种推荐: 鲑 165%, 海水温水 140%, 温水淡水 110-125%, 罗非/对虾 100-115%。

15.3 鱼池 DO 计算 v1.1

鱼池出水 DO = DO_sat × satTarget% − deltaDO_peak
deltaDO_peak = o2FishPeakHourly / flowM3h × 1000 (mg/L)

v1.1 用小时峰值评估 DO 降幅 (而不是 v1.0 的日均), 反映投喂后代谢尖峰。

15.4 CO₂ 脱气塔 v1.1 新增独立设备

CO₂ 脱气塔在 v1.1 作为独立设备建模。v1.0 没有此模型 (仅在增氧节点扣 70% 写死, 不存在独立的塔/风机/CAPEX/OPEX)。v1.1 用户首次看到此设备会有独立的设备投资和运行成本——这是新增的成本项, 不是"修正"。

15.4.1 设计参数

参数v1.1 默认依据 / 工业实测
G:L 气液比3国内中型项目 2-4, 鲑鱼场 5-7
风压2.5 kPa (可调)Delta Cooling Towers 实测 0.87-1.12 kPa
风机效率0.60罗茨风机 0.5-0.7, 离心风机 0.65-0.8
单次脱气率65%G:L=3 浅塔实际值

默认参数与 Summerfelt (2009) 文献区间 (5-15 kW per 1000 m³/h) 和 Delta Cooling Towers 工业项目实测对齐。

15.4.2 风机功率公式

风量 (m³/h) = 循环流量 × G:L
风机功率 (kW) = 风量/3600 × 风压(Pa) / 风机效率 / 1000

举例 (三文鱼成鱼 Q=10155 m³/h):

10155 × 3 / 3600 × 2500 / 0.6 / 1000 = 35 kW

全部 8 鱼种 CO₂ 风机合计约 600 kW, 落入 Summerfelt (2009) 文献区间。

15.4.3 稳态 CO₂ 浓度

CO₂_tank (mg/L) = CO₂_Daily (kg/d) × 10⁶ / (流量 (m³/h) × 24 × η_co2)

单位推导: CO₂_Daily [kg/d] × 10⁶ [mg/kg] = mg/d ÷ (Q × 24) [m³/d ÷ 1000=L/d? 注意] 等价等式: CO₂_tank (mg/L) = CO₂_Daily(kg/d) × 1000 (g/kg) × 1000 (mg/g) / [流量(m³/h) × 24 (h/d) × 1000 (L/m³) × η_co2] 化简: CO₂_tank (mg/L) = CO₂_Daily(kg/d) × 10⁶ / [流量(m³/h) × 24,000 × η_co2] 即代码中的: co2Tank = co2Daily × 10^6 / (flowM3h × 24 × 1000 × eta_co2)

海水阈值 15 mg/L, 淡水阈值 20 mg/L。脱气塔关闭时 η_co2 改为 co2FallbackEff (默认 5%)。

📌 v1.5+ 公式订正
旧手册写的 × 10 是单位换算错误 (kg/d → mg/L 必须 × 10⁶, 即 1 kg = 10⁶ mg)。 代码一直是对的, 仅手册书写错误。

15.5 氧锥拓扑 mainline vs bypass v1.7 新增

v1.6 之前 iRAS 在 UI 上已经有"主流路径 / 旁路(节能)"两个选项, 但 calcPumpPowerBreakdown / aerationSpec / OPEX 完全没区分 — 是"UI 假承诺"。v1.7 通过全网调研 (Linde SOLVOX cone / Pentair AES Speece / PR Aqua PPC / Global Seafood Advocate 2019) 确认: 国际主流大型 RAS 项目 "more often plumbed in side stream configuration"

15.5.1 两种拓扑的物理对比

对比项mainline (主流串联)bypass (旁路高压)
主流路径全部循环水过氧锥主流不过锥, 旁路 5-30% 流量过高压锥后回流
主泵附加压头+ 2-15 m (锥子压损)0 (主泵不过锥)
氧锥规格流量= 全循环流量 Q= Q × bypassRatio (小流量)
氧锥操作压力低压 10-21 psi (LHO)高压 1-3 bar (SOLVOX / PPC)
典型 DO 出锥130-170% 饱和度200-400% 饱和度
额外设备旁路泵 P-602 + 混合点 M-601
代表设备Pentair AES LHOLinde SOLVOX cone / PR Aqua PPC
适用规模中小型项目 (Q < 5000 m³/h)大型项目 (Q > 5000 m³/h)

15.5.2 v1.7 字段定义

字段默认钳制范围工程对照
o2ConeTopologymainlinemainline / bypass
o2ConeHeadLoss_m3 m[0, 15]LHO 2-3 / Speece 5-10 / SOLVOX 110 至 15
o2BypassRatio (%)10[5, 30]国际主流 8-15%
o2BypassPumpHead_m20 m[5, 40]玻璃钢 20 / SOLVOX 25-30
o2BypassPumpEta (%)60[30, 85]高压泵典型

15.5.3 主泵扬程公式

mainline: H_pump_total = H_main + o2ConeHeadLoss
bypass: H_pump_total = H_main (主泵不过锥, 旁路泵独立加压)

15.5.4 旁路泵 P-602 设计

Q_bypass = Q × o2BypassRatio
P_bypass (kW) = Q_bypass × ρ × g × o2BypassPumpHead_m / (3.6 × 10⁶ × o2BypassPumpEta)

P-602 加入设备清单 / CAPEX / OPEX / 热平衡, 1+1 备用模块化 (与主泵 P-101 同口径)。

15.5.5 经济对比 (1000 t 三文鱼实测)

mainline (12 m) bypass (9 m + 旁路 20 m) 主泵 P-101 610.4 kW 488.3 kW 旁路泵 P-602 0 92.2 kW ─────────────────────────────────────────────── 综合泵功率 610.4 kW 580.5 kW (-29.8 kW) 氧锥规格 Q = 10155 m³/h Q × 0.1 = 1015 m³/h 氧锥 CAPEX 355 万元 36 万元 (缩 10×) 设备 CAPEX 4,454 万 3,842 万 (-612 万) 工程总投资 9,800 万 8,454 万 (-1,346 万) 综合成本(元/kg) 36.70 35.41 (-1.29 元/kg) 全投资 IRR 14.90% 17.31% (+2.41%)

15.6 air 模式增氧 — K-101 鱼池增氧鼓风机 v1.7 新增

v1.7 之前 air 模式 (空气曝气) 只有 K-302 (BF 曝气风机), 缺鱼池主增氧风机。 温水鱼种 (罗非鱼 / 鳗鱼 / 对虾) 的大规模 RAS 设计无法用工具自动生成 — 而这是 80%+ 中国 RAS 项目的实际场景。

v1.7 加入 K-101 鱼池增氧鼓风机:

参数典型值说明
设备类型罗茨风机 + 微孔盘管曝气SOTE 18-22% (Wagner 2010)
风量基准o2FishPeakHourly × safety / SOTE按鱼池峰值氧需反推
风压15-25 kPa含曝气盘水深 + 管阻
风机效率60%罗茨风机典型
备用1+1 备鱼池主增氧故障即鱼窒息, 必须双机轮换
P&ID tagK-101与 K-302 (BF 风机) 区分
风量 Q_air (Nm³/min) = o2FishPeakHourly × safety / (60 × SOTE × 0.28)
功率 P_K101 (kW) = Q_air/60 × P_kPa × 1000 / 60% / 1000

K-101 加入 modulePerSpec (与 K-302 同口径 1+1 备), P&ID / 设备清单 / 报告章节 / 财务测算 全链路同步。

15.7 常见故障

故障原因措施
DO 偏低液氧用完 / 增氧塔堵塞检查液氧; 清洗增氧塔
CO₂ 偏高脱气塔风量不足增加风量; 清洗填料
液氧蒸发过快储罐真空层失效联系厂商检修
微孔盘气泡变大盘面堵塞 / 老化酸洗或更换
bypass 模式 DO 不均匀M-601 混合点设计不当 / 旁路比过高检查混合管段长度; 旁路比降到 10-15%
K-101 跳闸 (air 模式)1+1 备用未启动 / 双机故障立即手动切换备用; 鱼池含氧 15-30 min 内告急

第 16 章热负荷与温控 v1.2 完全重写

📌 v1.2 重写 v1.1 用单点能量平衡估算热泵负荷, 忽略车间空气节点。v1.2 引入双节点 Picard 迭代同时求解空气热平衡 + 空气湿平衡。本章给出基础理论框架, 详细求解过程见第 29 章。

16.1 项目级常量管理 (v1.2)

v1.1 把热平衡参数(气温/源水温/U 值/换气率)放在每个阶段的 proc 字段里, 用户必须在每个阶段重复设置。 v1.2 把所有热平衡参数挪到顶部全局配置区 globalThermalConfig(项目级常量影响所有阶段)。具体参数表见第 28 章和附录 B。

16.2 补水温差

补水量 = V_total × 日换水率 (m³/d) [v1.1: V_total 口径]
热负荷 = V × 1000 × 4.18 × ΔT / 3600 (kWh/d)

16.3 板式换热回收

换热后温度 = T_source + η_HX × (T_target − T_source)
η_HX = 75% (iRAS 默认)

16.4 热泵补齐 v1.1 制热/制冷分离

📌 v1.1 制热 COP_heat 默认 4.0, 制冷 COP_cool 默认 3.0 (v1.0 一刀切 3.5)。热带降温场景下 v1.0 电费低估 30-40%。

16.5 双节点稳态求解概述 v1.2

v1.1 假设池热散到一个无限大"室外", 实际上车间空气会被池水加热到比室外高 5-12°C, 导致池散热小于"对室外气温散热"; 同时湿度积累让蒸发率下降。v1.2 引入双节点稳态:

节点 1 (水体热平衡):
Q_metabolism + Q_equipment − Q_pool_loss − Q_evap × λ = 热泵供热/制冷 (Q_HP)
其中:
Q_pool_loss = UA_pool × (T_w − T_room)
Q_evap × λ = 蒸发潜热散热, 计入水体(蒸发夺走水体显热)

节点 2 (空气热平衡):
UA_pool × (T_w − T_room) =
UA_envelope × (T_room − T_air)
+ UA_floor × (T_room − T_ground)
+ UA_vent × (T_room − T_air)
− Q_internal_air

节点 3 (空气湿平衡):
m_evap × cover = m_air × ACH × (W_room − W_air_outdoor)

迭代变量: T_room (室温), W_room (室内空气含湿量)

详细模型、求解算法和输出字段见第 29 章。

📌 v1.5+ 公式订正
蒸发潜热 Q_evap × λ 从空气节点移到水体节点。物理: 蒸发夺走的是水体显热 (水分子带走能量), 空气节点只接收 潜热增量 (W_room 上升), 不直接接收 λ 项。 旧手册写在空气节点是 v1.2 推导笔误, 代码一直把 λ 项放在水节点的能量平衡里 (正确)。

16.6 热泵分类型 v1.2

类型COP_heatCOP_cool源温适用
空气源 air3.5-4.02.5-3.0T_air常规, 气温 > -10°C
水源 water4.5-5.53.5-4.5T_source有大量地下水/河水
地源 ground4.0-5.03.5-4.5T_ground (≈年均)有埋管空间
海水开式 seawater_open4.5-5.53.5-4.5T_seawater沿海/海上养殖, 钛板换热

16.7 工程注意

第 17 章水泵与管道设计

17.1 泵功率公式

P (kW) = 9.81 × Q (m³/h) × H (m) / (3600 × η)

17.2 iRAS 6 路泵 v1.1 默认

流量默认扬程 m默认效率说明
主循环泵循环流量 Q1268%全流量, 功率最大
反硝化支路泵Q × 反硝化比1060%分流到缺氧反应器
UV 支路泵Q × UV 比1060%
AOP 支路泵Q × AOP 比1060%
蛋分支路泵Q × 蛋分比660%蛋分塔近距离
RDF 反冲泵Q × 反冲比855%间歇运行

主循环泵功率通常占全场电费 30-40%, 是最大单一用电设备。变频泵可节能 15-25%。

17.3 管径选择

管道类型推荐流速 m/s说明
吸水管0.8-1.2流速太高→气蚀
压水管1.5-2.5流速太高→水锤
排污管1.0-2.0需自清洗流速防沉积
溢流管0.5-1.0重力流
管径 D = √(4Q / πv) Q=m³/s, v=m/s, D=m

17.4 水头损失

组成典型值说明
静水头0.5-2 m鱼池到处理间高差
管道摩阻2-4 m管长和管径决定
设备阻力3-5 mRDF+BF+UV 各约 1-2m
安全余量1-2 m10-20%
合计8-15 miRAS 默认 12m

第 18 章电气与控制系统

18.1 在线水质监测

参数传感器精度价格位置
DO荧光法±0.1 mg/L5,000-15,000鱼池+BF 出水
pH玻璃电极±0.012,000-5,000鱼池+BF
温度Pt100±0.1°C500-1,000鱼池+源水
TAN离子选择±0.1 mg/L20,000-50,000鱼池 (可选)
ORP铂金电极±5 mV3,000-8,000蛋分 (臭氧控制)
液位超声波±1 cm1,000-3,000鱼池/集水井

18.2 PLC 控制逻辑

18.3 报警系统

报警项阈值响应
DO 低< 5 一级 / < 3 紧急启动备用风机/氧气瓶
pH 低< 6.5加 NaHCO₃; 减投饵
TAN 高> 2.0 mg/L减投饵; 检查 BF
温度异常偏离 ±2°C检查热泵
液位低低于设定补水; 检查泄漏
停电市电中断柴油发电机 30 秒启动
泵故障电流异常切换备用泵

必备: ① 柴油发电机 (全场 100% 负荷); ② UPS (PLC + 传感器); ③ 短信报警; ④ 每月演练应急切换。

第 19 章生物安全与消毒策略

RAS 的封闭特性天然具有高生物安全性, 但一旦病原体进入会因循环水快速传播。"防进来"比"治病"更重要。新建项目应在投产前建立完整的生物安全 SOP。

19.1 进水消毒

19.2 人员与工具

19.3 隔离与检疫

19.4 日常消毒

对象消毒剂浓度/方法频率
工具碘伏100 ppm 浸泡 10 min每次使用后
脚池次氯酸钠200 ppm 每周换液每天
车间地面过氧化氢3% 喷洒每周
鱼池 (空池)高锰酸钾50 ppm 浸泡 24h批次间
进水UV60 mJ/cm²持续

第 20 章系统启动与调试

新建 RAS 满负荷需要 8-12 周。急于满负荷投鱼是新手最常犯错误——BF 处理能力是逐步建立的, 超过 BF 能力 = TAN 暴涨 = 大量死鱼。

20.1 启动流程

阶段时间操作关键指标
设备调试第 0-1 周空载运行; 试压设备无故障 48h
BF 培菌第 1-4 周加菌剂; 投 NH₄Cl 至 TAN 2-4TAN 下降, NO₂ 出现
硝化成熟第 4-6 周继续培菌TAN<0.5 且 NO₂<0.3 持续 3 天
低负荷放鱼第 6-8 周放容量 20-30%TAN<1.0 稳定
加负荷第 8-12 周每周加 15-20%指标安全
满负荷第 12+ 周正常生产TAN<1, NO₂<0.3, DO>6

20.2 常见启动问题

问题原因措施
TAN 不降温度太低/pH 太低升温 25-30°C; 补碱度
NO₂ 持续高Nitrobacter 未建立正常! 耐心等待
放鱼 TAN 暴涨投鱼太多/投饵太猛减投饵 50%; 部分换水
pH 快速下降碱度耗尽加 NaHCO₃
鱼不吃食应激检查水质; 减投饵
第 四 部 分
经济评估

第 21 章运行成本

21.1 成本构成

成本项计算基础典型占比
饲料日投饵量 × 饲料单价50-65%
电费 (水泵)主泵+支路泵 × 24h8-15%
电费 (BF 风机)BF 粗孔曝气 × 24h3-5%
电费 (UV)UV 功率 × 24h1-3%
电费 (热泵)热泵电耗 × 电价0-15%
液氧o2Fish / 吸收率 × 单价5-12%
臭氧电费臭氧发生器 × 24h1-3%
CO₂ 风机电费脱气风机 × 24h2-5%
甲醇反硝化 × 单价1-3%
NaHCO₃碱度投加 / 0.595 × 单价1-3%
杂项电费照明/控制/备用1-2%
CAPEX 折旧 v1.1设备/8年 + 土建/20年 / 36515-30%

21.2 峰值 vs 日均双口径

峰值成本通常是日均的 1.5-2.5 倍。设备必须能承受峰值, 但日常运行接近日均。

📌 v1.1 OPEX 中的化学品/臭氧消耗在年均工况下按 feedRatio (= feedAvg/feedMax) 折算。

21.3 CAPEX 折旧 v1.1 新增

v1.1 引入双年限折旧:

年折旧 = 设备 / 设备年限 + 土建 / 土建年限 (元/年)
日折旧 = 年折旧 / 365 (自动加到 totalAvg/totalMax)

21.4 综合成本对比 v1.3 含折旧 + 工况

1000 t/yr 三文鱼基准, 4 工况下综合成本对比(v1.3):

工况气温/源水OPEX折旧综合 元/kg对照文献
青岛冬-7°C/4°C22.413.635.96
海南冬17°C/22°C24.113.637.70
海南夏32°C/27°C26.913.640.47
挪威 Bergen 冬0°C/4°C22.713.636.30Atlantic Sapphire 60-70(含人工/苗)

其他鱼种(v1.1 数据, 单一工况):

鱼种OPEX 元/kg折旧综合文献
大菱鲆 500t31.916.248.146.32 (海水鱼体系)
罗非鱼 1000t19.99.529.325-30
石斑鱼 300t27.412.740.160-100 (不含苗)
加州鲈 500t27.912.740.630-50
鳜鱼 300t35.916.051.980-100
鳗鲡 300t44.615.259.860-100 (含苗)
对虾 100t25.716.241.925-35

注: iRAS 不计算人工/苗种/管理费等"软成本", 综合成本通常比文献偏低 20-30%, 工程实际中另行加计。

第 22 章投资概算

22.1 设备直购成本

设备计价单位默认单价推荐范围
转鼓微滤机元/(m³/h)500300-800
MBBR 生物滤池元/m³4,5003,000-6,000
UV 消毒器元/kW10,0008,000-15,000
蛋白分离器元/m³25,00015,000-40,000
臭氧发生器元/(kg/d)100,00080,000-150,000
AOP 反应池元/m³35,00020,000-50,000
反硝化反应器元/m³5,5004,000-8,000
脱气塔元/m²20,00015,000-30,000
主循环泵+变频元/kW3,5002,500-5,000
支路泵元/kW3,0002,000-4,000
罗茨风机元/kW4,0003,000-6,000
液氧增氧塔元/(m³/h)350200-500
液氧储罐 (3t)元/套100,00080,000-120,000
热泵 (空气源)元/kW2,0001,500-3,000
板式换热器元/m²1,200800-2,000

22.2 工程总投资

工程总投资 = 设备直购 × 放大系数
放大系数 = 2.2 (默认, 范围 1.8-2.5)

放大系数涵盖: 管道阀件 (15-20%)、电气仪表 (10-15%)、自控 (5-10%)、安装调试 (10-15%)、土建厂房 (30-40%)、设计监理 (5-8%)。

概算精度约 ±40%, 实际造价以厂商报价为准。iRAS 投资概算用于项目可行性初判, 不能替代正式造价。

22.3 氧锥拓扑对 CAPEX 的影响 v1.7 新增

v1.7 引入氧锥拓扑 mainline / bypass 区别建模, 对 CAPEX 有显著影响 (1000 t 大西洋三文鱼实测):

设备项mainlinebypass差异说明
液氧增氧塔 A-601355 万 (Q=10155 m³/h)36 万 (Q×10% = 1015 m³/h)-319 万规模缩 10× (按小流量定价)
主泵 P-101213 万 (610 kW)171 万 (488 kW)-42 万主泵不过锥, 扬程降低
旁路泵 P-602028 万 (92 kW)+28 万v1.7 新增独立设备
液氧储罐 V-901≈ 相同≈ 相同~0站级 CAPEX 按日耗推算, 拓扑不影响日耗
设备 CAPEX 小计4,454 万3,842 万-612 万
工程总投资 (× Lang 2.2)9,800 万8,454 万-1,346 万 (-14%)含放大系数
📝 拓扑选择建议 1000 t 以上大型三文鱼 / 鳟鱼项目建议默认采用 bypass 拓扑, 节省 CAPEX 600+ 万元 + 电费 25+ 万元/年。中小型项目 (Q < 5000 m³/h) 用 mainline 更简单, 因为 旁路高压锥的单价规模效应需要项目规模达到一定阈值才显现。
第 五 部 分
设计案例

第 23 章案例: 三文鱼 1000 t/年 v1.7 数据 (mainline / bypass 对比)

本案例使用 iRAS v1.7 默认的大西洋三文鱼参数, 演示完整计算并对比两种氧锥拓扑。气候选用挪威 Bergen 冬季工况 (T_air=0°C, T_source=4°C, T_target=15°C), 这是工程实际三文鱼项目最具代表性的气候。

23.1 阶段参数

参数苗种幼鱼中鱼成鱼
体重 (g)0.1→5050→500500→20002000→5000
月数6864
FCR1.01.11.21.20
蛋白 (%)50454240
feedRate (%)3.01.51.00.95
peakFactor1.51.51.51.5
水温 (°C)12141515
密度 (kg/m³)30506070
盐度 (‰)001230
损耗 (%)15843
turnover4.03.02.53.0
o2FishFactor0.200.210.230.23
recDepth_m v1.31.53.05.07.0
subSystemCount v1.31234

23.2 关键计算结果 v1.7 实测

项目苗种幼鱼中鱼成鱼全场
存塘 (t)341132237413
鱼池水体 V_tank (m³)1118362,1933,3856,525
工艺水体 V_process (m³)~89~469~1,015~1,798~3,371
系统总水体 V_total (m³)~200~1,305~3,208~5,183~9,896
循环流量 Q (m³/h)4242,4655,48110,155
日峰值饲料 feedMax (kg/d)~36~321~1,073~2,251~3,681
液氧 (kg/d)1188298~7721,169
o2DailyTotal (kg/d, 全场)1,113
o2PeakHourly (kg/h, 全场)69.56
A-601 峰值供氧 (kg/h)41.20
CO₂ 风机 K-201 (kW)74191169309
模块数1234
每套池数1222
反推池径 D (m)109.410.58.8
📝 注 上表"全场"列是 4 阶段累加值, 但成鱼期 (idx=3) 通常是设计基准 (最大流量、最大氧需), 单阶段数据中成鱼期最有工程意义。"o2DailyTotal 1113 kg/d / A-601 41.20 kg/h / peakHourly 69.56 kg/h" 是本案例三个关键指标, 与 SUMMARY 参考 100% 吻合

23.3 综合成本 (Bergen 冬季 + bypass 拓扑, v1.7 默认)

下表是 v1.7 工具实际输出值 (jsdom 跑测,默认气候即挪威 Bergen 冬季, 默认价格电 0.7 元/kWh / 液氧 1.2 元/kg / 饲料按各阶段单价):

项目mainline (主流串联)bypass (旁路高压)差异
CAPEX 设备直购4,455 万元3,843 万元-612 万元
CAPEX 工程总投资 (× Lang 2.2)9,800 万元8,454 万元-1,346 万元
年折旧 (设备 8 年 / 土建 20 年)824 万元/年711 万元/年-113 万元/年
年 OPEX (不含折旧)2,846 万元/年2,830 万元/年-16 万元/年
年综合成本 (含折旧)3,670 万元/年3,541 万元/年-129 万元/年
OPEX (元/kg)28.46 元/kg28.30 元/kg-0.16
折旧 (元/kg)8.24 元/kg7.11 元/kg-1.13
综合成本 (元/kg 鱼, 含折旧)36.70 元/kg35.41 元/kg-1.29
对照文献Atlantic Sapphire (Florida) 60-70 元/kg (含人工/苗/管理)
Nordic Aqua 宁波 (中国境内) Q2 2024 实测生产成本 41 元/kg ≈ EUR 5.23/kg
偏差说明iRAS 不含人工 (~6-8 元/kg) / 苗种 (~5 元/kg) / 管理与营销 (~3-5 元/kg)。
iRAS 36.70 + 软成本 ~16 元/kg ≈ 53 元/kg,与挪威 Salmon Evolution 目标"与海水网箱打平"(NOK 50-60/kg ≈ 35-42 元/kg) 量级一致。
📝 重要更正 (vs v1.7 早期手册) 早期版本误把 OPEX (28 元/kg) 写成综合成本 (22 元/kg)。 真实情况: v1.7 跑 1000 t 三文鱼 mainline 默认气候: OPEX = 28.46 元/kg, 折旧 = 8.24 元/kg, 综合成本 = 36.70 元/kg。 bypass 拓扑可省 1.29 元/kg, 主要来自 CAPEX 节省 → 折旧降低。这个综合成本与 v1.3 历史四工况基线 (Bergen 36.30 元/kg) 完美吻合 ✓

23.4 工况敏感性分析 (沿用 v1.3 4 工况基线)

工况OPEX 元/kg综合 元/kg差异原因
青岛冬22.435.96基准 (v1.3 数据)
海南冬24.137.70仍需制冷 (气温 > 水温)
海南夏26.940.4732 → 15°C 制冷温差大
挪威 Bergen 冬22.736.30气温温和 + 源水冷
📝 注 上表是 v1.3 数据 (mainline 拓扑)。挪威 Bergen 冬 36.30 元/kg 与 v1.7 jsdom 实测 36.70 元/kg 完美吻合 (偏差 +1%, 在迭代版本计算精度内)。bypass 拓扑可再降约 1.29 元/kg, 见 23.3。SUMMARY 路线图已规划 v1.8 用 bypass + v1.7 字段重新生成 4 工况基线 JSON。

关键观察: 三文鱼养殖目标水温 12-15°C, 海南夏季气温 32°C, 制冷温差 17-20°C, 制冷电费暴增。不建议在海南养三文鱼。挪威成本接近青岛, 长期运行更优。

23.5 mainline vs bypass 拓扑详细对比 v1.7 新增

1000 t 大西洋三文鱼成鱼期 (Q=10155 m³/h, V_total ≈ 5183 m³, 4 模块) 在 Bergen 冬季工况下的两种拓扑全工况对比:

对比项mainline (主流串联)bypass (旁路高压)差异
动力消耗
主泵 P-101 功率610.4 kW488.3 kW-122.1 kW
旁路泵 P-602 功率0 (无)92.2 kW+92.2 kW
综合泵功率610.4 kW580.5 kW-29.8 kW (节能)
年电费节省 (8000 h × 0.7 元/kWh)~17 万元/年
设备规格
液氧增氧塔 A-601 流量10155 m³/h (全流量)1015 m³/h (10% 旁路)缩 10×
主泵扬程 H_main12 m + 3 m (锥子) = 15 m12 m (主泵不过锥)-3 m
旁路泵扬程 H_bypass20 m+20 m (高压锥)
CAPEX (设备直购)
主泵 P-101213 万元171 万元-42 万
液氧增氧塔 A-601355 万元36 万元-319 万
旁路泵 P-602028 万元+28 万
液氧储罐 V-901 (站级)~ 相同~ 相同~0
设备 CAPEX 小计4,454 万元3,842 万元-612 万
工程总投资 (× Lang 2.2)9,800 万元8,454 万元-1,346 万 (-13.7%)
财务指标 (售价 60 元/kg, 30%自有/70%贷款 @ 5%/8年, 折现率 8%, 项目 15 年)
全投资 IRR14.90%17.31%+2.41%
自有资金 IRR~22%25.23%+3%
全投资 NPV (万元)5,7797,031+1,252
静态回收期 (年)6.195.57-0.62
动态回收期 (年)8.637.42-1.20
BEP (% 满产)62.3%57.9%-4.3%
📌 v1.7 工程结论
1000 t 大西洋三文鱼项目采用 bypass 拓扑相比 mainline: 综合成本省 1.29 元/kg + 工程 CAPEX 节省 1346 万 + 全投资 IRR 提升 2.4 个百分点 + 静态回收期缩短 0.62 年。 这与国际主流大型三文鱼 RAS 多采用 side stream 拓扑的工程实际一致 (Linde SOLVOX cone / PR Aqua PPC / Global Seafood Advocate 2019 "more often plumbed in side stream configuration")。

第 24 章案例: 大菱鲆 500 t/年

📌 v1.1 修正 大菱鲆密度 v1.1 修正 (15/30/40 → 25/50/60 kg/m³), 与莱州东方海洋等循环水项目实测一致。

24.1 阶段参数

参数苗种中鱼成鱼
体重 (g)3→5050→250250→600
月数486
FCR0.91.11.3
蛋白 (%)555048
feedRate (%)3.01.50.9
水温 (°C)161616
密度 (v1.1)255060
盐度 (‰)282828
损耗 (%)1563
turnover3.52.52.0
o2FishFactor0.210.210.21
recDepth_m v1.30.81.01.0

大菱鲆是底栖鱼, 水深保持在 1m 以内, 主要靠底面积容纳生物量。

24.2 关键结果

项目数值对照
全场存塘809 t
鱼池水体 V_tank5,121 m³v1.0 配置 8,036 m³ (-36%)
CAPEX 工程总投资9,626 万元v1.0 14,192 万元 (-32%)
年折旧809 万元
年综合成本2,403 万元
综合成本48.1 元/kg文献 46.32 元/kg (+4%)

结论: v1.1 大菱鲆综合成本与国家海水鱼产业体系 2019 调查数据偏差仅 +4%, 建模精度达到工程概算可用水平。

第 25 章案例: 南美白对虾 100 t/年

25.1 清水 RAS vs Biofloc 对比

对比维度清水 RASBiofloc (BFT)
水处理核心微滤机 + 生物滤池靠异养菌絮团原位处理 TAN
有无生物滤池有 (MBBR)无 (最大区别)
水色清澈褐色 (絮团 300-800 mg/L)
碳氮比无要求C:N ≥ 15:1 关键
DO 需求按 RAS 标准极高 (异养菌+硝化菌+虾)
密度2-15 kg/m³3-8 kg/m³
投资
运行难度中等高 (絮团管理, 翻塘风险)
iRAS 支持

25.2 阶段参数 v1.1

参数苗种养成前期养成后期
体重 (g)0.01→11→1010→30
月数11.51.5
FCR1.01.31.6
蛋白 (%)454238
feedRate (%)1053
水温 (°C)303030
密度 (kg/m³)2815
盐度 (‰)151515
损耗 (%)20108
o2BaseFactor + Q100.31 + 1.3同左同左
o2FishFactor (30°C)0.400.400.40
recDepth_m v1.31.52.02.5
📌 v1.1 对虾 Q10 设为 1.3 (而非鱼类 2.0), 反映高温(30°C)代谢饱和。o2FishFactor 0.40 与 Boyd 2018 实测一致。

25.3 关键结果

项目数值
全场存塘6.25 t (稳态)
鱼池水体 V_tank~558 m³
CAPEX 工程总投资1,922 万元
年折旧162 万元
年综合成本419 万元
综合成本41.9 元/kg
对照中国南方 RAS 实测 25-35 元/kg, 偏高 +20% (规模效应)

对虾 100t 项目偏高的原因是规模效应—— iRAS 不考虑 CAPEX 0.7 次方缩放。规模 500 t/yr 后综合成本会降到 30 元/kg 左右。

第 六 部 分
工程化与平台功能 (v1.3 新增)

第 26 章多模块并联与 N+1 备用v1.3 新章

26.1 为什么大型项目要做模块化

对于年产量 > 1000 t/yr 的大型 RAS 项目, 工程上不会建一个超大的循环系统, 而是拆成多个模块并联。原因:

典型工业项目模块数(参考 Atlantic Sapphire / Stolt Sea Farm / 国内大型项目):

规模典型模块数单模块产能
< 500 t/yr1 套全场 1 模块
500-1500 t/yr2-3 套200-500 t/yr/模块
1500-5000 t/yr3-6 套500-1000 t/yr/模块
> 5000 t/yr5+ 套1000+ t/yr/模块

26.2 iRAS 模块化逻辑

v1.3.x 中, 每个阶段的模块数完全由用户决定。在阶段卡片头部"🔧 套数" subSystemCount输入(默认 1, 范围 1-10):

numModules = stage.subSystemCount
单模块流量 = 阶段总流量 / numModules
单模块设备容量 = 阶段总设备容量 / numModules

不同阶段可以有不同模块数。比如:

这种"阶段递增"的配置反映了大型 RAS 的实际工程: 苗种期对单一系统依赖性高, 成鱼期产量大需要分散风险。

26.3 模块内 N+1 备用配置

每个模块内部按以下规则配置备用:

设备单机容量上限配置规则备用
主循环泵 P-10175 kWceil(模块流量功率 / 75)+1 台备用
BF 风机 K-30230 kWceil(BF 风量功率 / 30)+1 台备用
RDF 微滤 F-201600 m³/hceil(模块流量 / 600)+1 台备用
UV D-4018 kW2 套并联各设计 75% 余量1 台坏 → 另 1 台加大流量
热泵 H-801200 kWceil(热泵装机 / 200)+1 台备用
支路泵(UV/蛋分/反硝化/AOP)30 kWceil(支路功率 / 30)+1 台备用
加药计量泵1 kWceil(加药功率 / 1)+1 台备用
O3 发生器 G-90150 kWceil(O3 功率 / 50)+1 台备用

26.4 全场共用站级设备

不按模块拆分, 全场共用 1 套(站级设备):

设备配置说明
液氧储罐 V-9012 罐(主+备)双罐设计, 切换不停产
应急柴油机1 台 100% 全场负荷30 秒启动 + 自动切换
变压器2 套(高/低压双套分列)故障切换
UPS1 套(PLC + 传感器 + 关键阀门)10-30 分钟续航

26.5 模块化对成本的影响

模块化的设备成本变化: 单台机组容量变小 → 单价(元/kW)略升, 但数量增加 → 总成本基本不变(±5%)。 真正的成本影响来自备用机组数: 4 个模块各配 N+1 备用, 比 1 个大模块配 N+1 备用多备 3 台。 典型规模(1000 t/yr, 4 模块)备用成本约占 CAPEX 5-10%。

📌 工程权衡
模块数越多, 故障隔离越好, 但备用成本越高。一般原则: 500 t/yr 以下不拆; 500-1500 拆 2-3 个; > 2000 拆 3-5 个。具体数字看可靠性要求(三文鱼必须高可靠, 罗非可降低)和投资预算。

第 27 章鱼池配置与池深统一v1.3 新章

27.1 鱼池形状与尺寸的工程逻辑

RAS 鱼池设计常见形状: 圆形(Cornell 双流式)/ 八角 / 方形带圆角。圆形最常见, 因为:

圆形池的关键尺寸:

27.2 反推算法

给定阶段总水体 V_tank、每套池数 poolsPerSet、套数 numModules, 反推单池直径:

总池数 = poolsPerSet × numModules
单池容积 V = V_tank / 总池数
反推直径 D = √(4V / πH)
高径比 = H / D

27.3 工程合理性自动告警

iRAS 在反推后自动检查工程合理性, 给出建议:

条件告警内容
D < 4 m反推单池直径偏小, 建议减少每套池数或减少并联套数
D > 30 m反推单池直径偏大, 建议增加每套池数或增加并联套数
D > 20 m 且 H < 3 m大池配浅池(H/D < 0.15), 流场易死角, 建议加深至鱼种推荐值
poolsPerSet = 1每套仅 1 池, 无主备/检修轮替, 工程惯例 ≥ 2

27.4 鱼种推荐池深 (recDepth_m)

v1.3 在 SPECIES_DB 各鱼种各阶段加入推荐池深字段。这反映不同鱼种的生物学特性和工程惯例:

鱼种苗种期幼鱼期中鱼期成鱼期生物学/工程理由
三文鱼1.53.05.07.0立体水柱泳鱼, 鲑科大池深, 商业 7-12m
大菱鲆0.81.01.0底栖鱼, 浅池单层底面积, 不超过 1.2m
石斑鱼1.22.53.0半底栖, 中等深度
加州鲈/鳜1.52.53.0淡水温水肉食鱼, 立体水柱
罗非鱼1.02.03.0耐高密度
鳗鲡0.81.52.0底栖, 中等浅
对虾1.52.02.5底栖+水柱中等

用户可在阶段卡片头部"🐟 每套 X 池 × Y m 深"调整 Y 值, 默认填入鱼种推荐值。

27.5 池深统一数据源 (v1.3.x)

v1.3 早期版本中, 池深有两个独立输入:

  1. 阶段卡片头的 stage.poolDepth_m: 用于池径反推
  2. 热平衡区的 stage.proc.poolDepth: 用于估算池表面积 → 车间面积 → 散热

这种双输入容易让用户混淆: 调了一处, 另一处没调, 结果池形和散热模型用了两个不同的池深。 v1.3.x 统一为一个数据源:

统一池深取值优先级 (calcThermalLoad 和 applyModularization 共用):
1. stage.poolDepth_m (用户在阶段卡片头部输入)
2. stage.recDepth_m (鱼种推荐, 来自 SPECIES_DB)
3. proc.poolDepth (老字段兼容, 老方案 JSON 加载用)
4. 1.5 m (最终默认)

同时删除了热平衡区的"池水深"输入, UI 上只在阶段卡片头部出现一次。同步更新加载老方案 JSON 的逻辑: 若老方案缺 recDepth_m, 加载时从 SPECIES_DB 按阶段名/index 自动补充, 保证池径反推和热平衡用同一个推荐值。

27.6 实例: 三文鱼成鱼 1000 t/yr

给定参数:

反推:

结果: 8 个 Φ8.8m × 7m 圆池, 每池 423 m³。无告警。

第 28 章全局气候与围护建模v1.2 新章

28.1 为什么要做项目级常量管理

v1.1 把热平衡参数放在每个阶段的 proc 字段里。同一个项目有 4 个阶段, 用户必须在 4 个阶段重复设置:

这些是项目级常量(同一个厂房无论养什么阶段都是同一份)。 v1.2 把它们挪到顶部全局配置区 globalThermalConfig, 一处设置, 所有阶段使用。

28.2 地区预设

v1.2 提供常见地区的气候预设, 选择后自动填入气温/源水温/RH 设计值:

地区冬季 T_air冬季源水夏季 T_air夏季源水年均RH 冬/夏
青岛(沿海)-7°C4°C32°C24°C13°C/14°C60%/78%
海南17°C22°C33°C30°C25°C/26°C75%/80%
挪威 Bergen0°C4°C17°C13°C8°C/8°C80%/75%
法罗群岛2°C6°C13°C11°C7°C/9°C85%/82%
北方内陆-15°C2°C30°C20°C8°C/12°C50%/65%
华南10°C15°C33°C26°C22°C/23°C72%/82%

用户选地区后仍可微调任何参数。

28.3 围护预设

预设U_wallU_roofU_floorU_window说明
工业标准0.80.60.82.550-80mm 聚氨酯保温, 双层玻璃
北欧高保温0.250.200.301.2200mm 岩棉, 三层 Low-E
简易厂房2.01.51.55.0钢板单层, 单层玻璃

28.4 4 工况对比设计

v1.2 引入3 工况对比(可扩展为更多), 每次计算同时跑:

工况用途
设计冬季 (design_winter)冬季最冷月用电峰值, 决定热泵装机
设计夏季 (design_summer)夏季最热月用电峰值, 决定制冷机装机
全年均值 (annual)年化电费, 用于成本核算和投资决策

当前选中的工况(thermalScenario)用于主显示, 其他工况结果在阶段卡片底部以对比表呈现。 热泵装机容量取冬/夏极值的较大者(装机定容)。

28.5 池盖类型(蒸发率系数)

池盖直接影响蒸发率, 是热平衡的关键变量:

池盖类型系数典型蒸发率适用场景
无盖 (none)1.0开放水面全蒸发室外或简易车间
部分覆盖 (partial)0.40约 40% 蒸发常规工业 RAS, 部分孔/进料口开
密闭 (full)0.10约 10% 蒸发北欧/严寒地区, 节能高保温

蒸发是潜热(2257 kJ/kg @100°C, 2502 kJ/kg @0°C), 是冬季最大的散热项之一。密闭池盖能省 20-40% 的暖气电费。

28.6 通风换气与热回收

通风换气次数 ACH(1/h)是车间空气质量的关键参数:

通风换气产生显热 + 潜热散热(冬季 T_air 比 T_room 低很多)。引入热回收可以回收 70-85% 的显热(板式或转轮换热器):

通风显热散热 = 0.34 × V_room × ACH × (1 − ventHeatRecovery) × (T_room − T_air)
单位: W (kW = 上式 / 1000)
0.34 = ρ_air × c_p,air / 3600 = 1.2 × 1005 / 3600 (W/(m³·K·ACH))

北欧项目通常 ventHeatRecovery = 0.7-0.85, 国内项目多为 0(无热回收)。

第 29 章双节点稳态热平衡求解v1.2 新章

29.1 v1.1 的简化模型与局限

v1.1 的热平衡是单点能量平衡: 给定 T_w(目标水温)和 T_air(室外气温), 直接算补水加热量 + 池散热量 + 鱼代谢热 + 设备产热, 平衡后给热泵装机。 这个模型忽略了车间空气节点, 默认池的热散到一个无限大的"室外"。问题在于:

29.2 双节点稳态方程 v1.6 公式订正

v1.2 引入双节点稳态求解器, 联立水节点能量、空气节点显热、空气节点湿度三个平衡:

节点 1: 水体热平衡 (kWh/d)
Q_metabolism + Q_equipment + Q_makeup − Q_pool_loss − Q_evap × λ = Q_HP

其中: Q_pool_loss = UA_pool × (T_w − T_room) × 24/1000
Q_evap × λ = 蒸发潜热散热, 计入水体(水分子带走显热)
Q_HP = 热泵供热(负)或制冷(正)负荷
节点 2: 空气显热平衡 (W) — 仅含显热
UA_pool × (T_w − T_room) + Q_air_gain
= UA_other × (T_room − T_air)
+ UA_floor × (T_room − T_ground)
+ UA_vent × (T_room − T_air)

其中 UA_other = U_wall × A_wall + U_roof × A_building + U_window × A_window
UA_pool = poolU × poolArea (poolU 室内 5, 室外 18 W/m²·K)
UA_floor = U_floor × A_building (地板对的是 T_ground, 不是 T_air)
UA_vent = 0.34 × V_room × ACH × (1 − ventHR)
Q_air_gain = 水泵/UV 进空气部分的产热 [W]
节点 3: 空气湿平衡 (kg/s)
m_evap × cover × A_water = m_air × ACH × (W_room − W_air_outdoor)

其中 m_evap = β × (W_sat(T_w) − W_room) (kg/(m²·s), 蒸发可能为负=凝结)
m_air = ρ_air × V_room (kg)
W_sat(T) = 用 Antoine 方程或 Magnus 公式估饱和含湿量
📌 v1.6 公式订正(29 章)
旧手册把 + Q_evap × λ 写在空气节点(节点 1)能量平衡里 — 这是 v1.2 推导笔误。 物理: 蒸发夺走的是水体显热(水分子蒸发吸热), 空气节点只接收水蒸气的湿度增量(W_room↑), 不直接接收 λ 项。代码 (calcThermalLoad) 一直把 λ 项放在水节点的能量平衡里, 物理正确。 v1.6 把手册 16.5 节和 29.2 节统一改为与代码一致的写法 — 水节点含 Q_evap × λ, 空气节点只有显热

耦合: T_room 影响 W_sat, W_room 影响 m_evap。Picard 迭代求解(实际代码用闭式解):

  1. 初值 T_room = T_air + 5, W_room = W_air_outdoor
  2. 用当前 T_room 算 W_sat, 用当前 W_room 算 m_evap
  3. 从节点 3(湿)解新 W_room(线性, 直接求)
  4. 从节点 2(空气显热)解新 T_room(线性, 闭式; 不含 m_evap×λ 项)
  5. 检查 |T_room_new − T_room_old| < 0.01°C → 收敛, 否则回 2
  6. 典型 1 次闭式解即可收敛(v1.2 闭式解已替代迭代,强通风/高湿度场景仍可用 Picard)

29.3 蒸发凝结模型

蒸发率 m_evap 用Lewis 关系简化:

m_evap [kg/(m²·s)] = β × (W_sat(T_w) − W_room)
β = h_c × Lewis_factor / c_p_air
β 典型值 ≈ 2-5 × 10^-4 kg/(m²·s) (室内静水, 弱风)

当 W_sat(T_w) > W_room: 蒸发(从池水夺热, 减少水温), m_evap > 0 当 W_sat(T_w) < W_room: 凝结(空气中水蒸气在冷池面凝结, 给池水放热), m_evap < 0

凝结情况典型: 海南夏季水温 12-15°C(冷水鱼)远低于室内空气 32°C/RH 80%, 大量凝结。 凝结水自动加到补水量减项, 反映节水效益:

V_condense_recovery = |m_evap| × A_water × 86400 / 1000 (m³/d)
V_makeup_actual = V_makeup − V_condense_recovery

29.4 求解输出

双节点求解器输出:

输出字段含义典型值范围
T_room车间空气温度T_air + 5 ~ +12°C
RH_room车间相对湿度55-80%
W_room车间含湿量0.008-0.020 kg/kg
m_evap_kg_h蒸发率(负值=凝结)±5 ~ ±50 kg/h
UA_envelope围护总 UA(无含地板)500-2000 W/K
Q_envelope围护总散热5-50 kW
Q_pool池散热5-30 kW
Q_vent通风散热2-15 kW
Q_makeup补水加热5-30 kW
Q_internal内部产热(鱼+设备)10-100 kW
heatPumpKWhPerDay热泵日电耗100-3000 kWh/d
heatPumpAvgKw热泵平均功率5-200 kW
cop当前工况 COP2.5-5.5

29.5 RH 警告与除湿建议

当车间 RH > 85% 时, 建议加除湿设备(否则结构件腐蚀、电气问题、霉菌)。 iRAS 在 RH 接近饱和时给出醒目提示, 建议:

29.6 热泵分类型

v1.2 拆出 4 种热泵, 各自的 COP 和源温模型:

类型COP_heatCOP_cool源温适用
空气源 air3.5-4.02.5-3.0T_air常规, 气温 > -10°C
水源 water4.5-5.53.5-4.5T_source有大量地下水/河水
地源 ground4.0-5.03.5-4.5T_ground (≈年均)有埋管空间
海水开式 seawater_open4.5-5.53.5-4.5T_seawater沿海/海上养殖, 钛板换热

不同热泵类型的源温和 COP 在 globalThermalConfig 中可独立设置。

第 30 章P&ID 工艺图自动生成v1.3.x 新章

30.1 P&ID 是什么

P&ID(Piping & Instrumentation Diagram, 工艺管道仪表流程图)是化工/环保工程的核心设计文档, 描述工艺单元、管道连接、流向、关键仪表和控制逻辑。一份完整的 P&ID 应包含:

30.2 iRAS 的 P&ID 生成方案

v1.3.x 提供独立 P&ID 查看器 pid.html, 主页面"📐 输出 P&ID"按钮一键打开。 数据流:

  1. 主页面计算完成后, "输出 P&ID"按钮把当前方案数据序列化为 JSON, 存到 localStorage
  2. 新窗口打开 pid.html, DOMContentLoaded 时自动读取 localStorage 数据并渲染
  3. 主页面再次点击时, 通过 postMessage 通知已打开的窗口刷新

30.3 工艺图布局

v1.3.x P&ID 采用两层布局:

主流层(上半部 y=560):

T-101 鱼池 → F-201 微滤 → B-301 MBBR → C-501 CO₂ 脱气 → S-101 集水池

P-101 主泵

A-601 增氧锥 → 回 T-101 (高架回水管)

旁路(从主流 T 接头取水):

公用工程层(下半部 y=900):

30.4 LOX 跟随阶段切换 (v1.3.x)

v1.3 早期版本中, LOX 储罐显示的氧耗是全场总和 (summary.totalLOXKgD), 不论用户切换到哪个阶段都不变。 v1.3.x 修复为当前阶段氧耗 (stage.load.o2DailyTotal), 切换阶段时同步更新。

右上角元数据"全场氧耗"保持不变(那本来就是全场总和的语义, 用于储罐选型)。

30.5 设备表自动生成

P&ID 底部自动生成设备表, 含 tag/中文名/规格/数量列表:

tag名称规格示例(成鱼期)
T-101鱼池8 池 × Φ8.8 × 7m, 423 m³/池
F-201微滤机Q 11170 m³/h, η_TSS 75%
B-301MBBR 生物滤池V 438 m³, η_TAN 70%, 填充率 50%
C-501CO₂ 脱气塔G:L 3, η 65%
S-101集水池200 m³
P-101主循环泵488 kW @ 12m H, 68% η, VFD
A-601增氧锥41.2 kg O₂/h, η 85%
S-701蛋白分离器(旁路)α 20%, V 74 m³
D-401UV 消毒(旁路)α 10%, 38 kW @ 40 mJ/cm²
R-801反硝化反应器(旁路)α 10%, V 129 m³
V-901LOX 储罐1113 kg/d (成鱼期), 3 t 储罐
X-901AOP 高级氧化(公用工程)α 5%, V 17 m³

30.6 主泵 VFD 联锁示例

主泵 P-101 由鱼池 DO-101 通过 PID 闭环控制 VFD 频率, 维持鱼池 DO ≥ 鱼种最低饱和度 × DO_sat。 联锁线在 P&ID 中表示为虚线箭头从 DO-101 指向 P-101_VFD

类似的联锁包括:

30.7 P&ID 导出

pid.html 支持SVG 直接打印(浏览器 Ctrl+P → 保存 PDF) 或SVG 文件下载(右键保存)。 SVG 矢量格式, 缩放无损, 可在 AutoCAD/Visio 中导入用作工程图基础。

第 31 章设备清单导出v1.3.x 新章

31.1 设备清单的工程用途

设备清单是工程设计的核心交付物之一, 用于:

31.2 iRAS 设备清单结构

v1.3.x 提供独立设备清单页面 equipment.html, 主页面"📋 导出设备清单"按钮一键打开。结构:

  1. 项目元信息: 鱼种 / 年产量 / 保存时间
  2. 各阶段独立清单(苗种期 / 幼鱼期 / 中鱼期 / 成鱼期)
  3. 每个阶段下:
    • 多模块并联说明(N 个模块)
    • 鱼池配置(数量 × 直径 × 深度 × 单池容积)
    • 主流设备(T-101 / F-201 / B-301 / C-501 / S-101 / P-101 / A-601)
    • 旁路设备(S-701 / D-401 / R-801)
    • 公用工程(V-901 / G-901 / X-901 / K-302 / H-801 / P-202)
    • 分阶段小计金额
  4. 全场合计(所有阶段加和)

31.3 N+1 备用展开

清单中每个设备按"主用 + 备用"展开。例如成鱼期 4 模块, 每模块主泵需要 P-101 75 kW × 2 台, 加 1 台备用:

模块号设备规格数量
模块 1P-101 主泵(主)75 kW × 488/4=122 / 75 ≈ 2 台2
模块 1P-101 主泵(备用)同上1
模块 2P-101 主泵(主)同上2
模块 2P-101 主泵(备用)同上1
…(其他模块)…
全场P-101 主泵 合计75 kW12 台(8 主 + 4 备)

31.4 单价与金额

设备单价在主页面"单价配置"区设置, equipment.html 自动读取并计算:

设备小计 = 数量 × 单台容量 × 单价(元/容量)
全场设备直购 = 各设备小计加和
工程总投资 = 设备直购 × 放大系数(默认 2.2)

放大系数 2.2 涵盖管道阀件 / 电气仪表 / 自控 / 安装调试 / 土建厂房 / 设计监理(详见第 22 章)。

31.5 单列静设备说明

某些设备虽然容积大、价值高, 但不需要 N+1 备用, 而是单列 + 旁路阀设计。例如:

设备不需备用的原因故障应对
MBBR 生物滤池 B-301填料缓冲数小时, 故障时减投饵旁路阀让水绕过, 投菌剂时再开
反硝化反应器 R-801NO₃ 累积慢, 可承受数日停机旁路阀让水绕过
蛋白分离器 S-701DOM 积累慢旁路阀
CO₂ 脱气塔 C-501风机有 N+1, 塔体本身没动风机切换
液氧储罐 V-901主+备双罐设计切换

设备清单上这些标识为"1 套 (单列 + 旁路阀)", 不展开成多份。

31.6 数据传递机制 (v1.3.x)

主页面和 equipment.html 之间的数据通过 localStorage 共享(iRAS_equipmentExport key), 跨 tab 同步:

  1. 主页面点"导出设备清单" → buildPIDExport() 生成完整数据 → 写入 localStorage
  2. 新窗口 window.open('equipment.html') 打开, DOMContentLoaded 读取 localStorage
  3. 主页面再次点击 → 写入新数据, postMessage 通知已打开窗口刷新
  4. 读取后立即清除 key, 防止刷新页面看到陈旧数据

v1.3.x 把 P&ID 和设备清单的 key 拆成两份(iRAS_pidExportiRAS_equipmentExport), 让两个独立窗口可以同时开着, 各自显示最新数据互不干扰。

第 32 章4 工况成本对比基线v1.3 新章

32.1 为什么要看 4 个工况

RAS 项目选址决定了气候/水温, 不同地区的成本差异巨大。v1.3 提供 4 个典型工况的标准化对比, 帮助:

32.2 标准基线(1000 t/yr 三文鱼)

使用 v1.3 默认参数, 4 工况下综合成本(含 CAPEX 折旧, 不含人工/苗种):

工况气温 / 源水OPEX折旧综合 元/kg
青岛冬季 设计冬-7°C / 4°C22.413.635.96
海南冬季 设计冬17°C / 22°C24.113.637.70
海南夏季 设计夏32°C / 27°C26.913.640.47
挪威 Bergen 设计冬0°C / 4°C22.713.636.30

32.3 关键观察

  1. 挪威成本接近青岛: 挪威冬季气温比青岛温和(0°C vs -7°C), 但源水温度都是 4°C 左右, 水温维持成本相近。挪威夏天比青岛凉(避免制冷), 长期成本可能更优。
  2. 海南夏季最贵: 三文鱼养殖目标水温 12-15°C, 海南夏季气温 32°C, 制冷温差 17-20°C, 制冷电费暴增。 不建议在海南养三文鱼。
  3. 海南冬季并不便宜: 冬季气温 17°C 高于水温 12°C, 仍需要制冷(虽然量小)。比青岛冬季略贵约 5%。
  4. OPEX 差异主要在热泵: 4 工况 OPEX 差异 4.5 元/kg, 其中热泵成本占 3-3.5 元/kg。

32.4 与文献的对照

文献来源项目综合成本 元/kgiRAS 偏差
Atlantic Sapphire (2023)迈阿密 RAS 三文鱼60-70(含人工/苗)iRAS 不含软成本, 偏低 25-30
Stolt Sea Farm欧洲大菱鲆40-50iRAS 大菱鲆 48.1, 接近
国家海水鱼体系 2019大菱鲆 500 t46.32iRAS 偏 +4%
国内罗非 RAS规模 1000 t25-30iRAS 29.3, 接近
南方对虾 RAS100 t25-35iRAS 41.9, 偏高(规模效应未建模)

整体精度: iRAS 与文献偏差 ±10-20%, 对于工程概算用途(±40% 精度)足够。

第 33 章投资财务评价v1.4 新章

33.1 为什么需要财务评价

iRAS v1.3 之前的输出停在工程概算: CAPEX 万元 / 年 OPEX 万元 / 综合成本元/kg。 对于设计院做可研报告、业主做投资决策、银行做贷款评审来说, 这些数据只是起点—— 真正的关键问题是:

v1.3 之前用户必须自己导出数据到 Excel, 手算这些指标, 一份可研里财务评价部分通常占 30% 篇幅, 设计院做一份要花 1-2 周。v1.4 加入独立的投资财务评价页面 (finance.html), 把这部分自动化, 把 1-2 周压缩到几小时。

33.2 数据流与系统架构

finance.html 是个独立页面, 与主页面通过 localStorage 数据共享:

主页面 calculateAll → lastResults.summary 写入财务汇总

buildPIDExport → JSON (含 summary.finance + perStage[i].cost/capex)

localStorage['iRAS_financeExport']

finance.html 读取 → 用户调财务参数 → 输出 8 张表 + KPI + 图

导出 Excel (.xlsx, 7 工作表) 或打印 PDF

这个架构有两个优点:

  1. 解耦: 工程评估和财务评价是两个独立工作流, 工程师管主页, 财务管 finance.html, 同一份方案可以两人协作
  2. 不污染主页: 财务参数改了不会触发主页重算, 用户能在 finance.html 里快速做敏感性分析

33.3 用户输入参数 (17 个)

用户在 finance.html 顶部填写 17 个财务参数, 分 3 组:

参数默认值说明
📈 销售与产能
📈售价 (元/kg)按鱼种 (鲑 60 / 大菱鲆 80 / 罗非 22 / ...)出厂价
📈项目年限 (年)15评价期, 5-30 可调
📈第 1 年达产率 (%)30BF 培菌 + 鱼苗培育期
📈第 2 年达产率 (%)70逐步加负荷
📈第 3 年起达产率 (%)100满产
📈建设期 (年)1含调试期
💵 投资与筹资
💵流动资金 (% × 年 OPEX)15启动期周转金, 末年回收
💵自有资金比例 (%)30剩余 70% 银行贷款
💵贷款利率 (% 年)5国内项目贷款典型
💵还款方式等额本息或等额本金
💵还款年限 (年)8通常短于项目年限
📊 评价参数
📊折现率 (% 年)8用于 NPV, 国内项目典型 8-10%
📊所得税率 (%)25企业所得税
📊销售税金附加 (%)0农产品免征, 其他 6-13%
📊修理费率 (% × 设备)2工程惯例
📊管理费率 (% × 销售)3典型工业项目
📊人工成本 (万元/年)按规模估30-360 万元/年, 5 万元/人/年

33.4 现金流模型

财务评价的核心是逐年现金流。每年有 4 类现金流入流出:

每年净现金流分解 (t = 0..N):

1. 经营性现金流 (Operating CF)
= 销售收入 (revenue × ramp[t])
− 销售税金 (revenue × salesTax)
− 经营成本 (variable × ramp + fixed)
− 修理费 (capex_eq × repairRate)
− 管理费 (revenue × adminRate)
− 人工成本
− 所得税 (max(0, profit) × taxRate)

2. 投资性现金流 (Investing CF)
t=0: − fixedAsset (建设期一次性投入)
t=1: − workingCapital (流动资金)
t=N: + workingCapital + 残值 (末年回收 + 残值)
其中残值 = 设备投资 × 5% + 土建投资 × 20% (v1.5+ 修正, 拆分计算)

3. 融资性现金流 (Financing CF) — 仅"自有资金"视角才考虑
t=0: + loan (收到贷款)
t=1..loanYears: − interest − principal (还本付息)

全投资 CF = (1) + (2)
自有 CF = 全投资 CF + (3)
📌 v1.5+ 残值公式订正
旧版残值 = fixedAsset × 5% 是经验值, 但设备和土建残值率不同: 设备(8 年折旧周期, 通常 10-15 年报废) 残值 ≈ 5% (机电二手价值); 土建(20 年折旧周期, 通常 30-50 年使用) 残值 ≈ 20% (建筑物结构残值)。 v1.5 起拆分计算: 残值 = totalCapexEquip × 5% + totalCapexCivil × 20%。 对千吨级三文鱼项目, 旧公式低估末年残值约 8-12%, 即 NPV 偏低 0.3-0.5 个百分点(对回收期/IRR 影响小)。

33.5 IRR 求解算法

IRR (Internal Rate of Return) 是让 NPV = 0 的折现率:

NPV(IRR) = Σ_t CF[t] / (1 + IRR)^t = 0

这是非线性方程, iRAS 用 Newton-Raphson 迭代求解:

IRR_(n+1) = IRR_n − f(IRR_n) / f'(IRR_n)
其中 f(r) = NPV(r), f'(r) = − Σ_t [t × CF[t] / (1+r)^(t+1)]

起点 r=0.10, 收敛阈值 |Δr| < 1e-7, 最多 100 次迭代

当 Newton-Raphson 发散(典型场景: 多重 IRR 解或符号变化异常),自动切换到二分法兜底: 在 [-0.5, 5.0] 区间求 NPV=0 的根。两种算法配合, 能稳定处理常见的 RAS 项目现金流。

33.6 输出: 8 张标准财务表

遵循国家发改委《建设项目可行性研究报告编制规范》要求:

名称主要内容
1投资估算设备 + 土建/安装/电气/自控/设计 + 流动资金 = 项目总投资
2资金来源运用 + 还款计划自有/贷款比例 + 逐年利息/本金
3总成本费用估算 (满产年)经营成本 + 修理 + 管理 + 人工 + 折旧 + 利息
4销售收入和税金逐年: 达产率 → 产量 × 单价 → 销售收入 / 税金
5损益表逐年: 销售收入 − 各项成本 − 税 = 净利润
6现金流量表 (全投资)逐年净 CF + 累计 CF + 累计折现 NPV → IRR
7现金流量表 (自有资金)含贷款流入 + 还本付息, 算股东 IRR
8敏感性分析5 因子各 ±20% 对 IRR 影响, 按摆动排序

33.7 关键指标 (6 个 KPI)

指标定义判断标准
全投资 IRR (税后)项目所有现金流(不含贷款)的内部收益率> 折现率 (8%) → 财务可行
自有资金 IRR股东视角 IRR (项目 CF + 贷款流入 − 还本付息)有杠杆放大效应
NPV (净现值)按折现率折算到 t=0 的累计 CF 现值> 0 → 项目可行
静态投资回收期累计净 CF 首次 ≥ 0 的年份(线性插值)RAS 项目典型 5-8 年
动态投资回收期累计折现 CF 首次 ≥ 0 的年份静态值 + 折现影响
盈亏平衡点 (BEP)满产年净利润 = 0 时的产量比例< 75% 满产 → 安全裕度大

33.8 敏感性分析 (5 因子)

对项目 IRR 进行单因子敏感性分析, 5 个常用因子:

因子变化幅度典型 IRR 摆动 (1000 t 三文鱼)
售价±20%±5-7% (最敏感)
CAPEX±20%±3-4%
产量±10%±1-2%
饲料成本±20% × share_feed±1-1.5%
电费±20% × share_elec±0.4-0.6%
📌 v1.5+ 敏感性 share 订正
旧版"share_feed = 55%, share_elec = 20%"是硬编码经验值, 仅适用于三文鱼 1000 t 项目。 v1.5 改为从实际 OPEX 占比动态读取: share_feed = feedCost / annualOpex。 对虾项目 share_feed 可能仅 35%, 罗非鱼项目可能 65%, 不同鱼种敏感性表会自动调整。 注意 share 是说明 项目特定 的敏感性, 不是计算公式 (计算还是 OPEX × ±20%, 只是用 share 解释"为什么饲料 ±20% 只让 IRR 摆动 ±1.5%而不是 ±5%")。

: 折现率不影响 IRR (IRR 是让 NPV=0 的折现率, 与折现率参数无关), 故不在敏感性表中。 但折现率会影响 NPV——折现率越高, NPV 越低。

结果按 IRR 摆动从大到小排序, 自动绘制 Tornado 图 (横向条形图)。 设计院通常拿这张图给业主看"这个项目的最大风险因子是什么"——RAS 项目永远是售价

33.9 实例: 1000 t 三文鱼 (青岛冬季工况)

用 v1.3 默认参数计算的方案数据 (CAPEX 14681 万 / 年 OPEX 2236 万 / 年折旧 1360 万), 售价 60 元/kg, 30% 自有 + 70% 贷款 (5%, 8 年), 折现率 8%, 项目 15 年:

关键指标

指标数值解读
全投资 IRR14.0%> 折现率 8%, 财务可行 ✓
自有资金 IRR20.5%贷款杠杆放大 6.5%
NPV @ 8%6878 万元正值, 创造 6878 万元股东价值
静态投资回收期6.3 年RAS 项目典型水平
动态投资回收期9.0 年含折现影响
盈亏平衡产量63.4% 满产实际产量打 7 折也能保本, 安全裕度大

敏感性 (按 IRR 摆动排序)

因子低值 IRR基准 IRR高值 IRR摆动
售价 ±20%8.0%14.0%19.1%11.1%
CAPEX ±20%17.6%14.0%11.3%6.3%
产量 ±10%12.3%14.0%15.6%3.3%
饲料成本 ±20%15.2%14.0%12.7%2.5%
电费 ±20%14.4%14.0%13.5%0.9%
📌 工程解读
售价是 RAS 项目的第一风险因子。三文鱼售价从 60 跌到 48 元/kg, IRR 立刻降到 8% (折现率), 项目变得勉强可行。 所以 RAS 投资决策的关键是锁定下游销售渠道和价格——很多 RAS 项目失败不是工程出问题, 而是建成后市场价格波动让财务模型崩溃。

33.10 鱼种横向对比

用同样的 1000 t 规模, 不同鱼种和价格:

鱼种售价 元/kg全投资 IRR自有 IRRNPV 万元回收期结论
三文鱼6014.0%20.5%68786.3 年主流, 可行
三文鱼 (高端)8022.2%37.1%181504.5 年有溢价能力, 优秀
大菱鲆8038.0%73.2%288733.0 年毛利高 (国内底栖鱼供需)
罗非鱼22−2.5%−6.1%−6798不回本RAS 不适合罗非, 池塘已够便宜

这张表能直接告诉投资人: 不要用 RAS 养罗非鱼 — 不是技术问题, 是经济问题。 RAS 适合高单价 (> 40 元/kg) + 集约化要求高 (海水/北方/特殊品种)的鱼种。

33.11 Excel 导出

finance.html 用 SheetJS 库 (CDN 加载) 实现一键导出 Excel, 输出 7 个工作表:

  1. 1. 投资估算
  2. 2. 资金来源 + 还款计划
  3. 3. 损益表 (逐年)
  4. 6. 全投资现金流
  5. 7. 自有资金现金流
  6. 8. 敏感性分析
  7. 关键指标 (汇总)

设计院拿到 Excel 后可以:

📌 离线兜底
SheetJS 库通过 CDN 加载, 网络不通时 Excel 导出失败, finance.html 会弹窗提示用户改用打印 (Ctrl+P → 保存为 PDF)。

33.12 局限性与未来扩展

v1.4 财务评价的当前局限:

未来 v1.5+ 可能扩展的方向:

但所有扩展之前的前提是: v1.4 当前版本被实际用户用过, 反馈了真实问题。设计院、咨询公司、业主用一个月, 告诉我们哪里反人类、哪里数字对不上、哪里需要补字段——这才是优先级最高的下一步。

第 34 章可研报告生成器v1.5 新章

34.1 为什么需要可研报告生成器

v1.4 之前的 iRAS 输出层级是工程数据(P&ID + 设备清单 + 财务评价 Excel), 用户拿到这些数据后还需要手动整合到 Word 可研报告里。 这是 RAS 项目从"概念"到"立项"过程中的最大手工劳动:

v1.5 加入可研报告生成器 (report.html), 把 iRAS 的工程和财务数据自动套到 16 章扩展版可研模板中, 输出 .docx 文档。 设计院拿到后只需补充项目身份、场址自然条件、土建详细设计等"必须人工填写"的内容即可。 预期能把可研撰写时间从 1-2 个月压缩到 1-2 周

34.2 关键产品决策

v1.5 在设计阶段做了 10 个产品决策, 每个决策都直接影响实现方式:

  1. 范围: 完整版 16 章扩展版(适合大型项目+敏感地区)
  2. 模板版本: 政府立项发改委版(向下兼容银行/农业/私募版)
  3. 章节结构: 16 章扩展版(把环境影响和社会效益独立成章)
  4. 项目元数据: 标准 20 字段(平衡填表门槛和报告完整性)
  5. Word 模板获取: 不替设计院排版, 只输出结构化内容(每个院都有自己的模板, 替他们决定版式反而帮倒忙)
  6. 文字模板: 模板 + 条件分支(不用 LLM。可研是法律文件, LLM 幻觉风险与收益不成比例)
  7. 模板引擎: nunjucks(Jinja2 JS 移植, 语法标准, 后续可移植)
  8. docx 生成: docx.js (npm: docx)(浏览器侧懒加载)
  9. 模板存储: 内嵌 <script type="text/njk-template">(保持 iRAS 单文件部署传统)
  10. 元数据 UI: report.html 内弹出表单(主页面保持简洁, 只有真正生成报告的用户才填表)
📌 决策 5 — 不替设计院排版
这是 v1.5 最重要的产品决策。商业级可研生成器(如某些 SaaS)倾向于做"成品 Word", 但每个设计院都有自己的红头模板, 版式、字体、目录格式各不相同。iRAS 替他们决定版式反而帮倒忙——他们还是要把内容复制到自己的模板里。 v1.5 输出"结构化内容", 用户拿到后复制到自己院的模板里, 是更尊重设计院工作习惯的做法。
📌 决策 6 — 不用 LLM
大模型可以生成流畅的文字, 但 v1.5 拒绝了这条路线。理由:
1. 可研是法律文件, 一句错误的政策引用可能让项目立项失败
2. 同一项目两次生成结果不一样, 业主审核会困惑
3. 难以追溯每句话的来源
4. 依赖网络和 API 配额, 离线不能用
模板 + 条件分支虽然文字略呆板, 但确定性可追溯性是可研的核心要求。

34.3 16 章覆盖范围与自动化率

v1.5 覆盖发改委可研编制规范的扩展 16 章。各章自动化率(数据来自 iRAS 自动填充的比例):

名称iRAS 自动化率主要数据来源
1总论85%项目元数据(20 字段)+ 工程财务汇总
2项目背景与建设必要性50%鱼种条件分支 + 用户填地区描述
3市场分析与建设规模60%鱼种市场概述 + 规模配置
4建设场址与建设条件30%大量 placeholder, 用户填具体场址
5工艺技术方案95%iRAS 工艺数据完整覆盖
6主要设备方案100%iRAS 设备清单完整覆盖
7土建/总图/公用工程50%车间面积/电力负荷自动 + 详细土建用户填
8节能节水70%用水量/电耗自动 + 节能措施模板
9环境影响评价60%排放估算自动 + 现状用户填
10劳动安全卫生与消防40%危险源识别 + 应急预案模板, 不替代安评
11项目实施与组织机构65%进度甘特表 + 人员配置(规模分支)
12投资估算与资金筹措95%iRAS 投资数据完整覆盖
13财务评价100%iRAS finance.html 数据完整覆盖(8 张标准财务表 + NPV/IRR + 敏感性)
14国民经济 / 社会效益评价50%就业/税收/产业链估算(简化版社会效益)
15风险分析70%4 类风险识别 + 敏感性 + 风险等级条件分支
16结论与建议80%整体结论模板 + 关键指标条件分支

加权平均自动化率约 70%, 用户需要补充的主要是: 项目身份(20 字段)、场址自然条件、土建详细设计、消防/环评/安评等专项内容(后三者必须由有资质机构出具)。

34.4 数据流与系统架构

主页面 calculateAll → lastResults.summary 写入财务汇总

buildPIDExport → JSON (含 farm + workflow + capex + opex + finance)

localStorage['iRAS_reportExport']

report.html 读取 → 8 命名空间 context 拼装

20 字段元数据表单(用户填) → context.project

nunjucks.renderString(template, context) × 16 章

parseTextToDocxChildren(text, dx, images) → docx 对象树

SVG 图(工艺流程/Tornado/累计 CF)→ Canvas → PNG → ImageRun

docx.Packer.toBlob() → 用户下载 .docx

34.5 8 个数据命名空间

模板里所有变量按以下命名空间组织:

命名空间来源主要字段
{{ project }}用户在 report.html 表单填写项目名/业主/地址/编制单位/资质等 20 字段
{{ farm }}iRAS 主页面方案鱼种/规模/各阶段配置/池数/池径反推/总存塘等
{{ workflow }}iRAS perStage 工艺数据RDF/BF/UV/CO₂/增氧锥/蛋分/反硝化/AOP/集水池
{{ thermal }}iRAS 热平衡求解(预留)气候/围护/COP/工况(v1.5.x 完善)
{{ capex }}iRAS 投资估算设备直购/土建/流动资金/总投资
{{ opex }}iRAS 运行成本年 OPEX/综合成本/折旧/电耗/液氧
{{ finance }}finance.html 计算逻辑params(17 个用户参数)+ indicators(6 个 KPI)+ tables(8 张表)
{{ meta }}报告元信息生成时间/iRAS 版本/模板版本

34.6 嵌入 docx 的 3 张图

位置内容
图 5-1第 5 章 5.2 工艺流程6 个主流程方框 + 主循环 + 旁路单元(条件显示)+ 关键参数表
图 13-1第 13 章 13.6.1 全投资现金流实线全投资 + 虚线自有资金, 含 0 轴线、网格、节点标记
图 13-2第 13 章 13.11 敏感性分析5 因子 Tornado 图, 按 IRR 摆动从大到小排序

图表生成路径: 浏览器侧 SVG → Canvas → PNG → docx ImageRun。中文字体使用系统 SimHei/Microsoft YaHei, 高清 2× DPR 渲染。

34.7 输出文件示例

三文鱼 1000 t/年项目, 用 v1.5 生成的可研报告大致规格:

项目数值
文件大小约 200-400 KB(含 3 张 PNG 图)
章节数16 章 + 封面
段落数约 800 段
表格数约 50 个
列表项约 380 个
粗体强调约 600 处
页数(Word 默认排版)约 60-80 页

34.8 局限性与未来扩展

v1.5 当前局限:

未来 v1.5.x / v1.6+ 可能扩展:

但所有扩展前提是: v1.5 当前版本被实际用户使用过, 反馈了真实问题

附 录
参数表 · 命名规范 · 文献

附录 A: 完整符号表

本附录列出 iRAS 所有版本 (v1.0 → v1.7) 使用的全部符号。

符号含义单位默认值 / 范围
质量平衡 + 浓度
Δ每循环污染物增量mg/L
R总通过率
C_out鱼池出水浓度mg/L
C_in回水浓度mg/L
η_rdfRDF TSS 去除率%75
η_bioBF TAN 硝化率%70
η_bioTSSBF TSS 去除率%20
VTR₂₀BF @20°C 容积去除率g/(m³·d)500
VDR₂₀反硝化 @20°Cg/(m³·d)800
θ_nit / θ_den温度系数1.10 / 1.08
f_sal盐度修正max(1−0.01S, 0.3)
SOTE氧传质效率%粗孔 8% / 微孔 18%
tanCoefTAN 系数0.092 (淡) / 0.110 (海水肉食)
RQ呼吸商1.0
v1.1 物理框架
o2BaseFactor20°C 基准0.13-0.32 (按物种)
o2Q10温度敏感2.0 (鱼) / 1.3 (虾)
o2FishFactor生效值base × Q10^((T−20)/10)
o2DOMfactor异养菌耗氧0.10 (可调 0.05-0.20)
peakFactor日内峰均比1.5
safety*工程安全系数1.10-1.15
COP_heat / cool热泵4.0 / 3.0
η_HX板换回收率%75
co2BlowerPressureCO₂ 风压kPa2.5
co2BlowerEfficiencyCO₂ 风机效率%60
depYearsEquip设备折旧8
depYearsCivil土建折旧20
V_tank鱼池水体biomass / density
V_process工艺水体各设备水量合计
V_total系统总水体V_tank + V_process
v1.2 全局气候 + 围护
T_air_design_winter设计冬季气温°C-7 (青岛)
T_air_design_summer设计夏季气温°C32 (青岛)
T_air_annual年均气温°C13 (青岛)
T_source_design_winter设计冬季源水温°C4
T_source_design_summer设计夏季源水温°C24
T_source_annual年均源水温°C14
T_ground地温(年均近似)°C= T_air_annual
RH_winter / summer / annual3 工况相对湿度%60 / 78 / 70
U_wall墙保温系数W/(m²·K)0.8
U_roof屋顶保温系数W/(m²·K)0.6
U_floor地板保温系数W/(m²·K)0.8
U_window窗保温系数W/(m²·K)2.5
windowRatio窗墙比0.05
airChangeRate (ACH)通风换气次数1/h1.0
ventHeatRecovery排风热回收效率0 (国内) / 0.7-0.85 (北欧)
poolCover池盖类型none(1.0)/partial(0.4)/full(0.1)
buildingArea_user用户指定车间面积0(自动 = 1.6 × poolArea)
buildingHeight车间净高m6
thermalScenario当前工况design_winter / design_summer / annual
climateRegion地区预设china_qingdao / china_hainan / norway_bergen / …
envelopePreset围护预设industrial_standard / nordic_premium / cheap
hpType热泵类型air / water / ground / seawater_open
v1.2 双节点稳态求解
T_room车间空气温度(求解输出)°CT_air + 5~12
RH_room车间相对湿度(求解输出)%55-80
W_room车间空气含湿量(求解输出)kg/kg0.008-0.020
m_evap蒸发率(负值=凝结)kg/(m²·s)±2-50 ×10⁻⁵
UA_envelope围护总 UA(墙+顶+窗)W/K500-2000
UA_pool池散热 UAW/KpoolU × poolArea
poolU池水对空气换热系数W/(m²·K)5(室内)/ 18(室外)
v1.3 模块化 + 鱼池配置
subSystemCount阶段并联套数(用户输入)1-10, 默认 1
numModules实际模块数 = subSystemCount≥ 1
poolsPerSet每套池数(用户输入)1-10, 默认 1
poolDepth_m单池深度(用户输入, 统一数据源)m0.5-8, 默认 = recDepth_m
recDepth_m鱼种推荐池深(SPECIES_DB)m0.8-7.0(按鱼种和阶段)
v1.3.x 新增
thermalSharingFactor共用车间散热系数0.5-1.0, 默认 0.75
V_sump集水池容积clamp(Q/30, 20, 200)
v1.4 投资财务评价 (finance.html)
price (元/kg)售价元/kg按鱼种 (鲑 60, 大菱鲆 80, 罗非 22, 等)
years项目评价年限5-30, 默认 15
ramp1/2/3逐年达产率%30 / 70 / 100 (默认)
buildYears建设期0.5-3, 默认 1
workingCap流动资金 (% × 年 OPEX)%0-50, 默认 15
ownPct自有资金比例%0-100, 默认 30
loanRate贷款利率% 年0-15, 默认 5
loanMethod还款方式等额本息 / 等额本金
loanYears还款年限1-20, 默认 8
discount折现率 (用于 NPV)% 年0-20, 默认 8
taxRate所得税率%0-35, 默认 25
salesTax销售税金附加%0-20, 默认 0 (农产品免征)
repairRate修理费率% × 设备/年0-10, 默认 2
adminRate管理费率% × 销售收入0-20, 默认 3
laborCost人工成本元/年按规模估, 30-360 万元/年
NPV净现值= Σ CF[t]/(1+r)^t, r=discount
IRR内部收益率%令 NPV=0 的折现率, Newton-Raphson 求解
BEP盈亏平衡产量% 满产满产年净利润=0 时的产量比例
v1.7 氧锥拓扑 + 系统排水
o2ConeTopology氧锥工艺拓扑mainline / bypass, 默认 mainline
o2ConeHeadLoss_m主路时氧锥附加压头m3, 钳 [0, 15]
o2BypassRatio旁路占比 (bypass 模式)%10, 钳 [5, 30]
o2BypassPumpHead_m旁路泵 P-602 扬程m20, 钳 [5, 40]
o2BypassPumpEta旁路泵效率%60, 钳 [30, 85]
aopO3Dose_gm3AOP 臭氧投加剂量g O₃/m³0.5, 钳 [0.1, 3.0] (v1.7 暴露字段, UI v1.8)
sim.discharge.V_m3_d系统排水日流量m³/d= V_total × exchangeDaily (v1.7 重构)
sim.discharge.TSS_kg_d系统排水 TSS 日量kg/d= load.tssDaily (质量守恒, 含污泥)
sim.discharge.TAN_mg_L系统排水 TAN 浓度mg/L= C_tank (溶解态)
sim.discharge.NO3_mg_L系统排水 NO₃ 浓度mg/L= C_tank (溶解态)

附录 B: 全局热配置参数表(完整)

iRAS v1.3 全局热配置(globalThermalConfig)字段一览:

字段名类型默认值UI 位置说明
climateRegion下拉china_qingdao地区预设选择后自动填 T_air/T_source/RH
thermalScenario下拉design_winter当前工况3 个工况切换
envelopePreset下拉industrial_standard围护预设选择后自动填 U 值
hpType下拉air热泵类型4 种
T_air_design_winter数值-7气温区°C
T_air_design_summer数值32气温区°C
T_air_annual数值13气温区°C
T_source_design_winter数值4源水温区°C
T_source_design_summer数值24源水温区°C
T_source_annual数值14源水温区°C
RH_winter数值60湿度区%
RH_summer数值78湿度区%
RH_annual数值70湿度区%
U_wall数值0.8围护区W/(m²·K)
U_roof数值0.6围护区W/(m²·K)
U_floor数值0.8围护区W/(m²·K)
U_window数值2.5围护区W/(m²·K)
windowRatio数值0.05围护区窗墙比
airChangeRate数值1.0围护区ACH (1/h)
ventHeatRecovery数值0围护区0-0.85
poolCover下拉partial围护区none / partial / full
buildingArea_user数值0(自动)围护区0 = 自动估算
buildingHeight数值6围护区m
thermalSharingFactor数值0.75围护区0.5-1.0, 共用车间散热系数

附录 C: 设备 tag 命名总表

iRAS 完全采用 ISA-5.1 规范 tag 命名(字母+数字), 字母含义:

首字母含义示例
TTank 储罐 / 池T-101 鱼池
FFilter 过滤器F-201 微滤
BBiofilter / Bio-reactor 生物反应器B-301 MBBR
KKompressor 压缩机/风机K-101 鱼池增氧风机 v1.7 / K-201 CO₂ 风机 / K-302 BF 风机
DDisinfection 消毒D-401 UV
CColumn 塔/柱C-501 CO₂ 脱气塔
SSeparator / Sump 分离 / 集水S-101 集水池, S-701 蛋分
AAerator / Absorber 增氧 / 吸收A-601 增氧锥
RReactor 反应器R-801 反硝化
XeXtra / 高级处理X-901 AOP
GGenerator 发生器G-901 O₃
VVessel 大型容器 / 储罐V-901 液氧站级储罐
HHeat 热相关H-801 热泵
PPump 泵P-101 主循环泵 / P-602 氧锥旁路泵 v1.7
MMixer 混合点 / 管段配件 v1.7M-601 bypass 模式氧锥旁路混合点 (T 接头)

位号编号规则

仪表位号

tag 前缀类型示例
DO溶解氧DO-101 鱼池 DO 探头
pHpH 计pH-102 鱼池 pH
TT温度变送器TT-103 鱼池温度
LT液位变送器LT-104 鱼池液位
PDS压差开关PDS-201 RDF 压差
AT分析变送器AT-302 BF 出口 TAN
FT流量变送器FT-101 主管流量

附录 D: 鱼种参数表

iRAS 内置 9 个物种 (含 custom), 这里展示主要 4 种, 完整见软件。

品种阶段体重(g)FCR蛋白%feed%密度T°C盐‰recDepth_m
三文鱼苗种0.1→5061.0503.0301201.5
幼鱼50→50081.1451.5501403.0
中鱼500→200061.2421.06015125.0
成鱼2000→500041.20400.957015307.0
大菱鲆 v1.1苗种3→5040.9553.02516280.8
中鱼50→25081.1501.55016281.0
成鱼250→60061.3480.96016281.0
加州鲈苗种1→5020.9503.5252501.5
中鱼50→25041.1451.8402502.5
成鱼250→60041.3401.2502503.0
对虾苗种0.01→111.04510230151.5
养前1→101.51.3425830152.0
养后10→301.51.63831530152.5

完整 9 物种 (含罗非鱼/石斑鱼/鳗鲡/鳜鱼) 见 iRAS 软件内置数据库。recDepth_m 是 v1.3 新增的鱼种推荐池深字段。

附录 E: 参考文献

v1.0 引用

  1. Timmons, M.B., Ebeling, J.M. (2010). Recirculating Aquaculture, 2nd ed. Cayuga Aqua Ventures.
  2. Sharma, B., Ahlert, R.C. (1977). Nitrification and nitrogen removal. Water Research, 11(10).
  3. Benson, B.B., Krause, D. (1984). DO in freshwater/seawater. Limnol. Oceanogr., 29(3).
  4. USEPA (2006). UV Disinfection Guidance Manual (UVDGM). EPA 815-R-06-007.
  5. Summerfelt, S.T. et al. (2009). Ozonation and UV in RAS. Aquacultural Engineering, 40(1).
  6. Mateju, V. et al. (1992). Biological water denitrification. Enzyme and Microbial Tech., 14(3).
  7. Rusten, B. et al. (2006). Nitrification in MBBR. Aquacultural Engineering.
  8. Chen, S. et al. (2006). Nitrification kinetics in RAS. Aquacultural Engineering, 34(3).
  9. Nijhof, M., Bovendeur, J. (1990). Fixed film nitrification. Aquaculture, 87(2).
  10. Colt, J. (2006). Water quality in reuse systems. Aquacultural Engineering, 34(3).

v1.1 新增引用

  1. Boyd, C.E. (2018). Oxygen Demand of Aquaculture Feeds. Responsible Seafood Advocate.
  2. Davidson, J. et al. (2014). Atlantic salmon performance in RAS. USDA FWI.
  3. PR Aqua (2016). Low Head Oxygenator Design Standards.
  4. Vinci, B. et al. (2016). Saltwater RAS LHO. Aquacultural Engineering, 75:22.
  5. 国家海水鱼产业技术体系 (2019). 大菱鲆养殖产业调查数据. CARS-50.
  6. Stolt Sea Farm (2024). Annual Report - Land-based Turbot.
  7. Atlantic Sapphire (2023). Bluehouse RAS Technical Documentation.
  8. Delta Cooling Towers Inc. CO₂ Air Stripper Case Studies. www.deltacooling.com
  9. Suez Degremont. Water Handbook - CO₂ Degasification.
  10. 中国水产科学研究院黄海水产研究所. 大菱鲆工厂化循环水养殖技术研究. 2015.
  11. Hansa Engineering. CO₂ Stripping Column Internals. RVT Process Equipment.

v1.2 新增引用 (热平衡 + 围护)

  1. ASHRAE Handbook (2021). HVAC Applications, Chapter 22: Animal Environments.
  2. Lewis, W.K. (1922). The evaporation of a liquid into a gas. ASME, 44.
  3. Antoine, C. (1888). Vapeurs: nouvelle relation entre les tensions et les températures.
  4. Magnus, G. (1844). Versuche über die Spannkräfte des Wasserdampfes. Annalen der Physik.
  5. SINTEF Energy Research (2020). Heat Pumps for Aquaculture RAS Facilities.
  6. Mota, V.C. et al. (2019). The effects of carbon dioxide on growth performance, welfare, and health of Atlantic salmon post-smolt in RAS. Aquaculture, 498: 578-586.

v1.3 新增引用 (模块化 + 鱼池设计)

  1. Cornell University. Twin Drain Tank Design Guidelines (Cornell dual-drain RAS).
  2. Timmons, M.B. & Summerfelt, S.T. (2007). Pilot-scale evaluation of dual drain rectangular and circular tanks. Aquacultural Engineering.
  3. Atlantic Sapphire (2023). Annual Report - Modular RAS Architecture.
  4. AKVA Group (2022). Land-based Aquaculture: Standard Module Design.
  5. Billund Aquaculture. Modular RAS Design Standards (Denmark).
  6. Stolt Sea Farm (2024). Sustainability Report - Production Capacity Distribution.

v1.3.x 新增引用 (P&ID + 工程化)

  1. ISA-5.1 (2009). Instrumentation Symbols and Identification.
  2. ANSI/ISA-5.4. Instrument Loop Diagrams.
  3. ISA-S5.5. Graphic Symbols for Process Displays.
  4. API RP 14E. Design and Installation of Offshore Production Platform Piping Systems.

v1.4 新增引用 (投资财务评价)

  1. 国家发展和改革委员会、住房和城乡建设部 (2006). 建设项目经济评价方法与参数 (第三版). 中国计划出版社.
  2. 国家发改委 (2002). 投资项目可行性研究指南.
  3. 住建部 (2017). 中央政府投资工程建设项目可行性研究报告编制深度规定.
  4. 财政部 (2024). 企业会计准则第 4 号——固定资产折旧.
  5. Brealey, R.A., Myers, S.C., Allen, F. (2020). Principles of Corporate Finance, 13th ed. McGraw-Hill (NPV/IRR 理论).
  6. SheetJS Community Edition (2024). xlsx.full.min.js — Spreadsheet Data Toolkit. https://sheetjs.com

v1.5 新增引用 (可研报告生成器)

  1. 国家发改委 (2017). 政府投资项目可行性研究报告编制深度规定 (2017 年发改投资规〔2017〕1638 号).
  2. 建设部 (2007). 投资项目可行性研究指南 (试用版).
  3. nunjucks (Mozilla, 2014-2024). A rich powerful templating language for JavaScript. https://mozilla.github.io/nunjucks/
  4. docx.js (Dolan Miu, 2018-2024). Easily generate .docx files with JS/TS. https://docx.js.org/

v1.6 新增引用 (站级设备 + 物理校正)

  1. Mateju, V. et al. (1992). Biological water denitrification - A review. Enzyme and Microbial Technology, 14(3): 170-183. (反硝化碱度回收 2.5-3.0 g CaCO₃/g N 工程实测)
  2. Henze, M. et al. (1991). Activated Sludge Model No. 1. IAWPRC Scientific and Technical Report No. 1.
  3. Crittenden, J.C., Trussell, R.R., Hand, D.W. et al. (2012). MWH's Water Treatment: Principles and Design, 3rd ed. John Wiley & Sons. (反硝化化学计量)
  4. Mokhayeri, Y. et al. (2008). Kinetic characterization of methanol-utilizing denitrifying biomass. Water Environment Research, 80(5).
  5. Suhr, K.I., Pedersen, P.B. (2014). Methanol dosing in RAS denitrification. Aquacultural Engineering, 60: 47-54.
  6. Tsukuda, S. et al. (2015). Methanol-driven denitrification in marine RAS. Aquaculture Research, 46(11).
  7. 赵越, 刘鹰 (2018). 大西洋鲑 MBBR 设计参数研究. 大连海洋大学学报, 33(4): 543-549.
  8. Drønen, K. (2021). Atlantic salmon land-based RAS in Norway: V_process/V_total benchmark.

v1.7 新增引用 (氧锥拓扑 + 系统排水重构)

  1. Linde Gas (2020). SOLVOX cone 60 / 110 oxygenation system. Technical datasheet. (氧锥操作压力 3.0-3.8 bar, side stream 拓扑)
  2. Pentair AES (2021). Speece Cone Oxygenation: Engineering Design Manual. (玻璃钢锥操作压力 10-21 psi)
  3. PR Aqua (2019). PPC (Pressurized Production Column) High Pressure Oxygen Contactor Design Standards. (高压氧合柱 3-4 bar)
  4. Global Seafood Advocate (2019). Oxygenation in land-based RAS: cone vs LHO. "more often plumbed in side stream configuration" (大型 RAS 旁路拓扑工程实际)
  5. Sharrer, M.J. et al. (2010). Process and economic evaluation of microscreen drum filter backwash water recycling. Aquacultural Engineering, 43(1): 9-19. (RDF 反冲水 95-98% 清液回流)
  6. Wagner, M. et al. (2010). Standardized oxygen transfer testing for fine-bubble aeration in aquaculture. Aquacultural Engineering, 43(2): 33-41. (微孔盘曝气 SOTE 18-22%)
  7. Schroeder, J.P. et al. (2010). Comparison of ozone and ultraviolet light for the disinfection of recirculating water in marine RAS. Aquacultural Engineering, 42(1): 1-10. (AOP 臭氧投加剂量经验范围)
  8. Gonçalves, A.A., Gagnon, G.A. (2011). Ozone application in recirculating aquaculture system: an overview. Ozone: Science & Engineering, 33(5): 345-367.

附录 F: 免责声明

本工具(iRAS)和本手册的计算结果仅供工程概算参考。

— 本手册全文完 —

© 2026 海南登登科技咨询有限公司 · 孙程
iras.cn · github.com/chengsunmail/iras · 6881509@qq.com
本手册整合 iRAS v1.0 → v1.5 全部内容 · 完整版