基于EMD的二氧化硫紫外差分吸收光谱降噪重构方法与流程

本发明属于气体浓度检测领域，涉及差分光学吸收光谱，尤其是一种基于EMD的二氧化硫紫外差分吸收光谱降噪重构方法。

背景技术：

随着我国火电厂等工业中的脱硫建设大力进行，“十二五”期间我国SO2的排放总量降低了6％。在“十三五”期间我国SO₂的排放总量继续巩固SO₂减排的态势，计划实现10％～15％的减排目标，实现火电厂等工业中SO₂的超低排放。这对基于紫外差分吸收光谱技术的SO₂浓度的检测提出了检测精度进一步提高和检测下限进一步降低的要求，其本质是需要更高信噪比的低浓度SO₂紫外差分吸收光谱处理方法。

二氧化硫(SO₂)和氮氧化物(NOx)是大气污染的主要成分,其中固定污染源的燃烧排放是主要来源之一。烟气连续监测系统(Continuous Emission Monitoring System，简称为CEMS)是环保监测的重要系统，其对污染源排放的烟气进行连续、实时地监测，其功能包含实时测量各个污染源排放烟气中的SO₂，NOx的排放量。差分吸收光谱技术(DOAS)以原理和结构简单、响应速度快、精度高等优点，被作为烟气连续监测系统(CEMS)采用的重要方法得到广泛应用。目前，这种应用主要在280～320nm的紫外波长段对SO₂的紫外差分吸收光谱信号进行特征提取，并通过一定的反演算法完成实时SO₂浓度的反演计算，但在数据采集后因缺乏较为统一的标准化处理，以及缺少相应的差分后重构计算，在精度上还存在一定欠缺。

经检索，发现两篇与本专利内容相关的文献，其中，《量子电学学报》2001年第S1期提供了一篇《差分光学吸收光谱(DOAS)技术在烟气监测中的应用研究》，该文将差分光学吸收光谱(Differential Optical Absorption Spectroscoy；DOAS)技术中的浓度反演方法运用到烟气中NO、NO₂、NH₃和SO₂的监测.利用它们在200～250nm附近具有不同强度吸收的特性，测量烟气的吸收光谱，用DOAS方法反演出NO、NO₂、NH3和SO₂分子的浓度，这种方法消除了烟气中烟尘、水汽和其它成分的影响，使在线实时测量成为可能，在实验室内对不同气体分子浓度样气进行了测量和浓度反演，经过现场对比测量证明这种方法是正确的.

2006年陈晓虎的一篇硕士论文提供了一篇《燃烧排放烟气监测的差分吸收光谱技术理论分析与实验研究》，该文是针对目前对燃烧排放烟气浓度测量方法多为离线分析，且测量的稳定性、准确性较差的现状，提出采用差分吸收光谱技术(DOAS)测量烟气浓度的方法，实现烟气浓度的在线测量。首先从吸收截面的测量入手，对烟气中主要的污染气体SO₂、NO、NO₂、NH₃的吸收特性进行详细研究，发现这些气体在局部区域具有准周期吸收特征，这为引入新的测量方法提供了契机；其次，对决定DOAS测量精度的被测气体吸收截面的诸多问题进行了深入研究，重点研究了压力、温度对被测气体吸收截面的影响，并针对被测气体的特征提出了相应的修正措施；再次，在实际测量中，要求所使用的光谱仪必须与以吸收截面为基础的参考光谱相匹配，为了满足匹配这一条件，文中将图像复原技术引入到光谱复原，使用维纳滤波和小波消噪实现了光谱的复原，将复原后的光谱与所使用光谱仪的仪器函数卷积，即可实现；最后，在对DOAS测量原理详细分析的基础上，对传统的最小二乘—差分吸收光谱法进行了详细的阐述。

经对比，本案的技术方案较上述两篇文献，其特点是在数据采集标准化处理以及针对实际测量数据的后期重构方面。

技术实现要素：

本发明的目的在于克服现有技术不足，提供一种针对SO₂紫外差分吸收光谱进行数据匹配滤波处理，自适应地完成实测差分吸收光谱特征的基于EMD的二氧化硫紫外差分吸收光谱降噪重构方法。

本发明采用的技术方案是：

一种基于EMD的二氧化硫紫外差分吸收光谱降噪重构方法，其步骤为：

⑴分别对二氧化硫紫外波段的标准吸收截面和实时采集光谱进行处理，提取差分吸收截面S_Tdacs和差分吸收光谱S_Ada；

⑵对S_Tdacs和S_Ada进行标准化，获得S’_Tdacs和S’_Ada；

⑶对S’_Tdacs和S_Ada进行7维EMD分解，获得对应的本征模函数，对应差分吸收截面的7个本征模函数IMF记为IMF(j)_dacs,对应差分吸收光谱的7个本征模函数IMF记为IMF(j)_da，其中j＝1，2，3……7；

