多普勒激光雷达弧形扫描风速与风向预测区间的闭合统计表达式评估

《Remote Sensing》：A Closed-Form Statistical Expression for Evaluating Wind Speed and Direction Prediction Intervals from Doppler Lidar Arc Scans

　　

低仰角多普勒测风激光雷达扫描，即弧形扫描，将传统的垂直廓线测量扩展至水平离站位置。该技术旨在实现对仪器周边多个远距离位置的测量。弧形扫描方法已广泛应用于风能领域，同时也是环境监测中危害评估的有力工具。该方法引入了传统垂直廓线扫描所不具备的特定挑战。有限的扫描角度限制了可用于可靠矢量提取的风向观测数量。因此，在业务化配置中，建立一种可靠的方法来估算反演风速和风向的不确定度区间具有重要意义。在此，研究人员推导出一种闭合统计表达式，用于评估风速和风向的预测区间。由于快速更新的业务化场景（如实时扩散模拟）产生的扫描次数有限，该框架专门设计为在数学上保持稳健，且仅利用对角方差项即可计算，从而绕开了对数值不稳定的互协方差矩阵的需求。该表达式通过激光雷达与声学风速计的测量实验进行了验证。风与扫描弧的夹角成为影响风向和风速反演不确定度的主要参数。据此得出了有关扫描参数和选址考虑的结论。

**论文解读：基于多普勒激光雷达弧形扫描的风速与风向预测区间闭合统计表达式研究**

**一、研究背景与问题提出**

可靠的离站水平距离风场测量对于诸多环境和安全应用至关重要，包括风能资源评估、化学品或放射性物质泄漏的大气扩散模拟、航空安全以及城市微气象研究。传统原位测量手段，如安装在气象塔上的声学风速计，虽精度较高，但空间覆盖有限、规模化部署成本高昂，且在偏远或海上区域往往难以实施。多普勒测风激光雷达通过探测大气中随风运动的气溶胶粒子的径向速度，并借助合适的扫描几何结构重构风矢量，为上述问题提供了一种强有力的替代方案。速度方位显示（VAD）和 Doppler 光束摆动（DBS）等扫描技术可重构仪器上方的垂直风廓线，而弧形扫描则通过固定较低仰角的方位角扫描，实现对特定水平距离区域水平风场的测量。尽管弧形扫描在风能等领域应用日益广泛，但对其反演风矢量的不确定度进行严格定量评估的工具仍然匮乏，现有研究多依赖经验对比或理想化模拟，缺乏解析模型。这一方法论空白严重制约了弧形扫描技术在需要快速、多点位风场评估的业务化应用（如危险物质扩散风险实时评估）中的可靠性。

**二、研究目标与总体结论**

为填补上述空白，本研究旨在开发并验证一个闭合解析框架，用于估算由弧形扫描反演所得风速和风向的预测区间（PIs）。研究人员基于广义最小二乘法和一阶误差传播理论，推导出风速和风向不确定度的解析表达式。研究发现，反演不确定度与线视向（LOS）径向速度的标准差、扫描几何结构（波束间距、弧跨度）、平均风速以及扫描弧与风向的相对夹角直接相关。通过部署在距离声学风速计数百米外的激光雷达实地测量实验，验证了该理论框架的有效性。实地数据证实了一个关键的几何权衡关系：当扫描弧与风向对齐时（相对夹角较小），风速反演的不确定度最小；而当扫描弧与风向垂直时（相对夹角较大），风向反演的不确定度最小。这一结论为环境监测中仪器最优选址和扫描配置提供了量化的、可操作的指导。论文发表于《Remote Sensing》。

**三、主要关键技术方法**

本研究采用的主要技术方法包括：1. **广义最小二乘（GLS）反演与协方差分析**：建立弧形扫描径向速度与水平风分量的线性关系，推导风分量协方差矩阵的解析表达式。2. **Delta方法（一阶误差传播）**：将笛卡尔风分量的协方差矩阵变换为风速和风向的方差表达式。3. **对角协方差近似正则化方法**：针对快速扫描场景下全协方差矩阵秩亏（不可逆）的问题，提出仅保留对角方差项的近似策略，以实现稳定的实时反演和闭合解析表达。4. **实地对比验证实验**：在以色列Sde-Boker校区开展为期五个月（2022年1月至5月）的测量实验，将Halo Photonics StreamLine XR多普勒激光雷达弧形扫描反演的10分钟平均风矢量，与375米外25米高处安装的RM Young 81000声学风速计观测数据进行系统对比。数据按风与扫描弧相对夹角（Θ）分区间进行评估。

**四、研究结果与分析**

