一种基于FTIR光谱的城市大气碳排放分布检测方法与流程

一种基于ftir光谱的城市大气碳排放分布检测方法
技术领域
1.本发明涉及co2碳排放技术领域，尤其涉及一种基于ftir光谱的城市大气碳排放分布检测方法。


背景技术：

2.大气中的温室气体主要包括二氧化碳(co2)、甲烷(ch4)、全氟碳化物(pfcs)、氢氟烃(hfcs)、氧化亚氮(n2o)等，这些气体大量吸收地表受热后向外发射的长波辐射，进而对地球大气产生温室效应。自工业化时代以来，人类生产活动导致全球温室气体浓度急剧升高，进而引起全球气候发生改变，二氧化碳排放主要来自人为的工业活动、交通和居民生活，由于城市的碳排放总量由多个功能分区排放累加得到，因此需要获得城市内部不同功能分区的二氧化碳排放详细状况。
3.一般来说，不同城市功能分区的二氧化碳排放量主要来自统计数据，原位测量技术最早是采用化学检测技术来测量大气中的碳排放，地基遥感和卫星遥感主要是利用光学方法和技术对碳排放进行监测。对大气中温室气体的监测主要是原位测量，具体有气相色谱法、奥氏气体分析法等，但是这些监测技术受人为因素的干扰，精度较差，且时效性差，然而，由于缺乏城市内部的空间信息，使得对于估算高空间分辨率尺度的功能分区二氧化碳排放带来一定困难，因此需要一种基于ftir光谱的城市大气碳排放分布检测方法。


技术实现要素：

