工业化凡纳滨对虾(Litopenaeus vannamei)养殖中三氯异氰尿酸消毒的优化：一种系统动力学(System Dynamics)方法

《Aquaculture Reports》：Optimizing trichloroisocyanuric acid disinfection in industrialized Litopenaeus vannamei culture: A system dynamics approach

　　

摘要：三氯异氰尿酸(Trichloroisocyanuric acid, TCCA)作为含氯消毒剂广泛用于对虾养殖业，但不当投加量会导致病原控制失效或对养殖生物产生毒性。本研究建立定量框架以优化工业化凡纳滨对虾(Litopenaeus vannamei)养殖中的消毒过程，整合病原灭活、消毒剂衰减及对虾安全性评估。实验室实验量化了TCCA对致病性弧菌(Vibrio spp.)的灭活动力学，评估了不同水质条件下余氯衰减规律及对仔虾的急性毒性。TCCA在0.3–0.5 mg/L可在3 h内实现弧菌96%–99%灭活，但因余氯下降尤其于低浓度下后期出现细菌复萌。余氯符合一级衰减动力学，受温度、总悬浮物(TSS)、化学需氧量(COD)、pH和溶解氧(DO)显著影响。仔虾96 h半数致死浓度(LC50)为1.319 mg/L，据此计算安全浓度(SC)为0.1319 mg/L。将上述数据集集成至耦合弧菌生长、消毒剂衰减及杀菌效能的系统动力学模型中，可动态模拟病原抑制及消毒后反弹，模型预测与独立实验观测吻合良好。本研究为平衡弧菌防控与对虾安全提供了实用决策支持工具，支持可持续水产养殖的精准消毒管理。

论文解读：工业化凡纳滨对虾养殖中三氯异氰尿酸(TCCA)消毒优化的系统动力学研究

一、研究背景与意义

凡纳滨对虾(Litopenaeus vannamei)集约化养殖系统中高密度环境易导致致病性弧菌(Vibrio spp.，如副溶血弧菌V. parahaemolyticus、哈维氏弧菌V. harveyi等)暴发引起"弧菌病"，造成重大经济损失。三氯异氰尿酸(Trichloroisocyanuric acid, TCCA)因广谱杀菌、成本低被广泛用于水体消毒，但目前投加多依赖经验，缺乏兼顾杀菌效果、余氯衰减及对虾毒性的动态量化模型——剂量不足致弧菌耐药与复发，过量则致对虾中毒及微生态破坏。现有研究多孤立关注杀菌效力、氯衰减或生物毒性某一端点，未整合为交互动态过程。因此研究人员开展本项研究，通过实验测定TCCA对弧菌灭活动力学、受环境影响余氯一级衰减参数及凡纳滨对虾急性毒性LC₅₀，构建耦合弧菌生长—消毒剂衰减—杀菌作用的系统动力学模型，为工业化养殖提供精准消毒决策依据。该论文发表于《Aquaculture Reports》。

二、主要关键技术方法

研究人员选用凡纳滨对虾养殖池常见病原性弧菌(副溶血弧菌、溶藻弧菌V. alginolyticus、坎贝尔弧菌V. campbellii、哈维氏弧菌)及健康仔虾(体长2.84±0.12 cm)为材料；设0、0.1、0.3、0.5 mg/L TCCA处理组进行弧菌杀灭时序实验，采用Chick–Watson模型拟合灭活动力学；通过等对数间距设置TCCA浓度梯度(0.1–5.0 mg/L)开展96 h急性毒性实验，Probit回归计算LC₅₀及安全浓度(SC＝96 h LC₅₀×0.1)；调控温度(20–37℃)、pH(6–9)、COD(5–40 mg/L)、TSS(10–300 mg/L)、DO(5–7 mg/L)单因素变量，按时间点取样测余氯(DPD分光光度法)，一级动力学C_t=C₀e^?kt拟合衰减常数k并引入各因子修正函数；单因子实验测弧菌生长曲线并以修正Gompertz模型获取最大比生长速率(μ_max)和迟滞期(λ)；基于上述参数在Stella Architect平台构建含微生物生长、消毒剂灭活、消毒剂浓度衰减三个子模块的系统动力学模型，定义菌群丰度N(t)与有效TCCA浓度DC为时变状态变量，用常微分方程描述耦合竞争动力学，以相对误差(RE)和平均相对误差(ARE)验证模拟值与实测值吻合度。

三、研究结果

3.1. Bactericidal effects of TCCA against Vibrio spp.

