一种粒雪密度剖面模型的构建方法及粒雪密度的获取方法

本发明属于寒区积雪与冰川观测，更具体地说，特别涉及一种粒雪密度剖面模型的构建方法及粒雪密度的获取方法。

背景技术：

1、粒雪密度是冰川演化、寒区水文、冰盖物质平衡、极地遥感反演研究中的核心基础参数，粒雪垂直密度剖面刻画了积雪自表层向下压实致密化的演化规律，是冰川厚度估算、冰芯年代解析、积雪融水径流模拟不可或缺的输入条件。

2、当前获取粒雪密度剖面主要分为两类方式：第一类为现场冰芯钻孔取样实测，通过分层采集冰芯样本称重测算密度。该方式数据直观，但存在显著局限：钻孔作业成本高昂、野外作业条件恶劣，仅能获得离散深度点位数据，难以得到平滑连续的纵向剖面；同时钻孔数量有限，空间代表性不足，大范围冰川区域无法实现全覆盖监测。第二类为经验模型模拟方法，现有主流粒雪密度模型大多采用以深度为唯一自变量的单指数拟合模型，仅考虑自重压实带来的密度变化。

3、现有技术存在多项固有缺陷：

4、建模物理维度单一，忽略多要素耦合作用。粒雪致密化过程是自重压实、颗粒重结晶、孔隙演变协同作用的结果，积雪存续时长、粒雪温度会显著改变重结晶速率，孔隙结构直接调控密度上限。传统单变量模型仅引入深度参数，未耦合时间、温度、孔隙条件，模型物理机理不完备，在存在季节性冻融的高山冰川、近极地区域拟合误差大幅上升。

5、缺乏剖面分层差异化校正机制。粒雪垂直剖面自上而下压实速率存在明显分层特征：浅层粒雪受大气冻融扰动强烈、压实速率快；中层压实平稳过渡；深层孔隙逐步闭合，密度趋于稳定。统一的全局拟合模型无法适配三层差异化演化规律，普遍出现浅层密度波动失真、层位交界预测跳变、深层预测密度无限持续增长的失真问题。

6、模型参数优化手段简单，易出现过拟合、欠拟合。传统模型依靠最小二乘拟合系数，全局寻优能力不足，难以适配多参量耦合非线性模型；同时缺少分层误差补偿方案，进一步放大剖面预测偏差。

7、模型静态固化，自适应能力弱。现有经验模型参数一旦确定长期不变，无法响应年际气候波动、降雪条件变化带来的粒雪压实规律改变，长期使用预测精度持续衰减，需要反复重新建模，工作量大。

8、因此，亟需提出一种兼顾多物理场耦合机理、分层误差校正、智能参数寻优的粒雪密度剖面模型构建方案，并配套形成一套高效、连续、可动态更新的粒雪密度获取方法，解决现有建模手段精度不足、适用性差、观测成本高的技术难题。

技术实现思路

1、为了解现有粒雪密度建模大多仅采用深度单变量拟合，忽略多物理参量耦合作用，缺少剖面分层校正机制，模型拟合失真严重；传统方式依靠冰芯离散取样，无法获取连续密度剖面，测算效率低、普适性差的技术问题，本发明提供一种粒雪密度剖面模型的构建方法及粒雪密度的获取方法，实现寒区粒雪垂直剖面密度高精度、连续化测算。

2、一种粒雪密度剖面模型的构建方法，包括以下步骤：

3、s1、采集目标积雪/冰川区域多维度实测数据，构建原始粒雪观测数据集，对数据集开展标准化预处理，得到标准化粒雪样本数据集；

4、s2、基于粒雪压实演化机理，构建包含深度衰减因子、时间累积因子、温度修正因子、孔隙补偿因子的多物理场耦合粒雪密度剖面初始理论模型；

5、s3、划分数据集，采用改进自适应粒子群优化算法对初始理论模型参数迭代寻优，确定最优基础模型参数组；

