绿色债券与金融市场之间的跨分位数依赖关系：一种交叉分位数图方法

《International Journal of Financial Studies》：Cross-Quantile Dependence Between Green Bonds and Financial Markets: A Cross-Quantilogram Approach

　　

为了探究过冷水滴（SLD）条件下FENSAP-ICE数值模拟的显著系统性偏差，本文提出了一种结合数值先验和领域知识辅助特征的岭回归校正方法，用于冰形几何参数预测。FENSAP-ICE对八个几何冰形参数的预测作为数值先验被纳入机器学习模型，将学习目标从“从头预测”转变为“校正系统性偏差”。基于SLD结冰物理机制构建了六个工程辅助特征，为小样本建模提供了具有物理意义的可调维度。针对下表面结冰极限的零冰情况，设计了结合逻辑回归分类和岭回归的两阶段模型。通过留一交叉验证，在29组NACA0012翼型SLD结冰风洞实验数据上进行评估，改进系统将总冰面积的对称平均绝对百分比误差（sMAPE）从79.47%降至34.65%，下表面冰角从105.19%降至28.08%，上表面结冰极限从61.75%降至35.48%，平均冰厚从42.12%降至20.89%，均优于FENSAP-ICE预测。消融实验进一步揭示，仅引入数值先验即可将预测误差降低约10个百分点，是性能提升的主要驱动力。所提方法计算成本低、物理一致性高，为小样本条件下SLD冰形预测提供了实用的偏置校正框架。此外，为解决小样本交叉验证可能带来的乐观偏差，还采用了嵌套交叉验证和多个线性基线模型，以验证校正方法的鲁棒性和相对竞争力。

**论文解读：基于数值先验与领域知识辅助特征的SLD结冰数值模拟偏置校正方法**

**研究背景与问题**

飞机结冰是航空安全中的持续威胁，尤其在过冷水滴（Supercooled Large Droplet, SLD）条件下，水滴中值体积直径（Median Volumetric Diameter, MVD）超过50 μm，发生显著变形、破碎和溅射，导致结冰问题更为复杂和危险。Broeren等人研究表明，当冰层扩展到防冰保护区域之外时，升力损失可在数分钟内超过30%，凸显了准确预测SLD冰形的紧迫性。工程实践中，主流数值结冰代码FENSAP-ICE虽已广泛验证，但在SLD条件下存在系统性低估，特别是上、下表面结冰极限，直接采用会导致防冰区域设计不足、保护区域后方出现严重残余结冰。现有物理信息机器学习（Physics-Informed Machine Learning, PIML）方法主要通过嵌入控制方程约束或进行领域迁移来缓解实验数据稀缺，但尚未系统应用于校正SLD结冰模拟的偏置。为填补这一空白，本研究将FENSAP-ICE预测视为数值先验，采用岭回归学习其与结冰风洞测量之间的系统性偏差，从而以低计算成本校正关键冰形参数。

**研究内容与结论**

研究人员针对SLD条件下结冰数值模拟的系统性偏差，提出了一种改进的冰形预测系统，该系统基于岭回归和物理先验，在29组NACA0012翼型SLD结冰风洞实验数据上通过留一交叉验证（Leave-One-Out Cross-Validation, LOOCV）进行评估。主要结论如下：引入数值先验是性能提升的主要驱动力，消融实验表明可将平均冰厚预测误差降低约10个百分点；在此基础上构建的六个工程辅助特征可进一步降低局部冰形参数（如上表面最大厚度）误差约2个百分点，实现基线校正与局部精细化。与原始FENSAP-ICE预测相比，改进系统将总冰面积（A_total）的sMAPE从79.47%降至34.65%，下表面冰角（θ_l）从105.19%降至28.08%，上表面结冰极限（x_u）从61.75%降至35.48%，下表面结冰极限（x_l）从59.00%降至37.21%，平均冰厚（t_avg）从42.12%降至20.89%。通过嵌套交叉验证、与替代线性基线模型比较以及系数稳定性分析，确认了校正的鲁棒性。该方法通过校正结冰极限和总冰面积，降低了原始FENSAP-ICE预测中防冰区域设计不足的风险，且由于校正作为计算成本低的后续处理步骤，可直接嵌入现有结冰模拟工作流，适用于飞机防冰和除冰系统的初步设计。该论文发表在《International Journal of Financial Studies》。注意：原文标题与内容不符，但解读基于实际论文内容。

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

研究人员搭建了基于岭回归的偏置校正框架，核心方法包括：（1）数值先验引入：将FENSAP-ICE预测的八个几何参数作为输入特征，使模型学习从“预测”转变为“校正”；（2）领域知识辅助特征：基于SLD结冰物理机制构建六个无量纲工程特征，涵盖液滴惯性、热力学冻结模式、破碎趋势及上下表面不对称性；（3）两阶段建模策略：针对下表面结冰极限（x_l），先采用逻辑回归分类判断是否存在可测量冰层，再对含冰样本进行岭回归回归；（4）模型评估：采用留一交叉验证（LOOCV）进行参数选择与性能评估，并辅以嵌套交叉验证以量化乐观偏差。实验数据来源于中国空气动力研究与发展中心FL-16结冰风洞，共29组NACA0012翼型SLD测试条件，涉及攻角、液态水含量（Liquid Water Content, LWC）、静温、MVD和结冰时间等参数。

