含治疗失败与自然恢复随机SITRS霍乱模型：框架、结果与公共卫生启示

《Franklin Open》：A stochastic SITRS cholera model with treatment failure and self-recovery: Framework, results, and public health insights

　　

霍乱在资源有限环境中仍是持续的公共卫生威胁，治疗可及性不一、自然恢复常见且环境条件高度不确定。现有模型很少同时捕捉这三个因素。研究人员构建了随机SITRS（易感-感染-治疗-恢复-易感）霍乱模型，纳入自然自我恢复（γ2）、治疗失败（?）以及治疗与未治疗病例的差异死亡率（τ1≠τ2），并将环境变异性以乘性噪声形式嵌入。分析得到两项关键发现：其一，环境噪声降低暴发潜力（R0S<R0D），即随机性对疾病传播起天然但温和的制动作用；其二，敏感性分析识别接触率（α）与人口流入（π）为主要传播驱动因子，而治疗可及性（σ）是最有效的控制参数。模拟显示，在不确定条件下临床恢复比自然恢复提供更可预测的暴发管理。对公共卫生规划者而言，最优策略是接触降低（卫生运动、社交距离）与改善治疗可及性相结合——单一干预次优。本工作为治疗失败与自然恢复并存的真实不确定场景下的霍乱防控提供可操作指导。

研究背景方面，霍乱由霍乱弧菌（Vibrio cholerae）引起，在资源有限社区中若不及时治疗病死率超50%，而及时口服补液或静脉输液可降至1%以下。传统模型强调经污染水源的环境传播，但直接接触（照护、共餐）同为重要途径却常被低估。现有数学模型演进虽已纳入环境储库、空间异质性、疫苗接种、媒体影响与最优控制，但仍存在三方面局限：聚焦直接传播者缺乏临床现实性；多为确定性模型未能刻画环境不确定性；自我恢复与治疗失败等临床通路极少被统一整合。临床上轻症可不经治疗自发恢复（γ₂），治疗失败（?）因药物无效或疗程中断使患者回流感染池，且未治疗（τ₁）与治疗者（τ₂）病死率差异显著。本研究首次在SITRS框架内同步纳入自我恢复、治疗失败与差异死亡率，推导随机再生数（R₀^S）并比较临床与自然恢复路径，论文发表于《Franklin Open》。

关键技术方法上，研究人员基于同质混合假设建立四仓室常微分方程组（S-I-T-R），将感染仅归因为直接人际接触（率α），并把治疗者视为理想隔离无传播；随后引入独立标准布朗运动B_i(t)与噪声强度β_i将确定性系统扩展为It?意义下随机微分方程（SDE），噪声以乘性形式作用于各仓室。参数取自流行病学文献与假定值（见表1），通过正性全局解存在唯一性证明、确定性R₀^D与随机R₀^S推导、局部与全局稳定性分析、PRCC敏感性分析以及多情景数值模拟展开研究，样本并非实证队列而是参数化数值实验。

研究结果各节要点如下。第1节引言阐明直接传播设定下临床现实缺失与随机性不足，提出本模型新颖性。第2节模型构建给出确定性ODE与随机SDE全方程组及流程图，列明九项假设（同质混合、仅直接传播、治疗者隔离、常数参数、指数逗留时间等）及其可松弛方向，并专述直接传播聚焦的三点理由（拥挤场所相关性、临床通路对齐、补环境模型之不足）。第3节定性分析证明解的正性与有界性，由下一代矩阵得R₀^D=α/(μ+τ₁+σ+γ₂)，随机情形得R₀^S<R₀^D的解析条件，并证Disease-free平衡与Endemic平衡的稳定性。第4节数值模拟与敏感性分析显示：α与π为敏感正向驱动，σ为最强负向控制参数；提高σ比提高γ₂更能压缩峰值；随β_i增大暴发更受抑但轨迹波动加剧；临床恢复路径方差小于自然恢复路径。第5节结论指出单一干预次优，接触削减+治疗可及性组合最优，并承认忽略环境储库与治疗期间部分传播为局限。

讨论与结论翻译部分，研究人员总结：环境噪声作为温和自然制动使R₀^S<R₀^D，意味着确定性模型会高估暴发风险；接触率α与招募率π是主要传播驱动，治疗率σ是最具杠杆的控制旋钮；在资源约束下应并行推进卫生运动（降α）与ORT/静脉补液覆盖（升σ），而非独押自然恢复或单纯隔离。治疗失败（?）与自然恢复（γ₂）并存时，临床恢复提供比自然恢复更可预测的疫情管理。该SITRS随机框架为治疗失败与自我恢复共现、且环境高度不确定的现场霍乱防控提供了可转化数学依据。