一种基于HITRAN数据库的吸收光谱快速获取方法与流程

一种基于hitran数据库的吸收光谱快速获取方法
技术领域
1.本发明属于大气光谱计算领域，涉及一种用于hitran数据库的分子吸收光谱的高速计算方法。


背景技术：

2.基于hitran数据库的大气吸收光谱在许多领域都有重要应用，如可以用于大气透过率分析、痕量气体探测、气体温度遥感等。
3.在众多计算hitran吸收光谱的方法中，逐线积分法是公认的最准确的方法。采用该方法计算时，需要知道计算光谱范围内每条谱线的参数，由此得到每条谱线的展宽线型，然后再按波数积分，计算速度极慢。
4.为了解决计算速度问题，科学家们发展了带模式计算、k
‑
分布计算等快速计算方法。这些方法虽然部分解决了计算速度问题，但是其前提都是根据hitran数据库事先进行逐线积分计算，再通过复杂算法得到近似的结果。其计算精度和光谱分辨率都受到事先计算的数据库的限制，使用的灵活性不足。


技术实现要素：

5.本发明解决的技术问题是：克服现有逐线积分方法计算速度的不足，设计了一种分子吸收光谱的高速计算方法，该方法可以大幅提高不同温度和气压参数下的分子吸收光谱计算速度。
6.本发明的技术解决方案是：
7.一种基于hitran数据库的吸收光谱快速获取方法，包括如下步骤：
8.(1)将光谱吸收线的强度按指定的光谱间隔进行合并；
9.(2)求出各光谱吸收线的加权展宽参数；
10.(3)依据仿真精度按照不同的加权展宽参数对光谱吸收线进行分组；
11.(4)对每组光谱吸收线求其综合加权展宽线型函数；
12.(5)通过快速傅里叶变换计算每组光谱吸收线的线强与其综合加权展宽线型函数的卷积，得到每组光谱吸收线的光谱吸收系数；
13.(6)将各组光谱吸收线的光谱吸收系数相加得到最终的吸收光谱。
14.所述的将光谱吸收线的强度按指定的光谱间隔进行合并，具体为：从hitran数据库中提取出分子x的处于最小波数ν
min
至最大波数ν
max
间线强数组{s
i
}，对应波数数组{ν
i
}，以及对应的温度展宽参数数组{α
di
}和气压展宽参数数组{α
li
}，i＝1,2,3
…
n为线强及波数的序号，对于等间隔波数数组{ν
0k
}中的波长位置ν
0k
，从波数数组{ν
i
}中找到满足以下条件的波数ν并组成临时数组{ν
sj
}：
[0015][0016]
其中，j＝1,2,3
…
o表示满足上述条件的波数序号，δν0为指定的光谱间隔，k＝1,
2,3
…
m为等间隔光谱序号。
[0017]
所述的求出各光谱吸收线的加权展宽参数，具体为：
[0018]
从{s
i
}、{α
di
}和{α
li
}中提取出与{ν
sj
}对应的参数数组{s
sj
}、{α
dsj
}和{α
lsj
}，计算出对应ν
0k
的合并线强s
0k
、加权温度展宽参数α
d0k
、加权气压展宽参数α
l0k
以及综合加宽α
a0k
：
[0019]
s
0k
＝∑s
sj
[0020][0021][0022][0023]
所述的依据仿真精度按照不同的加权展宽参数对光谱吸收线分组，具体为：设定相对误差要求e％，将综合加宽α
a0k
分为p组，
[0024][0025]
式中符号表示向上取整，其中第l组的综合加宽α
a0l
都满足以下条件：
[0026][0027]
l＝1,2,3
…
p。
[0028]
所述的对每组光谱吸收线求其综合加权展宽线型函数，具体为：
[0029]
计算第l组综合加宽对应的合并线强{s
alm
}、加权温度展宽{α
adlm
}、加权气压展宽{α
allm
}，对应波数{ν
a0lm
}，式中角标m为第l组的光谱线序号m＝1,2,3
…
q；
[0030]
设置与{ν
0k
}等长的线强数组{s
0tk
}，并将该数组中各位置都设为0；
[0031]
根据波数数组{ν
a0lm
}中的波数位置将第l组合并线强{s
alm
}置入{s
0tk
}中的相应位置生成第l组线强数据{s
0tk
}，并由{ν
0k
}与{s
0tk
}组成线强函数s
k
(ν)；
[0032]
计算第l组综合加宽的加权温度展宽和加权气压展宽：
[0033][0034][0035]
由此得到温度展宽线型f
dl
(ν)为：
[0036]
气压展宽线型f
ll
(ν)为：
[0037]
并据此得到综合加权展宽线型函数f
l
(ν)为：式中表示傅里叶变换，表示逆傅里叶变换。
[0038]
所述的通过快速傅里叶变换计算每组光谱吸收线的线强与其综合加权展宽线型函数的卷积，得到每组光谱吸收线的光谱吸收系数，具体为：
[0039]
由傅里叶变换及逆傅里叶变换计算第l组光谱吸收系数σ
l
(ν)：
[0040][0041]
所述的最终吸收光谱系数σ(ν)＝∑σ
l
(ν)。
[0042]
一种处理装置，包括：
[0043]
存储器，用于存储计算机程序；
[0044]
处理器，用于从所述存储器调用并运行所述计算机程序，以执行上述的方法。
[0045]
一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序或指令，当所述计算机程序或指令被执行时，实现上述的方法。
[0046]
一种计算机程序产品，所述计算机程序产品包括指令，当所述计算机程序产品在计算机上运行时，使得计算机执行上述的方法。
[0047]
本发明与现有技术相比的优点在于：
[0048]
(1)与逐线积分每个光谱线使用不同的线型函数进行计算相比，本发明对所有光谱线采用同一综合线型展宽函数实现了吸收截面的快速计算。
[0049]
(2)为了提高计算精度，本发明根据综合加宽及误差要求对吸收光谱线进行分组，使每组光谱线的综合线型展宽函数与该组光谱线原对应的线型展宽函数相对误差满足精度要求；
[0050]
(3)本发明对分组后的光谱线分别与其对应的综合线型展宽函数进行卷积(由快速傅里叶变换及逆傅里叶变换实现)计算得到分组光谱吸收截面，再将各组光谱吸收截面进行累加得到吸终的吸收光谱。与逐线积分相比速度可以提高2～3个数量级，同时依据加权展宽参数对谱线分组从而将计算精度保持在一定的范围之内，最终实现了计算速度和精度的平衡。
附图说明
[0051]
图1为本发明方法的流程框图；
[0052]
图2为本发明实施例中根据逐线积分计算方法和本发明计算方法得到的吸收光谱对比图；
[0053]
图3为图2吸收光谱的局布放大对比图。
具体实施方式
[0054]
如图1所示，为本发明方法的主要流程，主要包括：第一步将光谱吸收线的线强按指定的光谱间隔进行合并，第二步求出各光谱吸收线的加权展宽参数，第三步依据仿真精度按照不同的加权展宽参数对光谱吸收线分组，第四步对每组光谱吸收线求其综合加权展宽线型函数，第五步通过快速傅里叶变换(fft)计算每组光谱吸收线的线强与其综合加权
展宽线型函数的卷积得到每组光谱吸收线的吸收系数，第六步将各组光谱吸收线的光谱吸收系数相加得到最终的吸收光谱。
[0055]
具体的，针对光谱研究中，需要计算某分子x从波数ν
min
至波数ν
max
间的间隔为δν0的波数数组{ν
0k
}的吸收光谱系数，可以采用如下方法进行处理。
[0056]
1、线强合并
[0057]
从hitran数据库中提取出分子x的处于最小波数ν
min
至最大波数ν
max
间线强数组{s
i
}、对应的波数数组{ν
i
}，以及对应的温度展宽参数数组{α
di
}和气压展宽参数数组{α
li
}。x为分子名称，ν为波数，s为吸收线强，i＝1,2,3
…
n，为线强及波数的序号，角标min表示最小，角标max表示最大。
[0058]
设等间隔波数数组为{ν
0k
}，δν0为指定的光谱间隔，则：
[0059]
ν
0k
＝v
min
(k
‑
1)δν0[0060]
式中，k＝1,2,3
…
m为等间隔光谱序号，且有ν
0m
≤ν
max
。
[0061]
对于波长位置ν
0k
，从波数数组{ν
i
}中找到满足以下条件的波数ν并组成临时数组{ν
sj
}：
[0062][0063]
角标s表示临时波数，j＝1,2,3
…
o表示满足上述条件的波数序号。
[0064]
从{s
i
}、{α
di
}和{α
li
}中提取出与{ν
sj
}对应的参数数组{s
sj
}、{α
dsj
}和{α
lsj
}，按照以下述公式计算出对应ν
0k
的合并线强s
0k
、加权温度展宽参数α
d0k
、加权气压展宽参数α
l0k
以及综合加宽α
a0k
：
[0065]
s
0k
＝∑s
sj
[0066][0067][0068]
另外再定义综合加宽α
a0k
，其计算公式为：
[0069][0070]
2、根据误差分组
[0071]
设相对误差要求为e％，则可以按照该误差将综合加宽α
a0k
分为p组，
[0072][0073]
上式中的符号表示向上取整运算。
[0074]
其中第l(l＝1,2,3
…
p)组的综合加宽α
a0l
都满足以下条件：
[0075][0076]
记对应第l组综合加宽α
a0l
的合并线强为{s
alm
}、加权温度展宽为{α
adlm
}、加权气压展宽为{α
allm
}，对应波数为{ν
a0lm
}，式中角标m为第l组的光谱线序号m＝1,2,3
…
q。
[0077]
3、计算第k组综合加宽的吸收光谱
[0078]
根据上述条件得到第l组对应合并线强{s
alm
}、加权温度展宽{α
adlm
}、加权气压展宽{α
allm
}，对应波数{ν
a0lm
}后，设置与{ν
0k
}等长的线强数组{s
0tk
}，并将该数组中各位置都设为0，角标t表示临时数组。
[0079]
根据波数数组{ν
a0lm
}中的波数位置将第l组合并线强{s
alm
}置入{s
0tk
}中的相应位置生成第l组线强数据{s
0tk
}。由{ν
0k
}与{s
0tk
}组成线强函数s
k
(ν)；
[0080]
计算第l组综合加宽的加权温度展宽和加权气压展宽：
[0081][0082][0083]
已知，温度展宽线型f
dl
(ν)为：
[0084]
已知，气压展宽线型f
ll
(ν)为：
[0085]
则线型函数f
l
(ν)为：
[0086]
式中表示傅里叶变换，表示逆傅里叶变换。
[0087]
由傅里叶变换及逆傅里叶变换计算光谱吸收系数σ
l
(ν)：
[0088][0089]
4、合并上述p组吸收光谱系数，计算最终吸收光谱系数：
[0090]
σ(ν)＝∑σ
l
(ν)。
[0091]
实施例
[0092]
若在光谱研究中，需要计算h2o分子第一同位素h
16
oh从波数ν
min
＝9000cm
‑1至波数ν
max
＝11000cm
‑1间的间隔为δν0＝0.001cm
‑1的波数数组{ν
0i
}的吸收光谱系数，可以应用本发明方法进行如下计算。
[0093]
1、线强合并
[0094]
从hitran数据库中提取出分子h
16
oh的处于波数ν
min
＝9000cm
‑1至波数ν
max
＝11000cm
‑1间线强数组{s
i
}及其对应波数数组{ν
i
}，并计算温度展宽数组{α
di
}和气压展宽数组{α
li
}，式中i＝1,2,3
…
n。
[0095]
以波长位置ν
0i
为例，从波数数组{ν
i
}中找到满足以下条件的波数并组成数组{ν
sj
}：
[0096]
ν
0i
‑
0.0005≤ν<ν
0i
0.0005
[0097]
从{s
i
}、{α
di
}和{α
li
}中提取出与{ν
sj
}对应的参数数组{s
sj
}、{α
dsj
}和{α
lsj
}，计算出对应ν
0i
的合并线强s
0i
、加权温度展宽α
d0i
和加权气压展宽α
l0i
，由此得到综合加宽α
a0i
。
[0098]
2、根据误差分组
[0099]
设相对误差要求为5％，则可以按照该误差将综合加宽α
a0i
分为20组，其中第k组的综合加宽α
a0
都满足以下条件：
[0100][0101]
对应第k组综合加宽的合并线强为{s
aki
}、加权温度展宽为{α
adki
}、加权气压展宽为{α
alki
}，对应波数为{ν
a0ki
}。
[0102]
3、计算第k组吸收光谱
[0103]
依据第k组的光谱参数{s
aki
}、{α
adki
}、{α
alki
}、{ν
a0ki
}，计算第k组光谱吸收系数。
[0104]
设置与{ν
0i
}等长的线强数组{s
0i
}并将该数组中各位置都设为0。根据波数数组{ν
a0ki
}中的波数位置将第k组合并线强{s
aki
}置入{s
0i
}中的相应位置生成第k组线强数据{s
0ki
}，则其线强函数为s
k
(ν)。按照以下公式计算第k组加权温度展宽和加权气压展宽：
[0105][0106][0107]
根据温度展宽和气压展宽参数计算线型函数f
k
(ν)，由傅里叶变换计算光谱吸收系数：
[0108][0109]
4、合并上述20组吸收光谱，计算最终吸收光谱：
[0110]
σ(ν)＝∑σ
k
(ν)。
[0111]
根据上述方法计算得到的吸收光谱见图2和图3(误差要求分别为5％和100％)，由图中可以看出，原逐线积分计算时间为383.667970秒，而本发明5％误差要求的计算时间为7.474211秒，本发明100％误差要求的计算时间为1.798916秒)。
[0112]
由此可以看出，采用本发明方法，5％误差要求(实际最大误差低于3％)的计算速度提高了约51倍，100误差要求(实际最大误差低于30％)的计算速度提高了213倍。
[0113]
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。