[发明专利]一种气溶胶颗粒物质量浓度的探测装置及探测方法在审

专利信息
申请号: 201710152960.3 申请日: 2017-03-15
公开(公告)号: CN106872324A 公开(公告)日: 2017-06-20
发明(设计)人: 何廷尧;饶志敏;华灯鑫 申请(专利权)人: 西安理工大学
主分类号: G01N15/06 分类号: G01N15/06
代理公司: 西安弘理专利事务所61214 代理人: 常娥
地址: 710048*** 国省代码: 陕西;61
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 气溶胶 颗粒 物质 浓度 探测 装置 方法
【权利要求书】:

1.一种气溶胶颗粒物质量浓度的探测装置,其特征在于,包括激光器(1)发射多个波段的激光脉冲,依次经由准直扩束系统(2)、反射镜(3)进入大气(4)中,产生的激光散射回波信号由望远镜(5)接收,望远镜(5)通过光纤(6)与分光系统(7)连接,分光系统(7)依次连接有光电探测器(8)、数据采集系统(9)和计算机(10)。

2.一种气溶胶颗粒物质量浓度的探测方法,其特征在于,具体包括以下步骤:

步骤1:激光器(1)发射出355nm、532nm和1064nm三个波段的激光脉冲,三个波段的激光脉冲经准直扩束系统(2)准直扩束后,由反射镜(3)反射进入大气(4)中;

步骤2:反射进入大气中的三个波段的激光脉冲与大气中的气溶胶颗粒物发射散射作用,产生的激光散射回波信号由望远镜(5)接收;

步骤3:由望远镜(5)接收的激光散射回波信号经光纤(6)引入分光系统(7)进行分光处理,得到355nm、532nm和1064nm三个波段的后向散射信号,分别由光电探测器(8)进行信号探测和光电转换;

步骤4:经光电转换之后产生的电信号由数据采集系统(9)进行数据采集,采集到的数据送入计算机(10)进行算法处理,得到单次散射的激光雷达功率方程;

步骤5,对步骤4得到的激光雷达功率方程通过反演算法得到气溶胶颗粒物质量浓度。

3.根据权利要求2所述的气溶胶颗粒物质量浓度的探测方法,其特征在于,步骤4中单次散射的激光雷达功率方程表述为消光系数σ(z,λ)与后向散射系数β(z,λ)的函数:

<mrow><mi>P</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mi>C</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mn>0</mn></msub><mo>&CenterDot;</mo><mfrac><mi>A</mi><msup><mi>z</mi><mn>2</mn></msup></mfrac><mo>&CenterDot;</mo><mi>Y</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&beta;</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mn>2</mn><msubsup><mo>&Integral;</mo><mn>0</mn><mi>z</mi></msubsup><mi>&sigma;</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><msub><mi>dz</mi><mn>1</mn></msub><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>

式中,P(z,λ)为距离z处的瞬时接收功率,C为激光雷达系统常数,P0激光发射功率,A望远镜接收面积,Y(z)系统几何重叠因子,σ(z,λ)为气溶胶消光系数,β(z,λ)为气溶胶后向散射系数。

4.根据权利要求3所述的气溶胶颗粒物质量浓度的探测方法,其特征在于,步骤5具体包括以下步骤:

步骤5.1:激光雷达功率方程由Fernald反演算法分别得到355nm、532nm和1064nm三个不同波段的气溶胶颗粒物消光系数:

<mrow><mi>&sigma;</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><msup><mi>z</mi><mn>2</mn></msup><msub><mi>T</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mn>3</mn><msub><mi>S</mi><mn>1</mn></msub><mo>/</mo><mn>4</mn><mi>&pi;</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>2</mn></mrow></msup><msub><mi>S</mi><mn>1</mn></msub></mrow><mrow><mi>C</mi><mi>E</mi></mrow></mfrac><mrow><mn>1</mn><mo>-</mo><mfrac><mrow><mn>2</mn><msub><mi>S</mi><mn>1</mn></msub></mrow><mrow><mi>C</mi><mi>E</mi></mrow></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>r</mi></msubsup><mi>P</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><msup><mi>z</mi><mn>2</mn></msup><msub><mi>T</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mn>3</mn><msub><mi>S</mi><mn>1</mn></msub><mo>/</mo><mn>4</mn><mi>&pi;</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>r</mi></mrow></mfrac><mo>-</mo><msub><mi>S</mi><mn>1</mn></msub><msub><mi>&beta;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>

式中,S1为大气气溶胶的雷达比,Tm(z)表示z处的大气分子的透射率,E表示激光雷达的发射脉冲能量;

