[论文解读] A Particle Filter for Stochastic Advection by Lie Transport (SALT): A case study for the damped and forced incompressible 2D Euler equation

Colin J. Cotter, Dan Crisan+3 more|arXiv (Cornell University)|Jul 27, 2019

Model Reduction and Neural Networks参考文献 22被引用 14

一句话总结

本文提出一种结合退火与抖动的粒子滤波方法，用于将稀疏、噪声较大的观测数据同化到通过随机对流输运（SALT）方法构建的阻尼且受迫的不可压缩二维欧拉方程的随机模型降维中。该方法仅使用约4.9×10⁴个自由度的降维SALT模型，便成功近似了具有约3.1×10⁶个自由度的高分辨率模拟所表示的真实状态，在普通计算硬件上实现了5个大尺度涡旋周转时间内的合理精度。

ABSTRACT

In this work, we combine a stochastic model reduction with a particle filter augmented with tempering and jittering, and apply the combined algorithm to a damped and forced incompressible 2D Euler dynamics defined on a simply connected bounded domain. We show that using the combined algorithm, we are able to assimilate data from a reference system state (the ``truth") modelled by a highly resolved numerical solution of the flow that has roughly $3.1 imes10^6$ degrees of freedom, into a stochastic system having two orders of magnitude less degrees of freedom, which is able to approximate the true state reasonably accurately for $5$ large scale eddy turnover times, using modest computational hardware. The model reduction is performed through the introduction of a stochastic advection by Lie transport (SALT) model as the signal on a coarser resolution. The SALT approach was introduced as a general theory using a geometric mechanics framework from Holm, Proc. Roy. Soc. A (2015). This work follows on the numerical implementation for SALT presented by Cotter et al, SIAM Multiscale Model. Sim. (2019) for the flow in consideration. The model reduction is substantial: The reduced SALT model has $4.9 imes 10^4$ degrees of freedom. Results from reliability tests on the assimilated system are also presented.

研究动机与目标

• 解决在地球流体动力学中使用粒子滤波进行高维、非线性数据同化所面临的挑战。
• 评估在基于随机对流输运（SALT）的模型降维框架下，结合退火与抖动的粒子滤波方法的性能。
• 评估当模型存在系统性偏差（即降维SALT模型与真实情况存在偏差）时，滤波算法的鲁棒性与准确性。
• 通过秩直方图分析与集合发散诊断，检验后验集合的可靠性。

提出的方法

• 信号通过阻尼且受迫的不可压缩二维欧拉方程的SALT降维版本建模，通过由空间相关布朗运动驱动的伊藤型随机微分方程引入随机性。
• 粒子滤波采用带顺序重要性重采样的自举滤波器，并通过退火稳定粒子退化现象，通过抖动保持集合多样性。
• 通过逐步增强观测对似然函数的影响，对似然函数应用退火，以改善粒子权重分布并减少退化。
• 在重采样后通过抖动扰动粒子位置，以防止集合坍塌并维持其发散性。
• 模型降维将自由度从约3.1×10⁶（真实情况）减少至约4.9×10⁴（SALT模型），从而实现在普通硬件上的高效计算。
• 观测在离散时间步长上通过包含289至1089个观测点的稀疏空间网格进行同化，观测噪声按λ=10进行缩放。

实验结果

研究问题

• RQ1 结合退火与抖动的粒子滤波能否有效将数据同化到高维地球流体系统的SALT降维模型中？
• RQ2 SALT降维模型在长时间尺度上对阻尼且受迫的二维欧拉方程真实动力学的近似程度如何？
• RQ3 SALT框架中的模型误差在多大程度上影响后验集合的可靠性与准确性？
• RQ4 在存在模型偏差与高维状态空间的情况下，退火与抖动在提升粒子滤波稳定性与性能方面发挥何种作用？

主要发现

• 该联合算法仅使用100个粒子和普通计算资源，成功在5000个精细分辨率时间步长（对应5个大尺度涡旋周转时间）内近似了真实状态。
• 尽管存在模型偏差，后验集合均值在大多数网格点上仍与真实状态保持合理接近，仅在[0.75, 0.5]和[0.5, 0.75]等特定点出现短暂偏差。
• 先验集合发散随时间增加，而后验集合发散保持稳定，表明重采样与滤波性能有效。
• 秩直方图分析显示存在一定程度的偏斜，尤其在观测点较少时（81 vs. 289），表明集合分布存在欠分散性，可能反映模型偏差。
• 通过抖动维持了集合多样性，通过退火防止了粒子退化，即使在模型不完美条件下也保持有效。
• 可靠性测试表明，尽管系统存在一定程度的偏差与欠分散性，后验集合仍具备足够的鲁棒性，适用于实际数据同化。