**研究结果**

**结果与分析**：基于29组NACA0012翼型SLD结冰风洞实验数据，通过LOOCV程序对改进系统进行全面评估。与FENSAP-ICE预测相比，改进系统在所有关键参数上均降低了sMAPE。对于总冰面积（A_total），sMAPE由79.47%降至34.65%；对于下表面冰角（θ_l），由105.19%降至28.08%；对于上表面结冰极限（x_u），由61.75%降至35.48%；对于下表面结冰极限（x_l），由59.00%降至37.21%；对于平均冰厚（t_avg），由42.12%降至20.89%。散点图显示预测值紧密围绕理想对角线，没有单调偏差或结构性偏置，证实岭回归框架结合数值先验和工程辅助特征有效校正了FENSAP-ICE的系统性偏差。全局参数（如A_total、t_avg）的预测鲁棒性最为突出，局部厚度参数（t_u、t_l）的预测精度也保持较高水平。结冰极限参数（x_u、x_l）的散点分布体现了差异化建模策略的效果，冰角参数（θ_u、θ_l）的散点结构反映了上下表面冰角形成机制的不同。

**消融实验**：沿三个核心维度（局部冰厚、结冰范围、总冰体积）进行消融研究，选择上表面最大厚度（t_u）、上表面结冰极限（x_u）和平均冰厚（t_avg）为代表参数，设计三种特征配置：M1（纯数据驱动，仅含原始条件参数）、M2（引入FENSAP数值先验特征）、M3（完整特征集，M2基础上增加六个工程特征）。对比M1与M2，引入FENSAP-ICE预测后，t_u的sMAPE降低5.71个百分点，x_u降低4.18个百分点，t_avg降低10个百分点，表明数值先验是主要性能驱动。采用M3完整特征集后，t_u的sMAPE进一步降至38.96%，x_u降至49.18%，t_avg改善至22.70%。工程特征对局部参数的调节作用更显著。

**扩展验证与鲁棒性分析**：进行了嵌套交叉验证以量化标准LOOCV评估的乐观偏差，结果表明嵌套评估与标准LOOCV的sMAPE基本一致，乐观偏差可忽略。将校正预测与Lasso、普通最小二乘（Ordinary Least Squares, OLS）和弹性网（ElasticNet）在相同交叉验证协议下进行比较，岭回归在多数目标参数上保持竞争力或更优。岭回归系数在29个交叉验证折中的稳定性分析显示，多数参数（尤其t_u、t_l、t_avg）的系数符号一致且变异系数（Coefficient of Variation, CV）远低于1，少数参数（如A_total、x_u、x_l的某些FENSAP先验系数、θ_u的LWC系数）的CV超过1，解释时应谨慎。

**结论与讨论**

研究结论部分翻译如下：针对过冷水滴条件下结冰数值模拟的系统性偏差，本文提出了一种基于岭回归和物理先验的改进冰形预测系统。通过29组NACA0012翼型结冰风洞实验数据的留一交叉验证评估，主要结论如下：（1）引入数值先验是性能提升的主要驱动力；消融实验表明，它将平均冰厚预测误差降低了约10个百分点。在此基础上构建的六个工程辅助特征可进一步降低局部冰形参数（如上表面最大厚度）的误差约2个百分点。两个组件分别实现了基线校正和局部精细化。（2）与原始FENSAP-ICE预测相比，改进系统将每个关键参数的sMAPE降低：总冰面积A_total从79.47%降至34.65%，下表面冰角θ_l从105.19%降至28.08%，上表面结冰极限x_u从61.75%降至35.48%，下表面结冰极限x_l从59.00%降至37.21%，平均冰厚t_avg从42.12%降至20.89%。（3）本研究受限于实验样本数量有限和速度条件覆盖范围；评估可能包含乐观偏差。未来工作将扩大数据规模，引入不确定性量化，并探索更直接捕获SLD特征微物理过程的特征。（4）通过嵌套交叉验证、与替代线性基线比较以及系数稳定性分析，确认了所提校正的鲁棒性。嵌套评估显示乐观偏差很小；校正对大多数目标参数仍具有竞争力；多数岭回归系数在各折中符号一致，少数系数幅度变化较大，解释时应谨慎。（5）通过校正结冰极限和总冰面积，所提方法降低了原始FENSAP-ICE预测中防冰区域设计不足的风险。由于校正是计算成本低的后续处理步骤，可直接嵌入现有结冰模拟工作流，因此非常适合飞机防冰和除冰系统的初步设计。