动态一致离散时间猴痘传播模型的构建及其在隔离政策评估中的应用

《BMC Public Health》：A dynamically consistent discrete-time model for monkeypox transmission with implications for quarantine policies

猴痘（Mpox）作为一种人畜共患传染病，对公共卫生构成重要挑战。为了解决传统连续时间模型在政策评估中可能出现的数值失真问题，研究人员构建了一种保持正性、有界性和平衡态等关键动力学性质的非标准有限差分（NSFD）离散时间模型。该研究明确计算了基本再生数 ?0，证明了其取决于人（?0h）和动物（?0a）的再生数，并分析指出单凭?0<1未必能确保疾病根除，而?0a>1足以导致人畜间疾病持续。该框架为基于证据的公共卫生决策和隔离策略评估提供了可靠工具。

　　

在公共卫生领域，数学模型是理解和预测传染病传播、评估干预措施效果的强大工具。猴痘（Mpox），一种与天花同属正痘病毒属的人畜共患病毒性疾病，近年来引起了全球范围的关注。尽管已有大量研究利用微分方程（连续时间模型）来模拟猴痘的传播动态，但在实际应用中，公共卫生数据的报告（如每日新增病例数）和干预措施的实施（如定期评估并调整隔离政策）往往发生在离散的时间点上。这使得连续时间模型在直接对接现实决策时，有时会显得有些“水土不服”。更为关键的是，当使用标准的数值方法（如欧拉法或龙格-库塔法）来求解这些连续模型时，可能会产生数值失真，例如出现不合理的负人口数或不稳定的解，从而误导政策判断。那么，能否构建一种既贴合公共卫生决策离散时间步特点，又能忠实反映疾病传播核心数学性质的模型框架呢？这正是本项研究试图回答的核心问题。

为此，研究团队在《BMC Public Health》上发表了一项研究，题为“A dynamically consistent discrete-time model for monkeypox transmission with implications for quarantine policies”。他们开发了一个离散时间的猴痘传播模型。其创新之处在于，他们并非简单地对连续模型进行离散化，而是精心设计了一种“非标准有限差分”（Nonstandard Finite Difference, NSFD）方案。这个方案的核心优势是能够严格保持对应连续模型的关键“动力学一致性”，包括解的正性（人口数不为负）、有界性（总人口不超过设定范围）以及平衡态（如疾病消亡或流行状态）的精确数学行为。这确保了模型在作为决策支持工具时，其数值模拟是稳健可靠的，不会因计算方法的缺陷而产生误导性结论。

为了开展这项研究，作者主要运用了以下关键方法：首先，构建了一个包含人和动物种群的仓室模型（Compartmental Model），并据此建立了非线性差分方程系统。其次，核心方法是设计了能保持动力学一致性的非标准有限差分（NSFD）离散化方案。再次，运用下一代矩阵（Next Generation Matrix）方法，在离散框架下严格推导并计算了模型的基本再生数?₀及其人与动物分量?₀^h和?₀^a。此外，还进行了系统的数学分析，包括平衡点存在性、局部及全局稳定性分析。最后，通过数值模拟，对比了NSFD方法与标准欧拉（Euler）法和龙格-库塔（Runge-Kutta）法的表现，并进行了参数敏感性分析以评估关键因素对传播的影响。

研究结果揭示了以下几个重要发现：

1. 1.
   模型框架与基本再生数：研究成功构建了一个动力学一致的离散时间猴痘传播模型。分析显示，系统总的基本再生数?₀可以明确表达为人类基本再生数?₀^h和动物基本再生数?₀^a的函数，具体形式为?₀= max{?₀^h, ?₀^a}。这直观地表明，疾病的传播潜力由人和动物宿主中更强的那个传播环节所主导。
2. 2.
   平衡态与稳定性分析：模型识别出三个重要的平衡态：无病平衡点（疾病在人和动物中均不存在）、仅动物地方性平衡点（仅在动物中持续流行）以及人和动物地方性平衡点（在人和动物中均持续流行）。稳定性分析表明，当?₀< 1时，无病平衡点是局部渐近稳定的。然而，研究有一个关键发现：仅凭?₀< 1可能不足以保证疾病被根除。这是因为即使总体?₀< 1，如果动物宿主内部的传播动力足够强，疾病仍可能在动物种群中持续存在，从而构成一个持续的传染源。
3. 3.
   疾病持续的条件：研究证明了一个强有力的充分条件：如果动物基本再生数?₀^a> 1，那么无论人类中的情况如何，疾病都将在人和动物种群中持续存在。这一发现凸显了动物宿主在猴痘流行病学中的核心作用。控制人类间的传播（降低?₀^h）若不伴随对动物源头的有效控制，可能无法实现消除疾病的目标。
4. 4.
   干预措施的影响：通过敏感性分析，研究量化了模型参数对疾病传播的影响。分析清晰地阐释了隔离措施和减少接触策略如何通过影响有效接触率，进而改变?₀^h和?₀^a，最终决定疾病是否能在人群中持续。这为评估不同隔离政策的潜在效果提供了定量的理论依据。
5. 5.
   数值方法的优越性：数值模拟验证了理论结果，并对比了所提出的NSFD方案与标准欧拉法和龙格-库塔法。结果表明，NSFD方案在不同时间步长下都能保持解的物理合理性和数值稳定性，而标准方法在较大步长时会出现违反正性或产生数值振荡等伪影，从而证明了NSFD框架在政策模拟应用中的可靠性和稳健性。

归纳研究的结论与讨论，本研究的核心贡献在于提出了一个为公共卫生决策“量身定制”的离散时间建模框架。它不仅从数学上保证了模型结果的可靠性（动力学一致性），而且其离散时间步的特性与真实的公共卫生数据上报和决策周期天然契合。研究得出的最重要结论是：在猴痘防控中，必须高度重视动物宿主的作用。单纯将人类群体中的基本再生数?₀^h控制在1以下，并不足以宣告胜利，如果动物源头（?₀^a）未被有效控制，疾病仍将持续威胁。这强烈提示，成功的防控策略需要人与动物部门联防联控，采取“One Health”（一体健康）的综合治理 approach。此外，该模型框架作为一个可靠的数字沙盘，能够用于模拟和比较不同隔离强度、隔离时长、接触追踪效率等策略的效果，从而为制定基于证据的、成本效益优化的公共卫生政策提供了有力的定量分析工具。这项研究将扎实的数学理论、稳健的计算方法与紧迫的公共卫生决策需求相结合，为应对猴痘等人畜共患传染病的挑战提供了新的方法论支持。