1. 济宁医学院公共卫生学院;
2. 济宁医学院口腔医学院, 济宁 272013
收稿日期:2021-03-26
基金项目:山东省高等学校科技计划项目(J16LD53)
Impact of air pollutants on influenza incidence in Jining city from 2009 to 2017 and its early warning
1. School of Public Health, Jining Medical University;
2. School of Stomatology, Jining Medical University, Jining 272013, China
流行性感冒(简称流感,influenza)是流感病毒引起的一种常见的急性呼吸道传染病,具有耐低温,传染性强,传播迅速等特点[1]。流感发病率受气候、大气污染物、受感人群等多种因素影响,是我国乃至全球高度关注的公共卫生问题[2]。20世纪以来,作为一种世界范围内流行的传染病,流感曾4次席卷全球,其中3次被认为是从中国开始爆发的[3]。每一次流感大爆发都给中国乃至世界造成不可磨灭的损失。目前普遍认为影响人体健康的大气污染物,主要包括PM2.5、PM10、SO2、NO2、CO、O3,这6种大气污染物常作为评价环境空气质量的参考指标,按照其对人体健康影响大小排序,依次为PM2.5、SO2、NOx[4]。已有研究表明大气污染物与人体呼吸系统疾病密切相关[5],易诱发鼻炎、咳嗽、哮喘等呼吸系统疾病[6]。本文通过对2014-2017年济宁市流感发病的流行病学特征进行分析,探讨大气污染物对流感发病的影响,为济宁市今后流感流行的预警预测提供依据。
1 资料与方法
1.1 资料来源
1.1.1 病例资料
通过中华人民共和国法定传染病网络直报系统,获得济宁市2014年1月1日-2017年12月31日确诊上报的流感病例,收集病例的年龄、性别、发病日期、疾病名称等信息。
1.1.2 污染物数据
通过中华人民共和国生态环境部数据中心收集济宁市2014年1月1日-2017年12月31日每日空气质量指数(air quality index,AQI)及每日大气污染物数据,包括每日PM2.5、PM10、SO2、NO2及CO的24h平均浓度,O3日最大8h平均浓度及AQI。
1.2 方法
每日大气污染物对每日流感发病影响分析:采用双向时间分层病例交叉设计与DLNM模型相结合[7],对济宁市每日流感病例发病数与每日PM2.5、PM10、SO2、NO2、CO和O3浓度之间的关系进行研究。对照选择及滞后时间设定同气温效应分析,将6种大气污染浓度分别作为直线变量引入模型,建立单一污染物模型,具体如下。
$
\begin{array}{l}
\log \left[ {{\rm{E}}\left( {{{\rm{Y}}_{\rm{t}}}} \right)} \right] = {\rm{ \mathsf{ α} }} + {{\rm{ \mathsf{ β} }}_{\rm{i}}}{{\rm{Z}}_{{\rm{i, t, l}}}} + {\rm{NS}}\left( {{\rm{do}}{{\rm{w}}_{\rm{t}}}} \right) \\
+ {\rm{NS}}\left( {{{\rm{T}}_{\rm{t}}}} \right) + {\rm{ \mathsf{ λ} }}\left( {{\rm{strat}}{{\rm{a}}_{\rm{t}}}} \right)
\end{array}
$
|
模型中Yt表示观察日t的流感新发病人数;E(Yt)表示观察日t当天的预期发病人数;α表示截距;Zi, t, l表示第i种污染物,由DLNM得到的滞后l天的浓度矩阵;βi表示第i种污染物浓度矩阵的回归系数;NS表示自然立方样条函数;dowt表示观察日t对应的露点温度;Tt表示观察日t对应的平均气温;stratat表示观察日t匹配的当年当月同一星期几产生的分层变量;λ表示stratat的回归系数。
根据我国《环境空气质量标准》(GB3095-2012)中的二类区域浓度限值,对各大气污染物浓度进行评价[8]。
1.3 统计学方法
采用Excel 2010建立每日流感发病数与大气污染物逐日数据库。R3.4.3软件完成所有数据分类汇总整理及变量间的简单相关分析。采用DLNM软件包建立每日流感发病数与不同滞后天大气污染物浓度交叉矩阵。所有统计分析均采用双侧检验,以P<0.05为有统计学意义。
2 结果
2.1 济宁市流感发病的流行病学特征
2014-2017年济宁市流感患病人数共1981例,其中男性1049例(52.95%),女性932例(47.05%)。0岁~患者193例(9.74%),3岁~患者550例(占27.76%),15岁~患者1238例(62.49%)。见表 1。
表 1
表 1 流感病例不同性别、年龄构成
变量 |
病例数 |
构成比/% |
性别 |
|
男 |
1049 |
52.95 |
女 |
932 |
47.05 |
年龄/岁 |
|
0~ |
193 |
9.74 |
3~ |
550 |
27.76 |
15~ |
1238 |
62.49 |
|
表 1 流感病例不同性别、年龄构成
|
2.2 大气污染物对流感发病的影响
2.2.1 2014-2017年济宁市大气污染物状况
2014-2017年济宁市AQI优、良、轻度污染、重度污染的天数分别为62d(3.4%)、1468d(81.5%)、214d(11.9%)、57d(3.2%),空气质量以良为主。
2014-2017年大气污染物PM2.5、PM10、SO2、CO、NO2的每日平均浓度分别是74μg/m3、131μg/m3、50μg/m3、1.2mg/m3、44mg/m3,超标总天数依次为518d、432d、20d、1d、48d。其中,SO2浓度均数逐年降低,达标率逐年上升,2017年达标率为100.0%。NO2的浓度均数变化不大,达标率较平稳。PM2.5和PM10达标率虽逐年上升,但在总体上仍低于其它污染物。CO达标率除2015年外,均为100.0%。O3的日最大8h平均浓度为106μg/m3,超标的天数为280d,年达标率78.9%~83.3%。详见表 2。
表 2
表 2 2014-2017年济宁市主要空气污染物年份评价
污染物 |
年份 |
浓度(μg/m3,CO为mg/m3) |
达标率/% |
min |
max |
x |
s |
SO2 |
2014 |
14 |
269 |
72 |
38 |
95.6 |
|
2015 |
14 |
166 |
57 |
27 |
99.2 |
|
2016 |
10 |
154 |
42 |
25 |
99.7 |
|
2017 |
6 |
81 |
27 |
13 |
100.0 |
NO2 |
2014 |
15 |
108 |
47 |
17 |
97.0 |
|
2015 |
9 |
106 |
44 |
17 |
95.3 |
|
2016 |
9 |
107 |
42 |
18 |
96.7 |
|
2017 |
12 |
114 |
41 |
17 |
98.1 |
PM10 |
2014 |
29 |
517 |
155 |
75 |
56.2 |
|
2015 |
34 |
594 |
140 |
69 |
65.2 |
|
2016 |
30 |
529 |
115 |
69 |
77.6 |
|
2017 |
27 |
758 |
113 |
63 |
82.7 |
PM2.5 |
2014 |
24 |
334 |
88 |
46 |
46.6 |
|
2015 |
19 |
368 |
81 |
45 |
58.6 |
|
2016 |
15 |
370 |
70 |
49 |
72.4 |
|
2017 |
14 |
259 |
57 |
37 |
80.5 |
CO |
2014 |
1.0 |
4.0 |
1.3 |
0.5 |
100.0 |
|
2015 |
1.0 |
5.0 |
1.3 |
0.5 |
99.7 |
|
2016 |
1.0 |
4.0 |
1.2 |
0.5 |
100.0 |
|
2017 |
0.0 |
4.0 |
1.1 |
0.5 |
100.0 |
O3 |
2014 |
9 |
255 |
108 |
58 |
79.5 |
|
2015 |
8 |
258 |
102 |
56 |
81.6 |
|
2016 |
9 |
223 |
100 |
52 |
83.3 |
|
2017 |
13 |
272 |
115 |
58 |
78.9 |
|
表 2 2014-2017年济宁市主要空气污染物年份评价
|
2.2.2 PM2.5对全人群和不同年龄段人群流感日发病人数的滞后效应
日平均浓度每增加10μg/m3,对全人群和15岁~人群流感日发病人数的最大RR值分别为1.03(1.00,1.05)、1.04(1.02,1.07),都出现在暴露后11d;对0岁~人群最大RR是1.06(1.03,1.09),出现在暴露后14d;3岁~人群RR值最大为1.03(1.02,1.07),出现在暴露后13d。详见表 3。
表 3
表 3 0~14d滞后PM2.5日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
0.99(0.98,1.01) |
0.98(0.95,1.02) |
1.00(0.98,1.03) |
0.99(0.97,1.02) |
lag1 |
1.00(0.98,1.02) |
0.93(0.90,0.97) |
0.98(0.96,1.01) |
1.02(0.99,1.05) |
lag2 |
0.99(0.97,1.01) |
1.05(1.01,1.08)* |
1.01(0.99,1.04) |
0.96(0.94,0.99) |
lag3 |
1.00(0.98,1.02) |
0.93(0.90,0.96) |
1.01(0.98,1.03) |
1.01(0.98,1.04) |
lag4 |
1.02(1.00,1.04)* |
1.01(0.97,1.04) |
1.02(0.99,1.04) |
1.03(1.00,1.05)* |
lag5 |
0.99(0.97,1.02) |
1.05(1.02,1.09)* |
1.00(0.97,1.02) |
0.98(0.96,1.01) |
lag6 |
1.00(0.98,1.02) |
0.99(0.96,1.03) |
1.02(0.99,1.04) |
0.99(0.96,1.02) |
lag7 |
0.98(0.96,1.00) |
0.97(0.94,1.01) |
0.99(0.96,1.02) |
0.97(0.95,1.00) |
lag8 |
1.00(0.98,1.02) |
1.00(0.96,1.03) |
0.99(0.97,1.02) |
1.01(0.98,1.03) |
lag9 |
1.01(0.99,1.03) |
1.00(0.97,1.04) |
1.02(0.99,1.05) |
1.01(0.98,1.03) |
lag10 |
0.99(0.96,1.01) |
0.96(0.93,1.00) |
1.00(0.97,1.03) |
0.98(0.96,1.01) |
lag11 |
1.03(1.00,1.05)* |
1.03(1.00,1.07)* |
0.99(0.96,1.02) |
1.04(1.02,1.07)* |
lag12 |
0.99(0.97,1.01) |
0.99(0.96,1.03) |
1.01(0.98,1.04) |
0.98(0.96,1.01) |
lag13 |
1.02(1.00,1.04)* |
1.00(0.97,1.04) |
1.03(1.00,1.06)* |
1.01(0.99,1.04) |
lag14 |
1.01(0.99,1.03) |
1.06(1.03,1.09)* |
1.01(0.99,1.04) |
1.00(0.98,1.02) |
注:*为P<0.05 |
|
表 3 0~14d滞后PM2.5日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
|
2.2.3 PM10对全人群和不同年龄段人群流感日发病人数的滞后效应
日平均浓度每增加10μg/m3,对全人群的最大RR值出现在暴露后4d为1.01(1.00,1.03)和11d为1.01(1.00,1.03),对0岁~、3岁~和15岁~人群的最大RR值分别出现在暴露后14d为1.03(1.02,1.06)、12d为1.02(1.00,1.04)和11d为1.02(1.00,1.04)。详见表 4。
表 4
表 4 0~14d滞后PM10日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
1.00(0.98,1.01) |
0.98(0.96,1.00) |
1.00(0.98,1.02) |
1.00(0.98,1.01) |
lag1 |
1.00(0.98,1.01) |
0.96(0.93,0.98) |
0.98(0.96,1.00) |
1.01(0.99,1.03) |
lag2 |
0.99(0.98,1.00) |
1.03(1.01,1.05)* |
1.01(0.99,1.02) |
0.97(0.96,0.99) |
lag3 |
1.00(0.99,1.02) |
0.95(0.93,0.97) |
1.01(0.99,1.03) |
1.01(0.99,1.03) |
lag4 |
1.01(1.00,1.03)* |
1.01(0.99,1.03) |
1.01(0.99,1.03) |
1.01(1.00,1.03)* |
lag5 |
1.00(0.99,1.01) |
1.03(1.01,1.05)* |
1.00(0.98,1.02) |
0.99(0.98,1.01) |
lag6 |
1.00(0.98,1.01) |
0.99(0.97,1.01) |
1.02(1.00,1.04)* |
0.99(0.97,1.01) |
lag7 |
0.99(0.97,1.00) |
0.99(0.97,1.01) |
0.98(0.96,1.00) |
0.99(0.97,1.01) |
lag8 |
1.00(0.99,1.02) |
0.98(0.96,1.01) |
1.00(0.98,1.02) |
1.01(0.99,1.02) |
lag9 |
1.00(0.99,1.02) |
1.00(0.98,1.02) |
1.01(0.99,1.02) |
1.00(0.98,1.01) |
lag10 |
1.00(0.98,1.01) |
0.98(0.96,1.01) |
1.00(0.98,1.02) |
1.00(0.98,1.01) |
lag11 |
1.01(1.00,1.03)* |
1.00(0.98,1.02) |
0.99(0.97,1.01) |
1.02(1.00,1.04)* |
lag12 |
1.00(0.99,1.02) |
1.01(0.99,1.03) |
1.02(1.00,1.03)* |
1.00(0.98,1.01) |
lag13 |
1.01(1.00,1.02)* |
0.99(0.97,1.01) |
1.01(0.99,1.03) |
1.01(1.00,1.03)* |
lag14 |
1.01(0.99,1.02) |
1.04(1.02,1.06)* |
1.01(0.99,1.02) |
1.00(0.99,1.02) |
注:*为P<0.05 |
|
表 4 0~14d滞后PM10日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
|
2.2.4 SO2对全人群和不同年龄段人群流感日发病人数的滞后效应
日平均浓度每增加10μg/m3,对全人群、0岁~、3岁~和15岁~人群最大RR分别是1.05(1.01,1.09)、1.04(0.99,1.10)、1.084(1.03,1.14)和1.06(1.02,1.11),暴露天数分别是11、8、2和4d。详见表 5。
表 5
表 5 0~14d滞后SO2日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
1.00(0.96,1.04) |
0.97(0.92,1.02) |
1.01(0.96,1.06) |
1.00(0.96,1.04) |
lag1 |
1.00(0.97,1.04) |
1.03(0.97,1.09) |
1.02(0.97,1.07) |
1.00(0.96,1.04) |
lag2 |
1.00(0.96,1.04) |
0.93(0.88,0.99) |
1.08(1.03,1.14)* |
0.98(0.94,1.02) |
lag3 |
1.00(0.97,1.04) |
0.97(0.92,1.03) |
0.98(0.93,1.02) |
1.01(0.97,1.06) |
lag4 |
1.04(1.01,1.08)* |
1.03(0.98,1.09) |
1.02(0.97,1.07) |
1.06(1.02,1.11)* |
lag5 |
0.98(0.94,1.02) |
0.99(0.93,1.05) |
0.98(0.93,1.03) |
0.97(0.92,1.01) |
lag6 |
1.01(0.97,1.05) |
1.02(0.97,1.08) |
1.03(0.98,1.09) |
1.01(0.96,1.05) |
lag7 |
0.98(0.94,1.02) |
0.94(0.89,0.99) |
1.00(0.95,1.05) |
0.97(0.93,1.02) |
lag8 |
0.99(0.96,1.03) |
1.04(0.99,1.10) |
1.00(0.96,1.05) |
0.98(0.94,1.02) |
lag9 |
1.01(0.98,1.05) |
0.95(0.90,1.00) |
1.06(1.02,1.11)* |
1.01(0.96,1.05) |
lag10 |
0.96(0.93,1.00) |
0.97(0.91,1.03) |
0.95(0.91,1.00) |
0.97(0.92,1.01) |
lag11 |
1.05(1.01,1.09)* |
1.01(0.95,1.07) |
1.04(0.99,1.09) |
1.06(1.02,1.11)* |
lag12 |
1.00(0.97,1.04) |
0.94(0.88,0.99) |
1.02(0.97,1.07) |
1.00(0.96,1.04) |
lag13 |
1.02(0.98,1.06) |
0.99(0.93,1.04) |
0.99(0.94,1.05) |
1.04(1.00,1.09)* |
lag14 |
1.00(0.97,1.04) |
1.01(0.96,1.06) |
1.03(0.98,1.08) |
0.99(0.95,1.03) |
注:*为P<0.05 |
|
表 5 0~14d滞后SO2日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
|
2.2.5 NO2对全人群和不同年龄段人群流感日发病人数的滞后效应
日平均浓度每增加10μg/m3,对全人群、3岁~和15岁~人群最大RR值分别为1.06(1.00,1.13)、1.13(1.05,1.22)和1.12(1.04,1.21),分别在暴露后11d、12d和11d。其中0岁~人群在暴露当日作用最为明显,RR最大值为1.10(1.00,1.19)。详见表 6。
表 6
表 6 0~14d滞后NO2日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
1.03(0.98,1.08) |
1.10(1.01,1.19)* |
1.00(0.94,1.07) |
1.02(0.96,1.09) |
lag1 |
0.99(0.94,1.05) |
0.85(0.78,0.94) |
0.93(0.88,1.00) |
1.06(0.99,1.13) |
lag2 |
1.01(0.95,1.06) |
1.08(0.98,1.18) |
1.05(0.99,1.11) |
0.95(0.89,1.02) |
lag3 |
1.02(0.97,1.08) |
0.89(0.80,0.98) |
1.04(0.98,1.11) |
1.05(0.98,1.13) |
lag4 |
1.05(0.99,1.12) |
1.04(0.94,1.15) |
1.09(1.02,1.16)* |
1.05(0.98,1.13) |
lag5 |
0.99(0.93,1.05) |
1.07(0.97,1.19) |
0.99(0.93,1.06) |
0.97(0.90,1.05) |
lag6 |
1.00(0.94,1.06) |
0.94(0.85,1.04) |
1.04(0.97,1.12) |
0.99(0.92,1.06) |
lag7 |
0.96(0.91,1.02) |
0.97(0.87,1.07) |
0.97(0.90,1.04) |
0.95(0.88,1.02) |
lag8 |
0.96(0.90,1.02) |
0.95(0.86,1.05) |
0.96(0.90,1.03) |
0.97(0.90,1.05) |
lag9 |
1.01(0.95,1.08) |
0.98(0.89,1.07) |
1.04(0.96,1.12) |
1.00(0.93,1.07) |
lag10 |
0.97(0.91,1.03) |
0.98(0.88,1.09) |
0.97(0.89,1.05) |
0.96(0.90,1.04) |
lag11 |
1.06(1.00,1.13)* |
1.02(0.92,1.13) |
0.96(0.89,1.03) |
1.12(1.04,1.21)* |
lag12 |
1.03(0.97,1.10) |
0.93(0.84,1.03) |
1.13(1.05,1.22)* |
1.01(0.93,1.08) |
lag13 |
1.03(0.96,1.10) |
0.99(0.90,1.09) |
1.03(0.95,1.12) |
1.05(0.97,1.13) |
lag14 |
1.02(0.96,1.07) |
1.01(0.93,1.10) |
1.08(1.00,1.16)* |
0.99(0.92,1.06) |
注:*为P<0.05 |
|
表 6 0~14d滞后NO2日平均浓度每增加10μg/m3对流感发病的RR(95%CI)
|
2.2.6 O3对全人群和不同年龄段人群流感日发病人数的滞后效应
日最大8h平均浓度每增加10μg/m3,对全人群、0岁~和3岁~人群最大RR分别是1.02(0.99,1.05)、1.15(1.08,1.22)和1.07(1.02,1.12),分别出现在暴露后9d、1d和6d。对15岁~人群最大RR出现在暴露后2d和9d,为1.02(0.99,1.06)。详见表 7。
表 7
表 7 0~14d滞后O3日最大8h平均浓度每增加10μg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
0.98(0.95,1.01) |
0.92(0.87,0.96) |
0.95(0.91,1.00) |
1.00(0.97,1.04) |
lag1 |
1.01(0.98,1.04) |
1.15(1.08,1.22)* |
1.02(0.97,1.07) |
0.98(0.94,1.02) |
lag2 |
1.01(0.98,1.04) |
0.92(0.87,0.98) |
1.00(0.96,1.04) |
1.02(0.99,1.06) |
lag3 |
0.99(0.96,1.02) |
0.99(0.94,1.04) |
0.95(0.91,1.00) |
1.00(0.97,1.04) |
lag4 |
1.01(0.98,1.04) |
1.00(0.95,1.05) |
1.06(1.01,1.11)* |
1.00(0.96,1.03) |
lag5 |
0.98(0.95,1.01) |
0.97(0.93,1.02) |
0.94(0.89,0.98) |
0.99(0.96,1.03) |
lag6 |
1.01(0.98,1.05) |
1.05(1.00,1.10)* |
1.07(1.02,1.12)* |
1.00(0.97,1.04) |
lag7 |
0.98(0.95,1.01) |
0.98(0.92,1.03) |
0.90(0.86,0.94) |
1.01(0.97,1.04) |
lag8 |
0.99(0.96,1.03) |
1.08(1.02,1.14)* |
1.01(0.97,1.06) |
0.97(0.94,1.01) |
lag9 |
1.02(0.99,1.05) |
0.95(0.90,1.01) |
1.05(1.00,1.10)* |
1.02(0.99,1.06) |
lag10 |
0.97(0.94,1.00) |
0.95(0.90,1.00) |
0.88(0.84,0.93) |
0.99(0.96,1.03) |
lag11 |
1.01(0.98,1.04) |
0.99(0.94,1.04) |
1.03(0.98,1.08) |
1.01(0.97,1.05) |
lag12 |
0.97(0.94,1.00) |
0.95(0.90,0.99) |
0.94(0.89,0.99) |
0.98(0.94,1.01) |
lag13 |
1.01(0.97,1.04) |
1.05(1.00,1.10)* |
1.02(0.97,1.07) |
1.00(0.96,1.03) |
lag14 |
0.99(0.96,1.02) |
1.01(0.96,1.06) |
0.89(0.85,0.93) |
1.02(0.99,1.05) |
注:*为P<0.05 |
|
表 7 0~14d滞后O3日最大8h平均浓度每增加10μg/m3对流感发病的RR(95%CI)
|
2.2.7 CO对全人群和不同年龄段人群流感日发病人数的滞后效应
日平均浓度每增加0.1mg/m3,对全人群、0岁~、3岁~和15岁~人群最大RR分别是1.01(0.99,1.04)、1.04(1.01,1.07)、1.02(1.00,1.05)和1.02(1.00,1.05),分别出现在暴露后11d、5d、9d和11d。详见表 8。
表 8
表 8 0~14d滞后CO日平均浓度每增加0.1mg/m3对流感发病的RR(95%CI)
滞后/d |
全人群 |
0岁~ |
3岁~ |
15岁~ |
lag0 |
1.01(0.99,1.03) |
0.98(0.95,1.01) |
1.01(0.99,1.03) |
1.01(0.99,1.03) |
lag1 |
1.00(0.98,1.02) |
0.96(0.93,0.99) |
0.98(0.96,1.01) |
1.01(0.99,1.04) |
lag2 |
0.99(0.97,1.01) |
1.03(1.00,1.06)* |
1.01(0.99,1.03) |
0.97(0.94,0.99) |
lag3 |
1.01(0.99,1.03) |
0.97(0.94,1.00) |
1.02(0.99,1.04) |
1.02(1.00,1.05)* |
lag4 |
1.01(0.99,1.03) |
1.00(0.97,1.03) |
1.02(0.99,1.04) |
1.01(0.99,1.03) |
lag5 |
1.01(0.99,1.03) |
1.04(1.01,1.07)* |
1.00(0.98,1.03) |
1.00(0.98,1.03) |
lag6 |
0.99(0.97,1.01) |
1.00(0.97,1.04) |
0.99(0.96,1.02) |
0.99(0.96,1.01) |
lag7 |
0.99(0.97,1.01) |
0.95(0.92,0.99) |
0.99(0.96,1.01) |
0.99(0.96,1.02) |
lag8 |
0.98(0.96,1.01) |
0.97(0.94,1.00) |
0.97(0.95,1.00) |
1.00(0.97,1.02) |
lag9 |
1.01(0.99,1.03) |
1.01(0.98,1.04) |
1.02(1.00,1.05)* |
1.00(0.97,1.02) |
lag10 |
1.00(0.98,1.02) |
1.02(0.98,1.05) |
0.98(0.95,1.01) |
1.00(0.98,1.03) |
lag11 |
1.01(0.99,1.04) |
1.00(0.97,1.03) |
0.99(0.96,1.02) |
1.02(1.00,1.05)* |
lag12 |
1.00(0.98,1.03) |
1.00(0.97,1.03) |
1.02(0.99,1.05) |
1.00(0.97,1.03) |
lag13 |
1.01(0.99,1.03) |
1.01(0.98,1.04) |
1.02(0.99,1.05) |
1.01(0.98,1.03) |
lag14 |
1.00(0.98,1.02) |
1.01(0.98,1.03) |
1.01(0.98,1.03) |
1.00(0.98,1.02) |
注:*为P<0.05 |
|
表 8 0~14d滞后CO日平均浓度每增加0.1mg/m3对流感发病的RR(95%CI)
|
3 讨论
山东省是以燃煤火力发电为主,煤炭的大量燃烧会造成大气污染物的大量排放[9]。相关研究指出,山东省城市PM2.5、PM10浓度最高的季节是冬季[10],与流感发病高峰为冬春季大致吻合。推测气候和大气污染物对流感日发病存在相关性影响。原因可能为PM2.5、PM10为颗粒物,可以作为流感病毒传播的载体[11],同时PM2.5、PM10可以阻挡太阳光及紫外线,增加了病毒的传播机会及在外界环境中存活的时间。2013年世卫组织首次指出PM2.5是主要的环境致癌物,长期吸入会对人体呼吸系统造成严重影响,且浓度愈高死亡风险愈高。SO2、NOx为刺激性气体,会影响呼吸系统,诱发咳嗽、哮喘等呼吸道症状,造成肺部结构的永久性损伤[12],当进入呼吸道以后,能够对呼吸道黏膜产生侵袭性,损害上皮细胞及肺泡的微结构,增加病毒感染的概率和人体的易感性。CO能与氧气竞争性结合血红蛋白,可增强机体的缺氧状态,降低了体液免疫和细胞免疫功能及应激性。因此,颗粒物和污染物可以降低人群的免疫机能,增大流感发病的风险。
本文结果显示,大气污染物对于流感发病的影响存在滞后效应。且随日平均浓度升高10μg/m3(CO为升高0.1mg/m3),6种大气污染物对流感滞后作用各不相同。PM2.5、PM10、CO和O3对0岁~人群滞后作用最大,分别出现在暴露后14d、14d、5d和1d。提示空气污染物浓度高时应减少0岁~年龄段人群雾霾天气外出次数和时间,加强防护措施(如戴口罩),降低流感患病风险。在暴露后相应天数父母应该加强对孩子身体状况的关注。SO2和NO2对3岁~人群滞后作用最大,分别出现在暴露后2d和12d。提示政府有关部门应加强环境保护,推广使用天然气和光电等清洁能源,制定有关汽车尾气排放的制度规定。火力发电厂应采取措施减少含硫含氮气体排放。其他污染物对不同年龄的作用效果不仅与污染物的种类有关系,而且与不同年龄段对该污染物的抵抗能力有关系。
综上所述,流感的发病规律受大气污染物等因素的影响,提示相关医疗卫生机构应该针对本地区人口构成及主要污染源的情况,提前做好相应的预防控制措施。有关部门应积极引导农民使用清洁燃料,降低大气污染物的排放,加强流感的宣传力度,提高农民对于流感的认识,提倡绿色能源的普及,工厂应合理安排工人工作休息时间。
利益冲突:所有作者均申明不存在利益冲突。