奶牛乳腺炎动态数学模型的最优控制与成本效益分析

《AppliedMath》：Optimal Control and Cost-Effectiveness Analysis of a Mathematical Model for Mastitis Dynamics in Dairy Cows

　　

乳腺炎仍是影响奶牛最重要的传染病之一，对动物健康、产奶量及牧场盈利能力造成显著影响。研究人员在既有非线性SIRS–P（易感-感染-恢复-易感并含环境病原体P分室）乳腺炎传播模型基础上，引入敏感性分析、最优控制与成本效益分析。模型包含易感牛、感染牛、恢复牛群体及环境病原体分室，兼顾直接牛-牛传播、经环境污染的间接传播、恢复、感染复发、补栏与淘汰。研究人员对基本再生数R0进行敏感性分析以识别对传播影响最强的参数；进而建立含时变控制的最优控制问题，控制项分别降低直接传播、降低环境污染、提升恢复率，应用庞特里亚金极大值原理（Pontryagin’s Maximum Principle）推导最优性系统，并采用前向-后向扫描法数值求解。通过数值模拟与ACER（平均成本效益比）、ICER（增量成本效益比）比较七种干预策略。结果表明完全整合策略S7（预防筛查+环境卫生+治疗）疾病降幅最大；但基于直接传播削减的策略在降低乳腺炎感染负担上资源分配最具成本效益。

研究背景与立项依据

乳腺炎是奶牛养殖中危害最严重的传染病之一，导致乳腺炎症、产奶量下降、奶质恶化、治疗费用上升与提前淘汰，兼具经济与动物福利双重负担。其流行病学由牛-牛接触、环境污染、挤奶卫生、病原菌脱落、泌乳阶段与管理实践共同驱动，且感染可复发、环境病原可持续储留。既有数学研究多聚焦传播动态与稳定性，缺少年内时变干预的经济学评价。研究人员在Zlatanovska等提出的SIRS–P框架上，补足敏感性分析、最优控制与ACER/ICER成本效益比较，以回答“何种干预组合既流行病学有效又经济高效”。论文发表于《AppliedMath》。

主要关键技术方法

研究人员采用四类方法：其一，对基本再生数R₀做归一化正向敏感性分析，量化参数相对变动对R₀的弹性；其二，在常微分SIRS–P系统中嵌入三个Lebesgue可测有界时变控制u₁（预防筛查降直接传播）、u₂（环境卫生降环境传播并增菌去除）、u₃（治疗升恢复率），构造含状态惩罚与二次控制代价的泛函；其三，用庞特里亚金极大值原理推导伴随系统与投影型最优控制表达式，借前向-后向扫描法在MATLAB离散时间网格数值求解；其四，基于模拟轨迹定义感染与环境负担减少量及总实施成本，计算ACER与ICER，按有效性排序剔除被支配与扩展支配策略，得到成本效益前沿。参数与初值沿用文献基线，无农场级队列校准。

研究结果

1. 2.
   Model Formulation：给出SIRS–P四分室（S,I,R,P）非线性ODE，含补栏Λ、自然淘汰μ、直接传播β₁SI/N、环境传播β₂SP/N、恢复γ、复发σ、脱菌η、环境菌衰减δ；明确各参数生物学意义，确立无控系统 endemic 持续条件。
2. 3.
   Sensitivity Analysis of the Basic Reproduction Number：推得R₀表达式并算归一化敏感性指数。Λ与β₁正向最强（β₁指数约0.5量级），β₂、ω（脱落）为正，δ、μ、γ为负，σ不显含于R₀故指数为零。结论：直接传播、补栏、环境菌负荷、恢复与淘汰是调控重点，过度淘汰不经济。
3. 4.
   Optimal Control Problem：引入u₁,u₂,u₃∈[0,u_i^max]，改直接项(1-u₁)β₁、环境项(1-u₂)β₂并加δ₁u₂去菌与γ₁u₃增恢复；泛函J积分加权I与P加二次代价。证存在性（凸闭有界控制集、状态有界、Arzelà–Ascoli弱收敛），由极大值原理得伴随方程终值λ_i(T)=0与u_i=clip(·,[0,u_i^max])投影式。
4. 5.
   Numerical Simulations：七策略S₁–S₇（单用u₁/u₂/u₃至全用三者）。无控时I与P维持endemic水平；S₁（仅u₁）大幅压低I与继发P；S₂仅降P对I有限；S₃升R降I；组合策略较单控更优，S₇同时抑I与P、保S最高。结论：多通路联合优于单点干预。
5. 6.
   Cost-Effectiveness Analysis：依感染减负排ACER/ICER，S₂,S₃,S₅严格支配，S₄扩展支配，前沿剩S₁,S₆,S₇；S₇感染减负最大但S₁感染ACER最低（每单位减负成本最小），S₇环境ACER亦最低。结论：S₇流行病学最优，S₁降感染负担最具性价比。
6. 7.
   Conclusions翻译：研究人员对SIRS–P乳腺炎模型做最优控制与成本效益分析，R₀敏感性指出补栏、直接传播、淘汰、恢复主导潜力；三时变控制下用庞特里亚金原理与前向-后向扫描得数值解；S₇整合预防筛查、环境卫生与治疗降幅最大，S₁仅削直接传播在感染ACER意义上最成本效益。实践应重早检、筛查、挤奶卫生、环境清消、病牛分隔与治疗。局限为参数来自文献非农场纵向校准，后续可加随机性、季节、分数阶与现场估计。

讨论总结

研究人员强调模型在定性指导而非场级预测的价值：敏感性分析将控制靶点收敛至直接接触阻断与恢复提升，最优控制显示时变强度优于恒强干预，成本效益分析分离“最有效”与“最划算”——完全整合策略适合追求疾病清零的管理目标，单靠预防筛查策略在预算受限时单位成本收益最高。文末指明未做参数辨识、未含随机与记忆效应，为后续场数据驱动与分数阶扩展留接口。