⑷利用式1计算每个对应的IMF(j)_dacs和IMF(j)_da的相关系数R(j)，为R(j)设定一个合适的T_R，筛选R(j)>T_R的IMF(j)_da，

⑸计算S’_Tdacs和S’_Ada的相关系数，记为Ra；

⑹利用筛选到的IMF(j)_da和Ra进行差分吸收光谱重构，获得最终输出的S_Rda，重构计算

如式2

S_Rda＝Ra·Σ_jIMF(j)_da (ifR(j)＞T_R) 式2。

而且，所述标准吸收截面是利用SO₂在295～305nm波长段内的差分吸收截面。

而且，所述实时采集光谱是利用SO₂在295～305nm波长段内的差分吸收光谱，然后通过计算机对光谱仪采集到的吸收光谱进行预处理后，得到SO₂的紫外差分吸收光谱。

本发明的积极效果和优点是：

本发明对实际采集光谱进行了重构，而为此本发明还进行了有针对性的标准化处理，结合EMD分解，经过本发明方法的处理后，重构的紫外差分吸收光谱的信噪比有明显提升，在低浓度区间不同浓度SO2的差分吸收光谱特征区分度得以增强，是测量精度能进一步提升。

附图说明

图1是本发明的SO₂紫外差分吸收光谱特征提取方法流程图；

图2是本发明中实时采集光谱的九个不同浓度的SO₂紫外差分吸收光谱图；

图3a是本发明中的标准差分吸收界面的EMD分解图；

图3b是浓度为0.00mg/m³的SO₂差分吸收光谱的EMD分解图；

图3c是浓度为57.43mg/m³的SO₂差分吸收光谱的EMD分解图；

图3d是浓度为85.14mg/m³的SO₂差分吸收光谱的EMD分解图；

图3e是浓度为113.71mg/m³的SO₂差分吸收光谱的EMD分解图；

图3f是浓度为127.43mg/m³的SO₂差分吸收光谱的EMD分解图；

图3g是浓度为171.71mg/m³的SO₂差分吸收光谱的EMD分解图；

图3h是浓度为286.29mg/m³的SO₂差分吸收光谱的EMD分解图；

图3i是浓度为434.57mg/m³的SO₂差分吸收光谱的EMD分解图；

图3j是浓度为573.43mg/m³的SO₂差分吸收光谱的EMD分解图；

图4是本发明重构后的SO₂紫外差分吸收光谱图。

具体实施方式

下面通过附图结合具体实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。

一种基于EMD的二氧化硫紫外差分吸收光谱降噪重构方法，见图1，其步骤为：

⑴分别对标准吸收截面和实时采集光谱进行处理，提取差分吸收截面S_Tdacs和差分吸收光谱S_Ada，标准吸收截面与实时采集光谱君是利用SO₂在295～305nm波长段内分别采集的差分吸收截面和差分吸收光谱，通过计算机对光谱仪采集到的吸收光谱进行预处理后，得到的九个不同浓度的SO₂紫外差分吸收光谱，光谱图见图2；

⑵对S_Tdacs和S_Ada进行标准化，获得S’_Tdacs和S’_Ada；

⑶对S’_Tdacs和S_Ada进行7维EMD分解，获得对应的IMFs，对应差分吸收截面的7个IMF记为IMF(j)_dacs,对应差分吸收光谱的7个IMF记为IMF(j)_da，其中j＝1，2，3……7，其中，差分吸收截面的EMD分解图和九个不同浓度的SO₂紫外差分吸收光谱的EMD分解图见图3a～图3j；

⑷利用式1计算每个对应的IMF(j)_dacs和IMF(j)_da的相关系数R(j)，为R(j)设定一个合适的T_R，筛选R(j)>T_R的IMF(j)_da，

⑸计算S’_Tdacs和S’_Ada的相关系数，记为Ra；

⑹利用筛选到的IMF(j)_da和Ra进行差分吸收光谱重构，获得最终输出的S_Rda，重构计算

如式2

S_Rda＝Ra·Σ_jIMF(j)_da (ifR(j)＞T_R) 式2。

利用标准差来对比评估重构前的原始紫外差分吸收光谱信号和利用发明方法获得的重构紫外差分吸收光谱信号的特征效果。以此验证发明方法对低浓度紫外差分吸收光谱信号信噪比的提高作用，其结果如下表1所示。

表1原始与重构紫外差分吸收光谱特征的标准差对比表

如图4和表1所示，可以看出，经过本发明方法的处理后，重构的紫外差分吸收光谱的信噪比有明显提升，在低浓度区间不同浓度SO₂的差分吸收光谱特征区分度得以增强。特别是0.00mg/m³对应的光谱数据，利用发明方法重构的紫外差分光谱的标准差是原始光谱信号标准差的0.5％。

尽管为说明目的公开了本发明的实施例和附图，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换、变化和修改都是可能的，因此，本发明的范围不局限于实施例和附图所公开的内容。