DENGUE TRANSMISSION INTERVENTION MODEL

登革热传播干预模拟器

温度驱动的人群 SEIR—伊蚊 SEI 确定性仓室模型,用于比较不同媒介防控、Wolbachia 与疫苗接种策略下的本地干预效果。

v1.1.6 服务检测中
适用说明:本工具用于比较不同媒介防控、Wolbachia 与疫苗接种策略在同一传播背景下的相对效果与实施时机。参数默认值依据公开文献[1–30]设定;应用于具体地区时,请依据当地血清学调查、伊蚊监测、病例报告与气象资料,按各参数下方“本地取值”说明调整,并以当地历史周病例数据核对曲线形态与量级。
预设干预
导出周序列(CSV)

参数拟合(本地化)

可选 · 以当地周报告病例序列估计雌蚊/人比、输入强度与报告率,并回填干预参数

可选

1

导入周报告病例数据

2

选择待估参数

主要结果

基线(无干预)与所选干预比较;病例数为模型报告病例数(= 发病数 × 报告率)。

基线累计报告病例—
干预累计报告病例—
干预峰值周报告病例—
累计病例相对减少—

—

周报告病例数

每周新发报告病例数(非累计)。藏青实线为基线(无干预);各干预(含当前所选与已保存干预)均以同一虚线显示,按颜色与图例区分;阴影为当前干预实施期。纵轴可切换线性或对数坐标。

计算中…

传播强度

左:基本再生数 R₀ 与有效再生数 Rₜ 的周变化(Rₜ = R₀ × √(Sₕ/Nₕ);Rₜ > 1 表示传播增长)。浅藏青实线为基线 R₀、藏青实线为基线 Rₜ,各干预 Rₜ 以同一虚线按颜色区分。右:R₀ 随日均气温变化(按当前雌蚊/人比)。

干预对比表

干预累计报告病例报告发病率(/10万)峰值周病例峰值周次较基线减少
尚无结果

AI 助手

依据本模型的参数设定与所列参考文献,解答模拟器使用、参数取值与本地化、结果解读及登革热流行病学与防控问题。

0/500
回答由大模型依据模型参数与文献生成,供专业参考,请结合当地监测资料研判;不提供个人诊疗建议。

模型结构

人群 SEIR 与伊蚊 SEI 两组仓室通过叮咬耦合;展开查看仓室示意图与方程。

登革热人群 SEIR—伊蚊 SEI 仓室模型结构示意图人群(宿主)伊蚊(媒介)Sₕ易感者Eₕ潜伏期感染者Iₕ感染者(有传染性)Rₕ免疫者λₕσ = 1/5.9 dγ = 1/5 d输入病例 ι(例/周)周报告病例= ρ · ΔC(发病累计)初始免疫 + 疫苗 VE×覆盖率K(t)成蚊容纳量(孳生地)Sᵥ易感伊蚊Eᵥ外潜伏期伊蚊Iᵥ感染性伊蚊羽化 μ(T)·KλᵥPDR(T)死亡 μ(T)死亡 μ(T)死亡 μ(T)K = m · M(T) · Nₕ成蚊密度↓ · 孳生地↓① 蚊→人:λₕ = a·b·Iᵥ / Nₕ人蚊接触↓ → a · Wolbachia → b② 人→蚊:λᵥ = a·pMI·Iₕ / Nₕ人蚊接触↓ → a · Wolbachia → pMI状态转移交叉传播橙色:随气温变化的参数(Caldwell 2021 性状函数)绿色:干预措施作用点 登革热人群 SEIR—伊蚊 SEI 仓室模型结构示意图人群(宿主)伊蚊(媒介)Sₕ易感者Eₕ潜伏期感染者Iₕ感染者Rₕ免疫者λₕ ①σ = 1/5.9 dγ = 1/5 d输入 ι→ 周报告 = ρ·ΔC初始免疫 + 疫苗K(t)成蚊容纳量Sᵥ易感伊蚊Eᵥ外潜伏期伊蚊Iᵥ感染性伊蚊羽化λᵥ ②PDR(T)K = m·M(T)·Nₕ死亡 μ(T)① 蚊→人:λₕ = a·b·Iᵥ / Nₕ② 人→蚊:λᵥ = a·pMI·Iₕ / Nₕ人蚊接触↓ → aWolbachia → b、pMI成蚊密度↓ · 孳生地↓ → K橙色:随气温变化的参数;虚线:交叉传播伊蚊各仓室均按 μ(T) 死亡;按 μ·K 羽化补充
模型方程

人群(Nₕ = Sₕ + Eₕ + Iₕ + Rₕ)

dSₕ/dt = −λₕ·Sₕ + μₕ(Nₕ − Sₕ)
dEₕ/dt =  λₕ·Sₕ − (σ + μₕ)Eₕ + ι
dIₕ/dt =  σ·Eₕ − (γ + μₕ)Iₕ
dRₕ/dt =  γ·Iₕ − μₕ·Rₕ
dC/dt  =  σ·Eₕ          周报告病例 = ρ·[C(t+7) − C(t)]
λₕ = a(T)·b(T)·Iᵥ/Nₕ

伊蚊(Nᵥ = Sᵥ + Eᵥ + Iᵥ;K(t) = m·M(T)·Nₕ,M(T) ∝ EFD·pEA·MDR/μ²,29 °C 时为 1)

dSᵥ/dt = μ(T)·K − λᵥ·Sᵥ − μ(T)·Sᵥ − δ·Sᵥ
dEᵥ/dt = λᵥ·Sᵥ − (PDR(T) + μ(T))·Eᵥ − δ·Eᵥ
dIᵥ/dt = PDR(T)·Eᵥ − μ(T)·Iᵥ − δ·Iᵥ
λᵥ = a(T)·pMI(T)·Iₕ/Nₕ        δ = μ(T)·max(Nᵥ/K − 1, 0)

再生数

R₀² = [a·b·m·M(T)·PDR / ((PDR + μ)·μ)] × [a·pMI / (γ + μₕ) · σ/(σ + μₕ)]
Rₜ  = R₀ · √(Sₕ/Nₕ)

干预:成蚊密度与孳生地措施使 K 乘以 (1 − 降幅);人蚊接触措施使 a 乘以 (1 − 降幅),均限于干预期;Wolbachia 使 b、pMI 乘以 (1 − 传病阻断效力 × 建立率)[30];疫苗于起始时将 VE × 覆盖率 比例的易感者移入 Rₕ。σ = 1/5.9 d[3],γ = 1/5 d,μₕ = 1/(75 × 365) d⁻¹;温度依赖性状取自 Caldwell 等[1],R₀(T) 形式参照 Mordecai 等[2]。

方法说明、适用范围与参考文献