[论文解读] Arctic Sea Ice State Estimation From Thermodynamic PDE Model

Shumon Koga, Miroslav Krstić|arXiv (Cornell University)|Jan 30, 2019

Arctic and Antarctic ice dynamics参考文献 24被引用 4

一句话总结

本文提出一种反推观测器，仅利用可测量的表层温度和冰层厚度数据，估计北极海冰中的三维温度分布。该方法相较于开环估计实现约十倍的更快收敛速度，具备理论稳定性保证，并通过包含盐度效应的热力学PDE模型数值仿真验证了性能。

ABSTRACT

Recent rapid loss of the Arctic sea ice motivates the study of the Arctic sea ice thickness. Global climate model that describes the ice's thickness evolution requires an accurate spatial temperature profile of the Arctic sea ice. However, measuring the complete temperature profile is not feasible within and throughout the Arctic icecap. Instead, measuring the ice's thickness is doable with the acquisition of data from submarine and satellite devices. In this paper, we develop a backstepping observer algorithm to estimate the temperature profile for the Arctic sea ice model via available measurements of sea ice thickness and sea ice surface temperature. The observer is designed in a rigorous manner to drive the temperature profile estimation error to zero, for a salinity-free sea ice model. Moreover, the proposed observer is used to estimate the temperature profile of the original sea ice model with salinity via numerical simulation. In comparison with the straightforward open-loop algorithm, the simulation results illustrate that our observer design achieves ten times faster convergence of the estimated temperature.

研究动机与目标

• 为解决由于热传感器有限，难以直接测量北极海冰空间分布温度分布的挑战。
• 开发一种仅使用可获取测量值（冰层厚度和表面温度）的状态估计方法，而非依赖冰-海洋界面处的温度梯度。
• 通过针对热力学PDE模型的严格反推设计，确保温度分布估计误差理论上收敛至零。
• 通过数值仿真在原始包含盐度的海冰模型上验证观测器性能，展示其相较于开环方法具有更优的收敛速度。

提出的方法

• 观测器采用适用于移动边界的抛物型PDE的反推方法设计，特别针对北极海冰热力学的Stefan问题公式。
• 观测器结构将冰层厚度和表面温度的边界测量作为输入，无需冰-海洋界面处的温度梯度。
• 通过解析推导观测器增益，确保在简化无盐度模型下估计误差的指数稳定性和收敛至零。
• 通过在真实物理参数和输入数据下对含盐度的原始模型进行数值仿真，将设计扩展至原始模型。
• 调参参数λ用于控制收敛速度与超调量之间的权衡，从而实现性能优化。
• 通过对比基于观测器的估计与开环估计，利用时间序列轮廓图和误差演化分析验证该方法。

实验结果

研究问题

• RQ1 能否设计一种反推观测器，仅使用表面温度和厚度测量数据，来估计北极海冰中的完整温度分布？
• RQ2 在缺乏温度梯度测量的情况下，所提出的观测器与开环估计方法相比，其收敛速度如何？
• RQ3 观测器增益参数（特别是λ）对收敛速度与估计超调量之间权衡的影响是什么？
• RQ4 当应用于包含盐度效应的原始热力学模型时，观测器设计在稳定性和准确性方面能保持到何种程度？

主要发现

• 所提出的反推观测器相较于开环估计，使估计温度分布的收敛速度提高约十倍。
• 观测器增益参数λ可有效调节收敛速度与超调量之间的权衡，λ = 5.0 × 10⁻⁶时实现性能平衡。
• 冰层厚度的估计误差初始为零，但因温度分布误差而增加，一旦温度估计收敛，误差即恢复至零，证实了精确温度估计的必要性。
• 数值仿真表明，简化模型下观测器的理论稳定性特性在应用于原始含盐度模型时仍可实际观测到。
• 即使在无法获取冰-海洋界面温度梯度的情况下，观测器仍能成功估计温度分布，而该参数在实际中难以测量。