步骤5.2:由多波段激光雷达探测到的气溶胶消光系数和后向散射系数信息,通过粒子谱反演算法,得到气溶胶颗粒物粒子谱的分布特征;

由米散射理论可知,气溶胶颗粒物的消光系数与其数密度有着直接的关系,波长为λ的气溶胶颗粒物光学参数与气溶胶颗粒物粒子谱之间满足第一类Fredholm积分方程:

<mrow><msub><mi>g</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>min</mi></msub><msub><mi>r</mi><mi>max</mi></msub></msubsup><mfrac><mn>3</mn><mrow><mn>4</mn><mi>r</mi></mrow></mfrac><msub><mi>Q</mi><mrow><mi>&beta;</mi><mo>/</mo><mi>e</mi><mi>x</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&lambda;</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mi>v</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>

式中,g表示光学参量,r气溶胶颗粒物粒子半径,m复折射率,λ波长,v(r)球形颗粒物粒子的体积分布函数,Qβ/ext(r,λ,m)是球形颗粒物粒子消光或者后向散射效率;Qβ/ext(r,λ,m)根据米散射理论计算得到,结合气溶胶颗粒物粒子谱n(z,r)得到消光系数σλ和后向散射系数βλ的函数表达式,

<mrow><msub><mi>&sigma;</mi><mi>&lambda;</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>min</mi></msub><msub><mi>r</mi><mi>max</mi></msub></msubsup><msub><mi>Q</mi><mrow><mi>&lambda;</mi><mo>,</mo><mi>e</mi><mi>x</mi><mi>t</mi></mrow></msub><mo>(</mo><mi>r</mi><mo>)</mo><msup><mi>&pi;r</mi><mn>2</mn></msup><mi>n</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mi>d</mi><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>

<mrow><msub><mi>&beta;</mi><mi>&lambda;</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>min</mi></msub><msub><mi>r</mi><mi>max</mi></msub></msubsup><msub><mi>Q</mi><mrow><mi>&lambda;</mi><mo>,</mo><mi>&beta;</mi></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><msup><mi>&pi;r</mi><mn>2</mn></msup><mi>n</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mi>d</mi><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>

式中,Qλ,ext和Qλ,β分别为气溶胶颗粒物消光效率和后向散射效率;

粒子体积分布函数v(r)通过加权因子Wj拟合样条曲线的方法得到,

<mrow><mi>v</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>W</mi><mi>j</mi></msub><msub><mi>B</mi><mi>j</mi></msub><mo>(</mo><mi>r</mi><mo>)</mo><mo>+</mo><msup><mi>&epsiv;</mi><mrow><mi>m</mi><mi>a</mi><mi>t</mi><mi>h</mi></mrow></msup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>

式中,Bj(r)为B样条曲线,εmath(r)为计算误差,K为样条函数的个数,因此通过加权因子Wj的求解获得气溶胶粒子谱分布;当计算误差εmath(r)趋于无限小时,得到球形粒子的体积分布,为:

<mrow><mi>v</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>W</mi><mi>j</mi></msub><msub><mi>B</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>

步骤5.3:将气溶胶颗粒物粒子消光效率与气溶胶颗粒物粒子谱分布特征做比值处理,得到气溶胶颗粒物质量消光效率(MEE);

颗粒物质量消光效率(MEE)表征的是对于任意给定的粒子谱分布和激光波长,单位体积内的消光系数与颗粒物质量浓度的比值,是粒子谱分布特性和消光效率的函数,表达式为

<mrow><mi>M</mi><mi>E</mi><mi>E</mi><mo>=</mo><mfrac><mrow><mi>&pi;</mi><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>min</mi></msub><msub><mi>r</mi><mi>max</mi></msub></msubsup><msup><mi>r</mi><mn>2</mn></msup><msub><mi>Q</mi><mrow><mi>e</mi><mi>x</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&lambda;</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mi>n</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mi>d</mi><mi>r</mi></mrow><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac><mi>&pi;</mi><mi>&rho;</mi><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>min</mi></msub><msub><mi>r</mi><mi>max</mi></msub></msubsup><msup><mi>r</mi><mn>3</mn></msup><mi>n</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mi>d</mi><mi>r</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>

式中,ρ为粒子标准密度;

步骤5.4:通过气溶胶消光系数与质量消光效率相结合获得颗粒物的质量浓度,如下式所示:

<mrow><mi>C</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow></mrow><mrow><mi>M</mi><mi>E</mi><mi>E</mi></mrow></mfrac><mi>g</mi><mo>/</mo><msup><mi>m</mi><mn>3</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow><mo>;</mo></mrow>

式中,σa(z,λ)为气溶胶消光系数。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安理工大学,未经西安理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201710152960.3/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top