| 主权项 |
一种微脉冲差分吸收激光雷达水汽时空分布反演方法，其特征在于：实现步骤如下：(1)测量水汽吸收峰探测波长λ<sub>1</sub>的线宽Δλ<sub>on</sub>和中心波长λ<sub>on</sub>，以及水汽非吸收峰探测波长λ<sub>2</sub>的精确线宽Δλ<sub>off</sub>和中心波长λ<sub>off</sub>；所述水汽吸收峰探测波长λ<sub>1</sub>具有宽线宽特征，中心波长λ<sub>on</sub>位于890‑980nm之间的某个水汽吸收峰上，线宽Δλ<sub>on</sub>大于2.5nm，能够覆盖水汽的多个吸收谱线；所述水汽非吸收峰探测波长λ<sub>2</sub>线宽Δλ<sub>off</sub>不限，为宽线宽或窄线宽的激光，中心波长λ<sub>off</sub>位于860‑890nm或者980‑1070nm范围内的水汽非吸收峰上；(2)利用大气分子光谱数据库HITRAN2008，即HITRAN Database提取水汽谱线线强，利用谱线卷积方法对HITRAN Database提取水汽谱线线强进行温度修正和谱线展宽修正，计算获取与测量同等条件，每个谱线不同高度的标准吸收截面：<img file="FDA0000697130980000011.GIF" wi="1276" he="80" />其中，i表示HITRAN Database提取不同波长的水汽谱线线强序列，T为绝对温度，ν为波数，S<sub>i</sub>(T)为温度和压力修正的线强函数，F<sub>iL</sub>(ν)为洛仑兹展宽修正函数，F<sub>iG</sub>(ν)多普勒展宽修正函数；δ<sub>i</sub>(ν,r)为修正后获取的标准吸收截面，其随高度不同而变化；(3)利用逐线积分方法，分别对所述步骤(1)测得水汽吸收峰探测波长线宽Δλ<sub>on</sub>和水汽非吸收峰探测波长线宽Δλ<sub>off</sub>谱线范围内，经步骤(2)修正后获取的标准吸收截面δ<sub>i</sub>(ν,r)进行逐线积分，分别获得不同高度水汽吸收峰探测光吸收截面δ(λ<sub>on</sub>,r)和非吸收峰探测光的吸收截面δ(λ<sub>off</sub>,r)：<img file="FDA0000697130980000012.GIF" wi="1243" he="126" />(4)由所述步骤(3)计算获得的不同高度水汽吸收峰探测光吸收截面δ(λ<sub>on</sub>,r)和非吸收峰探测光的吸收截面δ(λ<sub>off</sub>,r)，计算不同高度处的水汽差分吸收差分截面Δδ(r)：Δδ(r)＝δ(λ<sub>on</sub>,r)‑δ(λ<sub>off</sub>,r)(5)列出米散射激光雷达方程，气溶胶在水汽吸收峰和非吸收峰探测波长处对应的米散射激光雷达方程为：<img file="FDA0000697130980000013.GIF" wi="959" he="106" />P(r)是激光雷达接收到距离其r处的后向散射信号功率(W)，P<sub>t</sub>是激光发射功率(W)，k为激光雷达系统常数(W·km<sup>3</sup>·S<sub>r</sub>)，C(r)激光雷达几何校正因子，距离r是时间的函数，β(r)是距离r处的总后向散射系数(km<sup>‑1</sup>·S<sub>r</sub><sup>‑1</sup>)，其中β(r)＝β<sub>α</sub>(r)+β<sub>m</sub>(r)，β<sub>α</sub>(r)和β<sub>m</sub>(r)分别是大气气溶胶和大气分子在距离r处的后向散射系数，α(r)是距离r处的总消光系数(km<sup>‑1</sup>)，α(r)＝α<sub>α</sub>(r)+α<sub>m</sub>(r)，α<sub>α</sub>(r)和α<sub>m</sub>(r)分别表示大气气溶胶和大气分子在距离r处的消光系数；(6)对激光雷达过渡区数据进行几何因子校正，使得激光雷达探测盲区小于40m；对激光雷达方程进行距离修正，方程两边同时乘以地面到被测气溶胶粒子群高度的平方r<sup>2</sup><img file="FDA0000697130980000021.GIF" wi="1476" he="102" />(7)确定标定高度r<sub>c</sub>、大气分子消光系数α<sub>m</sub>(r)、气溶胶消光系数标定值α<sub>α</sub>(r<sub>c</sub>)、大气分子消光后向散射比S<sub>m</sub>、大气气溶胶消光后向散射比S<sub>α</sub>；根据Fernald提出激光雷达求解方法，建立求解水汽吸收峰探测波长对应中心波长λ<sub>on</sub>和水汽非吸收峰探测波长对应中心波长λ<sub>off</sub>的气溶胶消光系数和后向散射系数的方法；所述标定高度r<sub>c</sub>通过选取几乎不含大气气溶胶粒子的清洁大气层所在的高度来确定，即P(r)r<sup>2</sup>/β<sub>m</sub>的最小值对应的高度r选择为r<sub>c</sub>所述大气气溶胶消光后向散射比S<sub>α</sub>，依赖于探测激光波长、大气气溶胶粒子的尺度谱分布和折射指数，数值范围在0到90之间；S<sub>α</sub>＝α<sub>α</sub>(r)/β<sub>α</sub>(r) (3)所述大气分子消光后向散射比S<sub>m</sub>由大气分子瑞利散射理论计算得到：S<sub>m</sub>＝α<sub>m</sub>(r)/β<sub>m</sub>(r)＝8π/3 (4)所述大气分子消光系数α<sub>m</sub>(r)和大气分子后向散射系数β<sub>m</sub>(r)，通过公式获得：<img file="FDA0000697130980000022.GIF" wi="1158" he="246" />式(5)中，λ为探测激光波长，所述气溶胶消光系数标定值α<sub>α</sub>(r<sub>c</sub>)，先由大气气溶胶散射比获得标定高度r<sub>c</sub>处的气溶胶后向散射系数边界值β<sub>α</sub>(r<sub>c</sub>)，根据公式(3)计算气溶胶消光系数标定值α<sub>a</sub>(r<sub>c</sub>)，大气气溶胶散射比关系式为：1+β<sub>α</sub>(r<sub>c</sub>)/β<sub>m</sub>(r<sub>c</sub>)＝1.08 (6)根据Fernald提出激光雷达求解方法，所述大气气溶胶消光系数，由标定高度r<sub>c</sub>处以下的大气气溶胶后向积分消光系数方程反演获取，大气气溶胶后向积分消光系数方程：<img file="FDA0000697130980000031.GIF" wi="1864" he="345" />所述大气气溶胶后向散射系数根据所述公式(3)大气气溶胶消光后向散射比的定义进行反演计算；(8)根据所述步骤(7)中所述大气气溶胶消光系数反演方法，反演获得水汽吸收峰探测波长对应中心波长λ<sub>on</sub>处的气溶胶消光系数α<sub>α</sub>(λ<sub>on</sub>,r)和水汽非吸收峰探测波长对应中心波长λ<sub>off</sub>处的气溶胶消光系数α<sub>α</sub>(λ<sub>off</sub>,r)，再根据所述步骤(7)大气气溶胶消光后向散射比的定义分别获取水汽吸收峰和非吸收峰探测波长对应中心波长处的气溶胶后向散射系数β<sub>α</sub>(λ<sub>on</sub>,r)和β<sub>α</sub>(λ<sub>off</sub>,r)；(9)列出水汽差分吸收激光雷达方程：<img file="FDA0000697130980000032.GIF" wi="1936" he="140" />其中P(λ<sub>on</sub>,r)和P(λ<sub>off</sub>,r)分别为所述水汽吸收峰探测波长和水汽非吸收峰探测波长的探测光在高度r处的大气回波强度，由微脉冲激光雷达实际测量获得；α<sub>α</sub>(λ<sub>on</sub>,r)、α<sub>α</sub>(λ<sub>off</sub>,r)、β<sub>α</sub>(λ<sub>on</sub>,r)和β<sub>α</sub>(λ<sub>off</sub>,r)分别为所述探测波长对应中心波长处的气溶胶消光系数和后向散射系数；(10)根据步骤(9)中公式(8)，由所述步骤(4)中计算得到的不同高度处的水汽差分吸收差分截面Δδ(r)的计算结果和所述步骤(7)对α<sub>α</sub>(λ<sub>on</sub>,r)、α<sub>α</sub>(λ<sub>off</sub>,r)、β<sub>α</sub>(λ<sub>on</sub>,r)和β<sub>α</sub>(λ<sub>off</sub>,r)的反演结果，对微脉冲激光雷达实测数据进行数据反演计算获得水汽浓度廓线；(11)根据所述步骤(6)微脉冲激光雷达探测盲区小于40m，可假设微脉冲激光雷达自动温湿度记录仪所在高度与微脉冲激光雷达有效数据高度间的水汽均匀，利用激光雷达自动温湿度记录仪测量的水汽浓度值对所述步骤(10)反演得到水汽浓度廓线进行比对和标定，实现对水汽时空分布的反演和标定。 |
