2022, 39(4):362-366, 381.doi:10.11836/JEOM21578
噪声是工业生产中最常见的职业病危害因素之一,噪声导致的听力问题日趋严峻。据统计,我国工人职业噪声性听力损失的患病率约为21.3%[1]。近年来,职业性噪声聋成为仅次于职业性尘肺病的第二大职业病[2],且每年以超过10%的速度递增[3]。此外,研究发现,随着工业的发展,工业噪声类型发生巨大变化,非稳态噪声已取代稳态噪声成为工业中最常见的噪声类型,目前大多数岗位接触到的噪声属于非稳态噪声[4]。非稳态噪声(又称复杂噪声)是由在稳态噪声背景下叠加瞬态高能量脉冲性噪声组成[5],是除稳态噪声类型以外所有噪声类型的总和,包含典型的脉冲噪声。不同于稳态噪声,非稳态噪声具有复杂的时域结构,相对于稳态噪声的能量分布是正态的,而非稳态噪声的能量分布是非正态的,具有时变的特点[6]。
现有的噪声暴露测量方法主要适用于稳态噪声,仅用噪声的能量指标——等效A声级(equivalent continuous A-weighted sound pressure level, LAeq)来评估噪声暴露量。研究表明,应用稳态噪声的测量方法可能会低估非稳态噪声实际导致的听力损失[7]。峰度是能间接反映噪声时域结构的敏感指标[8]。国内外研究团队通过动物实验[9] 和人群调查[10-12] 发现,峰度可以作为噪声能量的辅助参量,用于非稳态噪声暴露的测量,峰度调整后的噪声能量指标比单一的能量指标能更有效反映非稳态噪声暴露对听力造成的损失。本文旨在全面地介绍基于峰度调整的职业性非稳态噪声暴露的测量方法,为非稳态噪声暴露的测量方法和相应标准制定提供思路。
现有噪声测量和评估标准中噪声的能量指标为强度,以8 h等效A声级(normalization of equivalent continuous A-weighted sound pressure level to a normal 8 h working day, LAeq,8 h)表示。LAeq,8 h不仅能反映环境噪声水平,还考虑到了噪声持续时间产生的影响。其测量时间比1 min LAeq更长,能显示噪声水平的变化,可真实反映劳动者实际接触噪声水平[13]。
根据“等能量学说”,噪声能量不仅与噪声强度有关,还与接噪时间密切相关,噪声强度越大,接噪时间越长,噪声能量就越大,造成听力损失的风险也越高[14]。累积噪声暴露量(cumulative noise exposure, CNE)是结合LAeq,8 h和接噪时间(工龄)的综合性指标,是能更全面反映噪声暴露的能量指标。其计算公式如下:
$$ {V_\rm{CNE}}=L_{\rm{Aeq,8\;h}} + {{10}} {\rm{lg}t} $$ | (1) |
(1)式中,VCNE是累积噪声暴露量,t为接噪时间(通常以年为单位)。式中常数10表明CNE符合“等能量学说”中的3 dB交换率,即接噪时间增加一倍(即2t),相当于LAeq,8 h增加3 dB。CNE与噪声性听力损失之间存在显著的剂量-反应关系,且与听力损失之间的相关性比LAeq,8 h更好[15]。
噪声的时域结构包括脉冲的峰值分布、持续时间分布、间隔时间分布等,若定量地分析这些变量来描述时域结构,在实际测量中可操作性差。于是,Erdreich[8]将这些必要变量归结为一个简单且易于计算的参数,提出用噪声峰度描述非稳态噪声的脉冲性,作为非稳态噪声时域结构的特异性指标,更具实用性。峰度本是描述随机变量概率分布相对于正态分布的统计量,计算公式如下:
$$ \beta=\dfrac{{\dfrac{1}{n}\sum\limits_{i=1}^n {{{\left( {x_{i} - \bar x} \right)}^4}} }}{{{{\left[ {\dfrac{1}{n}\sum\limits_{i=1}^n {{{\left( {x_{i} - \bar x} \right)}^2}} } \right]}^2}}} $$ | (2) |
(2)式中,β为峰度,xi是第i个值,
噪声能量与噪声强度和接噪时间均密切相关,所以峰度调整能量有两种调整方案,从CNE的公式也可看出,一种是调整接噪时间t,一种是调整噪声强度LAeq,8 h。
有学者曾提出峰度调整接噪时间的方案[6],根据CNE的公式,将公式中常数项10修改成
$${V_{{\text{CNE'}}}} = {L_{{\text{Aeq}},{\text{8}}\;{\text{h}}}} + \frac{{{(\text{ln}}\beta + 1.9)}}{{{\text{lg}}2}} \times {\text{lg}}t$$ | (3) |
(3)式中,β为所测噪声的峰度,VCNE'为峰度调整后的CNE,可用于衡量非稳态噪声的暴露量,评价非稳态噪声导致的听力损失。
Goley等[9]曾提出峰度调整噪声强度LAeq,8 h的方案,直接给LAeq,8 h增加含有峰度的调整项
$$ L{{'}}_{\rm{Aeq,8\;h}}=L_{\rm{Aeq,8\;h}} + \lambda {\text{lg(}}\beta /{{3)}}$$ | (4) |
(4)式中,β为所测噪声的峰度值,
当噪声类型为非稳态噪声时,峰度大于3,则
$$ V_{\text{PTS}}=b_{\text{0}} + {b_{{\text{1}}}{L_{\text {Aeq},8\;{\rm{h}}}}} + {{\epsilon }} $$ | (5) |
$$ V_{\text{PTS}}=b_{\text{0}} + {b_{\text{1}}{L_{\text {Aeq},8\;{\rm{h}}}}} + b_{\text{1}} {λ\text{lg}{(\beta /3)}} + {{\epsilon }} $$ | (6) |
$$ V_{\text{PTS}}=b_{\text{0}} + {b_{\text{A}}{L_{\text {Aeq},8\;{\rm{h}}}}} + b_{\text{B}} {\text{lg}}\left( {{{\beta /3}}} \right) + {{\epsilon }} $$ | (7) |
其中,ε为随机误差,且ε~N(0,σ2)。
稳态噪声测量仍依照现有的稳态噪声测量方法,本文介绍的峰度调整方法是针对工业中的非稳态噪声,所以适用对象为接触非稳态噪声的劳动者。稳态噪声的峰度值等于3,然而在实际的工业生产中,稳态噪声环境少,峰度为3~10的准稳态噪声环境多,其能量分布接近正态分布。稳态/准稳态噪声的典型工种有造纸行业的制浆工、纺织行业的纺纱工等,非稳态噪声的典型工种有家具厂的枪钉工、金属制造业的焊工、锻工、冲压工等[12]。
Qiu等[20] 在进行动物实验时发现,当LAeq达到很高水平时[例如大于95 dB(A)]时非稳态噪声与稳态噪声造成的听力损失无明显差异。Zhang等[7]也发现,随着噪声强度和工龄的增加,非稳态噪声与稳态噪声造成听力损失的差异也逐渐减弱,当噪声强度大于94 dB(A)而工龄超过10年时,峰度小于75的非稳态噪声导致的听力损失与稳态噪声所造成的趋向一致。此外,研究人员还发现,当噪声的LAeq小于70 dB(A)时,噪声峰度的大小与听力损失无关。这表明在噪声能量很高或很低的情况下,噪声能量对听力损失的作用远大于时域结构,峰度作为辅助参量参与评估非稳态噪声所致听力损失有一定的适用范围,只有在一定的强度范围内才能起到显著作用,即峰度作用阈值。目前关于峰度作用阈值的人群研究较少,要进一步证明基于峰度调整的非稳态噪声测量方法的有效性,还需要开展大规模人群调查研究峰度作用阈值。
稳态噪声只需测量LAeq,8 h,使用2型或以上、具有A计权、“S(慢)”档的声级计或个人噪声剂量计即可。非稳态噪声则不仅需要测量LAeq,8 h,还需要测量峰度。许多文献[7,10,12] 报道应用了特定的噪声剂量计记录声波,该仪器具有频率计权A、时间计权F,使用“预极化1/4英寸自由场型测试电容传声器”(1英寸=2.54 cm),录音精度为48 kHz/32 bit,灵敏度高(2.24 mV·Pa−1)。噪声强度测量范围为40~141 dB(A),频率范围为20 Hz~20 kHz。根据谢红卫[22] 和Zhang[12] 等学者的报道,他们采用能记录声波的特定噪声剂量计先录下劳动者整个班次接触到的噪声,之后由专业人员通过MATLAB软件编程分析噪声波形,从而计算得到LAeq,8 h和峰度值。这种方法不能直接读数,需要特定软件对声波进行分析,在实际测量时仍存在一些不便。已有研发团队将峰度测量和计算的软件程序内嵌在噪声剂量计中,在实测LAeq,8 h的基础上,自动计算峰度调整噪声能量指标,直接给出相应数值(
计算峰度的程序中,时间窗口设为40~60 s,窗口移动时不重叠。每个时间窗口可以得到一个峰度值,将这些峰度的平均值作为该劳动者接触噪声的最终峰度值[16]。Hamernik等[16]发现,采用40 s或更长的时间窗口时,得到的峰度达到一个稳定的值。Tian等[23]推荐60 s作为听力损失相关峰度计算的时间窗口长度,因为其兼具了较好的计算性能和听力损失评估准确度。在噪声采样频率为48 kHz的条件下,综合考虑了计算时效性、噪声动态特性的跟踪能力,经过不同窗口长度对不同噪声暴露的计算对比,最终选定计算窗口长度为40~60 s。
为正确选择测量仪器、测量方法和测量时间等,在测量前需进行现场调查。通过现场勘察和问询的方式,了解工作场所的面积、空间、工艺区划,调查噪声设备布局、工艺流程的划分、各生产程序的噪声特征、噪声变化规律,以及各岗位接触劳动者的数量、工作路线、工作方式、停留时间、工作班次等。现场调查需确保在正常生产状况下进行,并覆盖服务范围内全部工作场所,调查时间至少覆盖1个工作日。
在现场调查的同时,对接触噪声的劳动者开展工作日写实,应包括写实对象及其所在的岗位、工作日内从事的工作内容及活动的工作场所和接触噪声的时间及频率;每个作业岗位选择1~2名劳动者作为写实对象;对多条生产状况相同生产线的同类岗位,选择有代表性的1~2条生产线上的劳动者进行写实;对生产随意性大,每个工作日工作量和工作内容不稳定的岗位,应选择接触噪声强度最大、接触时间最长的劳动者进行工作日写实;对生产连续、稳定的作业岗位,或每个工作日生产状况相同的岗位,按1个工作日写实;对周期性生产作业的岗位,按生产周期写实。
在现场调查和工作日写实的基础上,制定噪声检测方案,绘制现场采样/测量点设置示意图。
稳态噪声需要测量劳动者接触噪声的LAeq,8 h。固定的工作岗位选用声级计,流动的工作岗位优先选用个体噪声剂量计,或对不同的工作地点使用声级计分别测量,并计算LAeq。非稳态噪声不仅要测量LAeq,8 h,还要测量噪声的峰度,需要用到能记录声波的特定噪声剂量计。测量前应将仪器充好电,并准备好传声器的防风罩,以备在风速较大处使用。正式测量前根据仪器校正要求完成测量仪器的校正。
正式测量之前,需对每个岗位噪声暴露情况进行预测量,判定噪声性质是属于稳态噪声[测量时间内声压级波动 < 3 dB(A)]还是非稳态噪声,预测量时间为一个完整工艺流程。稳态噪声与非稳态噪声采用不一样的测量方法。预测量可以在现场调查时完成。
通过预测量确定为稳态噪声的,测量LAeq,8 h时,使用2型或以上的声级计并设置为A计权、“S(慢)”档,取值为声级(LpA)或LAeq。将声级计放置在三角架上,或手持声级计,但应保持测试者与传声器的间距>0.5 m。传声器放置在劳动者耳部高度(站姿为1.50 m,坐姿为1.10 m),指向声源方向。每个岗位测量3次,取平均值。流动岗位使用个体噪声剂量计。
如“4仪器要求”中介绍,峰度可以通过录下噪声波形然后用MATLAB软件分析计算得到,也可以通过内置峰度测量和计算的噪声剂量计直接读出数值,因此非稳态噪声的测量也有相应的两种测量方法——“声波分析法”和“仪器直读法”。
(1)声波分析法。该法采用能记录声波的特定噪声剂量计,夹在劳动者的衣领处,固定在肩部,保持在劳动者耳部高度。录音时将传声器竖起,充分暴露于噪声,先录下劳动者整个班次接触到的噪声。录音时应在正常生产情况下进行。若风速超过3 m·s−1时,传声器应戴防风罩,应尽量避免电磁场的干扰。录音完成之后,由专业人员通过MATLAB软件编程分析噪声波形,从而计算得到LAeq,8 h和峰度值,并计算峰度调整噪声能量指标(CNE'或
(2)仪器直读法。该法采用内置峰度测量和计算的噪声剂量计,夹在劳动者的衣领处,固定在肩部,保持在劳动者耳部高度。测量时将传声器竖起,充分暴露于噪声。峰度是通过对每个时间窗口内的峰度取平均值而得出,峰度计的时间窗口为40~60 s,为尽可能测量出真实噪声水平,非稳态噪声应至少测量30 min,但需涵盖一个完整的作业周期。测量应在正常生产情况下进行。若风速超过3 m·s−1时,传声器应戴防风罩,应尽量避免电磁场的干扰。测量完成之后,通过仪器直接读出峰度值、实测LAeq,8 h值,或直接读出经峰度调整后的
两种方法没有本质差别,均可以要求记录或分析劳动者整个班次接触的噪声。为了方便实际操作,仪器直读法要求涵盖一个完整的作业周期,至少测量30 min以上。
测量记录应包括以下内容:测量日期、测量时间、气象条件(温度、相对湿度、风速)、仪器校准、测量地点(单位、厂矿名称、车间和具体测量岗位)、被测仪器设备型号和参数、测量仪器型号、噪声接触时长和测量时长、测量数据、峰度调整后数据、计算公式和过程、测量人员和陪同人员签名等。
现有的噪声测量标准是基于“等能量学说”建立的,“等能量学说”的原理是噪声对听力的损害仅与噪声能量成正比,即不管何种类型的噪声,只要其能量相等,造成的听力损害都应该相同[6]。而大量动物实验[18-21,24-25] 和流行病学调查[10-12] 证明了,在相等噪声能量的条件下,非稳态噪声的暴露比稳态噪声能产生更严重的听觉系统损害,这表明“等能量学说”并不适用于非稳态噪声。因此,现有的基于“等能量学说”的噪声暴露检测、评估标准也不适用于非稳态噪声。现有标准的问题在于并未考虑到非稳态噪声的时域结构对听力的影响,而单单用能量来描述噪声暴露,故并不能完全反映非稳态噪声导致的听力损失。非稳态噪声暴露的描述应该结合噪声的能量指标和时域结构指标,即本文中介绍的使用峰度调整噪声能量。
根据“等能量学说”,峰度调整可以有两种方案,一种是调整接噪时间t,另一种调整噪声强度(LAeq,8 h)。峰度调整接噪时间的方案已在人群研究中得到验证。Zhao等[11] 分别收集了非稳态噪声暴露和稳态噪声暴露工人的噪声和听力数据,发现调整前两组工人的剂量-反应关系曲线存在显著差异,进行峰度调整后两曲线几乎重合。同样,Zhang等[12] 研究发现在进行峰度调整前,非稳态噪声暴露者听力损失患病率比稳态噪声暴露者平均高7.63%,而在峰度调整CNE后仅仅只高1.12%。这表明峰度调整后使得不同类型的噪声在等能量条件下产生相等的效应,非稳态噪声性听力损失被低估的情况得到改善。Xie等[10] 还分别建立PTS关于CNE和CNE'的多元线性回归方程,发现峰度调整使得方程的决定系数升高,且CNE'对回归的相对贡献增大,这也说明了峰度调整CNE在人群中的有效性。未来仍需开展样本量更大的人群研究,用更丰富的方法验证这一调整方案。
另一种调整方案——峰度调整LAeq,8 h,相关的研究较少。Goley等[9] 曾开展动物实验验证这一调整效果,发现峰度调整后显著提升了LAeq与听力损失之间的相关性,大大减小了对非稳态噪声性听力损失的低估。调整后的指标
[1] |
ZHOU J, SHI Z, ZHOU L, et al. Occupational noise-induced hearing loss in China: a systematic review and meta-analysis[J]. BMJ Open, 2020, 10(9): e039576. DOI: 10.1136/bmjopen-2020-039576 |
[2] |
中华人民共和国国家卫生健康委员会. 2019年我国卫生健康事业发展统计公报[EB/OL]. [2021-05-27]. http://www.nhc.gov.cn/guihuaxxs/s10748/202006/ebfe31f24cc145b198dd730603ec4442.shtml. National Health Commission of the People’s Republic of China. Statistical bulletin on health development in China in 2019[EB/OL]. [2021-05-27]. http://www.nhc.gov.cn/guihuaxxs/s10748/202006/ebfe31f24cc145b198dd730603ec4442.shtml. |
[3] |
聂武, 胡伟江. 船舶工业企业噪声危害调查[J]. 中国工业医学杂志, 2016, 29(3): 167-170. NIE W, HU W J. Survey on noise hazard in shipbuilding industry[J]. Chin J Ind Med, 2016, 29(3): 167-170. |
[4] |
CHEN Y, ZHANG M, QIU W, et al. Prevalence and determinants of noise-induced hearing loss among workers in the automotive industry in China: A pilot study[J]. J Occup Health, 2019, 61(5): 387-397. DOI: 10.1002/1348-9585.12066 |
[5] |
SUTER A H. Occupational hearing loss from non-Gaussian noise[J]. Semin Hear, 2017, 38(3): 225-262. DOI: 10.1055/s-0037-1603726 |
[6] |
邱伟, 张美辨, 徐维超, 等. 峰度在评估复杂噪声所引起听力损失中的应用[J]. 中华耳科学杂志, 2016, 14(6): 701-707. DOI: 10.3969/j.issn.1672-2922.2016.06.002QIU W, ZHANG M B, XU W C, et al. The application of the kurtosis metric in evaluating hearing trauma from complex noise exposures[J]. Chin J Otol, 2016, 14(6): 701-707. DOI: 10.3969/j.issn.1672-2922.2016.06.002 |
[7] |
ZHANG M, XIE H, ZHOU J, et al. New metrics needed in the evaluation of hearing hazard associated with industrial noise exposure[J]. Ear Hear, 2021, 42(2): 290-300. DOI: 10.1097/AUD.0000000000000942 |
[8] |
ERDREICH J. A distribution based definition of impulse noise[J]. J Acoust Soc Am, 1986, 79(4): 990-998. DOI: 10.1121/1.393698 |
[9] |
GOLEY G S, SONG W J, KIM J H. Kurtosis corrected sound pressure level as a noise metric for risk assessment of occupational noises[J]. J Acoust Soc Am, 2011, 129(3): 1475-1481. DOI: 10.1121/1.3533691 |
[10] |
XIE H W, QIU W, HEYER N J, et al. The use of the kurtosis-adjusted cumulative noise exposure metric in evaluating the hearing loss risk for complex noise[J]. Ear Hear, 2016, 37(3): 312-323. DOI: 10.1097/AUD.0000000000000251 |
[11] |
ZHAO Y M, QIU W, ZENG L, et al. Application of the kurtosis statistic to the evaluation of the risk of hearing loss in workers exposed to high-level complex noise[J]. Ear Hear, 2010, 31(4): 527-532. DOI: 10.1097/AUD.0b013e3181d94e68 |
[12] |
ZHANG M, QIU W, XIE H, et al. Applying kurtosis as an indirect metric of noise temporal structure in the assessment of hearing loss associated with occupational complex noise exposure[J]. Ear Hear, 2021, 42(6): 1782-1796. DOI: 10.1097/AUD.0000000000001068 |
[13] |
陈辉, 蒋健敏, 张美辨, 等. 工业噪声测量指标研究进展[J]. 浙江预防医学, 2016, 28(4): 376-379. CHEN H, JIANG J M, ZHANG M B, et al. Research progress of industrial noise measurement index[J]. Zhejiang J Prev Med, 2016, 28(4): 376-379. |
[14] |
SHI Z, ZHOU J, HUANG Y, et al. Occupational hearing loss associated with non-Gaussian noise: a systematic review and meta-analysis[J]. Ear Hear, 2021, 42(6): 1472-1484. DOI: 10.1097/AUD.0000000000001060 |
[15] |
谢红卫, 张美辨, 全长健. 累计噪声暴露量与人听力损失的剂量反应关系研究[J]. 浙江预防医学, 2014, 26(4): 340-344,348. XIE H W, ZHANG M B, QUAN C J. A study on dose-response relationship between cumulative noise exposure and human hearing loss[J]. Zhejiang J Prev Med, 2014, 26(4): 340-344,348. |
[16] |
HAMERNIK R P, QIU W, DAVIS B. The effects of the amplitude distribution of equal energy exposures on noise-induced hearing loss: the kurtosis metric[J]. J Acoust Soc Am, 2003, 114(1): 386-395. DOI: 10.1121/1.1582446 |
[17] |
LEI S F, AHROON W A, HAMERNIK R P. The application of frequency and time domain kurtosis to the assessment of hazardous noise exposures[J]. J Acoust Soc Am, 1994, 96(3): 1435-1444. DOI: 10.1121/1.410287 |
[18] |
HAMERNIK R P, QIU W. Energy-independent factors influencing noise-induced hearing loss in the chinchilla model[J]. J Acoust Soc Am, 2001, 110(6): 3163-3168. DOI: 10.1121/1.1414707 |
[19] |
QIU W, DAVIS B, HAMERNIK R P. Hearing loss from interrupted, intermittent, and time varying Gaussian noise exposures: the applicability of the equal energy hypothesis[J]. J Acoust Soc Am, 2007, 121(3): 1613-1620. DOI: 10.1121/1.2434692 |
[20] |
QIU W, HAMERNIK R P, DAVIS B. The kurtosis metric as an adjunct to energy in the prediction of trauma from continuous, nonGaussian noise exposures[J]. J Acoust Soc Am, 2006, 120(6): 3901-3906. DOI: 10.1121/1.2372455 |
[21] |
QIU W, HAMERNIK R P, DAVIS R I. The value of a kurtosis metric in estimating the hazard to hearing of complex industrial noise exposures[J]. J Acoust Soc Am, 2013, 133(5): 2856-2866. DOI: 10.1121/1.4799813 |
[22] |
谢红卫, 唐仕川, 周莉芳, 等. 非稳态噪声累积暴露量与听力损失的关系[J]. 环境与职业医学, 2015, 32(1): 56-60. XIE H W, TANG S C, ZHOU L F, et al. Relationship between cumulative noise exposure and human hearing loss caused by non-stationary noise[J]. J Environ Occup Med, 2015, 32(1): 56-60. |
[23] |
TIAN Y, DING W, ZHANG M, et al. Analysis of correlation between window duration for kurtosis computation and accuracy of noise-induced hearing loss prediction[J]. J Acoust Soc Am, 2021, 149(4): 2367. DOI: 10.1121/10.0003954 |
[24] |
DAVIS R I, QIU W, HAMERNIK R P. Role of the kurtosis statistic in evaluating complex noise exposures for the protection of hearing[J]. Ear Hear, 2009, 30(5): 628-634. DOI: 10.1097/AUD.0b013e3181b527a8 |
[25] |
DAVIS R I, QIU W, HEYER N J, et al. The use of the kurtosis metric in the evaluation of occupational hearing loss in workers in China: implications for hearing risk assessment[J]. Noise Health, 2012, 14(61): 330-342. DOI: 10.4103/1463-1741.104903 |
[基金项目] 浙江省重点研发项目(2015C03039);2018年度浙江省151人才工程培养项目(无编号);2016年度浙江省卫生创新人才培养项目(无编号);职业健康标准前期研究项目(20210102);浙江省医药卫生基金(2019KY057);浙江省医药卫生基金(2021KY120)
[作者简介]
[收稿日期] 2021-12-06
引用格式
胡勇,
施志豪,
高向景, 等.
基于峰度调整的职业性非稳态噪声暴露测量方法[J].环境与职业医学,
2022, 39(4): 362-366, 381.
doi:10.11836/JEOM21578.
HU Yong , SHI Zhihao , GAO Xiangjing , XIN Jiarui , ZHOU Lifang , ZHANG Meibian . Measuring method of occupational non-Gaussian noise exposure based on kurtosis adjustment.Journal of Environmental & Occupational Medicine, 2022, 39(4): 362-366, 381. doi:10.11836/JEOM21578.