CN103852443A - 免标定调制光谱的实现方法 - Google Patents
免标定调制光谱的实现方法 Download PDFInfo
- Publication number
- CN103852443A CN103852443A CN201410101717.5A CN201410101717A CN103852443A CN 103852443 A CN103852443 A CN 103852443A CN 201410101717 A CN201410101717 A CN 201410101717A CN 103852443 A CN103852443 A CN 103852443A
- Authority
- CN
- China
- Prior art keywords
- signal
- harmonic
- harmonic signal
- frequency
- component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000000051 modifying Effects 0.000 title claims abstract description 48
- 238000001228 spectrum Methods 0.000 title claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 13
- 230000000875 corresponding Effects 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 8
- 239000000969 carrier Substances 0.000 claims description 6
- 238000000862 absorption spectrum Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000003595 spectral Effects 0.000 description 9
- 238000010521 absorption reaction Methods 0.000 description 7
- 239000007789 gas Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000000243 solution Substances 0.000 description 4
- 238000004611 spectroscopical analysis Methods 0.000 description 4
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 238000002156 mixing Methods 0.000 description 3
- 230000003287 optical Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000001809 detectable Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- QGZKDVFQNNGYKY-UHFFFAOYSA-N ammonia Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwLjAnIGhlaWdodD0nMzAwLjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PScxMzguMCcgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiM0Mjg0RjQnID5OPC90ZXh0Pgo8dGV4dCB4PScxNjUuNicgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiM0Mjg0RjQnID5IPC90ZXh0Pgo8dGV4dCB4PScxOTMuMicgeT0nMTg2LjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiM0Mjg0RjQnID4zPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NS4wJyBoZWlnaHQ9Jzg1LjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PSczNS4wJyB5PSc1My42JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojNDI4NEY0JyA+TjwvdGV4dD4KPHRleHQgeD0nNTEuMCcgeT0nNTMuNicgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzQyODRGNCcgPkg8L3RleHQ+Cjx0ZXh0IHg9JzY3LjAnIHk9JzYyLjknIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTVweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiM0Mjg0RjQnID4zPC90ZXh0Pgo8L3N2Zz4K N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- MYMOFIZGZYHOMD-UHFFFAOYSA-N oxygen Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwLjAnIGhlaWdodD0nMzAwLjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8cGF0aCBjbGFzcz0nYm9uZC0wIGF0b20tMCBhdG9tLTEnIGQ9J00gMTA1LjgsMTcwLjUgTCAxOTUuNywxNzAuNScgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMCBhdG9tLTAgYXRvbS0xJyBkPSdNIDEwNS44LDEyOS41IEwgMTk1LjcsMTI5LjUnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+Cjx0ZXh0IHg9JzM1LjcnIHk9JzE3MC4wJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPHRleHQgeD0nMjQwLjMnIHk9JzE3MC4wJyBjbGFzcz0nYXRvbS0xJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NS4wJyBoZWlnaHQ9Jzg1LjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8cGF0aCBjbGFzcz0nYm9uZC0wIGF0b20tMCBhdG9tLTEnIGQ9J00gMjEuNCw0Ny44IEwgNjIuOCw0Ny44JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0wIGF0b20tMCBhdG9tLTEnIGQ9J00gMjEuNCwzNi4yIEwgNjIuOCwzNi4yJyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8dGV4dCB4PSc2LjEnIHk9JzUzLjYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjNweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8dGV4dCB4PSc2NC4wJyB5PSc1My42JyBjbGFzcz0nYXRvbS0xJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000630 rising Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Abstract
本发明基于可调谐二极管激光吸收光谱，采用适应恶劣环境的波长调制技术，提供了一种免标定调制光谱的实现方法，该方法包含数字锁相信号解调步骤、调制光谱仿真步骤和测量参数提取步骤。其具体内容是，利用数字锁相算法对光电信号解调得到其谐波信号；通过准确的理论计算得到谐波信号的仿真数据库；结合谐波信号和仿真数据库提取流场速度、温度和组分浓度。
Description
免标定调制光谱的实现方法
技术领域
[0001] 本发明属于光学测量技术领域，涉及信号处理与数值仿真等方法，是针对气相分子的吸收光谱测量开发的通用、可靠测量方法。
背景技术
[0002] 光学测量技术在科研实验、环境监测和燃烧场诊断等领域获得了越来越广泛的应用。吸收光谱在理论上需要同时测量吸收前和吸收后的光强信息才能获得分子的吸收信息。在实际应用中，往往采用参考光束或函数拟合等方法确定吸收前的光强，即光谱基线。这些方法难以获得准确的光谱基线信息。此外，调制方法虽然能明显的提高测量信号的信噪比，但是调制光谱的信号强度与调制深度、谱线线型和强度调制等多方面的因素相关，因此调制光谱的定量测量也一直存在困难。
[0003] 美国斯坦福大学K.H.Lyle, R.K.Hanson等人，曾利用调制光谱来增强航空发动机进气道氧气吸收光谱测量的信噪比。由于未在理论上解决调制光谱强度与线型函数间的定量对应关系，此处仅用光谱峰的中心位置来测量速度，浓度的测量则是通过直接吸收光谱实现的。此后，G.B.Rieker从理论上分析了波长调制光谱测量的原理，给出了谐波信号的表达式，着重分析了测量可能的误差来源。研究显示利用If信号归一化2f信号之后，借助准确的理论仿真可实现任意压力、浓度和温度条件下的测量。验证了水分子和二氧化碳浓度和燃烧场温度的测量结果，显示波长调制光谱较之直接吸收光谱信噪比提高约4倍。该文献并未提出详细的实现方法，仅讨论了调制光谱测量的可行性。
[0004] 南京理工大学李宁等人，采用仿真的方法得到谐波信号理论值，利用二次谐波峰谷比来消除初始光强和探测器增益对信号强度的影响。因为波峰和波谷在时间上不是同时出现的，当测量过程中初始光强波动时，波峰、波谷对应的初始光强不相同，不能消除其对信号的影响。另外，该方法必需扫描得到完整的谐波信号才能提取波峰和波谷信息，对于谱线加宽严重或者出现重叠时，往往难以实现。清华大学热能工程系车璐等人，研究了利用奇数次谐波信号重建吸收谱线型。通过对吸收线型进行傅里叶级数和泰勒级数分析，得到线型函数与奇数次谐波的关系，并在实验室条件下验证了氨气的浓度测量。该方法中，谐波信号的强度仍与初始光强相关。初始光强虽然可以预先测量，但是在干扰严重的环境下，诸多非吸收因素导致的衰减致使初始光强在测量过程中出现明显的变化，此时测量结果将出现严重的偏差。因此该方法并未真正实现免标定测量，只能在实验室条件下应用。
[0005] 本发明基于可调谐半导体激光吸收光谱(tunable diode laser absorptionspectroscopy, TDLAS)技术，提供一种免标定调制光谱测量的实现方法，不仅利用了调制光谱的优势，而且从理论上解决了光谱基线提取的问题，可实现真正意义上的免标定测量。
发明内容
[0006] 本发明提供了一种免标定调制光谱的测量方法，实现基于可调谐半导体激光吸收光谱技术的真正意义上的免标定测量。包括:[0007] 步骤一、数字锁相信号解调步骤；
[0008] 步骤二、调制光谱仿真步骤；
[0009] 步骤三、测量参数提取步骤；
[0010] 其特征在于:
[0011] 步骤一的数字锁相信号解调步骤包括:(1)通过数据采集设备读取探测光电信号、采样率和数据长度L ; (2)根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号；(3)将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号；(4)根据采样频率、低通截止频率和滤波阶次，创建低通滤波算法；(5)将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，并分别计算R量；
[0012] 步骤二的调制光谱仿真步骤包括:(I)根据一定条件下的谱线参数进行数值积分计算，得到线型函数；(2)根据实测条件下的激光参数和计算得到的线型函数，计算谐波信号；(3)通过计算不同谱线参数对应的谐波信号，得到仿真数据库；
[0013] 步骤三的测量参数提取步骤包括:(1)提取谐波信号的多普勒频移量，并计算流场速度；(2)通过提取谐波信号的最大值来获得谐波信号的高度值；(3)通过迭代算法，对比仿真谐波信号的高度值和实验谐波信号的高度值，求解流场的温度和组分浓度。
[0014] 本发明利用独特的数据处理方法和准确的数值仿真，解决了调制光谱定量测量的问题。由于算法功能齐全，简化了测量硬件的需求，具有如下优点:
[0015] (I)数字锁相信号解调部分不仅替代了锁相放大器，而且拓展了其功能。数字锁相算法实现了载波与参考信号混频和低通滤波等锁相放大器的全部功能，还有硬件不能实现的功能。数字锁相算法不受调制频率的限制，原则上可以实现任意频率的信号解调；算法可以无限制的扩展通道数量，解调更高阶次的谐波信号。同时利用零相位数字滤波器算法解决了低通滤波带来的相位延迟问题。
[0016] (2)实现了波长调制光谱的免标定测量。调制光谱仿真部分，利用一次谐波归信号一化二次谐波信号，消除初始光强对谐波信号强度的影响，在原理上实现了免标定测量。该部分从吸收光谱的基本方程出发，通过准确的理论模型计算得到谐波信号，为实验测量提供对比值。
[0017] (3)创建了一种独特的算法，解决三维空间搜索解的问题。通过迭代来求解温度、浓度与谐波信号高度比之间的对应关系，降低了仿真计算量。
附图说明
[0018] 图1为本发明的免标定调制光谱的实现方法的流程图；
[0019] 图2为本发明的线型函数和谐波信号仿真示意图；
[0020] 图3为本发明的仿真谐波信号及高度比与温度的关系示意图；
[0021] 图4为本发明的谐波信号多普勒频移提取算法示意图；
[0022] 图5为本发明的温度、浓度迭代求解算法流程图。
具体实施方案[0023] 结合附图和实施案例对本发明的免标定调制光谱的实现方法作进一步详细描述。图1给出了本发明的流程图。
[0024] 本发明免标定调制光谱的实现方法，实施步骤如下:
[0025] 步骤一、数字锁相信号解调
[0026] 免标定调制光谱的实现方法的第一步，是利用数字锁相算法处理探测信号，共分为5个任务，实施方案如下:
[0027] ①、采集测量信号
[0028] 测量信号采集，由数据采集设备(如数据采集卡)读取探测光电信号、采样率和数据长度L。探测光电信号是经过调制的高频载波信号，包含了被测流场的信息。
[0029] ②、生成解调信号
[0030] 该任务根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号。参考信号的采样率和数据长度与探测信号相同，信号频率等于调制频率。
[0031] ③、计算混频信号
[0032] 数字混频计算，相当于锁相放大器中的混频过程，即将载波信号与参考信号相乘。为了求解一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，需要将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号。
[0033] ④、生成低通滤波器
[0034] 根据采用频率、低通截止频率和滤波阶次，创建低通滤波算法，即低通滤波器。本方法采用巴特沃斯(Butterworth)滤波器，因为其幅度响应曲线在通带内具有最平坦的特性，且在通带和阻带内幅度响应曲线是单调变化的。
[0035] ⑤、计算谐波信号
[0036] 将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)。并分别计算R量。引入采用零相位数字滤波算法，将消除谐波信号引滤波而产生的相位延迟，具体做法是将混频信号滤波运算、反向、再次滤波运算、再次反向。
[0037] 步骤二、调制光谱仿真
[0038] ①、计算线型函数
[0039] 气相分子的吸收光谱线型函数通常以Voigt线型函数描述，Voigt线型函数是Gauss线型函数和Lorentz线型函数的卷积,没有解析解,只能通过数值积分计算。
[0040]②、仿真谐波信号
[0041] 谐波信号仿真是数值仿真部分的核心，其目的是根据理论计算特定压力、浓度、温度条件下，某一谱线的谐波信号。谐波信号不仅与被测环境有关，还与谱线参数和激光参数有关，建立准确的理论模型，是温度与浓度测量的基础。
[0042] 激光穿过距离为L的均匀介质时，光强随传播距离呈指数衰减，可由Beer-Lambert 定律表不:
[0043] In —] = -KL = -PSu(Tmv)ZL ( I )[0044] 其中It和Ici分别是激光入射和透射光强，kv是光谱吸收系数，P是气体总压，X是待测气体组分浓度，Φ是线型函数，V是激光频率，Svi(T)是所用谱线Vi在温度T时的谱线强度。通常通过调制DFB激光器注入电流以改变激光的频率，同时激光器发射光强也随之改变，即方程(I)中的Itl为非定值。
[0045] 在扫描谱线的同时以高频信号调制激光波长，测量其谐波信号，即波长调制光谱。激光器的注入电流受到正弦调制，激光频率可以写成:
[0046]
[0047] 其中V是激光器中心频率，a是调制深度。Beer-Lambert定律改写成:
[0049] 其中IO是激光平均功率，10是归一化光强调制幅度，Ψ是频率调制和强度调制的相位差。方程(3)的指数部分反映了分子吸收信息，将其进行Fourier展开，可以表示为:
[0051] 透射光强受到调制和分子吸收共同作用，一次谐波(Rlf)、二次谐波(R2f)信号表示为:
[0054] 某条件下的线型函数和2f/lf信号如图2所示。
[0055] ③、构建仿真数据库
[0056] 在激光参数和谱线参数确定的条件下，谐波信号是温度和浓度的函数。在调制光谱中，温度是通过两条谱线谐波信号的强度之比来确定的，浓度则由已知温度条件的谐波信号强度确定。仿真数据库正是谐波信号强度比与温度和浓度对应关系的数据集，如图3所示。由于完整的数据库是三维的(谐波信号高度比、温度和浓度)，为了降低数据库的数据量和提供求解精度，每次仿真得到谐波信号高度比、温度，通过迭代仿真得到最终准确的温度和浓度。
[0057] 步骤三、测量参数提取步骤
[0058] ①、计算流场速度
[0059] 流场速度计算的任务是通过时间-频率对应关系将谐波信号的时间轴变换成频率轴，然后提取多普勒频移，最后根据光线交叉角度计算气流速度。为了准确获取谐波信号间的频率差，先将谐波信号按照最大值和两侧极小值在O〜I之间做规范化处理；对半高以上的信号进行等高插值，获得信号强度相同点的频率值；将上升沿和下降沿的等高点的频率差求平均，即是整个谐波信号的多普勒频移。多普勒频移提取过程见图4。
[0060] ②、信号高度提取
[0061] 调制光谱中，温度和浓度信息和谐波信号的高度相关，求解温度和浓度之前需要通过合适的算法准确提取谐波信号的高度。由于经过了低通滤波，谐波信号较光滑，谐波信号高度取最大值即可。
[0062] ③、温度、浓度求解
[0063] 在激光参数和谱线参数确定的条件下，谐波信号是温度和浓度的函数。在调制光谱中，温度是通过两条谱线谐波信号的强度之比来确定的，浓度则由已知温度条的谐波信号强度确定。如果建立谐波信号高度比，关于温度和浓度的三维数据库，则由测量信号可在数据库中查找温度和浓度值。但是这样导致数据库数据量的剧增，因为准确求解需要温度和浓度的递增量远小于实验误差。求解的算法是，设定初始浓度，建立谐波信号高度比和温度的对应关系，查询高度比得到温度；由温度谐波信号高度求浓度，并修正初始浓度，
[0064] 迭代该过程直至浓度和温度收敛。温度和浓度求解流程见图5。
Claims (1)
1.一种免标定调制光谱的实现方法，包括: 步骤一、数字锁相信号解调步骤； 步骤二、调制光谱仿真步骤； 步骤三、测量参数提取步骤； 其特征在于: 步骤一的数字锁相信号解调步骤包括: (1)通过数据采集设备读取探测光电信号、采样率和数据长度L ； (2)根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号； (3)将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号； (4)根据采样频率、低通截止频率和滤波阶次，创建低通滤波算法； (5)将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，并分别计算R量； 步骤二的调制光谱仿真步骤包括: (1)根据一定条件下的谱线参数进行数值积分计算，得到线型函数； (2)根据实测条件下的激光参数和计算得到的线型函数，计算谐波信号； (3)通过计算不同谱线参数对应的谐波信号，得到仿真数据库； 步骤三的测量参数提取步骤包括: (1)提取谐波信号的多普勒频移量，并计算流场速度； (2)通过提取谐波信号的最大值来获得谐波信号的高度值； (3)通过迭代算法，对比仿真谐波信号的高度值和实验谐波信号的高度值，求解流场的温度和组分浓度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410101717.5A CN103852443A (zh) | 2014-03-19 | 2014-03-19 | 免标定调制光谱的实现方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410101717.5A CN103852443A (zh) | 2014-03-19 | 2014-03-19 | 免标定调制光谱的实现方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103852443A true CN103852443A (zh) | 2014-06-11 |
Family
ID=50860356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410101717.5A Pending CN103852443A (zh) | 2014-03-19 | 2014-03-19 | 免标定调制光谱的实现方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103852443A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105548072A (zh) * | 2015-12-04 | 2016-05-04 | 东南大学 | 一种同时测量高温气体二维瞬态温度场和浓度场的方法 |
| CN106770024A (zh) * | 2017-01-19 | 2017-05-31 | 东南大学 | 基于免基线波长扫描直接吸收光谱的气体浓度测量方法 |
| CN106969800A (zh) * | 2017-05-03 | 2017-07-21 | 东南大学 | 一种利用单条谱线同时测量气体温度和浓度的装置和方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07280729A (ja) * | 1994-04-08 | 1995-10-27 | Mitsubishi Heavy Ind Ltd | ガス濃度測定装置 |
| WO2008112955A1 (en) * | 2007-03-14 | 2008-09-18 | Spectrasensors, Inc. | Pressure-invariant trace gas detection |
| CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
| CN103185704A (zh) * | 2011-12-27 | 2013-07-03 | 株式会社堀场制作所 | 气体计测装置和气体计测装置的波长调制幅度的设定方法 |
| CN103852444A (zh) * | 2014-03-19 | 2014-06-11 | 中国人民解放军装备学院 | 免标定调制光谱测量系统 |
-
2014
- 2014-03-19 CN CN201410101717.5A patent/CN103852443A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07280729A (ja) * | 1994-04-08 | 1995-10-27 | Mitsubishi Heavy Ind Ltd | ガス濃度測定装置 |
| WO2008112955A1 (en) * | 2007-03-14 | 2008-09-18 | Spectrasensors, Inc. | Pressure-invariant trace gas detection |
| CN103185704A (zh) * | 2011-12-27 | 2013-07-03 | 株式会社堀场制作所 | 气体计测装置和气体计测装置的波长调制幅度的设定方法 |
| CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
| CN103852444A (zh) * | 2014-03-19 | 2014-06-11 | 中国人民解放军装备学院 | 免标定调制光谱测量系统 |
Non-Patent Citations (3)
| Title |
|---|
| 李晗等: "可调谐二极管激光吸收光谱二次谐波信号的模拟与分析", 《光谱学与光谱分析》, 30 April 2013 (2013-04-30) * |
| 潘虎等: "免标定波长调制光谱在气体温度和浓度测量中的应用", 《红外与激光工程》, vol. 43, no. 3, 31 March 2015 (2015-03-31) * |
| 齐汝宾等: "准连续激光调制吸收谱的多谐波分析", 《光谱学与光谱分析》, 30 March 2012 (2012-03-30) * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105548072A (zh) * | 2015-12-04 | 2016-05-04 | 东南大学 | 一种同时测量高温气体二维瞬态温度场和浓度场的方法 |
| CN105548072B (zh) * | 2015-12-04 | 2019-06-21 | 东南大学 | 一种同时测量高温气体二维瞬态温度场和浓度场的方法 |
| CN106770024A (zh) * | 2017-01-19 | 2017-05-31 | 东南大学 | 基于免基线波长扫描直接吸收光谱的气体浓度测量方法 |
| CN106969800A (zh) * | 2017-05-03 | 2017-07-21 | 东南大学 | 一种利用单条谱线同时测量气体温度和浓度的装置和方法 |
| CN106969800B (zh) * | 2017-05-03 | 2019-10-11 | 东南大学 | 一种利用单条谱线同时测量气体温度和浓度的装置和方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102680020B (zh) | 一种基于波长调制光谱技术的气体参数在线测量方法 | |
| CN103852444A (zh) | 免标定调制光谱测量系统 | |
| CN102680428B (zh) | 一种基于一次谐波信号的气体温度和浓度在线测量方法 | |
| CN101105446B (zh) | 差分吸收光谱空气质量检测系统 | |
| CN103604774B (zh) | 基于非线性调谐提高激光气体分析灵敏度的方法和装置 | |
| CN102590138A (zh) | 一种基于激光吸收光谱的气体浓度在线测量方法 | |
| CN103323115A (zh) | 基于波长调制的气体吸收谱线线宽和线型系数的测量方法 | |
| CN103868884B (zh) | 一种基于调制系数的气体吸收率在线测量方法 | |
| CN103868885A (zh) | 基于复合多次谐波的气体浓度在线测量方法 | |
| CN103852443A (zh) | 免标定调制光谱的实现方法 | |
| CN104126106A (zh) | 一种用于气体感测的免校准扫描波长调制光谱的方法 | |
| CN109696415B (zh) | 一种基于快速傅里叶变换的气体吸收率在线测量方法 | |
| CN101320060A (zh) | 快速相位计 | |
| CN103886199B (zh) | 调制光谱信号的谐波小波分析方法 | |
| CN103645154A (zh) | 一种利用太赫兹光谱信号提取材料光学常数的方法 | |
| CN105259548A (zh) | 一种用于fmcw绝对距离测量技术中色散失配校正方法 | |
| CN104251820A (zh) | 一种基于Duffing系统的光声池微弱光声信号检测方法 | |
| MacMartin et al. | Frequency domain multimodel analysis of the response of Atlantic meridional overturning circulation to surface forcing | |
| CN102305940A (zh) | 流体因子提取方法 | |
| Chen et al. | Intercomparison between observed and simulated variability in global ocean heat content using empirical mode decomposition, part I: modulated annual cycle | |
| CN105629224A (zh) | 调频连续波雷达高精度测距方法 | |
| CN105954735A (zh) | 一种用于fmcw绝对距离测量技术中改进的高速色散失配校正方法 | |
| CN103196872A (zh) | 一种基于积分浊度计获取气溶胶粒子谱分布的方法 | |
| CN203872140U (zh) | 一种用于荧光信号解调的正交锁相放大器装置 | |
| CN207198034U (zh) | 太赫兹光谱分析系统及设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140611 |