发明名称 基于气溶胶扰动校正的差分吸收激光雷达CO<sub>2</sub>探测方法
摘要 本发明公开了一种基于气溶胶扰动校正的差分吸收激光雷达CO<sub>2</sub>探测方法,首先,在CO<sub>2</sub>浓度反演模型中考虑了气溶胶扰动影响;其次,进行1064<i>nm</i>和1572<i>nm</i>波长的米散射激光雷达的同时观测,获取1064<i>nm</i>和1572<i>nm</i>波长下气溶胶的消光系数和后向散射系数的线性经验关系;最后,进行常规双波长交替发射的差分吸收激光雷达CO<sub>2</sub>探测时,采用1064<i>nm</i>波长的米散射激光雷达连续探测,基于消光系数和后向散射系数的线性经验关系获得<i>on</i>波长和<i>off</i>波长下的消光系数和后向散射系数,结合CO<sub>2</sub>浓度反演模型获得CO<sub>2</sub>浓度反演结果。本发明反演结果更精确,在差分吸收激光雷达探测大气痕量气体领域具有很好的应用前景。
申请公布号 CN105510260A 申请公布日期 2016.04.20
申请号 CN201610031555.1 申请日期 2016.01.18
申请人 武汉大学 发明人 马昕;刘刚;龚威;朱忠敏;韩舸;相成志;梁艾琳;郑玉新;张腾
分类号 G01N21/31(2006.01)I;G01S7/497(2006.01)I 主分类号 G01N21/31(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 胡艳
主权项 一种基于气溶胶扰动校正的差分吸收激光雷达CO<sub>2</sub>探测方法,其特征是,包括:步骤1,考虑气溶胶扰动的影响,提出基于气溶胶扰动校正的CO<sub>2</sub>浓度反演模型:<math><![CDATA[<mrow><msub><mi>N</mi><mrow><msub><mi>co</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mn>1</mn></msub><msub><mi>r</mi><mn>2</mn></msub></msubsup><mo>&lsqb;</mo><msub><mi>&sigma;</mi><mrow><msub><mi>CO</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mrow><msub><mi>CO</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mi>d</mi><mi>r</mi></mrow></mfrac><mo>&lsqb;</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&Delta;</mi><mi>r</mi></mrow></mfrac><mi>ln</mi><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&Delta;</mi><mi>r</mi></mrow></mfrac><mi>ln</mi><mfrac><mrow><mi>&beta;</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mi>&beta;</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mi>&beta;</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>&beta;</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mn>1</mn></msub><msub><mi>r</mi><mn>2</mn></msub></msubsup><mo>&lsqb;</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mo>)</mo><mo>&rsqb;</mo><mi>d</mi><mi>r</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000908904560000011.GIF" wi="1892" he="143" /></maths>,其中,<img file="FDA0000908904560000013.GIF" wi="156" he="76" />为CO<sub>2</sub>在高度r处的浓度;Δr表示大气层厚度,Δr=r<sub>2</sub>‑r<sub>1</sub>,r<sub>1</sub>为较低高度,r<sub>2</sub>为较高高度,<img file="FDA0000908904560000012.GIF" wi="211" he="116" /><img file="FDA0000908904560000014.GIF" wi="234" he="64" />和<img file="FDA0000908904560000015.GIF" wi="240" he="64" />分别为高度r处on波长和off波长下CO<sub>2</sub>气体的吸收截面积;α(r,λ<sub>on</sub>)和α(r,λ<sub>off</sub>)分别为高度r处采用on波长和off波长激光雷达时气溶胶的消光系数;β(r<sub>1</sub>,λ<sub>on</sub>)和β(r<sub>1</sub>,λ<sub>off</sub>)分别为高度r<sub>1</sub>处采用on波长和off波长激光雷达时气溶胶的后向散射系数;β(r<sub>2</sub>,λ<sub>on</sub>)和β(r<sub>2</sub>,λ<sub>off</sub>)分别为高度r<sub>2</sub>处采用on波长和off波长激光雷达时气溶胶的后向散射系数;步骤2,进行1064nm波长的米散射激光雷达探测,获得米散射激光雷达在各高度的回波信号强度;同时,将CO<sub>2</sub>探测差分吸收激光雷达发射波长固定为off波长进行探测,获得1572nm波长米散射激光雷达在各高度的回波信号强度;利用Fernald激光雷达反演法反演气溶胶的消光系数,经线性拟合得两波长下消光系数的线性经验关系;步骤3,结合激光雷达比和步骤2反演的消光系数,获得气溶胶的后向散射系数,经线性拟合得两波长下后向散射系数的线性经验关系;步骤4,进行on和off双波长交替发射的差分吸收激光雷达CO<sub>2</sub>探测时,采用1064nm波长的米散射激光雷达连续探测,利用Fernald激光雷达反演法反演1064nm波长下气溶胶的消光系数,同时获得对应的后向散射系数;结合消光系数和后向散射系数的线性经验关系,获得1572nm波长下气溶胶的消光系数和后向散射系数,按时刻将1572nm波长下的消光系数和后向散射系数区分为on波长或off波长下的消光系数和后向散射系数;步骤5,将on波长和off波长下的消光系数和后向散射系数带入CO<sub>2</sub>浓度反演模型,获得CO<sub>2</sub>浓度反演结果。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学