6、s4、将粒雪剖面划分为浅层快速压实层、中层平稳过渡层、深层致密恒定层，针对不同层位分别采用对应的残差校正策略开展分层补偿修正，构建得到高精度粒雪密度剖面模型。

7、优选的，步骤s1中标准化预处理包括采用3σ准则剔除异常值，区分短时缺失、连续缺失数据采用不同方式补全，并通过极值归一化统一各物理量纲。

8、优选的，步骤s2中多物理场耦合粒雪密度剖面初始理论模型表达式为：

9、。

10、优选的，步骤s3所述改进自适应粒子群优化算法引入自适应惯性权重与动态学习因子，设置双重迭代终止条件，以实测密度与模型计算密度均方误差最小化为优化目标。

11、优选的，步骤s4中浅层采用小波降噪残差修正，中层采用分段线性渐变残差补偿，深层采用恒定阈值残差校准。

12、一种粒雪密度的获取方法，采用上述粒雪密度剖面模型实现，包括：采集目标点位实时基础参数，输入粒雪密度剖面模型运算得到粒雪密度结果，对输出密度开展精度核验；周期性新增实测数据对模型进行动态更新优化。

13、与现有技术相比，本发明具有如下有益效果：

14、本发明构建多物理场耦合初始理论模型，综合考虑深度、时间、温度、孔隙结构耦合影响，摆脱传统单变量拟合局限，模型物理机理更加完备，区域普适性更强。本发明对粒雪垂直剖面实施分层差异化残差校正，针对性解决浅层密度波动、层间拟合突变、深层密度极值失真缺陷，实现全深度剖面高精度拟合。本发明采用改进自适应粒子群算法完成参数全局寻优，有效抑制模型过拟合、欠拟合问题，提升模型参数校准精度与稳定性。本发明配套完整的粒雪密度获取流程，无需大量冰芯钻孔取样，能够输出连续纵向密度剖面，大幅降低野外观测工作量。本发明设置明确的数据精度判定标准，对模型输出结果自动核验，保障密度测算结果可靠可控。本发明支持模型周期性动态迭代更新，可适配冻融循环、积雪年限变化带来的粒雪演化规律改变，长期保持测算精度。本发明能够适用于极地冰川、高山永久积雪、季节性积雪等多种场景，适用范围广，能够为寒区水文、冰川演化研究提供可靠数据支撑。

1.一种粒雪密度剖面模型的构建方法及粒雪密度的获取方法，其特征在于，包括如下步骤：

2.根据权利要求1所述的一种粒雪密度剖面模型的构建方法，其特征在于，步骤s1中标准化预处理的具体步骤包括：

3.根据权利要求1所述的一种粒雪密度剖面模型的构建方法，其特征在于，步骤s2中多物理场耦合粒雪密度剖面初始理论模型的计算公式为：

4.根据权利要求1所述的一种粒雪密度剖面模型的构建方法，其特征在于，步骤s3中改进自适应粒子群优化算法的优化过程包括：

5.根据权利要求1所述的一种粒雪密度剖面模型的构建方法，其特征在于，步骤s4中分层残差校正的具体实施方式为：

6.一种粒雪密度的获取方法，其特征在于，应用权利要求1-5任一项所述构建的粒雪密度剖面模型执行，包括如下步骤：

7.根据权利要求6所述的一种粒雪密度的获取方法，其特征在于，步骤s10中各参数的采集精度控制标准为：

8.根据权利要求6所述的一种粒雪密度的获取方法，其特征在于，步骤s30中精度核验的具体方法为：

9.根据权利要求6所述的一种粒雪密度的获取方法，其特征在于，还包括模型动态自适应迭代优化步骤：

10.根据权利要求6所述的一种粒雪密度的获取方法，其特征在于，所述获取方法适用于极地冰川、高山永久积雪、山地季节性积雪、高原积雪的全深度粒雪密度连续监测；能够同步输出单点瞬时密度数据、纵向连续密度剖面、区域积雪密度空间分布特征参数，克服传统人工取样离散、分段拟合误差大、无法实现全域连续精准测算的技术短板。