4.本发明的目的是要提供一种基于ftir光谱技术的城市大气co2碳排放分布的检测系统。
5.为达到上述目的，本发明是按照以下技术方案实施的：
6.本发明包括以下步骤：
7.a采集城市大气中ftir光谱数据，对采集到的光谱数据进行预处理；
8.b通过对所述光谱数据采用极限学习算法进行光谱反演，获得光谱数据中碳排放的时间序列；
9.c对获得的碳排放的柱浓度进行误差分析，使用雅可比矩阵可以评估不同高度的碳排放浓度变化；
10.d对所述时间序列的分布变化提取碳排放气体浓度分布，进行温室气体的碳排放通量计算，分析通量计算误差。
11.进一步地，所述光谱数据的预处理包括校正,去噪处理和归一化处理。
12.进一步地，在步骤b中所述获得光谱数据中碳排放的时间序列的方法包括：
13.a加权欧式距离作为碳排放预测区间评估指标，筛选历史样本单元并优化极限学习机训练集；
14.b提出极限学习机参数混合寻优算法，通过精英保留策略遗传算法与分位数回归优化极限学习机模型隐层输入及输出权重与偏置参数，并采用训练后的模型预测碳排放区
间；
15.c基于时间序列变化的历史数据进行碳排放区间预测。
16.进一步地，在步骤d中所述提取碳排放气体浓度分布的方法包括
17.1)设定时间序列的长度为n，将时间序列分成m个时间段，每个时间段的长度为l(l《n)；
18.2)对于每个时间段，根据不同时间点的碳排放浓度值进行设定权重因子；
19.3)计算碳排放浓度平均值,将每个时间段的碳排放浓度平均值作为该时间段内碳排放气体浓度分布的代表值,假设第i个时间段的碳排放浓度平均值为mi，则有：
20.mi＝(∑(c1i*wi1)+∑(c2i*wi2)+...+∑(cnii*win))/∑(wi1+wi2+...+win)
21.其中，c1i，c2i，...，cnii是第i个时间段内每个时间点的碳排放浓度值，wi1，wi2，...，win是对应的权重因子。
22.进一步地，所述权重因子的计算方法包括
23.a计算该时间段内所有时间点的碳排放浓度值与该时间段平均碳排放浓度值的差值的绝对值之和：
24.s＝∑(|c1i-m1|+|c2i-m2|+...+|c12i-m12|)
25.b设定权重因子：
26.c如果s《10，则该时间段的权重因子为1；
27.如果10《＝s《20，则该时间段的权重因子为s/10；
28.如果s》＝20，则该时间段的权重因子为0。
29.进一步地，在步骤c中所述雅可比矩阵由最优估计方法来确定状态向量的变化值(δx)，具体关系式为：
[0030][0031]
其中j是雅可比矩阵，s
∈
是测量误差矩阵，假定是由光谱噪声值组成的对角矩阵，s0是先验协方差矩阵，其包含估计的先验误差值，r是光谱拟合残差，u表示当前状态向量最优估计值，u0表示先验状态向量。
[0032]
一种电子设备，包括：
[0033]
处理器；以及
[0034]
被安排成存储计算机可执行指令的存储器，所述可执行指令在被执行时使所述处理器执行所述权利要求1～6之任一所述方法。
[0035]
本发明的有益效果是：
[0036]
本发明是一种基于ftir光谱的城市大气碳排放分布检测方法，与现有技术相比，本发明具有以下技术效果：
[0037]
1本发明通过ftir光谱对碳排放分布检测能全面充分的反映监测区域的碳排放的现状和发展趋势，从而保证监测结果的准确性，为节能减排、生态环境保护等规划提供科学依据，对所述时间序列的分布变化提取碳排放气体浓度分布处理和分析，极大的缩短了数据处理时间，提高了碳监测工作效率，该方法不仅可以提高分析精度，同时具有较好的可解释性，可以直接应用于基于ftir光谱技术的城市大气co2碳排放分布的检测系统中。
附图说明
[0038]
图1为本发明一种基于ftir光谱的城市大气碳排放分布检测方法的步骤流程图；
[0039]
图2为本发明一种基于ftir光谱的城市大气碳排放分布检测方法的电子设备图；
具体实施方式
[0040]
下面以及具体实施例对本发明作进一步描述，在此发明的示意性实施例以及说明用来解释本发明，但并不作为对本发明的限定。
[0041]
本发明基于ftir光谱的城市大气碳排放分布检测方法包括以下步骤：
[0042]
如图1所示，在本实施例中，包括以下步骤：
[0043]
本发明包括以下步骤：
[0044]
a采集城市大气中ftir光谱数据，对采集到的光谱数据进行预处理；
[0045]
b通过对所述光谱数据采用极限学习算法进行光谱反演，获得光谱数据中碳排放的时间序列；
[0046]
在反演计算过程中，反演多采用缩放先验廓线的方法，垂直柱浓度和气体柱平均dmfs的关系是分别为：
[0047]
其中：pgr是地表气压，mh2o是水汽分子质量，m
dryair
是干燥空气分子质量，c
h2o
是水汽柱浓度。以大气中co2的柱浓度作为基底来计算目标气体柱平均干空气混合比浓度可以达到降低反演系统误差的作用，气体的柱平均干空气混合比(column-averageddry-airmolefraction，dmfs)由反演的气体的柱浓度计算，其关系式如下
[0048][0049]
c对获得的碳排放的柱浓度进行误差分析，使用雅可比矩阵可以评估不同高度的碳排放浓度变化；
[0050]
柱浓度是采用从地面到高空垂直柱中的碳排放、在0℃和1个标准大气压下的总层厚度来反映大气中碳排放二氧化碳的含量.
[0051]
d对所述时间序列的分布变化提取碳排放气体浓度分布，进行温室气体的碳排放通量计算，分析通量计算误差。
[0052]
在本实施例子中，所述光谱数据的预处理包括校正,去噪处理和归一化处理。
[0053]
在本实施例子中，在步骤b中所述获得光谱数据中碳排放的时间序列的方法包括：
[0054]
a加权欧式距离作为碳排放预测区间评估指标，筛选历史样本单元并优化极限学习机训练集；
[0055]
b提出极限学习机参数混合寻优算法，通过精英保留策略遗传算法与分位数回归优化极限学习机模型隐层输入及输出权重与偏置参数，并采用训练后的模型预测碳排放区间；
[0056]
极限学习机是一种单层前馈神经网络随机生成隐层输入权重和偏置，确定隐层输出权重具有n个输入层神经元、l个隐层神经元、m个输出层神经元的elm结构,对于n个离散训练样本输入数据为
x
xi＝[x
i1
，x
i2
，
…
，x
in
]
t
目标输出数据为ti＝[t
i1
，t
i2
，
…
，t
im
]
t
.
[0057]
c基于时间序列变化的历史数据进行碳排放区间预测。
[0058]
在本实施例子中，在步骤d中所述提取碳排放气体浓度分布的方法包括
[0059]
1)设定时间序列的长度为n，将时间序列分成m个时间段，每个时间段的长度为l(l《n)；
[0060]
2)对于每个时间段，根据不同时间点的碳排放浓度值进行设定权重因子；
[0061]
3)计算碳排放浓度平均值,将每个时间段的碳排放浓度平均值作为该时间段内碳排放气体浓度分布的代表值,假设第i个时间段的碳排放浓度平均值为mi，则有：
[0062]
mi＝(∑(c1i*wi1)+∑(c2i*wi2)+...+∑(cnii*win))/∑(wi1+wi2+...+win)
[0063]
其中，c1i，c2i，...，cnii是第i个时间段内每个时间点的碳排放浓度值，wi1，wi2，...，win是对应的权重因子。
[0064]
在本实施例子中，所述权重因子的计算方法包括
[0065]
a计算该时间段内所有时间点的碳排放浓度值与该时间段平均碳排放浓度值的差值的绝对值之和：
[0066]
s＝∑(|c1i-m1|+|c2i-m2|+...+|c12i-m12|)
[0067]
b设定权重因子：
[0068]
c如果s《10，则该时间段的权重因子为1；
[0069]
如果10《＝s《20，则该时间段的权重因子为s/10；
[0070]
如果s》＝20，则该时间段的权重因子为0。
[0071]
在本实施例子中，在步骤c中所述雅可比矩阵由最优估计方法来确定
[0072][0073]
在本实施例子中，雅可比矩阵用于确定假定状态向量的最佳变化并计算柱浓度反演的平均核，平均核表示不同高度处浓度变化对反演结果扰动的影响。雅可比矩阵由最优估计方法来确定状态向量的变化值(δx)，具体关系式为：
[0074][0075]
其中j是雅可比矩阵，s
∈
是测量误差矩阵，假定是由光谱噪声值组成的对角矩阵，s0是先验协方差矩阵，其包含估计的先验误差值，r是光谱拟合残差，u表示当前状态向量最优估计值，u0表示先验状态向量。
[0076]
图2是本技术的一个实施例电子设备的结构示意图。请参考图2，在硬件层面，该电子设备包括处理器，可选地还包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(random-accessmemory，ram)，也可能还包括非易失性存储器(non-volatilememory)，例如至少1个磁盘存储器等。当然，该电子设备还可能包括其他业务所需要的硬件。
[0077]
处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是isa(industrystandardarchitecture，工业标准体系结构)总线、pci(peripheralcomponentinterconnect，外设部件互连标准)总线或eisa(extendedindustr
ystandardarchitecture，扩展工业标准结构)总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图2中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
[0078]
存储器，用于存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。
[0079]
处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行，在逻辑层面上形成一种基于基于ftir光谱的城市大气碳排放分布检测装置。处理器，执行存储器所存放的程序，并具体用于执行前述任意一种基于基于ftir光谱的城市大气碳排放分布检测方法。
[0080]
上述如本技术图1所示实施例揭示的一种基于基于ftir光谱的城市大气碳排放分布检测方法可以应用于处理器中，或者由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器，包括中央处理器(centralprocessingunit，cpu)、网络处理器(networkprocessor，np)等；还可以是数字信号处理器(digitalsignalprocessor，dsp)、专用集成电路(applicationspecificintegratedcircuit，asic)、现场可编程门阵列(field－programmablegatearray，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本技术实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本技术实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
[0081]
该电子设备还可执行图1中一种基于基于ftir光谱的城市大气碳排放分布检测方法，并实现图1所示实施例的功能，本技术实施例在此不再赘述。
[0082]
本技术实施例还提出了一种计算机可读存储介质，该计算机可读存储介质存储一个或多个程序，该一个或多个程序包括指令，该指令当被包括多个应用程序的电子设备执行时，执行前述的任意一种基于基于ftir光谱的城市大气碳排放分布检测方法。
[0083]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0084]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0085]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0086]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0087]
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
[0088]
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。
[0089]
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitorymedia)，如调制的数据信号和载波。
[0090]
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0091]
本领域技术人员应明白，本技术的实施例可提供为方法、系统或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0092]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，包括：a采集城市大气中ftir光谱数据，对采集到的光谱数据进行预处理；b通过对所述光谱数据采用极限学习算法进行光谱反演，获得光谱数据中碳排放的时间序列；c对获得的碳排放的柱浓度进行误差分析，使用雅可比矩阵可以评估不同高度的碳排放浓度变化；d对所述时间序列的分布变化提取碳排放气体浓度分布，进行温室气体的碳排放通量计算，分析通量计算误差。2.根据权利要求1所述基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，所述光谱数据的预处理包括校正,去噪处理和归一化处理。3.根据权利要求1所述基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，在步骤b中所述获得光谱数据中碳排放的时间序列的方法包括：a加权欧式距离作为碳排放预测区间评估指标，筛选历史样本单元并优化极限学习机训练集；b提出极限学习机参数混合寻优算法，通过精英保留策略遗传算法与分位数回归优化极限学习机模型隐层输入及输出权重与偏置参数，并采用训练后的模型预测碳排放区间；c基于时间序列变化的历史数据进行碳排放区间预测。4.根据权利要求1所述基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，在步骤d中所述提取碳排放气体浓度分布的方法包括1)设定时间序列的长度为n，将时间序列分成m个时间段，每个时间段的长度为l(l<n)；2)对于每个时间段，根据不同时间点的碳排放浓度值进行设定权重因子；3)计算碳排放浓度平均值,将每个时间段的碳排放浓度平均值作为该时间段内碳排放气体浓度分布的代表值,假设第i个时间段的碳排放浓度平均值为mi，则有：mi＝(∑(c1i*wi1)+∑(c2i*wi2)+...+∑(cnii*win))/∑(wi1+wi2+...+win)其中，c1i，c2i，...，cnii是第i个时间段内每个时间点的碳排放浓度值，wi1，wi2，...，win是对应的权重因子。5.根据权利要求4所述基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，所述权重因子的计算方法包括a计算该时间段内所有时间点的碳排放浓度值与该时间段平均碳排放浓度值的差值的绝对值之和：s＝∑(|c1i-m1|+|c2i-m2|+...+|c12i-m12|)b设定权重因子：c如果s<10，则该时间段的权重因子为1；如果10<＝s<20，则该时间段的权重因子为s/10；如果s>＝20，则该时间段的权重因子为0。6.根据权利要求1所述基于ftir光谱的城市大气碳排放分布检测方法，其特征在于，在步骤c中所述雅可比矩阵由最优估计方法来确定状态向量的变化值(δx)，具体关系式为：
其中j是雅可比矩阵，s
∈
是测量误差矩阵，假定是由光谱噪声值组成的对角矩阵，s0是先验协方差矩阵，其包含估计的先验误差值，r是光谱拟合残差，u表示当前状态向量最优估计值，u0表示先验状态向量。7.一种电子设备，包括：处理器；以及被安排成存储计算机可执行指令的存储器，所述可执行指令在被执行时使所述处理器执行所述权利要求1～6之任一所述方法。

技术总结
本发明公开了一种基于FTIR光谱的城市大气碳排放分布检测方法，采集城市大气中FTIR光谱数据，对采集到的光谱数据进行预处理；通过对所述光谱数据采用极限学习算法进行光谱反演，获得光谱数据中碳排放的时间序列；对获得的碳排放的柱浓度进行误差分析；对所述时间序列的分布变化提取碳排放气体浓度分布，进行温室气体的碳排放通量计算，分析通量计算误差。本发明本发明通过FTIR光谱对碳排放分布检测能全面充分的反映监测区域的碳排放的现状和发展趋势，从而保证监测结果的准确性，对所述时间序列的分布变化提取碳排放气体浓度分布处理和分析，可以直接应用于基于FTIR光谱技术的城市大气CO2碳排放分布的检测系统中。的城市大气CO2碳排放分布的检测系统中。的城市大气CO2碳排放分布的检测系统中。


技术研发人员：张丽丽 包正义 张文豪 潘帅帅 王春梅 王更科 臧文乾 王栋
受保护的技术使用者：廊坊空间信息技术研发服务中心
技术研发日：2023.06.16
技术公布日：2023/9/12