基于EEMD-MFDFA的镇江港PM_%282.5%29演化特征分析.pdf

江苏大學学报（自然科学版）JOURNAL OF JIANGSU UNIVERSITY(Natural Science Edition)D0I:10.3969/j.issn.1671 7775.2023.06.0172023年11月第44卷第6 期Nov.2023Vol.44No.6开放科学（资源服务）标识码（OSID）：基于EEMD-MFDFA的镇江港PM2.5演化特征分析郑敏学，毛静，孙智灏，朱（1.江苏大学安全与环境工程学院，江苏镇江2 12 0 13；2.镇江港务集团有限公司，江苏镇江2 12 0 0 1）摘要：为定量分析镇江港PM2.5浓度时间序列的内在演化特征，以镇江港2 0 19 2 0 2 0 年的2 0 56 6 7组PM2.5时均浓度数据为基础，采用基于集合经验模态分解（EEMD）的多重分形去趋势波动分析（M FD FA）算法进行研究.结果表明：各监测点PM2.5日均浓度演化均具有显著的长期持续性特征，现阶段PM2.5浓度升高可能会导致未来一段时间内PM2.5浓度的持续升高，不同监测点PM2.5浓度的多重分形谱参数存在差异性，多重分形特征强的区域应侧重即时性防治，反之则应注重长期性监管；运用自组织临界性理论，发现15个监测点PM2.5浓度波动呈现幂律分布特征，符合自组织临界性特征，港口系统存在发生大气高污染的风险。关键词：PM2.5；镇江港；集合经验模态分解；多重分形去趋势波动分析；自组织临界理论中图分类号：X513引文格式：郑敏学，毛静，孙智灏，等：基于EEMD-MFDFA的镇江港PM2.s演化特征分析 J.江苏大学学报（自然科学版），2023,44(6):738-744.Evolution characteristics of PM2.5 in ZhenjiangPort based on EEMD-MFDFA龙，李志坚，居居珍?，黄爱国文献标志码：A文章编号：16 7 1-7 7 7 5(2 0 2 3)0 6-0 7 3 8-0 7ZHENG Minxue,MAO Jing,SUN Zhihao,ZHU Long,LI Zhjian?,JU Zhen?,HUANG Aiguo?(1.School of Safety and Environmental Engineering,Jiangsu University,Zhenjiang,Jiangsu 212013,China;2.Zhenjiang Port GroupCo.,Ltd.,Zhenjiang,Jiangsu 212001,China)Abstract:To quantitatively analyze the internal evolution characteristics of PM2.s time series in ZhenjiangPort,taking 205 667 PM2.5 hour concentration series in Zhenjiang port from 2019 to 2020 as basic data,the multifractal detrended fluctuation analysis(MFDFA)method based on ensemble empirical modedecomposition(EEMD)was used.The results show that the evolution of PM2.concentration at eachmonitoring point has significant long-term persistence characteristics.The rise of PM2.5 concentration atpresent stage may lead to the rise of PM2.5 concentration in the future.The multifractal spectrumparameters of PM2.s concentration at different monitoring points are different.Areas with strongmultifractal characteristics should be focused on immediate prevention and control,on the contrary,long-term supervision should be paid attention to.According to the theory of self-organized criticality,thefluctuation of PM2.s concentration at 15 monitoring points presents the characteristics of power-lawdistribution,which is in line with the self-organized criticality.收稿日期：2 0 2 2-0 4-0 5作者简介：郑敏学（19 6 7 一），男，江苏南京人，副教授（），主要从事环境监测、风险评价及安全管理等研究.毛静（19 9 6 一），女，浙江衢州人，硕士研究生（2 2 2 19 0 9 0 0 3 ），主要从事环境监测及安全管理研究。739第6 期郑敏学等：基于EEMD-MFDFA的镇江港PM2.5演化特征分析Key words:PM2.5;Zhenjiang Port;ensemble empirical mode decomposition;multifractal detrendedfluctuation analysis;theory of self-organized criticality镇江港地处江苏省镇江市西部，是长江三角洲地区重要的对外开放港口，港口地理坐标为北纬321110.25321127.25,东经119 3 7 56.53 1194016.20,2020年上半年货物吞吐量达17 3 0 5万t,身全国前十 1.港口在散货作业、场地堆存及车辆运输中所产生的PM2.5是构成区域大气污染的主要因素.随着港口铁矿石、煤炭等易起尘散货吞吐量的快速增长,PM2.5污染的控制难度进一步加大,极大影响职工健康、港区工作环境及周边生态环境的安全 2 1.中共中央、国务院2 0 2 1年印发的国家综合立体交通网规划纲要提出，要建立绿色化水平居于世界前列的国家综合立体交通网，港口作为交通网的重要组成部分，必须全面提升PM2.5治理水平，构建现代化绿色港口.镇江港共有15个监测点，分布密度高，覆盖各个生产场所，可以充分反映港口PM2.s的分布特点.研究镇江港的PM2.5分布规律及其演化特征，分析港口污染形成的机制，对提高PM2.5治理水平具有重要参考意义.PM2.s的演化过程既受人为、气象等因素的影响,同时也受到大气化学反应等机制作用的影响，因此PM2.5演化具有非线性、非平稳性等复杂特征.目前传统数理方法无法消除非平稳时间序列内在的自相关性和识别外在趋势所导致的伪相关现象 3 .多重分形去趋势波动分析（MFDFA)方法可直接从非线性复杂系统入手，整体上定量描述系统演化过程，获取其内在的变化规律 4.然而MFDFA算法的去趋势部分仍然存在不足，当序列太长或其趋势不符合多项式类型时，分析结果会不准确.为了克服该方法的不足，文献5将经验模态分解（EM D)算法作为去趋势模块嵌人到MFDFA算法当中,EMD在时间序列数据处理方面具有明显优势，但也存在“模态混合”的问题，即在不同时间尺度的信号可能被归类成同一个固有模态函数（IMF），或相同时间尺度的信号可能被归类成不同IME.为了克服上述缺点，文献 6 提出了集合经验模态分解(EEMD)算法,该算法可以更准确表示变量趋势.目前基于EEMD的MFDFA算法已用于数字高程模型和降水动力学等方面的研究 7-8 .为了探讨镇江港PM2.5浓度的演化特征，首次采用EEMD-MFDFA算法，分析港口PM2.5时间序列的多重分形特性,并基于自组织临界(SOC)理论，阐明港口不同监测点粉尘排放对PM2.5的非线性动力学影响，将为港口PM2.5治理提供科学依据1数据与方法1.1石研究数据来源研究中所运用的PM2.5时均浓度数据均采自20192020年镇江港的15个空气质量自动监测点.监测原理为光散射法，共收集数据为2 0 6 2 45组,经审查核验，其中有效数据为2 0 56 6 7 组.20192020年镇江港PM2.s时均浓度p（PM 2.5）监测点及空间分布如图1所示.人12141110151193730E119380E11938”30E119390E1193930E119400E1194030E图1220192020年PM2.5时均浓度监测点及空间分布1.2基于EEMD的MFDFA算法1.2.1MFDFA算法根据文献 9 ,MFDFA算法包括5个步骤1）运用具有时间序列特征的数据（k=1，2，,N)构建新的序列Y（i）,i为时间序列特征数据中的第i个数据，即Y(i)=Z(x-x),i=1,2,N,(1)=式中：为平均值.2）将序列Y(i)分成N,个独立子序列,s为时间尺度.为了避免尾部数据的丢失，对其进行逆向处理,得到2 N,个独立子序列.3）对每个子序列通过最小二乘法得到局部协方差函数F(u,s),即当u=1,2,N,时,F(,s)=Y(-1)s+i)-y.(i)3;(2)当u=N,+1,N,+2,2N,时,F(n,s)=Y(N-(-N.)+i)-y,(i),(3)S式中：y（i)为第u段数据上的第i点对应的趋势62.92375.m40.97740江苏大学学报（自然科学版）第44卷函数.4）所有2 N,个独立子序列的q阶波动函数通过以下等式获得：2Ns元2(F(w.).0,2NsexpZ In(F(u,s),=0.1.=5）广义Hurst 指数h(q)可通过计算函数F（s)与时间尺度s的双对数曲线斜率得到.通过Legendre变换,可得到奇异指数、多重分形谱函数f()和广义Hurst指数h(q).三者关系如下：=h(q)+qh(q),f()=q(-h(qg)+1.对于单分形时间序列,h（q)独立于;对于多重分形时间序列,h(q)随q的变化而变化.h(2）0.5,表明时间序列存在长期记忆性，且具有正效应，数值越大,表示持续性越强.=max-mmn为多重分形谱宽度，该值越大，多重分形特征越强，数据波动越剧烈0 1.AF(）=f（mn）-f（mx）为最小概率子集和最大概率子集的分形维数之差.f(）0 表示PM2.5浓度位于较低值的概率大于位于较高值的概率，1.2.2EEMD 算法EEMD算法本质是一种使用正态分布的白噪声序列辅助数据分析方法-12 根据文献13 ，EEMD算法包括3 个步骤.1）将正态分布的白噪声序列nm（t)添加到原始系列x(t)中,第次加人白噪声序列后函数为xm(t)=x(t)+nm(t).(7)2）使用EMD方法将加人白噪声序列的信号分解为K个固有模态函数 IMF分量 C篇(t)和残差r(t).3）在重复上述步骤过程中加入不同的白噪声序列.加入白噪声信号的次数为10 0 次，白噪声序列等级为0.2 级，可获得较好的分解效果 14-15,最终获得分解的相应IMF分量和残差.将分解得的IMF集合取平均值，用以抵消白噪声序列,得到EEMD分解后的IMF分量和残差如下：(8)1KTm(t).式中：K为加入高斯白噪声的总体数量；为添加白噪声的次数；j为对原序列进行EEMD分解得到的第j个分量;cm(t)和rm（t)分别为第m次分解得到的第j个IMF分量及其对应的剩余残差分量；c(t)为进行第K次分解得到的第j个IMF分量；r(t)为最后得到的残差。1.2.3基于EEMD的MFDFA算法为了优化传统MFDFA算法中的去趋势问题，将EEMD算法嵌入MFDFA算法的第3 步，用残差(4)替代原始序列,保持其他步骤不变 16 .在EEMD-MFDFA算法中,第3 步表述如下：当U=1,2,N,时,(w,s)=(-1)s+i)r();(9)当u=N,+1,N,+2,2N,时,(5)F(v,s)=IZY(N-(u-N.)s+i)-r(t).(6)(10)1.3自组织临界理论自组织临界理论可以阐释复杂系统中宏观方面出现幂律分布特征的动力机制，临界状态下的大气系统相对稳定，微小的干扰可能会导致连锁反应，使系统发生“崩塌”.系统达到自组织临界状态的最显著表征之一是幂律分布特征.若PM2.5时均浓度时间序列波动具有自临界组织特性，需满足以下关系式 1lgM(p A p o)-入lg Ap,式中：p为某时刻n的PM2.s时均浓度;p为PM2s时均浓度的波动值，Ap=pn+1-pn;入为标度指数;M为PM2.5时均浓度波动值大于某一个PM2.5时均浓度波动值po的次数.2结果与分析2.1空间分布特征图1中，临江区域总体PM2.s时均浓度（58.57g/m）高于非临江区域PM2.5时均浓度（52.6 7 g/m）.港口原始PM2.s时均浓度变化范围为10.0 223.9g/m.依据GB30952012环境空气质量标准，研究期间有18%的日均浓度超过国家二级1标准，空气质量等级为轻度污染及以上.14号监测K点浓度较高，日均浓度为6 2.9 2 g/m,10号监测点浓度偏低，日均浓度为40.9 7 g/m，最大浓度差为21.95g/m.14号监测点位于集装箱作业区域与码头生产区域运输线交汇点，在集装箱运输、码头货物作业时车辆来往频繁，产生了包含作业扬尘、道路扬尘及车辆尾气等多种尘源，从而造成PM2.5浓度持续偏高.10 号监测点靠近办公区，远离生产作(11)741第6 期郑敏学等：基于EEMD-MFDFA的镇江港PM2.5演化特征分析业区域，因而PM2.5浓度偏低.802.2基于EEMD-MFDFA算法的分析0对15个监测点在2 0 2 0 年的PM2.5日均浓度数-80据分别进行了EEMD分解，发现12 个监测点数据800序列长度均为3 6 6 个,其余3 个监测点数据序列长-80度不足3 6 6 个.为了保持PM2.5日均浓度数据连续30性,研究中仅对序列长度为3 6 6 个的12 个监测点PM2.5日均浓度数据进行分析.12 个监测点分布均匀,涵盖了整个港口，对其PM2.5日均浓度数据进行多重分形分析可以反映港口PM2.5日均浓度演化特征.1号监测点PM2.5日均浓度处于12 个监测点的中位（排序为第7).以1号监测点PM2.5日均浓度为例,计算得到的IMF分量和残差r(t)结果如图2 所示，其中加入白噪声次数为10 0 次，白噪声序列等级为0.2 级.通过图2 中的残差曲线，可以清晰地了解原始数据的总体变化趋势图3 为通过对2 0 2 0 年镇江港12 个监测点PM2.5日均浓度数据进行分解得到的残差曲线.由图3可知：港口总体PM2.5浓度变化呈U形分布；1月和12 月PM2.5浓度较高，6 9 月PM2.5浓度偏低，7月PM2.5浓度最低.120r6001206001206001206001206001206002020-01-01 2020-02-29 2020-04-29 2020-06-28 2020-08-27 2020-10-26 2020-12-25温度是影响PM2.5浓度变化的重要因素之一.通过Pearson相关分析方法，计算并统计各监测点PM2.5时均浓度与温度的相关系数，发现PM2.s时均浓度与温度呈负相关的关系，相关系数为-0.48 2-0.113.这是因为高温时有较强的气流活动，有利-30200-2020SHNI0-20209HNI0-20100-1090060302020-01-01 2020-02-29 2020-04-29 2020-06-28 2020-08-27 2020-10-26 2020-12-25图2 1号监测点PM2.5日均浓度的EEMD分解结果1401号700602号40201203号6001404号7005号1407008号1407002020-01-01 2020-02-29 2020-04-29 2020-06-28 2020-08-27 2020-10-26 2020-12-25日期(a)1-5,8号图3 12 个监测点PM2.5日均浓度数据的残差曲线日期日期(b)9-14号于污染物扩散.同时，镇江6 一9 月气温高、降水多，而1月和12 月气温低、降水少,温度、降雨等是导致港口PM2.5浓度序列的残差曲线呈U形变化的重要因素.运用EEMD-MFDFA算法，对港口12 个监测点PM2.5浓度进行多重分形分析，得到图4所示的9号10号11号12号13号14号742江苏大学学报（自然科学版）第44卷一1号Hurst 指数图.2.672.42.22.01.81.61.4L-6由图4可知,各监测点的h(q)随着阶数q值的增大而减小,这表明PM2.5浓度时间序列有多重分形特征.各监测点的h(2)均显著大于0.5,表明PM2.5日均浓度波动不服从经典马尔可夫函数的随机行为,具有显著的长期持续性.即过去排放的PM2.5浓度会对未来的PM2.5浓度造成影响.2、9 和14号监测点的h(2)相对较大,分别为2.3 1、2.3 0 和2.3 1,说明这3 个监测点所在区域的PM2.5浓度长期持续性较高.镇江港位于长江边，三面建有围墙与防风网,风速多为3 级及以下,3、4和5级风的PM2.5时均浓度数据分别为总监测数据的9 9.7 6%0.2 0%和0.0 4%，不利于粉尘向外部扩散，港口区域处于半封闭状态，半封闭状态会加剧PM2.5浓度长期持续性不变.PM2.5日均浓度的多重分形谱f(）和奇异指数的关系曲线如图5所示.由图5可知：所有监测点的多重分形谱f(）曲线都不对称,都呈现右偏抛物线形状；f(）0 时，表明镇江港PM2.5浓度在较低水平上波动的概率占主导地位.曲线具有明显的右截断和左长尾，归因于时间序列对大幅度局部浓度波动的敏感性.多重分形谱分析表明，大波动的影响在PM2.5日均浓度时间序列中起主导作用.1.51.00.5F0F-0.5F-1.00.91.21.5图5PM2.5日均浓度多重分形谱和奇异指数的关系曲线由图5还可知：变化范围为0.59 1.11;3、10和13 号监测点相对较大，分别为0.9 6、1.11和0.8 5,说明这3 个监测点所在区域PM2.5多重分形谱演化出较高波动性;2、12 和14号监测点相对较小，分别为0.6 1、0.59 和0.6 5,反映这3 个监测点所在区域PM2.s多重分形谱演化波动性较小.对比不同监测点的多重分形计算结果可知:323二89-4202图 4Hurst 指数图58号1.82.12.42.73.0号和13 号监测点位于堆场附近，装卸作业和转运作业多,PM2.s浓度多重分形特征较强，受作业时段影响，波动性较剧烈；10 号监测点位于门口与办公区域，远离作业场所，主要排放源为汽车扬尘和尾气，多重分形特征较强,对于这些区域的PM2.5污染应侧46号号号号号号二重即时性防治;2、9 和14号监测点靠近作业区及主要输运道路，生产作业及运输过程中的作业扬尘、道路扬尘和柴油车尾气是增大PM2.5浓度的重要因素，特别是频繁来往的集装箱作业车辆和散货转运车辆，造成该区域PM2.5浓度有较强的持续性，多重分形特征较弱，波动性较稳定，对于这些区域的PM2.5污染应侧重于长期性、全时段的防治.2.3自组织临界性为了进一步认识港口PM2.5时均浓度演化的动力机制,计算各监测点PM2.s浓度波动的累积频度统计分布,利用最小二乘法对PM2.5浓度进行线性回归,得到各监测点标度指数.图6 为各监测点PM2.5时均浓度波动值的累积分布.由图6 可知，各个监测点PM2.5时均浓度的标度指数（即图6 中曲线斜率）范围为2.193.05,拟合优度R均达到0.9 9,证明拟合程度较好,15个监测点的PM2.5浓度均符合幂律分布关系，变化趋势总体相同,可判断港口PM2.5系统为自组织临界系统.由于各个监测点所处区域的主要尘源和环境不同，导致不同区域PM2.5的幂律分布规律略有不同.3.22.82.4F2.0F1.61.22.82.4201.61.23.63.22.82.42.0F1.61.0图6 各监测点PM2.s时均浓度波动值的累积分布由于港口各监测点PM2.5日均浓度均呈现长期持续性发展,未来的PM2.5浓度受过去排放的影响，一6 号一7 号8号9号10号一11号一12 号一13 号14号-15号1.11.21gAp1.31.41.5743第6 期郑敏学等：基于EEMD-MFDFA的镇江港PM2.5演化特征分析港口生产作业和货物运输产生的粉尘浓度会进一步累积,易在局部空间内形成较严重的污染事件.若浓度持续增加,系统达到临界状态,少量的 PM2.5排放就会形成严重污染事件.港口各区域PM2.5浓度呈现幂律分布结构,是判断具有 SOC 特征的重要依据.因此自组织临界性是港口PM2.5演化的动力学机制之一，会导致港口未来发生高浓度污染事件的发生.3、10 和13 号监测点PM2.5浓度波动性大，即时性的排放较频繁.而2、9 和14号监测点PM2.5浓度呈长期持续性增强状态,预示未来PM2.5浓度将维持在较高浓度.大波动性与强持续性都会打破港口PM2.s的临界状态,因此加强局部高风险区域粉尘管理,有利于维持整个系统的平衡与稳定.通过对空间分布和多重分形特征的分析，辨识出港口内2、3、9 10、13 和14号等6 个高风险区域.对于2、9 和14号监测点所在区域应加强道路管理，减少车流量，积极推广使用国VI排放标准的柴油，确保道路清扫工作规范化，要求车辆加装催化型颗粒物捕集器.增加3、10 和13 号监测点所在区域的绿化系数,可有效降尘.在3 号和13 号监测点附近的装卸作业中，应降低落差高度，在落料处增设除尘设施，为降低散货物料运输过程中的PM2.5污染，转运过程宜采用密闭输送。3结 论1）各监测点的PM2.5日均浓度均具有长期持续性和多重分形特征,过去产生的PM2.5会影响未来一段时间的PM2.5浓度.作业扬尘、道路扬尘和柴油车尾气是导致PM2.5浓度长期持续性增强的主要因素，因此应注重污染的长期性防治.堆场装卸作业、物料转运作业和车辆行驶中产生的扬尘以及汽车尾气导致PM2.5日均浓度波动性增大，应注重即时性防治.2）多重分形全面地描述了港口PM2.5系统的复杂动态特性，基于EEMD的MFDFA算法表现出优越的性能,多重分形谱参数（、h（）和f(）可作为定量评价区域PM2.5浓度演化的指标.3）镇江港三面建有围墙与防风网，具有半封闭性,风速多为3 级及以下(PM2.5时均浓度数据占比为9 9.7 6%）.港口PM2.5系统已经达到自组织临界状态，存在发生大气高污染的风险，生产作业中局部微小的PM2.5污染会通过SOC机制形成较大污染.因此,加强局部高风险区域粉尘管理，有利于维持整个系统的平衡与稳定。参考文献(References)1张晓旭.镇江港：打造港口经济绿色样板逆势跑出产业强市加速度”J.大众投资指南，2 0 2 0（2 8）：3-5.ZHANG X X.Zhenjiang port:build a green model ofport economy and run out of the acceleration of strongindustrial city against the trend J.Public InvestmentGuide,2020(28):3-5.(in Chinese)2 董德明，杜山山,黄亚司，等长春市冬春季环境空气中PM2.5污染特征与来源解析 J吉林大学学报（理学版）,2 0 2 0,58(5):12 7 8-12 8 6.DONG D M,DU S S,HUANG Y S,et al.Pollution cha-racteristics and source apportionment of PM2.s in winterand spring in Changchun city J.Journal of Jilin Uni-versity(Science Edition),2020,58(5):1278-1286.(in Chinese)3谢焕丽，何红弟.香港港口PM2.5和PMio的多重分形特征 J.大气与环境光学学报,2 0 19,14（3）：17 9-190.XIE H L,HE H D.Multifractal property between PM2.5and PMio in Hongkong port J.Journal of Atmosphericand Environmental Optics,2019,14(3):179-190.(inChinese)4WANG Q Z.Multifractal characterization of air pollutedtime series in China J.Physica A:Statistical Mecha-nics and its Applications,D01:10.1016/j.physa.2018.09.065.5ZHANG C,WANG X F,CHEN S B,et al.A study ondaily PM2.5 concentrations in Hong Kong using the EMD-based MFDFA method J.Physica A:Statistical Me-chanics and its Applications,DOI:10.1016/j.physa.2019.121182.6 WU Z H,HUANG N E.Ensemble empirical mode de-composition:a noise-assisted data analysis method J.Advances in Adaptive Data Analysis,2009,1(1):1-41.7 曹建军,那嘉明,汤国安，等.DEM随机噪声误差去除方法研究：以黄土沟沿线地形剖面为例J.地理与地理信息科学,2 0 17,3 3(2):12-18,12 7.CAO J J,NA J M,TANG G A,et al.Method for noise re-moval of DEMs:a case study on terrain profile of loessshoulder line J.Geography and Geo-Information Sci-ence,2017,33(2):12-18,127.(in Chinese)8 AACBAZO M N,NCOBI G K,ALAMOU E,et al.As-sessing nonlinear dynamics and trends in precipitation byensemble empirical mode decomposition（EEM D)a n d744江苏大学学报（自然科学版）第44卷fractal approach in Benin Republic(West Africa)J:Complexity,D0I:10.1155/2021/3689397.9 COMEZ-COMEZ J,CARMONA-CABEZAS R,ARIZA-VILLAVERDE A B,et al.Multifractal detrended fluc-tuation analysis of temperature in Spain(19602019)J.Physica A:Statistical Mechanics and its Applica-tions,D0I:10.1016/j.physa.2021.126118.10 师冲,任燕,汤何胜,等.基于改进MF-DFA 和随机森林的液压阀故障诊断 J.重庆理工大学学报（自然科学),2 0 2 1,3 5(6):2 8 0-2 8 8.SHI C,REN Y,TANG H S,et al.Fault diagnosis of hy-draulic valve based on improved MF-DFA and randomforest J.Journal of Chongqing University of Technology(Natural Science),2021,35(6):280-288.（i n C h i-nese)11孙世政，韩宇，党晓园，等.嵌人式角位移传感器动态测量误差建模与补偿 J.重庆理工大学学报（自然科学),2 0 2 2,3 6(3):7 8-8 5.SUN S Z,HAN Y,DANG X Y,et al.Research on mo-deling and compensation for embedded angular displace-ment sensor dynamic measurement errors J.Journal ofChongqing University of Technology(Natural Science),2022,36(3):78-85.(in Chinese)12郭翠娟，李德冲，荣锋，等。基于最大相关波形延拓改进的EEMD方法 J重庆邮电大学学报（自然科学版),2 0 17,2 9(6):7 6 8 7 7 5.GUO C J,LI D C,RONG F,et al.An improved EEMDmethod based on maximal correlation wave form exten-sion J.Journal of Chongqing University of Posts andTelecommunications（Na t u r a l Sc i e n c e Ed i t i o n）,2 0 17,29(6):768-775.(in Chinese)13 林杜娟，刘春琼，黄红良，等.成都冬季PM2.5演化的EE-MD分解及自组织临界态 J.环境科学与技术，2019,42(8):133-141.DU J,LIU C Q,HUANG H L,et al.Multiscale EEMDdecomposition of PM2.s in Chengdu during winter andself-organized criticality of air pollution J.Environ-mental Science and Technology,2019,42(8):133-141.(in Chinese)14LIU D,XIAO Z H,HU X,et al.Feature extraction of ro-tor fault based on EEMD and curve code J.Measure-ment,D01;10.1016/j.measurement.2018.12.009.15封杨，黄筱调，陈捷，等。大型回转支承非平稳振动信号的EEMD-PCA降噪方法J南京工业大学学报（自然科学版）,2 0 15,3 7(3):6 1-6 6.FENG Y,HUANG X D,CHEN J,et al.EEMD-PCAbased denosing method for non-stationary vibration sig-nals of large-size slewing bearings J.Journal of Nan-jing Tech University(Natural Science Edition),2015,37(3):61-66.(in Chinese)16ZHANG X K,ZHANG G,QIU L,et al.A modified multi-fractal detrended fluctuation analysis（M FD FA）a p-proach for multifractal analysis of precipitation inDongting Lake Basin,China J.Water,DOI:10.3390/w11050891.17吴波,刘春琼，张娇，等.COVID-19期间区域大气高污染发生的非线性动力机制 J：中国环境科学,2 0 2 1，41(5):2028-2039.WU B,LIU C Q,ZHANG J,et al.Nonlinear dynamicmechanism of regional high air pollution during COVID-19 lockdown J.China Environmental Science,2021,41(5):2028-2039.(in Chinese)(责任编辑赵鸥)