摘要
目的: 以贝尔面瘫患者为研究对象，回顾性分析贝尔面瘫发病与发病时运气因素的相关性，旨在探讨运气理论的科学性，为临床运用运气理论防治贝尔面瘫提供依据;同时，为进一步明确深圳市的气候环境对贝尔面瘫发病的影响，充分体现中医学运气理论中“天人相应”的学术思想，本研究选择分布滞后非线性模型来分析气温、相对湿度、风速、水汽压、大气压、能见度、总云量、日照时间、日蒸发量、太阳辐射等气象因素对贝尔面瘫发病的影响，以获得更具价值的研究结果，为降低气象因素对人体健康的不利影响提供科学依据。 方法: 本研究为回顾性研究，研究第一部分根据病例筛选标准，收集2010年1月20日-2020年1月19日全部符合上述标准的住院病例，共纳入患者841例，采集所需病例信息并推算出患者准确的发病日期。根据天干地支推演五运六气，同时收集各运气时段的住院患者总人数，计算出同时段非贝尔面瘫患者的人数。对纳入的各运气时段住院患者中贝尔面瘫患者人数及非贝尔面瘫患者人数，采用卡方检验进行统计分析。同时，本研究从节气及季节角度对贝尔面瘫患者住院人数的占比，采用卡方检验进行统计分析。 研究第二部分根据病例筛选标准，收集2009年1月1日-2020年2月29日全部符合上述标准的住院病例，共纳入患者863例，采集所需病例信息并推算出患者准确的发病日期。气象资料根据深圳市国家气候观象台双龙自动气象观测站记录，由深圳市气象局提供，包括2009年1月1日-2020年2月29日共11年逐日的平均气温(℃)、相对湿度(％)、风速(0.1m/s)、水汽压(hpa)、大气压(hpa)、能见度(km)、总云量(％)、日照时间(h)、日蒸发量(mm)和太阳辐射(W/m2)等气象资料。应用R软件，选择分布滞后非线性模型，在深圳市的气象环境条件下分析上述气象因素对贝尔面瘫发病的影响效应。 结果: 1.基于运气学说的贝尔面瘫发病的探索性研究 1.1基于岁运太过不及理论对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，最高的是辛年，为1.01％，为水运不及之年;最低的是庚年，为0.44％，为金运太过之年。且差异有统计学意义(x2=33.864，P=0.000＜0.05)。 1.2基于主气六气理论对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，主气为初之气-厥阴风木时最高，为0.91％;主气为二之气-少阴君火时最低，为0.63％。但差异无统计学意义(x2=10.806，P=0.055＞0.05)。 1.3基于司天六气理论对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，太阴湿土司天时最高，为0.93％;少阴君火司天时最低，为0.42％。且差异有统计学意义(x2=13.316，P=0.021＜0.05)。 1.4基于在泉六气理论对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，少阴君火在泉时最高，为0.97％;少阳相火在泉时最低，为0.49％。且差异有统计学意义(x2=25.319，P=0.000＜0.05)。 2.基于节气、季节因素的贝尔面瘫发病的探索性研究 2.1基于二十四节气对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，节气为立春时最高，为1.15％;节气为小雪时最低，为0.44％。但差异无统计学意义(x2=34.811，P=0.054＞0.05)。 2.2基于季节(按节气划分)对贝尔面瘫患者的住院人数占比进行分析，结果发现，贝尔面瘫患者住院人数的占比，春季时最高，为0.77％;冬季时最低，为0.69％;但差异无统计学意义(x2=2.035，P=0.565＞0.05)。根据二十四节气，将四季划分为十二个季节进一步分析，结果发现，贝尔面瘫患者住院人数的占比，孟春时最高，为1.11％;孟冬时最低，为0.50％;且差异有统计学意义(x2=27.319，P=0.004＜0.05)。 3.在深圳市的气象环境条件下，基于分布滞后非线性模型分析气象因素对贝尔面瘫发病的影响 3.1在气温的总体效应分析中，随着气温的升高RR值呈逐渐减小的趋势。当日平均气温＜25℃时，会增加贝尔面瘫的发病风险;当日平均气温＞25℃时，会降低贝尔面瘫的发病风险;当日平均气温＞30℃时其与疾病的发生无统计学关联，提示气温对贝尔面瘫发病的影响是非线性的，不同的温度时贝尔面瘫的发病风险不同。不同日平均气温与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均气温在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.2在相对湿度的总体效应分析中，随着相对湿度的增加RR值呈逐渐减小的趋势。当日平均相对湿度＜76％时，会增加贝尔面瘫的发病风险;当日平均相对湿度＞78％时，会降低贝尔面瘫的发病风险，并且与疾病的发生均有统计学关联，提示相对湿度对贝尔面瘫发病的影响是非线性的，不同的相对湿度时贝尔面瘫的发病风险不同。不同日平均相对湿度与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均相对湿度在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.3在风速的总体效应分析中，随着日平均风速的增加，贝尔面瘫的发病风险也逐渐增加，从日平均风速大于2.6m/s开始，RR值增加较明显，提示风速对贝尔面瘫发病的影响是非线性的，随着风速的增加，贝尔面瘫的发病风险呈逐渐增加的趋势。不同日平均风速与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均风速在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.4在水汽压的总体效应分析中，随着水汽压的升高RR值呈逐渐增大的趋势。当日平均水汽压＜229hpa时，会降低贝尔面瘫的发病风险;当日平均水汽压＞231hpa时，会增加贝尔面瘫的发病风险，并且与疾病的发生均有统计学关联，提示水汽压对贝尔面瘫发病的影响是非线性的，不同的水汽压时贝尔面瘫的发病风险不同。不同日平均水汽压与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均水汽压在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.5在大气压的总体效应分析中，随着大气压的升高，RR值呈先逐渐增大后逐渐减小的趋势。当日平均大气压＜1005hpa时，会增加贝尔面瘫的发病风险;当日平均大气压＞1006hpa时，会降低贝尔面瘫的发病风险，其中日平均大气压在1010-1021hpa区间时，与疾病的发病有统计学关联，提示大气压对贝尔面瘫发病的影响是非线性的，不同的大气压时贝尔面瘫的发病风险不同。不同日平均大气压与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均大气压在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.6在能见度的总体效应分析中，随着能见度的增加，RR值呈先逐渐减小后逐渐增大的趋势。当日平均能见度＜17km或＞26km时，会增加贝尔面瘫的发病风险;当日平均能见度在18-25km之间时，会降低贝尔面瘫的发病风险，提示能见度对贝尔面瘫发病的影响是非线性的，不同的能见度时贝尔面瘫的发病风险不同。不同日平均能见度与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均能见度在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.7在总云量的总体效应分析中，随着总云量的增加，RR值呈先缓慢减小后逐渐增大的趋势。当日平均总云量＜66％或＞71％时，会增加贝尔面瘫的发病风险;当日平均总云量在67-69％之间时，会降低贝尔面瘫的发病风险，提示总云量对贝尔面瘫发病的影响是非线性的，不同的总云量时贝尔面瘫的发病风险不同。不同日平均总云量与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均总云量在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.8在日照时间的总体效应分析中，随着日照时间的延长，RR值呈先逐渐减小后逐渐增大再逐渐减小的趋势。当日照时间＜5h或＞10h时，会降低贝尔面瘫的发病风险;当日照时间在6-9h之间时，会增加贝尔面瘫的发病风险，提示日照时间对贝尔面瘫发病的影响是非线性的，不同的日照时间时贝尔面瘫的发病风险不同。不同日照时间与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日照时间在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.9在日蒸发量的总体效应分析中，随着日蒸发量的增加，RR值呈先逐渐增大后逐渐减小再逐渐增大的趋势。当日蒸发量＜1mm或在4-7mm区间时，会降低贝尔面瘫的发病风险;当日蒸发量在2-3mm区间或＞8mm时，会增加贝尔面瘫的发病风险，提示日蒸发量对贝尔面瘫发病的影响是非线性的，不同的日蒸发量时贝尔面瘫的发病风险不同。不同日蒸发量与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日蒸发量在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 3.10在太阳辐射的总体效应分析中，随着太阳辐射的增加，RR值呈缓慢减小的趋势。当日平均太阳辐射＜157W/m2时，会增加贝尔面瘫的发病风险;当日平均太阳辐射＞159W/m2时，会降低贝尔面瘫的发病风险，其中日平均太阳辐射在248-261W/m2区间时，与疾病的发病有统计学关联，提示太阳辐射对贝尔面瘫发病的影响是非线性的，不同的太阳辐射时贝尔面瘫的发病风险不同。不同日平均太阳辐射与贝尔面瘫发病关联作用的3D图和热力图均提示不同的日平均太阳辐射在不同的滞后天数下对贝尔面瘫发病的影响效应不同。 结论: 1.贝尔面瘫患者住院人数的占比，从岁运太过不及角度看，最高的是辛年，为水运不及之年;从六主气角度看，主气为初之气-厥阴风木时最高;从客气司天六气角度看，最高的是太阴湿土司天之年;从客气在泉六气角度看，少阴君火在泉时最高;从季节角度看，孟春(立春、雨水二节气)时最高。 2.在深圳市的气象环境条件下，各气象因素对贝尔面瘫发病的影响均是非线性的。除风速外，各气象因素对贝尔面瘫发病风险的影响，随着日均值的逐渐增加，呈先增加后降低或先降低后增加的趋势。而在风速的总体效应分析中，随着日平均风速的逐渐增加，贝尔面瘫的发病风险呈逐渐增加的趋势。此外，各个气象因素在不同的滞后天数下对贝尔面瘫发病的影响效应均是不同的。
NETL