[发明专利]一种基于四阶龙格‑库塔和模拟退火的等值盐密预测方法

摘要

本发明提供一种基于四阶龙格‑库塔和模拟退火的等值盐密预测方法，包括以下步骤(1)对经典预测模型进行参数动态化处理，得到以一阶欧拉方程表示的离散化的经典预测模型；(2)引入四阶龙格—库塔算法对步骤(1)得到的欧拉方程进行高阶求导迭代以抑制欧拉模型的离散误差；(3)借助模拟退火算法对模型参数进行全局最优估计，从而实现高精度的等值盐密度预测模型，消除了经典等值盐密度累计规律公式在零时刻预测结果为负值的错误，避免了新误差的带入；有效的抑制了预测误差，获得更为精确的预测结果；解决非线性高阶模型拟合的困难，提高拟合速度，快速获取参数最优估计值；解决了经典预测模型参数固定的弊端，实现动态参数估计。

主权项

一种基于四阶龙格‑库塔和模拟退火的等值盐密预测方法，其特征在于，包括以下步骤：(1)对经典预测模型进行参数动态化处理，得到以一阶欧拉方程表示的离散化的经典预测模型；(2)引入四阶龙格—库塔算法对步骤(1)得到的欧拉方程进行高阶求导迭代以抑制欧拉模型的离散误差；(3)借助模拟退火算法对模型参数进行全局最优估计，从而实现高精度的等值盐密度预测模型；所述步骤(1)的经典预测模型为ESDD＝A×(1‑K×exp(‑t/τ))   (1)其中，ESDD为等值盐密度，A为饱和等值盐密度，与绝缘子的结构型号和运行环境有关；K为常数；t为积污时间；τ为表征积污速率的常数；对ESDD进行参数动态化处理的方法为切线型离散化处理；对ESDD进行离散化处理过程为：对式(1)进行切线型离散化处理，可得式(2)如下：Sk＝A(1‑exp(‑k/τ))   (2)采用k时刻的一阶近似Sk+1-Sk=A(exp(-kτ)-exp(-k+1τ))≈A×exp(-kτ)×1τ---(3)]]>因此，k+1时刻ESDD可以表示为Sk+1=Sk+Aτ·exp(-kτ)---(4).]]>