TCCA各处理组初期(1–3 h)均快速降低活菌数，0.1 mg/L组6–12 h出现明显复萌，0.3与0.5 mg/L组3 h内灭活率达96%–99%且延缓反弹。Chick–Watson模型拟合显示V. campbellii(n=0.322)与V. harveyi(n=0.119)符合经典模型(R²>0.8)，V. parahaemolyticus与V. alginolyticus拟合较差(R²=0.685、0.633)且n为负值，故后两者在模型中采用基于实测数据的图形函数描述。

3.2. Acute toxicity of TCCA to juvenile L. vannamei

TCCA呈剂量与时间依赖性致毒，高浓度(≥0.707 mg/L)出现急游、侧卧、体色发白及死亡。24、48、72、96 h LC₅₀分别为3.428、2.265、1.912、1.319 mg/L，据此得96 h安全浓度(SC)为0.1319 mg/L。

3.3. Residual chlorine decay of TCCA and influencing factors

余氯衰减符合一级动力学(R²=0.711–0.975)。温度升高显著加速衰减(k由20℃时0.083 h^?1升至37℃时0.309 h^?1)；TSS与COD升高加速消耗(k由纯水0.075 h^?1分别升至300 mg/L TSS时0.234 h^?1和40 mg/L COD时0.286 h^?1)；pH升高减缓衰减(k由pH 6时0.183 h^?1降至pH 9时0.098 h^?1)；DO升高促进衰减(k由5 mg/L时0.035 h^?1升至7 mg/L时0.165 h^?1)。

3.4. Effects of environmental factors on Vibrio growth

修正Gompertz模型拟合良好(R²=0.90–1.00)。温度25–30℃、pH 7–8、适度盐度(约25‰)及适度有机物(COD)与TSS促进弧菌生长(μ_max升高，λ缩短)；DO影响具菌株特异性(兼性厌氧特性)。各因子校正函数被纳入系统动力学微生物生长子模块。

3.5. Development of the disinfectant reaction kinetics model

微生物生长用Logistic项(带环境承载力K及各环境因子校正的有效比生长速率μ_ef)描述，TCCA灭活用Kc=k_d·(C_dis)^z描述，余氯衰减用dC/dt=?k_ef·C_t、k_ef=k_d·fa(T)·fa(pH)·fa(TSS)·fa(COD)·fa(DO)描述，三者耦合成常微分方程组并在Stella Architect实现。

3.6. Model prediction results

设定典型养殖水质参数(30℃, pH 8, COD 15 mg/L, TSS 50 mg/L, DO 6.5 mg/L, 初始TCCA 0.3 mg/L)运行模型，模拟显示弧菌先降后随余氯衰减部分复萌趋于新平衡，余氯呈指数衰减。独立实验验证显示模拟与实测余氯及弧菌丰度时序趋势一致，ARE较低，模型预测性能可接受。

四、讨论与结论总结

讨论指出0.3–0.5 mg/L TCCA为推荐短期消毒范围——低于此(0.1 mg/L)杀菌不持续易复发，高于此逼近对虾安全阈值；弧菌种间敏感性差异可能与胞外多糖及生物膜有关。96 h SC(0.1319 mg/L)是毒理参考阈值而非推荐用量，生产中应结合余氯监测与水质的修正调整。高温、高TSS、高COD、高DO加速TCCA分解，碱性条件延长余氯稳定，与Arrhenius关系及HOCl/OCl^?解离平衡相符。修正Gompertz参数量化了环境对弧菌增殖的促进/抑制，支持模型中校正因子构建。模型局限含未单独参数化生物膜/胞外聚合物耐受亚群、未纳入底泥交换与昼夜波动等，可作后续扩展。

结论(Concluding remarks)：

工业化对虾养殖中TCCA消毒效力受病原菌种间敏感性差异影响，余氯持久性受水质强烈调控——升温、高悬浮物及高有机负荷加速TCCA源余氯衰减，pH 6–9范围内较高pH增强余氯稳定性，高溶解氧亦促进余氯消耗。同期环境因子调控弧菌生长：升温一般提高μ_max缩短λ，盐度与pH存在最适范围偏离则抑制，适中COD与TSS促进弧菌增殖，DO效应菌株特异。基于上述机制在Stella Architect构建耦合弧菌生长、消毒剂衰减及杀菌效应的系统动力学模型，模拟结果与实验观测时空趋势及量级吻合，ARE验证拟合优度。该模型可模拟多样水质与投药情景下弧菌种群动态，为凡纳滨对虾养殖中精准消毒投加、弧菌感染风险评估及病害防控策略优化提供量化决策支持工具。