**4.1 风速与风向反演验证**

通过对比激光雷达与声学风速计在不同Θ区间的观测数据，研究发现，随着Θ从平行向垂直方向增加，风速反演的相关系数（R2）呈中等程度下降，归一化均方根误差（NRMSE）相应增大。相反，风向反演的偏差虽无明显Θ依赖性，但其误差分布的四分位距（IQR）随Θ增加而减小，表明在近垂直对齐时风向反演的一致性更高。这一观测结果与理论预测的趋势定性吻合。

**4.2 风速和风向误差分析**

对观测误差分布的分析进一步证实了理论预测的几何权衡关系。随着Θ从平行（0-10°）向垂直（80-90°）转变，相对风速误差的概率分布逐渐展宽，而风向误差的概率分布则逐渐收窄。风速误差IQR随Θ增大而增大，风向误差IQR随Θ增大而减小，其趋势与理论模型预测（图3）定性一致。此外，二维直方图分析揭示了风速与风向误差之间存在耦合关系，在中等Θ值下最为显著，这源于同时从同一组LOS速度估计u和v分量时协方差矩阵非对角项的影响。

**4.3 对角协方差近似验证**

为评估忽略协方差矩阵非对角项的影响，研究利用40分钟观测窗口（确保全矩阵可逆）进行了对比分析。结果显示，对于相对风速不确定度，对角近似与全矩阵计算的中位差异仅为0.8%；对于风向不确定度，中位差异为-6°，但极端情况下（5-95百分位）差异可达-28.5°至7.5°。分析表明，相邻波束平均相关性（R）是导致差异的主要因素，高R值对应于差异较大的情况。进一步研究发现，R在强对流不稳定或极稳定大气条件下较高，分别与大尺度相干热力结构和水平二维相干结构有关。

**4.4 基于方差估计的预测区间**

预测区间（PI）的覆盖概率评估显示，无论是使用对角近似还是全协方差矩阵，风速和风向的经验覆盖概率（ECP）在k=1时略低于但接近标准正态分布参考值，而在k值较大时则显著低于参考值。分位数-分位数（Q-Q）图分析表明，两种不确定性估计公式均能较好地再现归一化误差分布的中心部分，但均呈现出比高斯分布更重的尾部特征。条件分析显示，风速不确定性覆盖受风速或湍流强度（TI）影响较小；而风向不确定性覆盖在强风条件下表现为过度覆盖，弱风条件下则覆盖不足。

**4.5 平均窗口时长对统计收敛和覆盖可靠性的影响**

通过评估10至50分钟不同观测窗口的ECP，研究发现10分钟短窗口下，方向依赖方差公式因样本量不足导致边界分散不足（风速覆盖约42%，风向约53%）。随着窗口延长至40-50分钟，覆盖概率向理论高斯期望值（68.3%）收敛（风速约61%，风向约66%）。采用空间均匀方差公式可有效缓解短窗口下因波束特定采样噪声导致的覆盖不足问题，将10分钟窗口的覆盖概率提升至约57%（风速）和约63%（风向）。固定时间尺度（40分钟）下的平稳段子采样分析显示，ECP随样本量增加向高斯期望系统收敛，但改善在15至20次扫描间趋于饱和，表明进一步增加扫描密度对方差估计的收益递减。

**五、总结与结论**

本研究推导并验证了用于估算多普勒激光雷达弧形扫描风速和风向预测区间的闭合解析框架。该框架揭示了风速与风向反演不确定度之间的几何权衡关系：风与扫描弧对齐时（小Θ角）最小化风速不确定度，垂直时（大Θ角）最小化风向不确定度。这一发现为仪器选址和扫描配置提供了量化指导。研究证实，在快速更新的业务化场景下，采用仅含对角方差项的正则化方法可有效实现稳定实时反演，且长时间平均下其经验覆盖概率与全协方差矩阵方法相比无显著差异。然而，在极端不稳定或稳定条件下，相邻波束径向速度的高相关性可能导致两种方法评估结果出现较大偏差。风向预测边界在强风条件下表现为过度覆盖，弱风下则略有不足。径向风速方差的统计收敛性对于预测边界的准确估计至关重要，短时间窗口下的收敛不足可通过应用均匀方差予以部分缓解。尽管当前单站点实地数据证实了模型对风-弧对齐依赖性的定性预测，但未来仍需在更复杂的扫描几何、地形类型和大气稳定度条件下进行系统性验证，并纳入相关的LOS方差，以全面确立该框架的稳健性和普适性。本研究为弧形扫描在需要快速、多点风场评估的业务化环境应用中的不确定度感知使用迈出了重要一步。